Stress analysis of plates with a circular hole reinforced by flange reinforcing member

This paper deals with the problem of stress analysis of plates with a circular hole reinforced by flange reinforcing member. The so called flange reinforcing member here means that the reinforcing member is built up by setting shapes or bars with any section shape on both sides of the plates along the edge of the hole. Two cases of external loads are considered. In one case the external loads are stressesσ_X^(∞),σ_Y^(∞),and τ_XY^(∞) acting at infinite point of the plate, and in the other the external loads are linear distributed normal stresses. The procedure of solving the problems mentioned above consists of three steps. Firstly, the reinforcing member is taken out from the plates and considered to be a circular bar being solved to determine its deformation under the action of radial force q₀(θ) and tangential force t₀(θ) which are forces acting upon each other between reinforcing member and plate. Secondly, the displacements of plate with a circular hole under the action of q₀(θ) and t₀(θ) and external loads are determined. Finally, forces q₀(θ) and t₀(θ) are obtained by the compatibility of deformations between reinforcing member and plate. Then the internal forces and displacements of reinforcing member and plate are deduced from q₀(θ) and t₀(θ) obtained.     

On forced vibration of anisotropic cylinders

Summary The dynamic response of a circular cylinder with thick walls of transverse curvilinear isotropy subjected to a uniformly distributed pressure varying periodically with time is analyzed by means of the Laplace transformation, and the exact solution is obtained in closed form. The previously obtained solutions for forced vibrations with isotropy, and free vibrations with transverse curvilinear isotropy are included as special cases of the general results reported here.Nomenclature t time - r, , z cylindrical coordinates - _ii components of normal strain - _ii components of normal stress - u radial displacement - c _ij elastic constant - mass density - c ² c ₁₁/ - ² c ₂₂/c ₁₁ - a, b inner, outer radius of the cylinder - , A, B constants - forced angular frequency - function defined by (9) - p, real, complex variables - constant defined by (14) - real number - , Lamé elastic constants - J (x) Bessel function of first kind - Y (x) Bessel function of second kind - I (x) modified Bessel function of first kind - K (x) modified Bessel function of second kind     

A direct method for constructing solutions of the Hamilton-Jacobi equation

Summary A new method that allows to establish the existence of global solutions of the classical Hamilton-Jacobi equation for a Hamiltonian H(y, x), periodic in (x₁,,x_d), is presented.

---

Sommario Si presenta un nuovo metodo che permette di Stabilire l'esistenza di soluzioni globali dell'equazione classica di Hamilton-Jacobi per Hamiltoniane H(y, x), periodiche in (x₁,,x_d).

     

The effect of background particulates on the delayed transition of a heated 9:1 ellipsoid

An experimental study was done to quantify the effects of a variety of background particulates on the delayed laminar-turbulent transition of a thermally stabilized boundary layer in water. A Laser-Doppler Velocimeter system was used to measure the location of boundary layer transition on a 50 mm diameter, 9:1 fineness ratio ellipsoid. The ellipsoid had a 0.15 m RMS surface finish. Boundary layer transition locations were determined for length Reynolds numbers ranging from 3.0 × 10⁶ to 7.5 × 106. The ellipsoid was tested in three different heating conditions in water seeded with particles of four distinct size ranges. For each level of boundary layer heating, measurements of transition were made for clean water and subsequently, water seeded with 12.5 m, 38.9 m, 85.5 m and 123.2 m particles, alternately. The three surface heating conditions tested were no heating, T = 10°C and T = 15°C where T is the difference between the inlet model heating water temperature, T _i, and free stream water temperature, T . The effects of particle concentration were studied for 85.5 m and 123.2 m particulates.The results of the study can be summarized as follows. The 12.5 m and 38.9 m particles has no measurable effect on transition for any of the test conditions. However, transition was significantly affected by the 85.5 m and 123.2 m particles. Above a length Reynolds number of 4 × 10⁶ the boundary layer transition location moved forward on the body due to the effect of the 85.5 m particles for all heating conditions. The largest percentage changes in transition location from clean water, were observed for 85.5 m particles seeded water.Transition measurements made with varied concentrations of background particulates indicated that the effect of the 85.5 m particles on the transition of the model reached a plateau between 2.65 particulates/ml concentration and 4.2 particles/ml. Measurements made with 123.3 m particles at concentrations up to 0.3 part/ml indicated no similar plateau.     

