Engineering Papers⌕ Search

Engineering topics

Liu, C. H.

Publications and source records attributed to Liu, C. H..

At least 91 records · Page 5

Numerical studies of motion and decay of vortex filaments

A computational code is developed for the integro-differential equations governing the motion of the centerlines of vortex filaments submerged in a background potential flow. These equations, which are derived from the method of matched asymptotic analysis, include the effect of decaying large-magnitude circumferential and axial velocity components in the vortical cores. Numerical examples are presented to assess the effect of large axial velocity and of nonsimilar initial profiles in vortical cores. The initial configurations of the filaments are chosen so as to fulfill the basic assumption of asymptotic analysis, which is the effective vortical core size is much smaller than all other length scales in the flowfield, e.g., the radius of curvature and interfilament distance. The computations are continued until the basic assumption is no longer valid, that is, when the merging or intersection of filaments have begun. Various types of local or global merging or intersection of filaments are classified and demonstrated by numerical examples.

Liu, C. H.↗

Implicit, vectorizable schemes for the flux-difference split, three-dimensional Navier-Stokes equations

Two hybrid upwind models are defined for solving the Euler equations. The algorithms both employ approximate factorization (AF) in crossplane and symmetric block Gauss-Seidel relaxation in the third direction. One approach adds an additional factorization step to lower the number of required grid point operations for inversion of the block tridiagonal matrices; however, the move permits only one third of the operations to be vectorized. Finite difference solutions are calculated on a C-H-type grid, in this case enveloping a slender, sharp-edged delta wing. Sample data are provided for the calculated vortex flow for Re of 10,000, at a 20.5 deg angle of attack, represented in a crossflow velocity vector plot and in a spanwise pressure coefficient distribution. The AF scheme, without additional factorization, when used with a grid covering 51 x 51 x 72 points provides a convergent solution with no time step lasting longer than 0.00001 sec.

Liu, C. H.↗

Initiation of breakdown in slender compressible vortices

The onset of vortex breakdown in compressible flows is investigated analytically for the case in which the flow is axially symmetric, the vortex is isolated, its axis is parallel to the main flow, and the vortex radius is small compared to the breakdown length. The conservation equations for mass, momentum, and energy are formulated and solved numerically using a finite-difference scheme, as described by Krause (1985); numerical results are presented in graphs and briefly characterized.

Krause, E.↗

Possibility of measuring gravity-wave momentum flux by single beam observation of MST radar

Vincent and Reid (1983) proposed a technique to measure gravity-wave momentum fluxes in the atmosphere by mesosphere-stratosphere-troposphere (MST) radars using two or more radar beams. Since the vertical momentum fluxes are assumed to be due to gravity waves, it appears possible to make use of the dispersion and polarization relations for gravity waves in extracting useful information from the radar data. In particular, for an oblique radar beam, information about both the vertical and the horizontal velocities associated with the waves are contained in the measured Doppler data. Therefore, it should be possible to extract both V sub Z and V sub h from a single beam observational configuration. A procedure is proposed to perform such an analysis. The basic assumptions are: the measured velocity fluctuations are due to gravity waves and a separable model gravity-wave spectrum of the Garrett-Munk type that is statistically homogeneous in the horizontal plane. Analytical expressions can be derived that relate the observed velocity fluctuations to the wave momentum flux at each range gate. In practice, the uncertainties related to the model parameters and measurement accuracy will affect the results. A MST radar configuration is considered.

Liu, C. H.↗

Effects of anisotropy on the frequency spectrum of gravity waves observed by MST radar

In the investigation of gravity waves using mesosphere-stratosphere-troposphere radar data, model gravity-wave spectra have been used. In these model spectra, one usually assumes azimuthal symmetry. The effect of spectral anisotropy on the observed spectrum is studied here. It is shown that for a general Garrett-Munk-type spectrum, the anisotropy does not affect the frequency spectrum observed by the vertically beamed radar. For the oblique beam, however, the observed frequency spectrum is changed. A general gravity wave spectrum including azimuthal anisotropy is considered.

