Engineering Papers⌕ Search

Engineering topics

Cheng, Y.

Publications and source records attributed to Cheng, Y..

At least 91 records · Page 5

Vertical Diffusivities of Active and Passive Tracers

The climate models that include a carbon-cycle need the vertical diffusivity of a passive tracer. Since an expression for the latter is not available, it has been common practice to identify it with that of salt. The identification is questionable since T, S are active, not passive tracers. We present the first derivation of the diffusivity of a passive tracer in terms of Ri (Richardson number) and Rq (density ratio, ratio of salinity over temperature z-gradients). The following results have emerged: (a) The passive tracer diffusivity is an algebraic function of Ri, Rq. (b) In doubly stable regimes (DS, partial derivative of T with respect to z > 0, partial derivative of S with respect to z < 0), the passive scalar diffusivity is nearly the same as that of salt/heat for any values of Rq < 0 and Ri > 0. (c) In DC regimes (diffusive convection, partial derivative of T with respect to z < 0, partial derivative of S with respect to z < 0, Rq > 1), the passive scalar diffusivity is larger than that of salt. At Ri = O(1), it can be more than twice as large. (d) In SF regimes (salt fingers, partial derivative of T with respect to z > 0, partial derivative of S with respect to z > 0, Rq < 1), the passive scalar diffusivity is smaller than that of salt. At Ri = O(1), it can be less than half of it. (e) The passive tracer diffusivity predicted at the location of NATRE (North Atlantic Tracer Release Experiment) is discussed. (f) Perhaps the most relevant conclusion is that the common identification of the tracer diffusivity with that of salt is valid only in DS regimes. In the Southern Ocean, where there is the largest CO2 absorption, the dominant regime is diffusive convection discussed in (c) above.

Canuto, V. M.↗

An Attempt to Derive the epsilon Equation from a Two-Point Closure

The goal of this paper is to derive the equation for the turbulence dissipation rate epsilon for a shear-driven flow. In 1961, Davydov used a one-point closure model to derive the epsilon equation from first principles but the final result contained undetermined terms and thus lacked predictive power. Both in 1987 and in 2001, attempts were made to derive the epsilon equation from first principles using a two-point closure, but their methods relied on a phenomenological assumption. The standard practice has thus been to employ a heuristic form of the equation that contains three empirical ingredients: two constants, c(sub 1 epsilon), and c(sub 2 epsilon), and a diffusion term D(sub epsilon) In this work, a two-point closure is employed, yielding the following results: 1) the empirical constants get replaced by c(sub 1), c(sub 2), which are now functions of Kappa and epsilon; 2) c(sub 1) and c(sub 2) are not independent because a general relation between the two that are valid for any Kappa and epsilon are derived; 3) c(sub 1), c(sub 2) become constant with values close to the empirical values c(sub 1 epsilon), c(sub epsilon 2), (i.e., homogenous flows); and 4) the empirical form of the diffusion term D(sub epsilon) is no longer needed because it gets substituted by the Kappa-epsilon dependence of c(sub 1), c(sub 2), which plays the role of the diffusion, together with the diffusion of the turbulent kinetic energy D(sub Kappa), which now enters the new equation (i.e., inhomogeneous flows). Thus, the three empirical ingredients c(sub 1 epsilon), c(sub epsilon 2), D (sub epsilon)are replaced by a single function c(sub 1)(Kappa, epsilon ) or c(sub 2)(Kappa, epsilon ), plus a D(sub Kappa)term. Three tests of the new equation for epsilon are presented: one concerning channel flow and two concerning the shear-driven planetary boundary layer (PBL).

Canuto, V. M.↗

Ocean Turbulence, III: New GISS Vertical Mixing Scheme