Lp-estimates for the nonlinear spatially homogeneous Boltzmann equation

This paper studies L^p-estimates for solutions of the nonlinear, spatially homogeneous Boltzmann equation. The molecular forces considered include inverse k^th-power forces with k > 5 and angular cut-off.The main conclusions are the following. Let f be the unique solution of the Boltzmann equation with f(v,t)(1 + ¦v²¦)^(s ₁ ^{+ /p)/2} L¹, when the initial value f ₀ satisfies f ₀(v) 0, f ₀(v) (1 + ¦v¦²)^(s ₁ ^{+ /p)/2} L¹, for some s₁ 2 + /p, and f ₀(v) (1 + ¦v¦²)^s/2 L^p. If s 2/p and 1 < p < , then f(v, t)(1 + ¦v¦²)^{(s s} ₁)^/2 L^p, t > 0. If s >2 and 3/(1+ ) < p < , thenf(v,t) (1 + ¦v¦²)^(s(s ₁ ^{+ 3/p))/2} L^p, t > 0. If s >2 + 2C₀/C₁ and 3/(l + ) < p < , then f(v,t)(1 + ¦v¦²)^s/2 L^p, t > 0. Here 1/p + 1/p = 1, x y = min (x, y), and C₀, C₁, 0 < 1, are positive constants related to the molecular forces under consideration; = (k – 5)/ (k – 1) for k^th-power forces.Some weaker conclusions follow when 1 < p 3/ (1 + ).In the proofs some previously known L-estimates are extended. The results for L^p, 1 < p < , are based on these L-estimates coupled with nonlinear interpolation.     

Thermal entrance region heat transfer for MHD laminar flow in parallel-plate channels with unequal wall temperatures

The problem of thermal entry heat transfer for Hartmann flow in parallel-plate channels with uniform but unequal wall temperatures considering viscous dissipation, Joule heating and axial conduction effects is approached by the eigenfunction expansion method. The series expansion coefficients for the nonorthogonal eigenfunctions are obtained by using a method for nonorthogonal series described by Kantorovich and Krylov [21]. Numerical results are obtained for the case with entrance condition parameter _o=1 and open circuit condition K=1. The parametric values of Ha=0, 2, 6, 10 and Br=0, –1 are considered for Hartmann and Brinkman numbers, respectively.

---

Zusammenfassung Das Problem der Wärmeübertragung im thermischen Einlauf einer Hartmannströmung im ebenen Spalt mit einheitlichen, aber ungleichen Wandtemperaturen wurde unter Berücksichtigung viskoser Dissipation, Joulescher Heizung und axialer Wärmeleitung mit Hilfe einer Entwicklung nach Eigenfunktionen behandelt. Die Koeffizienten der Entwicklung für nichtorthogonale Eigenfunctionen wurde nach einer Methode für nichtorthogonale Reihen nach Kantorovicz und Krylow [21] berechnet. Numerische Ergebnisse werden für den Eintrittsparameter _o=1 und die Bedingung für den offenen Stromkreis K=1 erhalten. Die Parameterwerte Ha=0, 2, 6, 10 und Br=0, –1 werden für die jeweiligen Werte der Hartmann- und der Brinckman-Zahl betrachtet.

---

Nomenclature a one-half of channel height - ¯B,B₀ magnetic field Induction vector and magnitude of applied magnetic field - Br Brinkman number, _f U_m ²/(k_c) - C_n,D_n coefficients in the series expansion of _e, see eq. 16 - c_p specific heat at constant pressure - ,E₀ electric field intensity vector and component - E_n,O_n even and odd eigenfunctions - Ha Hartmann number, (/_f)^1/2 B_o a - h₁,h₂ local heat transfer coefficients at lower and upper plates - ¯J,J_y electric current density vector and component - K external loading parameter, E_o/(B_o U_m) - k thermal conductivity - Nu₁, Nu₂ local Nusselt numbers, h₁,a/k and h₂a/k, respectively - P fluid pressure - Pe Peclet number, PrRe - Pr Prandtl number, C_{p f}/k - q₁,q₂ rates of heat transfer per unit area,^–k(T/Z)Z=–a'^–k(T/Z) Z=a respectively - Re Reynolds number, U_ma/u_f - T,T₀,T₁,T₂ fluid temperature, uniform entrance temperature, uniform but different lower and upper plate temperatures, respectively - T_b,T_m bulk temperature and (T₁+T₂)/2 - U,U_m,u axial, mean and dimensionless velocities, respectively - ¯V velocity vector - X,Z axial and transverse coordinates - x,z dimensionless coordinates - _n,_n even and odd eigenvalues - ,₀,_b dimensionless fluid, entrance and bulk temperatures, respectively - _c,_e,_f characteristic temperature difference (T₂-T_m), and dimensionless fluid temperatures, defined by eq. (10) - _e,_f magnetic permeability and viscosity of fluid - fluid density - electric conductivity - viscous dissipation function - (1-)/2     