Liu, C. H.↗

Doppler effects on velocity spectra observed by MST radars

Recently, wind data from mesophere-stratosphere-troposphere (MST) radars have been used to study the spectra of gravity waves in the atmosphere (Scheffler and Liu, 1985; VanZandt et al., 1985). Since MST radar measures the line-of-sight Doppler velocities, it senses the components of the wave-associated velocities along its beam directions. These components are related through the polarization relations which depend on the frequency and wave number of the wave. Therfore, the radar-observed velocity spectrum will be different from the original gravity-wave spectrum. Their relationship depends on the frequency and wave number of the wave as well as the propagation geometry. This relation can be used to interpret the observed data. It can also be used to test the assumption of gravity-wave spectrum (Scheffler and Liu, 1985). In deriving this relation, the background atmosphere has been assumed to be motionless. Obviously, the Doppler shift due to the background wind will change the shape of the gravity-wave power spectrum as well as its relation with the radar-observed spectrum. Here, researcher's investigate these changes.

Scheffler, A. O.↗

The first operation and results of the Chung-Li VHF radar

The Chung-Li Very High Frequency (VHF) radar is used in the dual-mode operations, applying Doppler beam-swinging as well as the spaced-antenna-drift method. The design of the VHF radar is examined. Results of performance tests are discussed.

Chao, J. K.↗

On the use of colour reflectivity plots to monitor the structure of the troposphere and stratosphere

The radar reflectivity, defined as the range squared corrected power of VHF radar echoes, can be used to monitor and study the temporal development of inversion layer, frontal boundaries and convective turbulence. From typical featurs of upward or downward motion of reflectivity structures, the advection/convection of cold and warm air can be predicted. High resolution color plots appear to be useful to trace and to study the life history of these structures, particularly their persistency, descent and ascent. These displays allow an immediate determination of the tropopause height as well as the determination of the tropopause structure. The life history of warm fronts, cold fronts, and occlusions can be traced, and these reflectivity plots allow detection of even very weak events which cannot be seen in the traditional meteorological data sets. The life history of convective turbulence, particular evolving from the planetary boundary layer, can be tracked quite easily. Its development into strong convection reaching the middle troposphere can be followed and predicted.

Rottger, J.↗

Initiation of breakdown in slender compressible vortices

The initiation of the breakdown process for axially symmetric compressible flows is investigated using a numerical solution of the conservation equations for mass, momentum, and energy. The vortex is isolated, with its axis parallel to the direction of the main stream, and the core radius is small compared to the breakdown length. Computations for several flowfields indicate that the breakdown of the solution is shifted further downstream with increasing Mach number until breakdown is no longer observed. In the subsonic case, the influence of the initial temperature distribution on the breakdown length of the solution is more pronounced than in the supersonic case, with heating of the core enhancing breakdown, and cooling delaying it. The breakdown of the solution is seen to always occur for nonvanishing axial velocity components.

Krause, E.↗

Admissible upstream conditions for slender compressible vortices

The influence of the compressibility on the flow in slender vortices is being studied. The dependence of the breakdown of the slender-vortex approximation on the upstream conditions is demonstrated for various Reynolds numbers and Mach numbers. Compatibility conditions, which have to be satisfied if the vortex is to remain slender, are discussed in detail. The general discussions are supplemented by several sample calculations.

Liu, C. H.↗

Smoothed Two-Dimensional Edges for Laminar Flow

New concept allows passive method for installing flaps, slats, iceprotection equipment, and other leading-edge devices on natural-laminar-flow (NLF) wings without causing loss of laminar flow. Two-dimensional roughness elements in laminar boundary layers strategically shaped to increase critical (allowable) height of roughness. Facilitates installation of leading-edge devices by practical manufacturing methods.

Holmes, B. J.↗

Numerical studies of the merging of vortices

