Efficiency of Mixing Forced by Unsteady Shear Flow

2009 ◽  
Vol 39 (5) ◽  
pp. 1150-1166 ◽  
Author(s):  
Ryuichiro Inoue ◽  
William D. Smyth

Abstract The dependence of mixing efficiency on time-varying forcing is studied by direct numerical simulation (DNS) of Kelvin–Helmholtz (KH) instability. Time-dependent forcing fields are designed to reproduce a wavelike oscillation by solving the equations of motion in a tilted coordinate frame and allowing the tilt angle to vary in time. Mixing efficiency Γc is defined as the ratio of potential energy gain to dissipation, both averaged over one forcing cycle and first examined via parameters characterizing waves: the minimum Richardson number Rimin and the normalized frequency of the forcing ω/N. The effect of Reynolds number Re0 and the initial random disturbance amplitude b are also examined. In the experiments presented, Γc varies between 0.21 and 0.36 and is controlled by the timing of two events: the emergence of KH billows and the arrival of the deceleration of the mean shear by the wavelike forcing. Here, Γc is higher than a canonical value of 0.2 when the deceleration phase of the forcing suppresses the less efficient turbulence after breakdown of KH billows. However, when Rimin and ω/N are small, KH billows start to develop before Rimin is achieved. Therefore, the forcing accelerates the mean shear and thereby sustains turbulence after the breakdown of KH billows. The canonical value is then reproduced in the DNS. Although larger values of Re0 and b intensify the development of KH billows and modify Γc, this effect is less significant when forcing fields act to sustain turbulence. The time-averaged Thorpe scale and Ozmidov scale are also used to see how mixing is modified by forcing fields and compared with past microstructure measurements. It is found that DNS also corresponds to past observations if the forcing accelerates the mean shear to sustain turbulence.

2015 ◽  
Vol 45 (10) ◽  
pp. 2522-2543 ◽  
Author(s):  
Alberto Scotti

AbstractThis paper uses the energetics framework developed by Scotti and White to provide a critical assessment of the widely used Thorpe-scale method, which is used to estimate dissipation and mixing rates in stratified turbulent flows from density measurements along vertical profiles. This study shows that the relevant displacement scale in general is not the rms value of the Thorpe displacement. Rather, the displacement field must be Reynolds decomposed to separate the mean from the turbulent component, and it is the turbulent component that ought to be used to diagnose mixing and dissipation. In general, the energetics of mixing in an overall stably stratified flow involves potentially complex exchanges among the available potential energy and kinetic energy associated with the mean and turbulent components of the flow. The author considers two limiting cases: shear-driven mixing, where mixing comes at the expense of the mean kinetic energy of the flow, and convective-driven mixing, which taps the available potential energy of the mean flow to drive mixing. In shear-driven flows, the rms of the Thorpe displacement, known as the Thorpe scale is shown to be equivalent to the turbulent component of the displacement. In this case, the Thorpe scale approximates the Ozmidov scale, or, which is the same, the Thorpe scale is the appropriate scale to diagnose mixing and dissipation. However, when mixing is driven by the available potential energy of the mean flow (convective-driven mixing), this study shows that the Thorpe scale is (much) larger than the Ozmidov scale. Using the rms of the Thorpe displacement overestimates dissipation and mixing, since the amount of turbulent available potential energy (measured by the turbulent displacement) is only a fraction of the total available potential energy (measured by the Thorpe scale). Corrective measures are discussed that can be used to diagnose mixing from knowledge of the Thorpe displacement. In a companion paper, Mater et al. analyze field data and show that the Thorpe scale can indeed be much larger than the Ozmidov scale.


1966 ◽  
Vol 25 ◽  
pp. 373
Author(s):  
Y. Kozai

The motion of an artificial satellite around the Moon is much more complicated than that around the Earth, since the shape of the Moon is a triaxial ellipsoid and the effect of the Earth on the motion is very important even for a very close satellite.The differential equations of motion of the satellite are written in canonical form of three degrees of freedom with time depending Hamiltonian. By eliminating short-periodic terms depending on the mean longitude of the satellite and by assuming that the Earth is moving on the lunar equator, however, the equations are reduced to those of two degrees of freedom with an energy integral.Since the mean motion of the Earth around the Moon is more rapid than the secular motion of the argument of pericentre of the satellite by a factor of one order, the terms depending on the longitude of the Earth can be eliminated, and the degree of freedom is reduced to one.Then the motion can be discussed by drawing equi-energy curves in two-dimensional space. According to these figures satellites with high inclination have large possibilities of falling down to the lunar surface even if the initial eccentricities are very small.The principal properties of the motion are not changed even if plausible values ofJ3andJ4of the Moon are included.This paper has been published in Publ. astr. Soc.Japan15, 301, 1963.


1988 ◽  
Vol 66 (7) ◽  
pp. 576-579
Author(s):  
G. T. Karahalios ◽  
C. Sfetsos

A sphere executes small-amplitude linear and torsional oscillations in a fluid at rest. The equations of motion of the fluid are solved by the method of successive approximations. Outside the boundary layer, a steady secondary flow is induced in addition to the time-varying motion.


Author(s):  
X. Tong ◽  
B. Tabarrok

Abstract In this paper the global motion of a rigid body subject to small periodic torques, which has a fixed direction in the body-fixed coordinate frame, is investigated by means of Melnikov’s method. Deprit’s variables are introduced to transform the equations of motion into a form describing a slowly varying oscillator. Then the Melnikov method developed for the slowly varying oscillator is used to predict the transversal intersections of stable and unstable manifolds for the perturbed rigid body motion. It is shown that there exist transversal intersections of heteroclinic orbits for certain ranges of parameter values.