An example of a quasiconvex function that is not polyconvex in two dimensions

We study the different notions of convexity for the function f () = ||² (||² – 2 det ) where ^2×2, introduced by Dacorogna & Marcellini. We show that f is convex, polyconvex, quasiconvex, rank-one convex, if and only if ¦¦ 2/3 2, 1, 1+ (for some >0), 2/3, respectively.     

The cross sectional area difference method — a new technique for determination of particle concentration by laser doppler anemometry

A new technique for the determination of particle concentration from the signals of a laser Doppler anemometer (LDA) is described. It is based on a statistical relation between the number of Doppler periods, or the amplitude of the Doppler signals, and the particle concentration. The technique allows the mass flux of the dispersed phase of a two-phase flow to be obtained from the data set of a conventional one-dimensional (ID) LDA. The technique has been called the cross sectional area difference method. Simulations and first experimental results are presented and discussed.List of symbols a, b, c half-axes of measurement control volume (mcv) - a ₁, b ₁, c ₁ half-axes of detection volume - c _L velocity of light - d _m beam waist diameter - d _p particle diameter - d _pc diameter of the calibration particle - d _pmin minimum detectable particle diameter - e elementary charge - h Planck's constant - i number of particle size classes - k wavenumber - m visibility - m refractive index - n(d _p ) particle concentration - n(d _pi ) concentration of ith particle class - n vector of n(d _pi ) - q exponent of size dependence of G(d _p ) - v _x x-velocity component - x fringe spacing - y ₀, z ₀ coordinates of particle trajectory and cross sectional area - A cross sectional area of mcv - A matrix of A ₁ - a ₁ cross sectional area of detection volume - A ₁ difference of neighbouring cross sectional areas - C _A normalisation constant for linear graduation of amplitude - C _N normalisation constant for Doppler periods - C _scat non-size-dependent factor of G(d _p ) - C _x normalisation constant for nonlinear graduation of amplitude - F() power spectral density - G(d _p ) integral scattering function - H number of accumulated counts - H _max maximum number of accumulated counts - I amplitude of Doppler signal - I _max I for a particle passing through the origin of the mcv - I _s trigger level - K logarithmic amplitude ratio - K _max logarithmic amplitude ratio for I _s - K _x degree of linear class width of amplitude - K _A degree of nonlinear class width of amplitude - N number of Doppler periods - N _m number of Doppler periods required by signal validation - N _max N for a particle passing through the origin of the mcv - N ₀ fringe number inside mcv along x-axis - P _L laser power - S ₀ particle arrival rate - S ₁ trigger rate - S ₁ contribution to trigger rate coming from A ₁ - S ₁ vector of S _1i - S _1i contribution to trigger rate coming from ith class of distribution - _Q quantum efficiency - wavelength of laser light - off-axis angle - elevation angle - angular frequency - beam intersection angle - phase difference     

Interactions in general continua with microstructure

This paper studies the behavior of the one dimensional Broadwell model of a discrete three velocity gas on a bounded domain 0 x 1 with specularly reflective boundary condition at x = 0, 1. For smooth initial data we show that the initial boundary value problem possesses a unique smooth solution which tends as t to a free state consisting of traveling waves f _1e (x – ct), f _2e (x + ct), f _3e (x) where each f _ie is 2-periodic. The convergence is in the weak* topology of an appropriate Orlicz-Banach state space. No smallness assumptions are made on the data.In memory of Ronald J. DiPerna     

A Strongly Degenerate Quasilinear Equation: the Parabolic Case

We prove the existence and uniqueness of entropy solutions of the Neumann problem for the quasilinear parabolic equation u_t=÷ a(u, Du), where a(z,)=f(z,), and f is a convex function of with linear growth as ||||, satisfying other additional assumptions. In particular, this class includes the case where f(z,)=(z)(), >0, and is a convex function with linear growth as ||||.     