We have found a new way to express the solutions of the RSM (Reynolds Stress Model) equations that allows us to present the turbulent diffusivities for heat, salt and momentum in a way that is considerably simpler and thus easier to implement than in previous work. The RSM provides the dimensionless mixing efficiencies Gamma-alpha (alpha stands for heat, salt and momentum). However, to compute the diffusivities, one needs additional information, specifically, the dissipation Epsilon. Since a dynamic equation for the latter that includes the physical processes relevant to the ocean is still not available, one must resort to different sources of information outside the RSM to obtain a complete Mixing Scheme usable in OGCMs. As for the RSM results, we show that the Gamma-alpha s are functions of both Ri and Rq (Richardson number and density ratio representing double diffusion, DD); the Gamma-alpha are different for heat, salt and momentum; in the case of heat, the traditional value Gamma-h = 0.2 is valid only in the presence of strong shear (when DD is inoperative) while when shear subsides, NATRE data show that Gamma-h can be three times as large, a result that we reproduce. The salt Gamma-s is given in terms of Gamma-h. The momentum Gamma-m has thus far been guessed with different prescriptions while the RSM provides a well defined expression for Gamma-m(Ri,R-rho). Having tested Gamma-h, we then test the momentum Gamma-m by showing that the turbulent Prandtl number Gamma-m/Gamma-h vs. Ri reproduces the available data quite well. As for the dissipation epsilon, we use different representations, one for the mixed layer (ML), one for the thermocline and one for the ocean;s bottom. For the ML, we adopt a procedure analogous to the one successfully used in PB (planetary boundary layer) studies; for the thermocline, we employ an expression for the variable epsilon/N(exp 2) from studies of the internal gravity waves spectra which includes a latitude dependence; for the ocean bottom, we adopt the enhanced bottom diffusivity expression used by previous authors but with a state of the art internal tidal energy formulation and replace the fixed Gamma-alpha = 0.2 with the RSM result that brings into the problem the Ri, R-rho dependence of the Gamma-alpha; the unresolved bottom drag, which has thus far been either ignored or modeled with heuristic relations, is modeled using a formalism we previously developed and tested in PBL studies. We carried out several tests without an OGCM. Prandtl and flux Richardson numbers vs. Ri. The RSM model reproduces both types of data satisfactorily. DD and Mixing efficiency Gamma-h(Ri,Rq). The RSM model reproduces well the NATRE data. Bimodal epsilon-distribution. NATRE data show that epsilon (Ri < 1) approximately equals 10epsilon(Ri > 1), which our model reproduces. Heat to salt flux ratio. In the Ri much greater than 1 regime, the RSM predictions reproduce the data satisfactorily. NATRE mass diffusivity. The z-profile of the mass diffusivity reproduces well the measurements at NATRE. The local form of the mixing scheme is algebraic with one cubic equation to solve.

Canuto, V. M.↗

Autonomous landmark tracking orbit determination strategy

In this paper, an orbit determination strategy is described that is fully autonomous and relies on a computer-based crater detection and identification algorithm that is suitable for both automation of the ground based navigation system and autonomous spacecraft based navigation.

landmark tracking orbit determination↗

Optical landmark detection for spacecraft navigation

In this paper, a new crater detection algorithm is suggested. Experimental studies show that this new algorithm can achieve sub-pixel accuracy in position, its detection rate is better than 90% and its false alarm rate is less than 5%. These good characteristics indicate that it is an ideal crater detection algorithm for spacecraft optical navigation.

optical↗

Autonomous landmark based spacecraft navigation system

An autonomous landmark based spacecraft navigation scheme is presented. This new schme involves the following data processing steps: image selection and planning; landmark detection; preliminary image matching; crater matching; database management; and finally landmark based orbit determination.

autonomous↗

Limbed Excursion Mechanical Utility Rover: LEMUR II

The assembly inspection, and maintenance requirements of permanent installations in space demand robotic agents that provide a high level of operational flexibility relative to the mass and volume of the robotic system.

robot limb manipulation robot vision↗

Passive imaging based hazard avoidance for spacecraft safe landing

In this paper, a new and fast passive sensing based hazard avoidance approach for spacecraft safe landing is suggested. This approach contains two important algorithms: a texture based landing site selection algorithm and a landing site slope estimation algorithm.

spacecraft safe landing passive imaging homography↗

New Third-Order Moments for the CBL

Turbulent convection is an inherently non-local phenomenon and a primary condition for a successful treatment of the CBL (convective boundary layer) is a reliable model of non-locality. In the dynamic equations governing the convective flux, the turbulent kinetic energy, etc., non-locality is represented by the third-order moments, TOMs. Since the simplest form, the so-called down gradient approximation (DGA), severely underestimates the TOMs (up to an order of magnitude), a more physical model is needed. In 1994, an analytical model was presented which was derived directly from the dynamical equations for the TOMs. It considerably improved the DGA but was a bit cumbersome to use and, more importantly, it was based on the quasi-normal (QN) approximation for the fourth-order moments. Here, we present a new analytic expression for the TOMs which is structurally simpler than the 1994 expression and which avoids the QN approximation. The resulting fit to the LES data is superior to that of the 1994 model.

Canuto, V. M.↗