USRE41398E1  Design of viscoelastic coatings to reduce turbulent friction drag  Google Patents
Design of viscoelastic coatings to reduce turbulent friction drag Download PDFInfo
 Publication number
 USRE41398E1 USRE41398E1 US11054719 US5471905A USRE41398E1 US RE41398 E1 USRE41398 E1 US RE41398E1 US 11054719 US11054719 US 11054719 US 5471905 A US5471905 A US 5471905A US RE41398 E1 USRE41398 E1 US RE41398E1
 Authority
 US
 Grant status
 Grant
 Patent type
 Prior art keywords
 equation
 boundary
 viscoelastic
 surface
 coating
 Prior art date
 Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
 Active
Links
Images
Classifications

 F—MECHANICAL ENGINEERING; LIGHTING; HEATING; WEAPONS; BLASTING
 F15—FLUIDPRESSURE ACTUATORS; HYDRAULICS OR PNEUMATICS IN GENERAL
 F15D—FLUID DYNAMICS, i.e. METHODS OR MEANS FOR INFLUENCING THE FLOW OF GASES OR LIQUIDS
 F15D1/00—Influencing flow of fluids
 F15D1/10—Influencing flow of fluids around bodies of solid material
 F15D1/12—Influencing flow of fluids around bodies of solid material by influencing the boundary layer

 B—PERFORMING OPERATIONS; TRANSPORTING
 B63—SHIPS OR OTHER WATERBORNE VESSELS; RELATED EQUIPMENT
 B63B—SHIPS OR OTHER WATERBORNE VESSELS; EQUIPMENT FOR SHIPPING
 B63B1/00—Hydrodynamic or hydrostatic features of hulls or of hydrofoils
 B63B1/32—Other means for varying the inherent hydrodynamic characteristics of hulls
 B63B1/34—Other means for varying the inherent hydrodynamic characteristics of hulls by reducing surface friction

 B—PERFORMING OPERATIONS; TRANSPORTING
 B64—AIRCRAFT; AVIATION; COSMONAUTICS
 B64C—AEROPLANES; HELICOPTERS
 B64C21/00—Influencing airflow over aircraft surfaces by affecting boundarylayer flow
 B64C21/10—Influencing airflow over aircraft surfaces by affecting boundarylayer flow using other surface properties, e.g. roughness

 B—PERFORMING OPERATIONS; TRANSPORTING
 B64—AIRCRAFT; AVIATION; COSMONAUTICS
 B64C—AEROPLANES; HELICOPTERS
 B64C23/00—Influencing airflow over aircraft surfaces, not otherwise provided for

 Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSSSECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSSREFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
 Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
 Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
 Y02T50/00—Aeronautics or air transport
 Y02T50/10—Drag reduction
 Y02T50/16—Drag reduction by influencing airflow
 Y02T50/166—Drag reduction by influencing airflow by influencing the boundary layer

 Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSSSECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSSREFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
 Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
 Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
 Y02T70/00—Maritime or waterways transport
 Y02T70/10—Measures concerning design or construction of watercraft hulls
 Y02T70/12—Improving hydrodynamics of hull
 Y02T70/121—Reducing surface friction

 Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSSSECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSSREFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
 Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
 Y02T—CLIMATE CHANGE MITIGATION TECHNOLOGIES RELATED TO TRANSPORTATION
 Y02T70/00—Maritime or waterways transport
 Y02T70/10—Measures concerning design or construction of watercraft hulls
 Y02T70/12—Improving hydrodynamics of hull
 Y02T70/121—Reducing surface friction
 Y02T70/123—Hull coatings, e.g. biomimicry
