Engineering PapersSearch

Engineering topics

Durbin, P. A.

Publications and source records attributed to Durbin, P. A..

At least 19 records

Non-Equilibrium Turbulence Modeling for High Lift Aerodynamics

This phase is discussed in ('Non linear kappa - epsilon - upsilon(sup 2) modeling with application to high lift', Application of the kappa - epsilon -upsilon(sup 2) model to multi-component airfoils'). Further results are presented in 'Non-linear upsilon(sup 2) - f modeling with application to high-lift' The ADI solution method in the initial implementation was very slow to converge on multi-zone chimera meshes. I modified the INS implementation to use GMRES. This provided improved convergence and less need for user intervention in the solution process. There were some difficulties with implementation into the NASA compressible codes, due to their use of approximate factorization. The Helmholtz equation for f is not an evolution equation, so it is not of the form assumed by the approximate factorization method. Although The Kalitzin implementation involved a new solution algorithm ('An implementation of the upsilon(sup 2) - f model with application to transonic flows'). The algorithm involves introducing a relaxation term in the f-equation so that it can be factored. The factorization can be into a plane and a line, with GMRES used in the plane. The NASA code already evaluated coefficients in planes, so no additional memory is required except that associated the the GMRES algorithm. So the scope of this project has expanded via these interactions. . The high-lift work has dovetailed into turbine applications.

Durbin, P. A.

Non-linear k-epsilon-v(sup 2)(bar) modeling with application to high-lift

The k-epsilon-v(sup 2)(bar) model has been investigated to quantify its predictive performance on two high-lift configurations: 2D flow over a single-element aerofoil, involving closed-type separation; 3D flow over a prolate spheroid, involving open-type separation. A 'code-friendly' modification has been proposed which enhances the numerical stability, in particular, for explicit and uncoupled flow solvers. As a result of introducing Reynolds-number dependence into a coefficient of the s-equation, the skin-friction distribution for the by-pass transitional flow over a flat plate is better predicted. In order to improve deficiencies arising from the Boussinesq approximation, a nonlinear stress-strain constitutive relation was adopted, in which the only one free constant is calibrated on the basis of DNS data, and the Reynolds-stress anisotropy near the wall is fairly well represented.

Lien, F. S.

Application of the k-epsilon-v(exp 2) model to multi-component airfoils

Flow computations around two-element and three-element configurations are presented and compared to detailed experimental measurements. The k-epsilon-v(exp 2)(bar) model has been applied and the ability of the model to capture streamline curvature effects, wake-boundary layer confluence, and laminar/turbulent transition is discussed. The numerical results are compared to experimental datasets that include mean quantities (velocity and pressure coefficient) and turbulent quantities (Reynolds normal and shear stresses).

Iaccarino, G.

Prediction of the backflow and recovery regions in the backward facing step at various Reynolds numbers

A four-equation model of turbulence is applied to the numerical simulation of flows with massive separation induced by a sudden expansion. The model constants are a function of the flow parameters, and two different formulations for these functions are tested. The results are compared with experimental data for a high Reynolds-number case and with experimental and DNS data for a low Reynolds-number case. The computations prove that the recovery region downstream of the massive separation is properly modeled only for the high Re case. The problems in this case stem from the gradient diffusion hypothesis, which underestimates the turbulent diffusion.

Michelassi, V.

Turbulence modeling for non-equilibrium flow

The work performed during this year has involved further assessment and extension of the k-epsilon-v(exp 2) model, and initiation of work on scalar transport. The latter is introduced by the contribution of Y. Shabany to this volume. Flexible, computationally tractable models are needed for engineering CFD. As computational technology has progressed, the ability and need to use elaborate turbulence closure models has increased. The objective of our work is to explore and develop new analytical frameworks that might extend the applicability of the modeling techniques. In past years the development of a method for near-wall modeling was described. The method has been implemented into a CFD code and its viability has been demonstrated by various test cases. Further tests are reported herein. Non-equilibrium near-wall models are needed for some heat transfer applications. Scalar transport seems generally to be more sensitive to non-equilibrium effects than is momentum transport. For some applications turbulence anisotropy plays a role and an estimate of the full Reynolds stress tensor is needed. We have begun work on scalar transport per se, but in this brief I will only report on an extension of the k-epsilon-v(exp 2) model to predict the Reynolds stress tensor.

Durbin, P. A.

A Reynolds stress model for near-wall turbulence

The paper formulates a tensorially consistent near-wall second-order closure model. Redistributive terms in the Reynolds stress equations are modeled by an elliptic relaxation equation in order to represent strongly nonhomogeneous effects produced by the presence of walls; this replaces the quasi-homogeneous algebraic models that are usually employed, and avoids the need for ad hoc damping functions. The model is solved for channel flow and boundary layers with zero and adverse pressure gradients. Good predictions of Reynolds stress components, mean flow, skin friction, and displacement thickness are obtained in various comparisons to experimental and direct numerical simulation data. The model is also applied to a boundary layer flowing along a wall with a 90-deg, constant-radius, convex bend.

Durbin, P. A.

A transport equation for eddy viscosity

A transport equation for eddy viscosity is proposed for wall bounded turbulent flows. The proposed model reduces to a quasi-homogeneous form far from surfaces. Near to a surface, the nonhomogeneous effect of the wall is modeled by an elliptic relaxation model. All the model terms are expressed in local variables and are coordinate independent; the model is intended to be used in complex flows. Turbulent channel flow and turbulent boundary layer flows with/without pressure gradient are calculated using the present model. Comparisons between model calculations and direct numerical simulation or experimental data show good agreement.

Durbin, P. A.

Local anisotropy in strained turbulence at high Reynolds numbers

It is shown that the hypothesis of local isotropy is implausible in the presence of significant mean rates of strain. In fact, it appears that in uniform shear flow near equilibrium, local isotropy can never constitute a systematic approximation, even in the limit of infinite Reynolds number. An estimate of the level of mean strain rate for which local isotropy is formally a good approximation is provided.

Durbin, P. A.

Considerations on the moving contact-line singularity, with application to frictional drag on a slender drop

It has previously been shown that the no-slip boundary conditions leads to a singularity at a moving contact line and that this presumes some form of slip. Present considerations on the energetics of slip due to shear stress lead to a yield stress boundary condition. A model for the distortion of the liquid state near solid boundaries gives a physical basis for this boundary condition. The yield stress condition is illustrated by an analysis of a slender drop rolling down an incline. That analysis provides a formula for the frictional drag resisting the drop movement. With the present boundary condition, the length of the slip region becomes a property of the fluid flow.

Durbin, P. A.

Observations of the frequencies in a sphere wake and of drag increase by acoustic excitation

Vortex shedding and instability wave frequencies have been measured in the wakes of spheres in the Reynolds number range 500 less than Re less than 60,000. The effect of acoustic excitation was examined and an interaction between the two frequency modes was found at the lower Reynolds numbers; through this interaction, external forcing at the instability frequency could change the vortex shedding frequency. The development of the mean wake was manipulated by forcing near to the dominant shear layer instability frequency. With this forcing, the separated shear layer moved closer to the surface of the sphere and the reversed flow region of the wake was shortened. Concomitantly, the base pressure decreased and drag increased.

Kim, H. J.

Investigation of the flow between a pair of circular cylinders in the flopping regime

The wakes of a pair of circular cylinders are grossly unsteady when the cylinders are separated in a direction normal to the approaching flow by less than one cylinder diameter. The wakes flop randomly between two asymmetric states. The time-scale for the flopping is several orders of magnitude longer than the timescale of vortex shedding, and also several orders of magnitude longer than the timescale for instability of the separating shear layers. When a splitter plate is positioned suitably on the centerline of the cylinders, the flopping can be stopped and the flow made to assume either of the asymmetric states, or a symmetric steady state. For a range of plate positions a new, periodic oscillation occurs. Acoustic excitation can also destroy the flopping mean flow, replacing it by a symmetric flow.

Kim, H. J.

Free-streamline analysis of deformation and dislodging by wind force of drops on a surface

Free-streamline theory is used to analyze the deformation and dislodging by wind pressure of drops of liquid adhered by surface tension to a solid surface. The critical Weber number for droplets to be dislodged is determined as a function of advancing and receding contact angle. Graphical results for drop shape are in good agreement with observation.

Durbin, P. A.

Roll-up of vorticity in adverse-pressure-gradient boundary layers

It is shown how the unsteady, nonlinear critical-layer equation determines the evolution of instability waves in a weak adverse-pressure-gradient boundary layer. Numerical solutions show that the nonlinearity halts the growth of these inviscidly unstable waves. The stabilizing effect of nonlinearity, in the present case, can be described as a consequence of either the increase (toward zero) of the phase jump across the critical layer or the roll-up of the critical-layer disturbance vorticity.

Goldstein, M. E.

Shear flow instability generated by non-homogeneous external forcing

An experiment has been designed and conducted in order to ascertain whether instability waves can be generated by nonhomogeneous forcing, using a biconvex vane located outside the mixing layer whose oscillation was induced by an electromagnetic shaker through a linkage. The vane was oscillated at 20 Hz, and the resulting spectra were computed by a spectrum analyzer. The data are judged to provide an example of instability waves generated solely through nonhomogeneous forcing.

Durbin, P. A.

An anemometer for highly turbulent or recirculating flows

An anemometer which determines flow velocity by ionizing air and sensing the convective displacement of the ions is described. It is suited to measurement in low speed, highly unsteady gas flows. Comparisons to hot wire spectra suggest the corona anemometer has adequate frequency response to make it a useful tool for fluid dynamics measurement.

Durbin, P. A.

Corona anemometry for qualitative measurement of reversing surface flow with application to separation control by external excitation

An corona anemometer which detects gas flow by the displacement of an ion beam is described, and experiments are performed using the anemometer to investigate the active control of diffusor separation by periodic forcing. The apparatus is applied to the separated flow over a rearward facing ramp. An oscillating vane is attached to the surface near the separation point. It is suggested that the enhancement in turbulent energy produced by the oscillating vane is due to drastic modification of the wake shear flow, and not to vane-produced turbulence.

Durbin, P. A.

Asymptotic analysis of corona discharge from thin electrodes

The steady discharge of a high-voltage corona is analyzed as a singular perturbation problem. The small parameter is the ratio of the length of the ionization region to the total gap length. By this method, current versus voltage characteristics can be calculated analytically.

Durbin, P. A.