2019 ◽  
Vol 89 (5) ◽  
pp. 646
Author(s):  
В.М. Жданов

AbstractBarodiffusion in slow flows of a gas mixture is studied with an approximation using hydrodynamic equations of motion for the individual mixture components. It is shown that consideration of the viscous momentum transfer and the contribution of Knudsen layers for the mixture flowing in a channel has a considerable effect on the value of the barodiffusion factor. The relations are obtained for the mean diffusion fluxes of components and for the total flux of the mixture in a circular cylindrical capillary; these relations are valid for moderately small Knudsen numbers used for calculation of the diffusion baroeffect and separation effect when the gas mixture flows in a set of capillaries connecting two volumes. The modification of the relations for the barodiffusion factor (and for the diffusion slip coefficient cross-linked with it) allows interpreting the sign alteration of these effects observed experimentally for some gas mixtures at intermediate Knudsen numbers.


Author(s):  
S. E. Abd El-Bar

Under the influence of some different perturbations, we study the stability of collinear equilibrium points of the Restricted Three Body Problem. More precisely, the perturbations due to the triaxiality of the bigger primary and the oblateness of the smaller primary, in addition to the relativistic effects, are considered. Moreover, the total potential and the mean motion of the problem are obtained. The equations of motion are derived and linearized around the collinear points. For studying the stability of these points, the characteristic equation and its partial derivatives are derived. Two real and two imaginary roots of the characteristic equation are deduced from the plotted figures throughout the manuscript. In addition, the instability of the collinear points is stressed. Finally, we compute some selected roots corresponding to the eigenvalues which are based on some selected values of the perturbing parameters in the Tables 1, 2.


2020 ◽  
Author(s):  
Kian Boon Law ◽  
Kalaiarasu M Peariasamy ◽  
Balvinder Singh Gill ◽  
Sarbhan Singh Lakha Singh ◽  
Bala Murali Sundram ◽  
...  

Abstract The susceptible-infectious-removed (SIR) model offers the simplest framework to study transmission dynamics of COVID-19, however, it does not factor in its early depleting trend observed during a lockdown. We modified the SIR model to specifically simulate the early depleting transmission dynamics of COVID-19 to better predict its temporal trend in Malaysia. The classical SIR model was fitted to observed total (I total), active (I), and removed (R) cases of COVID-19 before lockdown to estimate the basic reproduction number. Next, the model was modified with a partial time-varying force of infection, given by a proportionally depleting transmission coefficient, βt, and a fractional term, z. The modified SIR model was then fitted to observed data over 6 weeks during the lockdown. Model fitting and projection were validated using the mean absolute percent error (MAPE). The transmission dynamics of COVID-19 was interrupted immediately by the lockdown. The modified SIR model projected the depleting temporal trends with lowest MAPE for I total, followed by I, I daily, and R. During lockdown, the dynamics of COVID-19 depleted at a rate of 4·7% each day with a decreased capacity of 40%. For 7–day and 14–day projections, the modified SIR model accurately predicted I total, I, and R. The depleting transmission dynamics for COVID-19 during lockdown can be accurately captured by time-varying SIR model. Projection generated based on observed data is useful for future planning and control of COVID-19.


Author(s):  
F. Zhang ◽  
B. J. Gilmore ◽  
A. Sinha

Abstract Tolerance allocation standards do not exist for mechanical systems whose response are time varying and are subjected to discontinuous forcing functions. Previous approaches based on optimization and numerical integration of the dynamic equations of motion encounter difficulty with determining sensitivities around the force discontinuity. The Alternating Frequency/Time approach is applied here to capture the effect of the discontinuity. The effective link length model is used to model the system and to account for the uncertainties in the link length, radial clearance and pin location. Since the effective link length model is applied, the equations of motion for the nominal system can be applied for the entire analysis. Optimization procedure is applied to the problem where the objective is to minimize the manufacturing costs and satisfy the constraints imposed on mechanical errors and design variables. Examples of tolerance allocation are presented for a single cylinder internal combustion engine.


Author(s):  
Yu Wang

Abstract A model is developed for analyzing mechanical systems with a pair of bodies with topological changes in their kinematic constraints. It is built upon the concept of Poincaré map rather than following the traditional methods of differential equations. The model provides a set of well-defined and naturally-discrete equations of motion and is capable of giving physical insights of dynamic characteristics of deadbeat convergence of multiple collisions and periodic or chaotic responses. The development of dynamic model and a local stability analysis are presented in Part 1, and the global analysis and numerical simulation are discussed in Part 2.


1964 ◽  
Vol 54 (1) ◽  
pp. 263-276
Author(s):  
J. E. Goldberg ◽  
J. L. Bogdanoff ◽  
D. R. Sharpe

abstract Stationary properties of the response of simple nonlinear systems to earthquake type (nonstationary) random disturbances are determined. The method used consists in generating member functions of the input processes, integrating the equations of motion using these as inputs to obtain output member functions, and then, from this ensemble of outputs, determining the statistical properties of interest. The method is flexible in that it is easy to apply to simple as well as complex linear and nonlinear systems; it yields a variety of statistical properties of engineering significance; it provides a method for approaching the design of aseismic structures; and more generally, shock resistant structures on a practical basis. In addition to determining and displaying input and corresponding output member functions, we determine individual as well as average velocity spectra, first passage time probabilities, and extreme value distributions of relative displacements. Comments are also offered on the utility of these results in design.


Sign in / Sign up

Export Citation Format

Share Document