On laminar flow through a uniformly porous pipe

Numerous investigations ([1] and [4–9]) have been made of laminar flow in a uniformly porous circular pipe with constant suction or injection applied at the wall. The object of this paper is to give a complete analysis of the numerical and theoretical solutions of this problem. It is shown that two solutions exist for all values of injection as well as the dual solutions for suction which had been noted by previous investigators. Analytical solutions are derived for large suction and injection; for large suction a viscous layer occurs at the wall while for large injection one solution has a viscous layer at the centre of the channel and the other has no viscous layer anywhere. Approximate analytic solutions are also given for small values of suction and injection.

Nomenclature

General r distance measured radially - z distance measured along axis of pipe - u velocity component in direction of z increasing - v velocity component in direction of r increasing - p pressure - density - coefficient of kinematic viscosity - a radius of pipe - V velocity of suction at the wall - r ²/a ² - R wall or suction Reynolds number, Va/ - f() similarity function defined in (6) - u ₀() eigensolution - U(0) a velocity at z=0 - K an arbitrary constant - B _K Bernoulli numbers Particular Section 5 perturbation parameter, –2/R - ² a constant, –K - x / - g(x) f()/ Section 6 perturbation parameter, –R/2 - ² a constant, –K - g() f() - g _c ()=g() near centre of pipe - * point where g()=0 Section 7 2/R - ² K - t (1–)/ - w(t, ) [1–f(t)]/ - ₀, ₁ constants - g() f()– ₀ - ₀/ - ₀ a constant - * point where f()=0     

On an analogue of the Euler-Cauchy polygon method for the numerical solution of ux y = f(x,y, u,ux, uy)

This paper develops, with an eye on the numerical applications, an analogue of the classical Euler-Cauchy polygon method (which is used in the solution of the ordinary differential equation dy/dx=f(x, y), y(x ₀)=y ₀) for the solution of the following characteristic boundary value problem for a hyperbolic partial differential equation u _xy =f(x, y, u, u _x , y _y ), u(x, y ₀)=(x), u(x ₀, y)=(y), where (x ₀)=(y ₀). The method presented here, which may be roughly described as a process of bilinear interpolation, has the advantage over previously proposed methods that only the tabulated values of the given functions (x) and (y) are required for its numerical application. Particular attention is devoted to the proof that a certain sequence of approximating functions, constructed in a specified way, actually converges to a solution of the boundary value problem under consideration. Known existence theorems are thus proved by a process which can actually be employed in numerical computation.

---

     

Berechnung des dynamischen Verhaltens von instationären Temperaturfeldern in zylindrischen Bauteilen

Zusammenfassung Im ersten Teil dieser Untersuchung wird zur Betrachtung des dynamischen Verhaltens instationärer Temperaturfelder in den Wandungen zylindrischer Rohre ein mathematisches Modell erstellt und mit Hilfe der Laplace-Transformation ausgewertet. Im einzelnen werden dabei die Übertragungsfunktionen der Rohrwandtemperaturen hergeleitet und für den Fall der Abweichung vom stationären Zustand unter dem Einfluß äußerer Störungen explizit dargestellt.Im zweiten Teil der Untersuchung wird das sich daraus ergebende dynamische Verhalten der Wandtemperatur fluiddurchströmter Rohre für einige Beispiele in Form von Ortskurven dargestellt.

---

Computation of the dynamic behaviour of unsteady-state temperature fields in cylindrical structures

In the first part of this paper a mathematical model is developed allowing the investigation of unsteady-state temperature fields in the walls of cylindrical pipes. Evaluation is done by means of Laplace-transformation. In particular the transfer function of the pipe wall temperature is derived, explicitly shown for the case of deviations from steady-state influenced by external disturbances. In the second part of this paper the resulting dynamical behaviour of the wall temperature of heat pipes containing a fluid is shown by means of Nyquist plots for several examples.

---

