Search NASA⌕ Search

SEARCH · Search NASA

Results for “Dissipative dynamics”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 325 records · Page 18

Large-scale computations in fluid mechanics; Proceedings of the Fifteenth Summer Seminar on Applied Mathematics, University of California, La Jolla, CA, June 27-July 8, 1983. Parts 1 & 2

Papers are presented on such topics as the use of semi-Lagrangian advective schemes in meteorological modeling; computation with high-resolution upwind schemes for hyperbolic equations; dynamics of flame propagation in a turbulent field; a modified finite element method for solving the incompressible Navier-Stokes equations; computational fusion magnetohydrodynamics; and a nonoscillatory shock capturing scheme using flux-limited dissipation. Consideration is also given to the use of spectral techniques in numerical weather prediction; numerical methods for the incorporation of mountains in atmospheric models; techniques for the numerical simulation of large-scale eddies in geophysical fluid dynamics; high-resolution TVD schemes using flux limiters; upwind-difference methods for aerodynamic problems governed by the Euler equations; and an MHD model of the earth's magnetosphere.

Engquist, B. E.↗

Accretion disk boundary layers in cataclysmic variables. 1: Optically thick boundary layers

We develop numerical models of accretions disks in cataclysmic variables (CVs), including and emphasizing the boundary layer region where the accretion disk meets the accreting white dwarf. We confine ourselves to solutions where the boundary layer region is vertically optically thick, and find that these solutions share several common features. The angular and radial velocities of the accreting material drop rapidly in a dynamical boundary layer, which has a radial width approximately 1%-3% of the white dwarf radius. The energy dissipated in this region diffuses through the inner part of the disk and is radiated from the disk surface in a thermal boundary layer, which has a radial width comparable to the disk thickness, approximately 5%-15% of the white dwarf radius. We examine the dependence of the boundary layer structure on the mass accretion rate, the white dwarf mass and rotation rate, and the viscosity parameter alpha. We delineate the boundary between optically thick and optically thin boundary layer solutions as a function of these parameters and suggest that by means of a careful comparison with observations it may be possible to estimate alpha in CVs. We derive an expression for the total boundary layer luminosities as a function of the parameters and show that it agrees well with the luminosites of our numerical solutions. Finally, we calcuate simple blackbody continuum spectra of the boundary layer and disk emission for our solutions and compare these to soft X-ray, EUV, and He II emission-line observations of CVs. We show that, through such comparisons, it may be possible to determine the rotation rates of the accreting stars in CVs, and perhaps also the white dwarf masses and the accretion rates. The spectra are quite insensitive to alpha, so the uncertainty in this parameter does not affect such comparisons.

Popham, Robert↗

Virtual Passive Controller for Robot Systems Using Joint Torque Sensors

This paper presents a control method based on virtual passive dynamic control that will stabilize a robot manipulator using joint torque sensors and a simple joint model. The method does not require joint position or velocity feedback for stabilization. The proposed control method is stable in the sense of Lyaponov. The control method was implemented on several joints of a laboratory robot. The controller showed good stability robustness to system parameter error and to the exclusion of nonlinear dynamic effects on the joints. The controller enhanced position tracking performance and, in the absence of position control, dissipated joint energy.

Aldridge, Hal A.↗

Wavenumber Vacillation in Weakly-Stratified Baroclinic Flows

The wavenumber selection in an unstable baroclinic flow has long been a subject of investigation in geophysical fluid dynamics. For a weakly unstable flow, the linearly most unstable wave is always the final realized wave. As the flow becomes more unstable, the dominant wave gradually shifts from the most unstable wave to a longer wave. However, when the supercriticality is reduced from the state with a dominant longwave, it may remain in that state even beyond the threshold where it first developed. Hysteresis in a quasigeostrophic system also displays a complex behavior: it occurs not only between the states of different dominant wavenumbers, but also between the states of identical dominant wavenumber but of different dynamic characteristics. In a recent numerical experiment using a high-resolution spectral Eady model with asymmetric Ekman dissipation, it is found that in a weakly stratified flow with strong meridional forcing the flow does not show a preference for any single wave. Instead, the dominant wave periodically varies among the linearly most unstable wave and the longer waves. The mechanism of wavenumber vacillation will be examined in terms of meridional heat flux, the wave PV gradient, and wave saturation. The importance of wave-wave and wave-mean flow interactions will also be investigated.

Chou, Shih-Hung↗

Microgravity Diode Laser Spectroscopy Measurements in a Reacting Vortex Ring