Abstract
Description
This application claims the benefit of No. 60/138,023, filed Jun. 8, 1999.
Since M. O. Kramer reported successful experimental results in 1957, there have been repeated attempts to reduce frictional drag in turbulent fluid flow over a surface by applying a passive compliant coating. Experimental results in this area have been mixed. Most investigators have reported a drag increase, while only a few have claimed drag reduction for turbulent flow. A number of theoretical studies have characterized the stability of the laminar boundary layer over a deforming surface and other studies have characterized the reaction of a coating to a fluctuating load. However, no rigorous analytical technique has been previously reported that has been used to successfully design a dragreducing coating for turbulent flow.
In the past, passive coatings were tested without specification and full characterization of critical physical parameters, such as the frequency dependent complex shear modulus, density, and thickness. In order to achieve and ensure drag reduction with a viscoelastic coating, a methodology is required for selecting appropriate material properties and for estimating anticipated drag reduction as a function of configuration and velocity.
Relevant background information for associated technical topics is available in the literature, and may be useful due to the technical complexity of this invention. A classical discussion of boundary layer theory, including formulation of NavierStokes and turbulent boundary layer equations, is provided in BoundaryLayer Theory, by Dr. Hermann Schlichting, published by McGraw Hill, New York, seventh edition, 1979. A discussion of structures and scales in turbulent flows can be found in Turbulence, 1975, McGraw Hill, written by J. O. Hinze, and in “Coherent Motions in the Turbulent Boundary Layer,” in Annual Review of Fluid Mechanics, 1991, volume 23, pp. 60139, written by Steven K. Robinson. Background on Reynolds stress types of turbulence models is found in the chapter, “Turbulent Flows: Model Equations and Solution Methodology,” written by Tom Gatski, and included in the Handbook of Computational Fluid Mechanics, published by Academic Press in 1996. Equations in fluid and solid mechanics are often expressed in indicial, or tensor, notation, for compactness. Chapter 2 in the text A First Course in Continuum Mechanics, by Y. C. Fung, PrenticeHall, Inc., Englewood Cliffs, N.J., 1977, provides a brief introduction into tensor notation for mechanics equations. An introduction to finite difference methods, which are used to solve the system of momentum and continuity equations for a turbulent fluid, is provided in the text, Computational Fluid Dynamics for Engineers, written by Klaus Hoffman, and published in 1989 by the Engineering Education System i Austin in Austin, Texas. Descriptions of measured and mathematically modeled physical properties of polymers are found in the text, Viscoelastic Properties of Polymers by J. D. Ferry, Wiley, New York, 1980, 3^{rd }edition. The article, “Loss Factor Height and Width Limits for Polymer Relaxation,” by Bruce Hartmann, Gilbert Lee, and John Lee, in the Journal of the Acoustical Society of America Vol. 95, No. 1, January 1994, discusses mathematical characterization of shear moduli for real viscoelastic, polymeric materials, including those approximated by the HavriliakNegami approach.
Recently in the international literature (K. S. Choi, X. Yang, B. R. Clayton, E. J. Glover, M. Atlar, B. N. Semonev, and V. M. Kulik, “Turbulent Drag Reduction Using Compliant Surfaces,” Proceedings of the Royal Society of London, A (1997) 453, pp. 22292240). Choi et al. reported experimental measurements of up to 7% turbulent friction drag reduction for an axisymmetric body coated with a viscoelastic material. These experiments were performed in the United Kingdom, using coatings designed and fabricated in Russia at the Institute of Thermophysics, Russian Academy of Sciences, Novosibirsk, by a team headed by B .N. Semenov. The basic design approach was outlined in “On Conditions of Modelling and Choice of Viscoelastic Coatings for Drag Reduction,” in Recent Developments in Turbulence Management, K. S. Choi, ed., 1991, pp. 241262, Dordrecht, Kluwer Publishers. The Novosibirsk design approach is semiempirical in nature, and does not take into account the full characterization of the complex shear modulus of the viscoelastic material, namely, the relaxation time of the material. The Novosibirsk design approach does take into account frequencydependent material properties. Furthermore, the Novosibirsk concept is valid only for a membranetype coating, such as a film which coats a foamrubber saturated with water or glycerine, and where only normal fluctuations of the surface are considered.
The structure of coatings intended for drag reduction has been addressed in the international literature, starting with the 1938 patent No. 669897, “An Apparatus for the Reduction of Friction Drag,” issued in Germany to Max O. Kramer. Kramer later received a patent in 1964, U.S. Pat. No. 3,161,385, and in 1971, U.S. Pat. No. 3,585,953 for coatings to extend laminar flow in a boundary layer. Soviet inventor's certificates, such as “A Damping Covering,” USSR patent 1413286, Publication 20.01.1974, Bulletin of the Inventions 14, by V. V. Babenko, L. F. Kozlov, and S. V. Pershin, “An Adjustable Damping Covering,” USSR patent 1597866, Publication 15.03.1978, Bulletin of the Inventions 110, by V. V. Babenko, L. F. Kozlov, and V. I. Korobov, and “A Damping Covering for Solid Bodies,” USSR patent 1802672, Publication 07.02.1981, Bulletin of the Inventions 15, by V. V. Babenko and N. F. Yurchenko, have also described the structure of dragreducing coatings comprised of viscoelastic materials. These inventor's certificates identified the threedimensional structure within a dragreducing coating, but do not address the methodology for choosing appropriate parameters of the viscoelastic materials to be used in the manufacture of such coatings. Structural features include multiple layers of materials, longitudinal, riblike inclusions of elastic, viscoelastic, or fluid materials, and heated elements. Viscoelastic coatings may be combined with other forms of structure, such as longitudinal riblets molded on or within the surface of the coating. As described in the international literature in publications such as “Secondary Flow Induced by Riblets,” written by D. B. Goldstein and T. C. Tuan, and published in the Journal of Fluid Mechanics, volume 363, May 25, 1998, pp. 115152, twodimensional, rigid riblets alone have been shown experimentally to reduce surface friction drag up to about 10%.
The present invention enables the design of a passive viscoelastic coating for the reduction of turbulent friction drag. Coatings with material properties designed using the methodology described in this invention have reduced friction drag by greater than 10%. The methodology of the present invention permits, as a first object of the invention, the specification of the frequency dependent complex shear modulus, the density, and the thickness of an isotropic viscoelastic material which will reduce turbulent friction drag relative to specific flow conditions over a rigid surface. Quantitative levels of drag reduction can be estimated. Mathematical detail is provided for the cases of turbulent flow over a rigid flat plate as well as a viscoelastic flat plate, where the invention accounts for both normal and longitudinal oscillations of the surface. A second object of the invention is the specification of material properties for a coating composed of multiple layers of isotropic viscoelastic materials. A third object of the invention is the specification of material properties for a coating composed of an anisotropic material. A fourth object of the invention is the minimization of edge effects for coatings of finite length. A fifth object of the invention is the stabilization of longitudinal vortices through combination of viscoelastic coating design with additional structure, such as riblets.
The methodology used herein to describe the interaction of a turbulent boundary layer (TBL) with a viscoelastic (VE) layer involves two tasks, 1) a fluids task, involving the calculation of turbulent boundary layer parameters, given boundary conditions for a rigid, elastic, or viscoelastic surface (herein referred to as the TBL problem), and 2) a materials task, involving the calculation of the response of a viscoelastic or elastic surface to a periodic forcing function which approximates the loading of the turbulent boundary layer. The invention focuses upon cation of amplitudes of layer. The invention focuses upon calculation of amplitudes of surface oscillations and velocities, and of the energy flux for a viscoelastic coating (hereinafter referred to as the VE problem). These two tasks are coupled by coefficients related to surface boundary conditions of energy absorption and surface oscillation amplitudes (hereinafter referred to as dynamic and kinematic boundary conditions, respectively). The TBL problem is first solved for a rigid surface, thus providing necessary input to describe the forcing function on the surface, and also providing baseline calculations of friction drag, for comparison. The VE problem is solved next, given a periodic forcing function that approximates the shear and pressure pulsations of a given boundary layer. Initial choices for material parameters are based on theoretical and empirical guidelines. Optimal material parameters are chosen, following a series of iterations, such that the following two criteria are met:

 1) the energy flux into the viscoelastic coating is maximum, and
 2) the amplitudes of surface oscillation are less than the thickness of the viscous sublayer of the turbulent boundary layer over the coating.