Formelzeichen a Temperaturleitzahl m²/sec - A, B, A^*, B^*, , ¯B Integrationskonstanten °C - ber, bei, ker, kei Kelvin-Funktionen - ber₁, bei₁, ker₁, kei₁ kelvin-Funklionerion - Bi Biot-Zahl - c spezifische Wärme kJ/kg K - F Übertragungsfunktion - i –1 (imaginäre Einheit) - I₀, K₀, I₁, K₁ modifizierte Bessel-Funktionen - N Nenner (Gl. (39)) - r Rohrradius m - R normierter Abstand von der Innenwand % - s (komplexe) Laplace-Variable 1/sec - t Zeitvariable sec - T Zeitkonstante sec - u Integrationsvariable (Gl. (15)) - Y₀₀, Y₁₀, Y₁₁ Hilfsfunktionen (Gl. (35)-(37)) - Wärmeübergangszahl kW/K m² - kleine Änderung - Laplace-Operator 1/m² - Umgebungstemperatur °C - Rohrwandtemperatur °C - Wärmeleitfähigkeit kW/K m - Dichte kg/m³ - (komplexe) Kennvariable (Gl. (11)) - Frequenz 1/sec - Variable (Gl. (45)) Indizes a Rohraußenwand - FDS Frischdampfsammelrohr - F Fluid - H Heizgas - i Rohrinnenwand - m Mittel - VD Verdampferrohr - W Rohrwand - 0 zum Zeitpunkt t=t₀ - -(Überstreichung) stationärer Zustand Herrn Prof. Dr.-Ing. R. Quack zum 65. Geburtstag gewidmet.     

Positively invariant regions for a problem in phase transitions

Positively invariant regions for the system v _t + p(W) _x = V _xx , W _t –V _x = W _xx are constructed where p < 0, w < , w > , p(w) = 0, w , > 0. Such a choice of p is motivated by the Maxwell construction for a van der Waals fluid. The method of an analysis is a modification of earlier ideas of Chueh, Conley, & Smoller [1]. The results given here provide independent L bounds on the solution (w, v).Dedicated to Professor James Serrin on the occasion of his sixtieth birthday     

The L1-Stability of Boundary Layers for Scalar Viscous Conservation Laws

Let v=v(x) be a non-trivial bounded steady solution of a viscous scalar conservation law u _t+f(u) _x =u _xx on a half-line R⁺, with a Dirichlet boundary condition. The semi-group of this IBVP is known to be contractive for the distance d(u, u)u–u₁ induced by L ¹(R⁺). We prove here that v is asymptotically stable with respect to d: if u ₀–vL ¹, then u(t)–v₁0 as t+. When v is a constant, we show that this property holds if and only if f(v)0. These results complement our study of the Cauchy problem [2].     

Symmetries and Global Solvability of the Isothermal Gas Dynamics Equations

We study the Cauchy problem associated with the system of two conservation laws arising in isothermal gas dynamics, in which the pressure and the density are related by the -law equation p() with =1. Our results complete those obtained earlier for >1. We prove the global existence and compactness of entropy solutions generated by the vanishing viscosity method. The proof relies on compensated compactness arguments and symmetry group analysis. Interestingly, we make use here of the fact that the isothermal gas dynamics system is invariant modulo a linear scaling of the density. This property enables us to reduce our problem to that with a small initial density.One symmetry group associated with the linear hyperbolic equations describing all entropies of the Euler equations gives rise to a fundamental solution with initial data imposed on the line =1. This is in contrast to the common approach (when >1) which prescribes initial data on the vacuum line =0. The entropies we construct here are weak entropies, i.e., they vanish when the density vanishes.Another feature of our proof lies in the reduction theorem, which makes use of the family of weak entropies to show that a Young measure must reduce to a Dirac mass. This step is based on new convergence results for regularized products of measures and functions of bounded variation.Acknowledgement P.G.L. and V.S. were supported by a grant from INTAS (01-868). The support and hospitality of the Isaac Newton Institute for Mathematical Sciences, University of Cambridge, where part of this research was performed during the Semester Program Nonlinear Hyperbolic Waves in Phase Dynamics and Astrophysics (January to July 2003) is also gratefully acknowledged. P.G.L. was also supported by the Centre National de la Recherche Scientifique (CNRS).     

Uniqueness and the maximum principle for quasilinear elliptic equations with quadratic growth conditions