The technique of Diode Laser Spectroscopy (DLS) with wavelength modulation is utilized to measure the concentration of methane in reacting vortex rings under microgravity conditions. From the measured concentration of methane, other major species such as water, carbon dioxide, nitrogen, and oxygen can be easily computed under the assumption of equilibrium chemistry with the method of Interactive Temperature with Assumed Chemistry (ITAC). The conserved scalar approach in modelling the coupling between fluid dynamics and combustion is utilized to represent the unknown variables in terms of the mixture fraction and scalar dissipation rate in conjunction with ITAC. Post-processing of the DLS measurements and the method of ITAC used in computing the species concentration are discussed. From the flame luminosity results, the increase in ring circulation appears to increase the fuel consumption rate inside the reacting vortex ring and the flame height for cases with similar fuel volumes. Preliminary results and application of ITAC show some potential capabilities of ITAC in DLS. The measured concentration of methane, and computed concentrations of water and carbon dioxide agree well with available results from numerical simulations.

Chen, Shin-Juh↗

Microgravity Diode Laser Spectroscopy Measurements in a Reacting Vortex Ring

The technique of Diode Laser Spectroscopy (DLS) with wavelength modulation is utilized to measure the concentration of methane in reacting vortex rings under microgravity conditions. From the measured concentration of methane, other major species such as water, carbon dioxide, nitrogen, and oxygen can be easily computed under the assumption of equilibrium chemistry with an iterative method called ITAC (Iterative Temperature with Assumed Chemistry). The conserved scalar approach in modelling the coupling between fluid dynamics and combustion is utilized to represent the unknown variables in terms of the mixture fraction and scalar dissipation rate in conjunction with ITAC. Post-processing of the DLS and the method used to compute the species concentration are discussed. From the flame luminosity results, ring circulation appears to increase the fuel consumption rate inside the reacting vortex ring and the flame height for cases with similar fuel volumes but different ring circulations. The concentrations of methane, water, and carbon dioxide agree well with available results from numerical simulations.

Chen, Shin-Juh↗

Model-Invariant Hybrid Computations of Separated Flows for RCA Standard Test Cases

NASA's Revolutionary Computational Aerosciences (RCA) subproject has identified several smooth-body separated flows as standard test cases to emphasize the challenge these flows present for computational methods and their importance to the aerospace community. Results of computations of two of these test cases, the NASA hump and the FAITH experiment, are presented. The computations were performed with the model-invariant hybrid LES-RANS formulation, implemented in the NASA code VULCAN-CFD. The model- invariant formulation employs gradual LES-RANS transitions and compensation for model variation to provide more accurate and efficient hybrid computations. Comparisons revealed that the LES-RANS transitions employed in these computations were sufficiently gradual that the compensating terms were unnecessary. Agreement with experiment was achieved only after reducing the turbulent viscosity to mitigate the effect of numerical dissipation. The stream-wise evolution of peak Reynolds shear stress was employed as a measure of turbulence dynamics in separated flows useful for evaluating computations.

Woodruff, Stephen↗

Updates on Modeling the Water Cycle with the NASA Ames Mars Global Climate Model

Global Circulation Models (GCMs) have made steady progress in simulating the current Mars water cycle. It is now widely recognized that clouds are a critical component that can significantly affect the nature of the simulated water cycle. Two processes in particular are key to implementing clouds in a GCM: the microphysical processes of formation and dissipation, and their radiative effects on heating/ cooling rates. Together, these processes alter the thermal structure, change the dynamics, and regulate inter-hemispheric transport. We have made considerable progress representing these processes in the NASA Ames GCM, particularly in the presence of radiatively active water ice clouds. We present the current state of our group's water cycle modeling efforts, show results from selected simulations, highlight some of the issues, and discuss avenues for further investigation.­

Global Circulation Models (GCMs)↗

Magnetotail Current Sheet Thinning and Magnetic Reconnection Dynamics in Global Modeling of Substorms