If the amplitude of oscillations exceeds the thickness of the viscous sublayer of the turbulent flow, then the oscillations effectively increase the roughness of the surface, thus leading to an increase in friction drag. Furthermore, as the phase speed of disturbances in the boundary layer exceeds the shear wave speed in the material, a resonant interaction with large amplitude waves occurs. These conditions are to be avoided. Moderate energy flux into the material, however, where energy is transformed into internal shear waves and eventually dissipated as heat, leads to a qualitative and quantitative change in the turbulent energy balance, with a consequent reduction in friction drag. Optimal physical properties of a viscoelastic material will vary with freestream velocity, position along a body, pressure gradient, and any other factors which influence the development of the boundary layer and the characteristics of local turbulent fluctuations.
By solving the TBL equations, the turbulent friction drag over a viscoelastic, elastic, or rigid surface can be quantitatively evaluated. In the case of a viscoelastic surface, where energy is absorbed and surface oscillations are nonzero, both dynamic and kinematic boundary conditions are specified. These boundary conditions are derived directly from the solution of the VE equations for energy flux and surface oscillation amplitudes, and then transferred into a dissipation boundary condition and Reynolds stress boundary conditions for solution of the TBL equations. Vertical oscillations influence the effective roughness of the surface, and the rootmean square (rms) value of the vertical oscillation amplitude is classified as the dynamic roughness. If the oscillation amplitudes are lower than the viscous sublayer thickness, it is appropriate to estimate Reynolds stresses as zero. The equations for a turbulent boundary layer describe turbulent diffusion as a gradient approximation, which accommodates the dynamic boundary condition, and nearwall functions are introduced to describe for different surfaces the redistribution of turbulent energy in the nearwall region.
The present invention will become more fully understood from the detailed description given below and the accompanying drawings, which are given by way of illustration only and thus are not limitative of the present invention, wherein:
The present invention identifies physical and geometric parameters of a viscoelastic coating that reduces turbulent friction drag under given flow conditions. Furthermore, the invention permits evaluation of the anticipated drag reduction effectiveness of a given material with known physical properties for a given body configuration and set of flow conditions. The methodology has been applied principally to the characterization of coatings for turbulent flow over flat plates and bodies of revolution, and can also be applied to more complex geometries having curvature and nonzero pressure gradients.
A fluid boundary layer is the very thin layer of fluid adjacent to a surface over which fluid is flowing. It is the region where frictional forces play a major role, and is where the flow adjusts from conditions at the surface to conditions in the freestream of the flow. The outer edge of the boundary layer is traditionally defined as that location where the ratio, β, of the mean velocity, U, to the freestream velocity, U_{28 }, is a constant which is approximately equal to 1:
The value of the constant, β, chosen will depend upon configuration and numerical stability concerns. For the case of a flat plate, an appropriate value for this constant is 0.9975.
A turbulent boundary layer is characterized by a spectrum of pressure and shear fluctuations, the frequency, phase speed, and amplitude characteristics of which are a function of such factors as freestream velocity, body configuration, surface conditions, and pressure gradient. With flow over a rigid surface, there is no motion of the surface. With an elastic or viscoelastic surface, the wall pressure and shear fluctuations act as a forcing function which can deform the surface, creating surface waves. With a viscoelastic surface, energy from the turbulent boundary layer may be absorbed and dissipated by the coating, thus necessitating proper specification of boundary conditions for both Reynolds stresses and energy absorption at the wall (i.e., kinematic and dynamic boundary conditions).
where τ_{w }is the shear stress at the wall and ρ is the density of the fluid. The interaction of the turbulent flow with a viscoelastic coating leads to the formation of a quasiperiodic surface wave. The motion and energy absorption of the coating (kinematic and dynamic boundary conditions) in turn affect the energy balance in the turbulent boundary layer and the value of the friction drag, the latter of which is the surface integral of the wall shear stress.
The methodology for this invention is schematically shown in FIG. 2. The present invention includes solutions for: 1) turbulent boundary layer (TBL) parameters, including friction drag over rigid, viscoelastic, or elastic plates; and, 2) energy absorption and oscillation amplitudes of a viscoelastic (VE) plate excited by a periodic load which approximates that of a turbulent boundary layer. These two parts of the solution are coupled by boundary, conditions, parts of the solution are coupled by boundary conditions, both for energy absorbed by the surface and for the amplitudes of surface motion, and are part of an overall methodology for selecting dragreducing coatings and for quantifying drag reduction for given flow conditions. Below, the two parts of the solution, as well as the methodology that couples them, are described.
Characterization of a Turbulent Boundary Layer Over Rigid, Elastic, or Viscoelastic Surfaces (TBL Problem)
In
General System of Equations of Continuity, Motion and Energy: Turbulent flow parameters are obtained through the solution of a system of equations of continuity, motion, and energy, with accompanying boundary conditions. These equations are developed from principles of conservation of mass, conservation of momentum (Newton's second law), and energy balance (as based on the first law of thermodynamics).
In Cartesian coordinates, the general equation of continuity for a compressible fluid having a density of ρ, and velocity components U, V, and W in the streamwise, normal, and transverse directions, is given by Equation 3, below.
Alternatively, Equation 3 can be written in terms of indicial notation (as in Equation 4), where x, y, and z are represented by x_{1}, x_{2}, and x_{3}, respectively, an where U, V, and W are represented by U_{1}, U_{2}, and U_{3}, respectively. It is implied that the index, i, can have a value of 1, 2, or 3, and that a repeated index of i indicates summation.
Equation (3) is further simplified for a fluid that is incompressible, i.e., where the density of the fluid is constant, the following applies:
The generalized equations of motion, termed the NavierStokes equations, are expressed in Cartesian coordinates for the case of an incompressible fluid with constant viscosity as:
In Equations 6:

 P is the mean pressure,
 g_{i }is the body force vector, due to external fields, such as gravity, which act on the element, and
 ν is the kinematic viscosity (assumed to be constant).
Note that the above expression, written in indicial notation, actually represents three equations, for the three components of velocity in the x, y, and z directions.
Turbulent velocity components may be described as the sum of the mean and fluctuating components, U_{i }and u′, respectively, where U_{1}, U_{2}, and U_{3}, are equivalent to U, V, and W and where u′_{1}, u′_{2}, and u′_{3 }are equivalent to u′, v′, and w′:
U_{i}=U_{i}+u_{i} (Equations 7a7b7c)
An overbar indicates timeaveraging:
U_{i}=Ū_{i} (Equations 8a8b8c)
Substituting Equations 7a7c into Equations 6a6c and timeaveraging yields the following system of three complex nonlinear secondorder partial differential equations of motion for turbulent flow:
In Equations 9a9c, the components
Closure of the generalized system of equations including the continuity equation (Equation 5) and the equations of motion (Equations 9a9c) for a turbulent flow requires seven additional equations to characterize the six Reynolds stresses,
There exist multiple approaches within the literature for developing additional equations for Reynolds stress terms in turbulent flow, but this invention adopts a Reynoldsstresstransporttype methodology. In this methodology, equations for Reynolds stresses take the following general form:
where P_{ij }is termed the production, Π_{ij }is termed the pressurestrain correlation tensor, J_{ijk }is termed the diffusive flux of the Reynolds stresses, and ε_{ij }is termed the dissipation tensor.
In the general case, equations for all six Reynolds stress terms, and for the energy dissipation rate must expressed. The equation for the isotropic dissipation rate, ε, is similar in structure to the equations for the transport of Reynolds stresses. Full mathematical expressions for the Reynolds stress and isotropic dissipation rate equations shall be expressed in the following section for the specific case of a twodimensional turbulent boundary layer.
In summary, the equations which are solved to determine turbulent boundary layer parameters include:

 the continuity equation, Equation 5,
 the three equations of motion, Equations 9a, 9b, 9c,
 six equations for Reynolds shear and normal stresses, Equations 12a12f, and
 the equation for isotropic dissipation rate, ε (Equation 16 below).
The methodology for the solution of turbulent flow parameters involves a finite difference approximation of the system of equations of motion and continuity, with accompanying boundary conditions.
Turbulent Boundary Layer Equations: Complete mathematical formulations are provided for the specific case of a turbulent boundary layer with a steady, twodimensional mean flow and a constant freestream velocity, U_{∞}. Twodimensional turbulent boundary layer equations, as termed in the literature, are derived from the general continuity equation (Equation 5) and equations of motion (Equations 9a9c), given the assumptions that:

 The mean transverse velocity, W, is zero.
 Gravitational forces can be neglected.
 The pressure gradient is approximately zero in the y direction.
 The mean velocity in the streamwise direction, U, is much greater than the mean velocity in the normal direction, V.
 The rate of change of parameters in the x direction is much smaller than the rate of change of parameters in the y direction.
The above assumptions permit simplification of the set of equations required to solve for turbulent boundary layer parameters to include the revised continuity equation, Equation (13):$\begin{array}{cc}\frac{\partial U}{\partial x}+\frac{\partial V}{\partial y}=0& \left(\mathrm{Equation}\text{\hspace{1em}}13\right)\end{array}$
and the equation of motion for the U velocity component (Equation 14):$\begin{array}{cc}U\frac{\partial U}{\partial x}+V\frac{\partial U}{\partial y}=\frac{1}{\rho}\frac{\partial \mathrm{??}}{\partial x}+v\frac{{\partial}^{2}U}{\partial {y}^{2}}\frac{\partial \stackrel{\_}{{u}^{\prime}{v}^{\prime}}}{\partial y}& \left(\mathrm{Equation}\text{\hspace{1em}}14\right)\end{array}$
Where transport equations for the six Reynolds stress components are required in the general case, the Reynolds shear stress components −
where P_{ij }is the production term, Π_{ij }is the pressurestrain correlation tensor, J_{ijk }is the diffusive flux of the Reynolds stresses, and ε_{ij }is the dissipation tensor. A fifth equation for ε is:
where the expression for viscous diffusion may alternatively be approximated as:
if required for numerical stability in solutions of viscoelastic, nonoscillating surfaces with limited grid points in the nearwall region.
In Equations (15a15d), the term P_{ij }may be expressed as:
In Equation (16), the term P_{Σ} may be expressed as:
The pressurestrain correlation tensor, Π_{ij}, which redistributes energy between different components of Reynolds stresses, may be expressed as:
where the π′_{ij,1 }terms represent nearwall redistribution of turbulent energy from the streamwise component to the normal and transverse components, the π′_{ij,2 }terms represent nearwall variation of the Reynolds stress tensor component production, and the π′_{ij,3 }terms represent nearwall redistribution of turbulent energy proportional to local vorticity:
Here
is a unique damping function for the nearwall region:
where:
and
Here, D_{ij }is a dissipation tensor:
is the gradient of turbulent and viscous diffusive flux of the Reynolds stresses in the boundary layer, where only one component remains in the boundarylayer representation:
where A is 6 in the equation for
except for Equation (16), where:
C_{i}=C_{ε} (Equation 30)
The dissipation tensor, ε_{ij}, is written as:
where ƒ_{s }characterizes flow in the nearwall region:
Equation (16) includes two functions, ƒ_{1 }and ƒ_{2}, which also introduce corrections for nearwall flows:
ƒ_{1}=1+0.8e^{−R} ^{ t } (Equation 34)
ƒ_{2}=1−0.2e^{−R} ^{ t } ^{2} (Equation 35)
Values of constants for flow over a flat plate are as shown in Table 1:
C_{1}  C_{2}  C_{ε1}  C_{a2}  C_{t}  C_{ε}  C_{1}′  C_{2}′  C_{3}′ 
1.34  0.8  1.45  1.9  0.12  0.15  0.36  0.45  0.036 
For the case of a twodimensional boundary layer, turbulent boundary layer parameters at different x and y locations are determined through solution of the continuity equation (Equation 13), the equation of motion in the xdirection (Equation 14), the transport equations for the Reynolds stresses
Boundary Conditions: Boundary conditions are values of parameters at the limits of the boundary layer, i.e., at the surface and the freestream. The freestream velocity is defined as U_{∞}. Boundary conditions at the surface are specified for Reynolds normal and shear stresses (kinematic boundary conditions), as well as for the isotropic dissipation rate (dynamic boundary condition). For an arbitrary geometry, the x and y coordinates of the surface must be specified. If the surface is a flat plate, the boundary will be along the line y=0.
Since oscillation amplitudes at the surface are small, linearized kinematic boundary conditions, where mean velocities at the surface are assumed to be zero, are appropriate. Boundary conditions for fluctuating velocity components at the surface of a flat plate are expressed as:
where ξ_{1 }and ξ_{2 }are the longitudinal and vertical surface displacement components, respectively, u* is the friction velocity (as previously defined), and Θ is the angle of the longitudinal axis relative to the mean flow in the x_{1}x_{3 }plane. With linearized boundary conditions, mean velocities at the wall are assumed to be zero. Surface displacements are approximated by the first mode of a Fourier series:
Here, α_{3 }is the wavenumber corresponding to the maximum turbulent energy in the boundary layer, and is given by:
where the energycarrying frequency, ω_{e}, is assumed to be:
and the phase speed corresponding to energycarrying disturbances in the boundary layer is assumed to be:
C≈0.8U_{∞} (Equation 42)
Since there is a range of frequencies which carry energy, as reported within the scientific literature, it is advantageous to also perform calculations for the case where:
In the absence of resonance, it is appropriate to timeaverage components of the Reynolds stress at the wall:
where ξ_{i} is the rms amplitude of the displacement. For a passive isotropic viscoelastic coating excited by a forced load, the response takes the form of a traveling wave, so that the phase shift between normal and longitudinal displacements, φ_{2}φ_{1}, will be approximately π/2, and the displacements, φ_{2}φ_{1} , will be approximately π/2, and the Reynolds shear stresses at the surface will be approximately zero. For anisotropic materials, the phase shift can be different, so that negative Reynolds shear stresses can be generated at the wall. For a rigid wall, there will be no motion at the wall, so that Reynolds shear and normal stresses shall be equal to zero.
The boundary condition for the isotropic dissipation rate is:
where the first term reflects viscous dissipation and the second reflects absorption of energy by the viscoelastic material. For a rigid surface, there is no energy absorption at the wall, so that the second term equals zero. The absorption of turbulent energy by the coating is equivalent to −
which is the diffusive flux of energy across the boundary, characterized using a gradient mechanism for turbulent diffusion. This expression of the dynamic boundary condition is compatible with the Reynolds stress transport methodology of turbulence closure.
Equations 13, 14, 15a15d, and 16 are solved for mean velocity components, Reynolds normal and shear stresses, and energy dissipation, given the kinematic and dynamic boundary conditions (Equations (44) through (48)) based on the solution of the viscoelasticity problem (as described in the following section). The problem is solved numerically, using finite difference approximations of the parabolic equations. Friction drag for a body with a viscoelastic coating is calculated as the integral of wall shear stresses, τ_{w}, over the surface of the body, where:
for a twodimensional body, and where μ=ρν is the dynamic viscosity. Comparison of results with those calculated for a rigid body of identical geometry under identical flow conditions leads to an estimation of anticipated friction drag reduction.
To reduce friction drag, it is necessary to minimize surface oscillation amplitudes, while maximizing the flux of turbulent energy from the flow into the coating, −
then the normal Reynolds stresses at the boundary (Equations (44)(46)) may be approximated as zero. For a coating which absorbs energy, with low levels of oscillation, shear stresses in the nearwall region of the boundary layer decrease, as does the production of turbulence in the boundary layer. For a coating that oscillates at amplitudes greater than that of the viscous sublayer, the surface can act as a dynamic roughness element and thereby enhance the level of turbulence generated within the boundary layer.
Response of a Viscoelastic Material to a Turbulent Boundary Layer (VE Problem)
The second part of the methodology determines the response of a viscoelastic material to a turbulent boundary layer (step 2 in FIG. 2). For a rigid surface, Reynolds stresses on the surface (Equations (44) to (47)) are zero, and the isotropic dissipation rate contains only the viscous term. However, for a viscoelastic material, the kinematic and dynamic boundary conditions are determined through solution of the twodimensional conservation of momentum equation for a viscoelastic material:
where ρ_{s }is the material density, ξ_{2 }and ξ_{2 }are the longitudinal where ρ_{s } is the material density, ξ _{1 } and ξ _{2 } are the longitudinal and normal displacements through the thickness of the coating, and σ_{ij }is the amplitude of the stress tensor. The stress tensor for a viscoelastic material is written for a KelvinVoigt type of material as:
σ_{ij}=λ(ω)ε^{s}δ_{ij}+2μ(ω)ε_{ij} ^{s} (Equations 52a52d)
where ε_{ij} ^{s }is the strain tensor:
and:
ε^{s}=ε_{ii} ^{s} (Equation 54)
λ(ω) is the frequencydependent Lame constant, which is defined in terms of the bulk modulus, K(ω), which can be reasonably approximated as the static bulk modulus, K_{0}, and the complex shear modulus, μ(ω):
Displacements, ξ_{i}, are approximated as periodic, in the form of Equation (39), and can be expressed as a function of potentials of longitudinal and transverse (shear) waves:
where ∇_{φ} is the gradient of φ and ∇×{right arrow over (Ψ)} is the curl of the vector {right arrow over (Ψ)}. Equation (54) can be rewritten as two decoupled equations for the two wave potentials:
Equations (57) and (58) can be solved for the potentials, φ and Ψ, and hence for displacements, velocities, and stresses through the thickness of the coating, if boundary conditions are specified. The coating is fixed at its base, so that the longitudinal and normal displacements are zero, and the shear stress and pressure load on the surface is known. Pressure and shear pulsations on the coatings are approximated as periodic functions, with a form similar to that of the displacements in Equation (39), but with the following magnitudes, respectively:
τ_{ω}=ρu*^{2} (Equation 59)
ρ_{rms}=K_{p}τ_{ω}K_{p}ρu*^{2} (Equation 60)
I_{ω} =ρμ _{φ} ^{2} (Equation 59 )
ρ_{rms} =K _{p} I _{(t)} =K _{p} ρμ* ^{2} (Equation 60 )
where K_{p }is the Kraichnan parameter, whose value is approximated as 2.5.
If shear pulsations are included, a phase shift between shear and pressure pulsations must also be introduced.
If calculations are performed for a unit load, then surface displacements under actual load will be:
Kinematic boundary conditions in Equations (44) to (47) for the turbulent boundary layer problem are rewritten in terms of output from the materials problem:
Dynamic boundary conditions are rewritten in the form:
where:
C_{k3}=C_{k2}K_{p}γ(ω)
and where γ(ω) is a dissipative function of the coating material.
The flux of turbulent fluctuating energy through the surface can be solved for directly, as:
−
but the nondimensionalized flux can also be approximated as a diffusive flux term, using the gradient diffusion approach:
where:
Equation 68 provides a basis for determining the value of the kinematic coefficient of turbulence diffusion, ε_{q}, on an absorbing surface. Substituting Equation 68 into Equation 66 yields the following expression for ε_{q}, defined as ε_{t}_{y=0}: 66 yields the following expression for {tilde over (ε)}_{q} , defined as {tilde over (ε)} _{t}_{y=0}:
Here, K^{+} _{max }is the maximum of turbulence kinetic energy, and K^{+} _{q }is the kinetic energy of the oscillating surface, both quantities nondimensionalized by U_{∞} ^{2}. y_{max} ^{+} is defined as the normal distance from the surface to the maximum of turbulence energy, nondimensionalized as follows:
and:
where:
Thus we can determine the dissipation rate at the wall based on (Equation 48).
Methodology to Choose Properties of a DragReducing Viscoelastic Material
A methodology to choose properties of a viscoelastic coating that reduces turbulent friction drag necessarily requires both of the previously described solutions for turbulent boundary layer parameters and response of a viscoelastic material.
For the case of twodimensional flow over a flat plate, the TBL problem is solved for a rigid plate, in order to determine the boundary layer thickness, δ, at a given freestream velocity, U_{∞}, and location. The boundary layer thickness is determined from the finite difference solution of the seven equations of continuity, motion in the xdirection, transport equations for Reynolds normal and shear stresses, and the equation for dissipation rate, assuming no motion at the wall. The external limit of the boundary layer is defined as that location where the ratio of the mean velocity to the freestream velocity is a constant, β, between 0.95 and 1.0.
The frequencydependent, complex shear modulus of a material, μ(ω), can be expressed in different mathematical forms, some of which approximate experimentally measured shear modulus data more accurately than others. A single relaxation time (SRT) material is one where the complex shear modulus is expressed using a single relaxation time, τ_{s}, and a single value for the dynamic shear modulus, μ_{2}. In and a single value for the dynamic shear modulus μ_{2} . In Equation (76), an SRT material would be represented for the case of N=1. A multiple relaxation time (MRT) material is one where N>1 in the representation for complex shear modulus in Equation (78).
The HavriliakNegami (HN) representation for the complex shear modulus, is given by Equation 74. This equation is more complex, but often is more suitable for describing real materials:
For a HN type of material, whose complex shear modulus is expressed in the form of Equation (77), μ_{∞} is the limiting highfrequency modulus and α_{HN }and β_{HN }are constants. For the type of polymeric materials used for dragreducing coatings, K(ω) is essentially constant, with a value of approximately 1×10^{8 }Pa.
It is recommended to first determine an optimal SRT type of material, and then to choose an HN type of material, whose properties can be created with available polymer chemistry.
An SRT material can be adequately characterized by the material thickness, H, the density, ρ_{s}, the static shear modulus, μ_{0}, the dynamic shear modulus, μ_{s}, and the relaxation time, τ. An appropriate density for the viscoelastic material, ρ_{s}, is within 10% that of water. For an SRT material, the initial guess for the static shear modulus of the material, μ_{0}, is:
μ_{0}=ρ_{s} C ^{2} (Equation 78)
based on the criterion that the speed of shear waves in the material is approximately the same as the phase speed of the energycarrying disturbances, C. This phase speed, C, is assumed to be 0.8 of the value of the freestream velocity, U_{∞} (Equation 42). If the convective velocity exceeds the shear wave velocity, an instability occurs, and large waves appear on the surface of the material, leading to an increase of drag for the coating.
An initial choice for thickness, H, for isotropic viscous materials, where
is:
and for isotropic, low viscosity materials, where
is:
The optimal desired thickness for a coating may be greater than practical for a given application. While isotropic coatings thinner than recommended in Equations 7980 can still be effective, anisotropic coatings that are stiffer in the normal dimension relative to the transverse and longitudinal dimensions can provide equivalent performance with significant reduction in thickness.
Given specified values of H, μ_{0}, and ρ_{s}, the VE problem, as expressed in Equations (57) and (58), is solved numerically for a matrix of values of τ_{s }and μ_{s }(i.e., for different values of the complex shear modulus) and for a range of wavenumbers. The wavenumber corresponding to the maximum turbulent energy in the boundary layer is:
where the frequency, ω_{e}, for maximum energycarrying disturbances is estimated by Equation (41). Calculations yield surface displacement amplitudes and the flux of turbulent fluctuating energy into the coating. The best combination of properties for a SRT material occurs where the surface displacement under actual load (Equation (61)) is less than the viscous sublayer thickness, and where the energy flux into the coating (Equation (57)) is at a maximum. Furthermore, it is desirable to maintain this criterion for a range of frequencies from approximately one decade below to one decade above the energycarrying frequency, ω_{e}.
Once a set of optimal values of τ_{s }and μ_{s }are determined for a set of specified values of H, μ_{0}, and ρ_{s}, the calculations are iterated using slightly different values of thickness, H, and static modulus, μ_{0}. From these calculations are chosen the optimal set of parameters for an SRT material (H, μ_{0}, ρ_{s}, τ_{s }and μ_{s}), given specified flow conditions and configuration.
The complex shear moduli of real polymeric materials, such as polyurethanes and silicones, which are candidates for viscoelastic coatings cannot be adequately described by the SRT representation. More complex MRT or HN representations of the shear modulus require multiple constants, and are less suitable for numerical parametric evaluation. Therefore, results for SRT materials are used to select candidate materials, such as described by the HN formulation (Equation (79)) which can be more readily fabricated in practice. As a guideline, it is desired to match the complex shear modulus curves of the target SRT material and the HN material (value and slope) over frequencies ranging from decade below to one decade above ω_{e}, with the most important matching being in the immediate vicinity of ω_{e}.
To design a multilayer isotropic coating, properties of a complex shear modulus, density, and thickness are specified for individual layers, and nonslip boundary conditions between layers are imposed. The properties of the upper layer are specified according to the methodology for a single layer, and the lower layers will have progressively lower static shear moduli, as optimized for lower freestream velocities. Thus, welldesigned multilayer coatings can reduce drag over a range of freestream velocities.
In the design of an anisotropic coating, the complex shear modulus has different values in the normal direction relative to the longitudinal and transverse directions (hereinafter termed transversely isotropic). If the viscoelastic material follows a singlerelaxation time model, then the static shear modulus, μ_{0}, the dynamic shear modulus, μ_{s}, and the relaxation time, τ_{s}, will differ with direction, as expressed in Equations (82) and (83). The static shear modulus in the normal direction, μ_{01}, will be greater than than in the longitudinaltransverse plane, μ_{02}. The complex shear modulus in the normal direction is expressed as:
while the shear modulus in the streamwise and transverse directions is expressed as:
For a viscoelastic, transversely isotropic material, surface oscillation amplitudes can be reduced relative to an isotropic material, while the level of energy flux into the material is increased. Thus, welldesigned anisotropic coatings will be significantly thinner than isotropic coatings associated with the same level of drag reduction.
Methodology to Choose Structure of a DragReducing Viscoelastic Coating
A further aspect of coating design is the choice of internal structure within the viscoelastic material. In practical applications of viscoelastic coatings, the coating will be finite in length, with leading, trailing, and side edges. The influence of the finite edges affects coating performance. Well posed edges can order and stabilize transverse and longitudinal vortical structures in the nearwall region of the flow and thereby delay the deformation of these vortical structures and enhance the stability of the flow. However, unstructured edges can accentuate the amplitude of oscillations of the viscoelastic material in this region. Local instabilities can degrade the performance of the coating, so that, even with a welldesigned material, the influence of the edges can lead to a drag increase. Hence, the coating is structured in the vicinity of finite edges. The thickness of the coating is decreased to minimize such oscillations, using techniques such as a rigid wedge underneath the coating, or other localized structure near an edge (FIG. 3). For large bodies, a continuous coating may be impractical or difficult to fabricate. An alternative design is a piecewise continuous coating, composed of finite segments of coating, where both the longitudinal and transverse edges of the coating system are organized to stabilize flow structures and to minimize adverse effects at the edges of each segment.
In addition to wellposed edges, viscoelastic coatings may be combined with surface structure to enhance the stabilization of longitudinal vortices along the length of the coating, and hence to increase the level of drag reduction through multiple physical mechanisms. Structure can include the placement of riblets on top of the viscoelastic coating, or the creation of socalled “inverse” riblets. In the latter case, a viscoelastic coating may be molded over ribs or ridges of rigid material, so that longitudinal riblet structures form when fluid flows over the viscoelastic surface.
The dimensions (scales) of the segments and the dimensions of the structures within the coating are selected as multiples of the transverse and longitudinal scales in the near wall turbulent flow. These scales vary with body speed, position along the body and when nonNewtonian additives, such as dilute aqueous solutions of highmolecular weight polymers, are present.
Claims (25)
C=0.8U_{∞}.
C_{k3}=C_{k2}K_{p}γ(ω)
−
ξ_{i}=a_{i1}e^{ia} ^{ e } ^{(x−Ct)}
σ_{ij}=λ(ω)ε^{s}δ_{ij}+2μ(ω)ε_{ij} ^{s}
μ_{s}=μ_{ω=∞}−μ_{0}
Priority Applications (2)
Application Number  Priority Date  Filing Date  Title 

US09546380 US6516652B1 (en)  19990608  20000410  Design of viscoelastic coatings to reduce turbulent friction drag 
US11054719 USRE41398E1 (en)  20000410  20050210  Design of viscoelastic coatings to reduce turbulent friction drag 
Applications Claiming Priority (1)
Application Number  Priority Date  Filing Date  Title 

US11054719 USRE41398E1 (en)  20000410  20050210  Design of viscoelastic coatings to reduce turbulent friction drag 
Related Parent Applications (1)
Application Number  Title  Priority Date  Filing Date  

US09546380 Reissue US6516652B1 (en)  19990608  20000410  Design of viscoelastic coatings to reduce turbulent friction drag 
Publications (1)
Publication Number  Publication Date 

USRE41398E1 true USRE41398E1 (en)  20100629 
Family
ID=24180174
Family Applications (2)
Application Number  Title  Priority Date  Filing Date 

US09546380 Active US6516652B1 (en)  19990608  20000410  Design of viscoelastic coatings to reduce turbulent friction drag 
US11054719 Active USRE41398E1 (en)  19990608  20050210  Design of viscoelastic coatings to reduce turbulent friction drag 
Family Applications Before (1)
Application Number  Title  Priority Date  Filing Date 

US09546380 Active US6516652B1 (en)  19990608  20000410  Design of viscoelastic coatings to reduce turbulent friction drag 
Country Status (5)
Country  Link 

US (2)  US6516652B1 (en) 
EP (1)  EP1272387A1 (en) 
JP (1)  JP2003530217A (en) 
RU (1)  RU2250175C2 (en) 
WO (1)  WO2001076934A1 (en) 
Cited By (1)
Publication number  Priority date  Publication date  Assignee  Title 

US20100100361A1 (en) *  20070228  20100422  Keio University  Numerical structureanalysis calculation system 
Families Citing this family (19)
Publication number  Priority date  Publication date  Assignee  Title 

US20030145980A1 (en) *  20020204  20030807  Sinha Sumon Kumar  System and method for using a flexible composite surface for pressuredrop free heat transfer enhancement and flow drag reduction 
US7206258B1 (en)  20050413  20070417  United States Of America As Represented By The Secretary Of The Navy  Dual response acoustical sensor system 
BE1016624A3 (en) *  20050602  20070306  Rompay Boudewijn Gabriel Van  Underwater surface for e.g. boat hull or propeller, includes macroscopic irregularities for improving hydrodynamics 
US8039055B2 (en) *  20060720  20111018  Cortana Corporation  Method to increase the efficiency of polymer drag reduction for marine and industrial applications 
DE102007016643A1 (en) *  20070405  20081009  Geislinger Gmbh  A frictional clamping connection and processes for their preparation 
JP5590442B2 (en) *  20100310  20140917  国立大学法人東京農工大学  Frictional resistance reduction device 
US8457939B2 (en)  20101230  20130604  Aerion Corporation  Generating inviscid and viscous fluidflow simulations over an aircraft surface using a fluidflow mesh 
US8437990B2 (en)  20110311  20130507  Aerion Corporation  Generating a simulated fluid flow over an aircraft surface using anisotropic diffusion 
US8538738B2 (en)  20110322  20130917  Aerion Corporation  Predicting transition from laminar to turbulent flow over a surface 
US8892408B2 (en)  20110323  20141118  Aerion Corporation  Generating inviscid and viscous fluid flow simulations over a surface using a quasisimultaneous technique 
US20130146217A1 (en) *  20111209  20130613  Nicholas Joseph Kray  Method of Applying Surface Riblets to an Aerodynamic Surface 
JP5916490B2 (en) *  20120409  20160511  中国塗料株式会社  Frictional resistance prediction method and ship bottom coating film, paint film performance evaluation method using the method, the coating film performance evaluation device 
CA2896758A1 (en)  20121231  20140703  Telvent Dtn Llc  Dynamic aircraft threat controller manager apparatuses, methods and systems 
EP2938538A4 (en)  20121231  20170111  Telvent Dtn Llc  Dynamic turbulence engine controller apparatuses, methods and systems 
US9845162B2 (en) *  20130503  20171219  The Boeing Company  Protective finish for wing tip devices 
CN103512844A (en) *  20131009  20140115  哈尔滨工程大学  Nonsmooth surface fluid friction resistance testing device and nonsmooth surface antidrag effect evaluating method 
WO2015095890A1 (en) *  20131222  20150625  Telvent Dtn Llc  Dynamic storm environment engine apparatuses, methods and systems 
JP2017527752A (en) *  20140728  20170921  ザ リージェンツ オブ ザ ユニバーシティー オブ コロラド，ア ボディー コーポレート  Phononic materials used to control the flow behavior 
RU178601U1 (en) *  20170914  20180411  Общество с ограниченной ответственностью "МЕТАЛЛ"  Compliant coating for the hull plating 
Citations (16)
Publication number  Priority date  Publication date  Assignee  Title 

US3435796A (en) *  19671113  19690401  Us Navy  Method and apparatus for drag reduction 
US3516376A (en) *  19680816  19700623  Tadeusz Kowalski  Structure for reducing the drag between a fluid and a solid body 
US4771799A (en) *  19871029  19880920  Conoco Inc.  Method for improving the performance of highly viscous concentrates of high molecular weight drag reducing polymers 
US4865271A (en) *  19870319  19890912  RollsRoyce Plc  Boundary layer devices 
US4932612A (en) *  19860225  19900612  Blackwelder Ron F  Method and apparatus for reducing turbulent skin friction 
US5020561A (en)  19900813  19910604  Atlantic Richfield Company  Drag reduction method for gas pipelines 
US5133519A (en)  19890421  19920728  Board Of Trustees Operating Michigan State University  Drag reduction method and surface 
US5342465A (en) *  19881209  19940830  Trw Inc.  Viscoelastic damping structures and related manufacturing method 
US5619433A (en)  19910917  19970408  General Physics International Engineering Simulation Inc.  Realtime analysis of power plant thermohydraulic phenomena 
US6024119A (en) *  19980420  20000215  The United States Of America As Represented By The Secretary Of The Navy  Flow control system having actuated elastomeric membrane 
US6287664B1 (en) *  19971114  20010911  William F. Pratt  Continuous wave composite viscoelastic elements and structures 
US6332593B1 (en) *  20000216  20011225  Brown University Research Foundation  Method and apparatus for reducing turbulent drag 
US6357374B1 (en) *  20000721  20020319  Cortana Corporation  Method and apparatus for increasing the effectiveness and efficiency of multiple boundary layer control techniques 
US20020173569A1 (en) *  20000515  20021121  Esa Karhu  Drag reducing composition 
US20050032029A1 (en) *  19991110  20050210  Trunk Frank J.  Method of multidimensional analysis of viscoelastic materials for stress, strain, and deformation 
US6892989B1 (en) *  20030529  20050517  The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration  Method for reducing the drag of bluntbased vehicles by adaptively increasing forebody roughness 
Family Cites Families (6)
Publication number  Priority date  Publication date  Assignee  Title 

DE669897C (en)  19390106  Versuchsanstalt Fuer Luftfahrt  Means for reducing the frictional resistance  
US3161385A (en)  19600615  19641215  Coleman Kramer Inc  Means and method for stabilizing laminar boundary layer flow 
US3585953A (en)  19690620  19710622  Max Otto Kramer  Means and method for stabilizing laminar boundary layer flow 
US5359574A (en) *  19930827  19941025  The United States Of America As Represented By The Secretary Of The Navy  Electromagnetically activated compliant wavywall 
US5901928A (en) *  19960614  19990511  Aptek, Inc.  Active turbulence control technique for drag reduction 
US5890681A (en) *  19970501  19990406  The United States Of America As Represented By The Secretary Of The Navy  Method for controlling microturbulence 
Patent Citations (16)
Publication number  Priority date  Publication date  Assignee  Title 

US3435796A (en) *  19671113  19690401  Us Navy  Method and apparatus for drag reduction 
US3516376A (en) *  19680816  19700623  Tadeusz Kowalski  Structure for reducing the drag between a fluid and a solid body 
US4932612A (en) *  19860225  19900612  Blackwelder Ron F  Method and apparatus for reducing turbulent skin friction 
US4865271A (en) *  19870319  19890912  RollsRoyce Plc  Boundary layer devices 
US4771799A (en) *  19871029  19880920  Conoco Inc.  Method for improving the performance of highly viscous concentrates of high molecular weight drag reducing polymers 
US5342465A (en) *  19881209  19940830  Trw Inc.  Viscoelastic damping structures and related manufacturing method 
US5133519A (en)  19890421  19920728  Board Of Trustees Operating Michigan State University  Drag reduction method and surface 
US5020561A (en)  19900813  19910604  Atlantic Richfield Company  Drag reduction method for gas pipelines 
US5619433A (en)  19910917  19970408  General Physics International Engineering Simulation Inc.  Realtime analysis of power plant thermohydraulic phenomena 
US6287664B1 (en) *  19971114  20010911  William F. Pratt  Continuous wave composite viscoelastic elements and structures 
US6024119A (en) *  19980420  20000215  The United States Of America As Represented By The Secretary Of The Navy  Flow control system having actuated elastomeric membrane 
US20050032029A1 (en) *  19991110  20050210  Trunk Frank J.  Method of multidimensional analysis of viscoelastic materials for stress, strain, and deformation 
US6332593B1 (en) *  20000216  20011225  Brown University Research Foundation  Method and apparatus for reducing turbulent drag 
US20020173569A1 (en) *  20000515  20021121  Esa Karhu  Drag reducing composition 
US6357374B1 (en) *  20000721  20020319  Cortana Corporation  Method and apparatus for increasing the effectiveness and efficiency of multiple boundary layer control techniques 
US6892989B1 (en) *  20030529  20050517  The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration  Method for reducing the drag of bluntbased vehicles by adaptively increasing forebody roughness 
NonPatent Citations (3)
Title 

A.D. Young, Boundary Layers, 1989, p. 12, AIAA Education Series, Washington, DC. 
P.A. Davidson, TurbulenceAn Introduction for Scientists and Engineers, 2004, pp. 130 and 131, Oxford University Press, New York. 
P.A. Davidson, Turbulence—An Introduction for Scientists and Engineers, 2004, pp. 130 and 131, Oxford University Press, New York. 
Cited By (2)
Publication number  Priority date  Publication date  Assignee  Title 

US20100100361A1 (en) *  20070228  20100422  Keio University  Numerical structureanalysis calculation system 
US8352219B2 (en) *  20070228  20130108  Keio University  Numerical structureanalysis calculation system 
Also Published As
Publication number  Publication date  Type 

US6516652B1 (en)  20030211  grant 
JP2003530217A (en)  20031014  application 
RU2250175C2 (en)  20050420  grant 
WO2001076934A1 (en)  20011018  application 
EP1272387A1 (en)  20030108  application 
RU2002129899A (en)  20040320  application 
Similar Documents
Publication  Publication Date  Title 

Luo et al.  A numerical simulation of unsteady flow in a twodimensional collapsible channel  
Cherukat et al.  Direct numerical simulation of a fully developed turbulent flow over a wavy wall  
Fung  An introduction to the theory of aeroelasticity  
Goldstein et al.  Modeling a noslip flow boundary with an external force field  
Goldstein  Scattering of acoustic waves into TollmienSchlichting waves by small streamwise variations in surface geometry  
Ahn et al.  Strongly coupled flow/structure interactions with a geometrically conservative ALE scheme on general hybrid meshes  
Liao et al.  Simulating flows with moving rigid boundary using immersedboundary method  
Lam et al.  On the condition of streak formation in a bounded turbulent flow  
Baarholm et al.  On fatigue damage accumulation from inline and crossflow vortexinduced vibrations on risers  
Celebi et al.  Nonlinear modeling of liquid sloshing in a moving rectangular tank  
Bourgault et al.  A finite element method study of Eulerian droplets impingement models  
Xue et al.  Constitutive model for quasi‐static deformation of metallic sandwich cores  
Michelin et al.  Resonance and propulsion performance of a heaving flexible wing  
Arnal et al.  Laminarturbulent transition prediction in threedimensional flows  
Wilson et al.  Unsteady RANS method for ship motions with application to roll for a surface combatant  
Sutherland  Finiteamplitude internal wavepacket dispersion and breaking  
Azcueta  Computation of turbulent freesurface flows around ships and floating bodies  
Monaghan et al.  SPH particle boundary forces for arbitrary boundaries  
Saric  Görtler vortices  
Zhu et al.  Mode coupling and flow energy harvesting by a flapping foil  
Aquelet et al.  Euler–Lagrange coupling with damping effects: Application to slamming problems  
Tsukahara et al.  DNS of turbulent Couette flow with emphasis on the largescale structure in the core region  
Sytine et al.  Convergence tests for the piecewise parabolic method and Navier–Stokes solutions for homogeneous compressible turbulence  
Staquet  Twodimensional secondary instabilities in a strongly stratified shear layer  
Librescu et al.  Implications of cubic physical/aerodynamic nonlinearities on the character of the flutter instability boundary 
Legal Events
Date  Code  Title  Description 

CC  Certificate of correction  
FPAY  Fee payment 
Year of fee payment: 12 