In this paper, we show that the maximum principle holds for quasilinear elliptic equations with quadratic growth under general structure conditions.Two typical particular cases of our results are the following. On one hand, we prove that the equation (1) {ie77-01} where {ie77-02} and {ie77-03} satisfies the maximum principle for solutions in H ¹()L(), i.e., that two solutions u ₁, u ₂H¹() L() of (1) such that u ₁u₂ on , satisfy u ₁u₂ in . This implies in particular the uniqueness of the solution of (1) in H ₀ ¹ ()L().On the other hand, we prove that the equation (2) {ie77-04} where fH^–1() and g(u)>0, g(0)=0, satisfies the maximum principle for solutions uH¹() such that g(u)¦Du|{²L¹(). Again this implies the uniqueness of the solution of (2) in the class uH ₀ ¹ () with g(u)¦Du|{²L¹().In both cases, the method of proof consists in making a certain change of function u=(v) in equation (1) or (2), and in proving that the transformed equation, which is of the form (3) {ie77-05}satisfies a certain structure condition, which using ((v₁ -v ₂)⁺)ⁿ for some n>0 as a test function, allows us to prove the maximum principle.     

Some remarks on differential equations of quadratic type

In this paper we study differential equations of the formx(t) + x(t)=f(x(t)), x(0)=x ₀ C HereC is a closed, bounded convex subset of a Banach spaceX,f(C) C, and it is often assumed thatf(x) is a quadratic map. We study the differential equation by using the general theory of nonexpansive maps and nonexpansive, non-linear semigroups, and we obtain sharp results in a number of cases of interest. We give a formula for the Lipschitz constant off: C C, and we derive a precise explicit formula for the Lipschitz constant whenf is quadratic,C is the unit simplex inR ⁿ, and thel ₁ norm is used. We give a new proof of a theorem about nonexpansive semigroups; and we show that if the Lipschitz constant off: CC is less than or equal to one, then lim_tf(x(t))–x(t)=0 and, if {x(t):t 0} is precompact, then lim_tx(t) exists. Iff¦C=L¦C, whereL is a bounded linear operator, we apply the nonlinear theory to prove that (under mild further conditions on C) lim_t f(x(t))–x(t)=0 and that lim_t x(t) exists if {x(t):t 0} is precompact. However, forn 3 we give examples of quadratic mapsf of the unit simplex ofR ⁿ into itself such that lim_t x(t) fails to exist for mostx ₀ C andx(t) may be periodic. Our theorems answer several questions recently raised by J. Herod in connection with so-called model Boltzmann equations.     

Thermodynamic activation functions of viscous flow of non-polar liquids

The thermodynamic activation function of viscous flow may be determined from the expression for the pre-exponential factor in the Eyring relationship (the viscosity coefficient), which is a function of density and relative permittivity, together with the thermal dependence of the viscosity coefficient. This method of determination is demonstrated for a series of n-alkanes C₆–C₂₀.List of symbols A, B, C parameters in empirical viscosity-temperature dependence - E molecular bond energies of liquid calculated from heat of evaporation in vacuum - E _v viscous flow activation energy - G activation Gibbs function of viscous flow - H _v heat of evaporation - H _vb heat of evaporation at boiling point - H viscous flow activation enthalpy - k Boltzmann constant - m weight of molecule - M relative molecular weight - N _A Avogadro constant - q total number of bonds in one mole of liquid - R gas constant - S _id entropy of one mole of compound in ideal gas state - S _sat entropy of one mole of compound in saturated vapour state - S activation entropy of viscous flow - T absolute temperature - T _c critical temperature - T _r reduced temperature - T _c critical volume - T _f free volume of liquid per molecule - V ₀ molar volume of molecules - V _m molar volume of liquid - relative permittivity - viscosity - density - _r reduced density     

The Conjugacy of Stochastic and Random Differential Equations and the Existence of Global Attractors

We consider stochastic differential equations in d-dimensional Euclidean space driven by an m-dimensional Wiener process, determined by the drift vector field f₀ and the diffusion vector fields f₁,...,f_m, and investigate the existence of global random attractors for the associated flows . For this purpose is decomposed into a stationary diffeomorphism given by the stochastic differential equation on the space of smooth flows on R^d driven by m independent stationary Ornstein Uhlenbeck processes z¹,...,z^m and the vector fields f₁,...,f_m, and a flow generated by the nonautonomous ordinary differential equation given by the vector field (_t/x)^–1[f₀(_t)+ _i=1 ¹ f_i(_t)z _t ⁱ ]. In this setting, attractors of are canonically related with attractors of . For , the problem of existence of attractors is then considered as a perturbation problem. Conditions on the vector fields are derived under which a Lyapunov function for the deterministic differential equation determined by the vector field f₀ is still a Lyapunov function for , yielding an attractor this way. The criterion is finally tested in various prominent examples.