Magnetotail current sheet thinning and magnetic reconnection are key elements of magnetospheric substorms. We utilized the global MHD model BATS-R-US with Adaptive Mesh Refinement developed at the University of Michigan to investigate the formation and dynamic evolution of the magnetotail thin current sheet. The BATSRUS adaptive grid structure allows resolving magnetotail regions with increased current density up to ion kinetic scales. We investigated dynamics of magnetotail current sheet thinning in response to southwards IMF turning. Gradual slow current sheet thinning during the early growth phase become exponentially fast during the last few minutes prior to nightside reconnection onset. The later stage of current sheet thinning is accompanied by earthward flows and rapid suppression of normal magnetic field component $B-z$. Current sheet thinning set the stage for near-earth magnetic reconnection. In collisionless magnetospheric plasma, the primary mechanism controlling the dissipation in the vicinity of the reconnection site is non-gyrotropic effects with spatial scales comparable with the particle Larmor radius. One of the major challenges in global MHD modeling of the magnetotail magnetic reconnection is to reproduce fast reconnection rates typically observed in smallscale kinetic simulations. Bursts of fast reconnection cause fast magnetic field reconfiguration typical for magnetospheric substorms. To incorporate nongyritropic effects in diffusion regions we developed an algorithm to search for magnetotail reconnection sites, specifically where the magnetic field components perpendicular to the local current direction approaches zero and form an X-type configuration. Spatial scales of the diffusion region and magnitude of the reconnection electric field are calculated self-consistently using MHD plasma and field parameters in the vicinity of the reconnection site. The location of the reconnection sites and spatial scales of the diffusion region are updated during the simulations. Such an approach allows quantifying the interaction between large-scale global magnetospheric dynamics and microphysical processes in diffusion regions localized near reconnection sites. To clarify the role of smallscale non-MHD effects in diffusion region on the global magnetospheric dynamic and to test different models of dissipation we perform simulations with steady southward IMF driving of the magnetosphere.

Kuznetsova, M. M.↗

Hydrodynamics of rotating stars and close binary interactions: Compressible ellipsoid models

We develop a new formalism to study the dynamics of fluid polytropes in three dimensions. The stars are modeled as compressible ellipsoids, and the hydrodynamic equations are reduced to a set of ordinary differential equations for the evolution of the principal axes and other global quantities. Both viscous dissipation and the gravitational radiation reaction are incorporated. We establish the validity of our approximations and demonstrate the simplicity and power of the method by rederiving a number of known results concerning the stability and dynamical oscillations of rapidly rotating polytropes. In particular, we present a generalization to compressible fluids of Chandrasekhar's classical results for the secular and dynamical instabilities of incompressible Maclaurin spheroids. We also present several applications of our method to astrophysical problems of great current interest, such as the tidal disruption of a star by a massive black hole, the coalescence of compact binaries driven by the emission of gravitational waves, and the development of instabilities in close binary systems.

Lai, Dong↗

A Lagrangian dynamic subgrid-scale model turbulence

A new formulation of the dynamic subgrid-scale model is tested in which the error associated with the Germano identity is minimized over flow pathlines rather than over directions of statistical homogeneity. This procedure allows the application of the dynamic model with averaging to flows in complex geometries that do not possess homogeneous directions. The characteristic Lagrangian time scale over which the averaging is performed is chosen such that the model is purely dissipative, guaranteeing numerical stability when coupled with the Smagorinsky model. The formulation is tested successfully in forced and decaying isotropic turbulence and in fully developed and transitional channel flow. In homogeneous flows, the results are similar to those of the volume-averaged dynamic model, while in channel flow, the predictions are superior to those of the plane-averaged dynamic model. The relationship between the averaged terms in the model and vortical structures (worms) that appear in the LES is investigated. Computational overhead is kept small (about 10 percent above the CPU requirements of the volume or plane-averaged dynamic model) by using an approximate scheme to advance the Lagrangian tracking through first-order Euler time integration and linear interpolation in space.

Meneveau, C.↗

Score-based deterministic density sampling

We propose a deterministic sampling framework using Score-Based Transport Modeling for sampling an unnormalized target density π given only its score ∇ log π. Our method approximates the Wasserstein gradient flow on KL($f_t$∥π) by learning the time-varying score ∇ log $f_t$ on the fly using score matching. While having the same marginal distribution as Langevin dynamics, our method produces smooth deterministic trajectories, resulting in monotone noise-free convergence. We prove that our method dissipates relative entropy at the same rate as the exact gradient flow, provided sufficient training. Numerical experiments validate our theoretical findings: our method converges at the optimal rate, has smooth trajectories, and is often more sample efficient than its stochastic counterpart. Experiments on high-dimensional image data show that our method produces high-quality generations in as few as 15 steps and exhibits natural exploratory behavior. The memory and runtime scale linearly in the sample size.

97 MATHEMATICS AND COMPUTING↗

Computation of turbulent flows