A newly improved method has been employed to study the roll-up and decay of trailing vortex sheets. This method is obtained by using a scheme formulated recently Liu et al., Sept. 1985 to improve the accuracy of the boundary data. The scheme is much more efficient than the previous one used by Weston and Liu (1982). Comparisons of the present numreical results with other numerical solutions of the Navier-Stokes equations are presented.

Weston, R. P.↗

Wall reflection of a viscous vortex ring

The behavior of a viscous axisymmetric vortex ring being reflected from a wall is investigated. The incompressible Navier-Stokes equations formulated in terms of the vorticity function and vector potential are numerically integrated by implicit finite difference methods. To specify the vector potential at a far boundary from the wall, the existing integral method used so far only for an unbounded domain is modified by a kind of image method. The trajectory of the vortex ring calcualted as a result closely resembles that observable from the experiment.

Sa, J. Y.↗

Interaction of decaying trailing vortices in ground shear

The drift of trailing vortices in a cross wind near to the ground is modeled by an unsteady, two-dimensional, rotational flowfield with a concentration of large vorticity in spots having finite total strength. The problem is analyzed by a combination of the method of matched asymptotic analyses for the decay of the vortical spots and the Euler solution for the unsteady rotational flow. A primary feature of the numerical scheme is that the grid size and time step depend only the length and velocity scales of the background flow and is independent of the effective core size of a vortical spot which can be much smaller than the grid size. Numerical results are presented to demonstrated the strong interaction between the trajectjories of the vortical spots and the redistribution of vorticity in the background flowfield.

Liu, C. H.↗

The speed of wave-wave interactions in the atmosphere

Resonant wave-wave interactions are considered. Studies of the interaction coefficient show that rapid transfer of wave action can take place in the disjoined parts of the spectrum for three processes, namely: elastic scattering, parametric subharmonic instability and induced diffusion. Of the three processes, the vertical shear plays a role in two. The vertical shear of a moderate scale interacts through elastic scattering to make the spectrum vertically symmetric. On the other hand, the vertical shear of a large scale interacts through induced diffusion and is responsible for diffusion in k sub z space. When interacting with a vertical shear, it is known that the vertical shear acts as a catalyst and is not involved in energy transfer. Consequently, in both elastic scattering and induced diffusion, the vertical shear does not gain or lose energy. Through parametric subharmonic instability the more energetic large-scale waves are feeding energy into moderate and small scale waves of an elevation angle of 60 deg or larger.

Yeh, K. C.↗

Vortex-dominated flow with viscous core structure

Recent theoretical studies of vortex-dominated flows are reviewed with special emphasis on those for which the viscous core structures play an important role. The problems to be described are: The interaction and merging of two-dimensional vortices and of curved vortex filaments, the roll-up and decay of trailing far wakes, and the initiation of vortex breakdown. The analysis utilizes finite-difference solutions of the Navier-Stokes equations complemented by asymptotic expansion techniques.

Liu, C. H.↗

An asymptotic theory for the interference of large aspect ratio swept wings and multiple propeller slipstreams

This paper presents an asymptotic method for the analysis of the interference of multiple tractor propeller slipstream with large aspect ratio swept wings. It is assumed that the height of the slipstream is of the order of the wing chord and its spanwise extent is of the order of the wing span. Three different flow regions are identified by employing different stretching transformations. Asymptotic expansions are made in each of the three regions, using the chord to span ratio as the small expansion parameter. The details of the nonuniform flow in the slipstream enter into the wing sectional analysis. In the outer limit, the wing shrinks to a swept lifting line, and the slipstream reduces to a thin sheet of jet carrying the momentum gain from the propeller. The curvature of this jet sheet results in a pressure difference which is represented by a vortex sheet. The governing equations are solved by discretization. Several examples are considered for which experimental data are available. Comparison of the present results with the experimental data as well as other numerical solutions showed generally good agreement.

Prabhu, R. K.↗