All of the work deals with a model that incorporates dynamical equations for the mean velocity field, for all components of the Reynolds stress tensor, and for the 'isotropic dissipation'. The trace of the Reynolds stress tensor provides a velocity scale for the turbulence. For completeness, a turbulence length or time scale is also needed; many workers impose the length scales, but in a generalized model it must be evolved by the turbulence itself. In order, the following special classes of flows are considered: (1) decay of isotropic turbulence, (2) return to isotropy in homogeneous turbulence without strain, (3) homogeneous turbulence with strain, and (4) inhomogeneous turbulence with strain and turbulent transport.

Reynolds, W. C.↗

Particle flux decrease-increase events at synchronous orbit and the temporal sequence of aurora during substorms

A systematic temporal correlation has been found between the energetic particle intensity variations measured at 6.6 earth radii and the development of large scale auroral features. The intensification and equatorward drift of eastwest oriented stable discrete homogeneous auroral arc systems coincide with the decrease in energetic particle intensity at 6.6 earth radii as the nightside magnetosphere develops into a more tail-like configuration. The subsequent major breakup of the aurora coincides with the recovery in particle intensity as the field returns to a more dipolar configuration. Since this prominent decrease-increase sequence must be related to the intensification or inward convection of the tail current plasma sheet configuration followed by its diversion or dissipation, the auroral correlation documented here closely links the auroral particle precipitation to the plasma sheet and tail current dynamics.

Erickson, K. N.↗

Calculated viscosity-distance dependence for some actively flowing lavas

The importance of viscosity as a gauge of the various energy and momentum dissipation regimes of lava flows has been realized for a long time. Nevertheless, despite its central role in lava dynamics and kinematics, it remains among the most difficult of flow physical properties to measure in situ during an eruption. Attempts at reconstructing the actual emplacement viscosities of lava flows from their solidified topographic form are difficult. Where data are available on the position of an advancing flow front as a function of time, it is possible to calculate the effective viscosity of the front as a function of distance from the vent, under the assumptions of a steady state regime. As an application and test of an equation given, relevant parameters from five recent flows on Mauna Loa and Kilauea were utilized to infer the dynamic structure of their aggregate flow front viscosity as they advanced, up to cessation. The observed form of the viscosity-distance relation for the five active Hawaiian flows examined appears to be exponential, with a rapid increase just before the flows stopped as one would expect.

Pieri, David↗

The self-consistent dynamic pole tide in non-global oceans

The dynamic pole tide is determined by solving Laplace tide equations which take into account the presence of continents in oceans, oceanic self-gravitation and loading, and mantle elasticity. Dynamical effects are found to be only mild. It is shown that the dynamical pole tide contributes about one day more to the Chandler period than a static pole tide would, and dissipates wobble energy at a very weak rate. It is noted that, depending on the wobble period predicted for an oceanless elastic earth, mantle anelasticity at low frequencies may nevertheless contribute negligibly to the Chandler period.

Dickman, S. R.↗

On rotational discontinuities in both two-fluid and hybrid models

Rotational discontinuities are studied in a two-fluid model that includes finite ion inertia dispersion and in a hybrid model in which the full ion dynamics is retained while the electrons are treated as a massless fluid. It is shown that as in previous dissipative MHD studies, a rotational discontinuity is unstable in both models and evolves to a 2-3 intermediate shock, a slow rarefaction wave, and other waves. In addition, it is shown that the so-called Walen relation, which holds exactly for rotational discontinuities, can also be well satisfied by tran-Alfvenic intermediate shocks. Thus intermediate shocks can be candidates for those observed structures that satisfy the Walen relation.

Wu, C. C.↗

Probability density distribution of velocity differences at high Reynolds numbers

Recent understanding of fine-scale turbulence structure in high Reynolds number flows is mostly based on Kolmogorov's original and revised models. The main finding of these models is that intrinsic characteristics of fine-scale fluctuations are universal ones at high Reynolds numbers, i.e., the functional behavior of any small-scale parameter is the same in all flows if the Reynolds number is high enough. The only large-scale quantity that directly affects small-scale fluctuations is the energy flux through a cascade. In dynamical equilibrium between large- and small-scale motions, this flux is equal to the mean rate of energy dissipation epsilon. The pdd of velocity difference is a very important characteristic for both the basic understanding of fully developed turbulence and engineering problems. Hence, it is important to test the findings: (1) the functional behavior of the tails of the probability density distribution (pdd) represented by P(delta(u)) is proportional to exp(-b(r) absolute value of delta(u)/sigma(sub delta(u))) and (2) the logarithmic decrement b(r) scales as b(r) is proportional to r(sup 0.15) when separation r lies in the inertial subrange in high Reynolds number laboratory shear flows.

Praskovsky, Alexander A.↗