Results for “dissipative dynamics” — Page 62

DOE OSTI2024

High-temperature structure, elasticity, and thermal expansion of ε-ZrH 1.8

Zirconium hydride is a promising candidate material for nuclear microreactor applications as a solid-state moderator component, owing to its favorable neutronics properties and good thermal stability over other metal hydrides. Here, in the present work, the crystal structure, thermal expansion, and elastic properties of the hydrogen-rich ε phase hydride were measured at elevated temperatures in the range 300–900 K. Samples were prepared by direct hydriding Zircaloy-4 metal – a nuclear-grade zirconium alloy. Room-temperature lattice parameters agree well with those reported from literature for unalloyed zirconium hydride and fall within an observed quadratic H-content dependence. The coefficients of thermal expansion, determined from lattice expansion and dilatometry, agree well within our work but were about 30 % lower than those reported by others for unalloyed hydrides. Density functional theory-based molecular dynamics simulations were used to compare with thermal expansion and elasticity measurements. Results showed lattice parameter temperature dependence and slope of thermal expansion align with those from measurements. Based on diffraction scans at select temperatures, ε phase remained stable in air up to at least 770 K. Likewise, dilatometry showed smooth thermal expansion up to the thermal decomposition temperature around 950 K. The precise decomposition temperature was not determined via diffraction due to sparse scanning. The complete elastic property measurements were gathered for ε-phase Ziracloy-4 hydride for the first time. Young's modulus was lower compared to the metal and δ hydride phases. High-temperature elasticity measurements were limited to <350 K due to acoustic dissipation effects.

36 MATERIALS SCIENCE↗

NASA NTRS1972

The origin of the moon.

Recent studies of the dynamics and thermodynamics of the early solar nebula have provided a basis for a theory of the origin of the moon which is consistent with presently developed views of the origin of the solar system. The hypothesis of inhomogeneous planetary accretion is extended. It is suggested that the anomalous properties of the moon, such as its enrichment in Ca, Al, Ti, and other refractories and its depletion in iron and volatiles can be explained if the bulk of the moon represents a high temperature condensate. It is proposed that the moon is composed chiefly of compounds that condense before iron and that the volatile content of the moon was brought in as a thin veneer after the solar nebula dissipated.

Anderson, D. L.↗

NASA NTRS1975

Further numerical calculations of the circulation of the atmosphere of Venus

Results of two-dimensional simulations of the deep circulation of Venus are presented. They prove that the high surface temperature can only be explained by the greenhouse effect, and that Goody and Robinson's (1966) dynamical model is not valid. Analytical radiative equilibrium solutions for a semi-grey atmosphere, both with and without an internal heat source, are presented. It is shown that the greenhouse effect is sufficient to produce the high surface temperature if the thermal optical depth is much greater than 100 and the ratio of solar to thermal optical depths does not exceed 0.005. A two-dimensional version of a three-dimensional model is used to test the validity of the new mechanism proposed by Gierasch (1975) to explain the 4-day circulation. Numerical experiments with large horizontal viscosities failed to show strong zonal velocities even for the case of large Prandtl numbers. It is observed that the dissipation of angular momentum introduced by the strong horizontal diffusion more than compensates for the upward transport of angular momentum due to the Hadley cell.

De Rivas, E. K.↗

DOE OSTI2024

Numerical implications of including drifts in SOLPS-ITER simulations of EAST

The inclusion of drifts in plasma edge codes like SOLPS-ITER is required to match simulation data with experimental profiles. However, this remains numerically challenging. In this paper, the effect of some numerical factors on the final plasma solution is investigated. This study is performed on three EAST simulations in the upper single null configuration: an attached purely deuterium case, an attached case with limited Ne-seeding, and a detached Ne-seeded case. The effects of the anomalous conductivity and anomalous thermo-electric coefficient on the plasma potential are investigated. Next, the effect of the employed grids is shown. In order to investigate these effects, accurate drift simulations are needed. Therefore, the employed time step and numerical parameters are discussed for the three studied simulations. For all presented simulations, it is verified that the restriction of the grid to the first flux surface tangent to the main chamber wall is sufficient within the context of non-extended simulations. This means that the main power dissipation takes place inside the simulated domain, and only a small fraction of the power is leaving the B2.5 grid through the grid boundary closest to the first wall. Finally, the effect of drifts on the asymmetry between the inner and outer target for EAST simulations is demonstrated.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2020

Dissipative Assembly of Macrocycles Comprising Multiple Transient Bonds

Abstract Dissipative assembly has great potential for the creation of new adaptive chemical systems. However, while molecular assembly at equilibrium is routinely used to prepare complex architectures from polyfunctional monomers, species formed out of equilibrium have, to this point, been structurally very simple. In most examples the fuel simply effects the formation of a single short‐lived covalent bond. Herein, we show that chemical fuels can assemble bifunctional components into macrocycles containing multiple transient bonds. Specifically, dicarboxylic acids give aqueous dianhydride macrocycles on treatment with a carbodiimide. The macrocycles are assembled efficiently as a consequence of both fuel‐dependent and fuel‐independent mechanisms; they undergo slower decomposition, building up as the fuel recycles the components, and are a favored product of the dynamic exchange of the anhydride bonds. These results create new possibilities for generating structurally sophisticated out‐of‐equilibrium species.

Hossain, Mohammad Mosharraf↗

DOE OSTI2022

Implications of inertial subrange scaling for stably stratified mixing

In this work, we investigate the effects of the turbulent dynamic range on active scalar mixing in stably stratified turbulence by adapting the theoretical passive scalar modelling and demonstrating their usefulness through consideration of the results of direct numerical simulations of statistically stationary homogeneous stratified and sheared turbulence. By analysis of inertial and inertial–convective subrange scalings, we show that the relationship between the active scalar and turbulence time scales is predicted by the ratio of the Kolmogorov and Obukhov–Corrsin constants, provided mean flow parameters permit the two subrange scalings to be appropriate approximations. We use the resulting relationship between time scales to parameterise an appropriate turbulent mixing coefficient $\varGamma \equiv \chi /\epsilon$, defined here as the ratio of available potential energy ($E_p$) and turbulent kinetic energy ($E_k$) dissipation rates. With the analysis presented here, we show that $\varGamma$ can be estimated by $E_p,E_k$ and a universal constant provided an appropriate Reynolds number is sufficiently high. This large Reynolds number regime appears here to occur at $ {{Re_b}} \equiv \epsilon / \nu N^{2} \gtrapprox 300$ where $\nu$ is the kinematic viscosity and $N$ is the characteristic buoyancy frequency. We propose a model framework for irreversible diapycnal mixing with robust theoretical parametrisation and asymptotic behaviour in this high-$ {{Re_b}}$ limit.

42 ENGINEERING↗

DOE OSTI2025

Entropy is an important design principle in the photosystem II supercomplex

Photosystem II (PSII) can achieve near-unity quantum efficiency of light harvesting in ideal conditions and can dissipate excess light energy as heat to prevent the formation of reactive oxygen species (ROS) under light stress. Understanding how this pigment–protein complex accomplishes these opposing goals is a topic of great interest that has so far been explored primarily through the lens of the system energetics. Despite PSII’s known flat energy landscape, a thorough consideration of the entropic effects on energy transfer in PSII is lacking. In this work, we aim to discern the free energetic design principles underlying the PSII energy transfer network. To accomplish this goal, we employ a structure-based rate matrix and compute the free energy terms in time following a specific initial excitation to discern how entropy and enthalpy drive ensemble system dynamics. We find that the interplay between the entropy and enthalpy components differ among each protein subunit, which allows each subunit to fulfill a unique role in the energy transfer network. This individuality ensures that PSII can accomplish efficient energy trapping in the reaction center (RC), effective nonphotochemical quenching (NPQ) in the periphery, and robust energy trapping in the other-monomer RC if the same-monomer RC is closed. We also show that entropy, in particular, is a dynamically tunable feature of the PSII free energy landscape accomplished through regulation of LHCII binding. These findings help rationalize natural photosynthesis and provide design principles for more efficient solar energy harvesting technologies.

59 BASIC BIOLOGICAL SCIENCES↗

NASA NTRS1989

Assessment of numerical techniques for unsteady flow calculations

The characteristics of unsteady flow motions have long been a serious concern in the study of various fluid dynamic and combustion problems. With the advancement of computer resources, numerical approaches to these problems appear to be feasible. The objective of this paper is to assess the accuracy of several numerical schemes for unsteady flow calculations. In the present study, Fourier error analysis is performed for various numerical schemes based on a two-dimensional wave equation. Four methods sieved from the error analysis are then adopted for further assessment. Model problems include unsteady quasi-one-dimensional inviscid flows, two-dimensional wave propagations, and unsteady two-dimensional inviscid flows. According to the comparison between numerical and exact solutions, although second-order upwind scheme captures the unsteady flow and wave motions quite well, it is relatively more dissipative than sixth-order central difference scheme. Among various numerical approaches tested in this paper, the best performed one is Runge-Kutta method for time integration and six-order central difference for spatial discretization.

Hsieh, Kwang-Chung↗

NASA NTRS2013

Carbon nanotube vacuum gauges with wide-dynamic range and processes thereof

A miniature thermal conductivity gauge employs a carbon single-walled-nanotube. The gauge operates on the principle of thermal exchange between the voltage-biased nanotube and the surrounding gas at low levels of power and low temperatures to measure vacuum across a wide dynamic range. The gauge includes two terminals, a source of constant voltage to the terminals, a single-walled carbon nanotube between the terminals, a calibration of measured conductance of the nanotube to magnitudes of surrounding vacuum and a current meter in electrical communication with the source of constant voltage. Employment of the nanotube for measuring vacuum includes calibrating the electrical conductance of the nanotube to magnitudes of vacuum, exposing the nanotube to a vacuum, applying a constant voltage across the nanotube, measuring the electrical conductance of the nanotube in the vacuum with the constant voltage applied and converting the measured electrical conductance to the corresponding calibrated magnitude of vacuum using the calibration. The nanotube may be suspended to minimize heat dissipation through the substrate, increasing sensitivity at even tower pressures.

Manohara, Harish↗

DOE OSTI2021

High-Fidelity Simulations of Gravity Currents Using a High-Order Finite-Difference Spectral Vanishing Viscosity Approach

This numerical work investigates the potential of a high-order finite-difference spectral vanishing viscosity approach to simulate gravity currents at high Reynolds numbers. The method introduces targeted numerical dissipation at small scales through altering the discretisation of the second derivatives of the viscous terms in the incompressible Navier-Stokes equations to mimic the spectral vanishing viscosity (SVV) operator, originally designed for the regularisation of spectral element method (SEM) solutions of pure advection problems. Using a sixth-order accurate finite-difference scheme, the adoption of the SVV method is straightforward and comes with a negligible additional computational cost. In order to assess the ability of this high-order finite-difference spectral vanishing viscosity approach, we performed large-eddy simulations (LES) of a gravity current in a channelised lock-exchange set-up with our SVV model and with the well-known explicit static and dynamic Smagorinsky sub-grid scale (SGS) models. The obtained data are compared with a direct numerical simulation (DNS) based on more than 800 million mesh nodes, and with experimental measurements. A framework for the energy budget is introduced to investigate the behaviour of the gravity current. First, it is found that the DNS is in good agreement with the experimental data for the evolution of the front location and velocity field as well as for the stirring and mixing inside the gravity current. Secondly, the LES performed with less than 0.4% of the total number of mesh nodes compared to the DNS, can reproduce the main features of the gravity currents, with the SVV model yielding slightly more accurate results. It is also found that the dynamic Smagorinsky model performs better than its static version. For the present study, the static and dynamic Smagorinsky models are 1.8 and 2.5 times more expensive than the SVV model, because the latter does not require the calculation of explicit SGS terms in the Navier-Stokes equations nor spatial filtering operations.

direct numerical simulation↗

NASA NTRS1988

Dynamical evolution of a solar coronal magnetic field arcade

Calculations of the long-term dynamical evolution of a solar coronal magnetic field arcade which is subjected to shearing photospheric flows are presented. The evolution is obtained by numerical solution of a subset of the resistive magnetohydrodynamic equations. For a simplified model of the bipolar magnetic field observed in the solar corona, it is found that photospheric flow produces a slow evolution of the magnetic field, with a buildup of magnetic energy. For certain photospheric shear profiles, the field configuration produced is linearly unstable to an ideal magnetohydrodynamic mode when the shear exceeds a critical value. The nonlinear evolution of this instability shows the spontaneous formation of current sheets. Reconnection of the magnetic field produces a rapid release of magnetic energy. The major fraction of the energy is dissipated resistively, while a small fraction is converted into kinetic energy of an ejected plasmoid. The relevance of these results to two-ribbon flares is discussed.

Mikic, Z.↗

NASA NTRS2005

Overview of Results from the Cassini Magnetospheric Imaging Instrument (MIMI) During the First Year of Operations

The MIMI investigation comprises three sensors covering the indicated energy ranges: the Ion and Neutral Camera (INCA) -- 7 keV/nuc 3 Mev) and protons (1.6 < E < 160 Mev) from the back end of the dual field-of-view telescope. The Saturn observation sequences began in January, 2004 and culminated in Saturn Orbit Insertion on July 1, 2004. The MIMI sensors observed substantial activity in interplanetary space for several months prior to SOI, including several interplanetary shocks associated with corotating interaction regions, numerous increases most likely originating from particle streams in the vicinity of the Saturnian bow shock and energetic neutral atoms (ENA) emanating from Saturn s magnetosphere. Results following SOI revealed: a dynamical magnetosphere with a day-night asymmetry and an 11-hour periodicity; several water-product ions (O+, OH+, H2O+), but little N+; inferred quantities of neutral gas sufficient to cause major losses in the trapped ions and electrons in the middle and inner magnetosphere; a Titan exosphere that is a copious source of ENA; INCA imaging through ENA has also revealed a previously unknown radiation belt residing inward of the D-ring that is most likely the result of double charge-exchange between the main radiation belt and the upper layers of Saturn s exosphere. Finally, there is ample evidence for the presence of substorm-like injections of plasma that subsequently corotates for a number of days before dissipating on the night-side magnetotail. The observations will be presented and discussed in the context of current theoretical models.

Krimigis, S. M.↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

DOE OSTI2026

Hydrogeological assessment of CO2 containment assurance and wellbore integrity at a Gulf Coast storage site

Abstract A large-scale carbon capture and storage (CCS) initiative on the Texas Gulf Coast serves as a premier demonstration of the U.S. Department of Energy’s CarbonSAFE program. Targeting deep saline formations, specifically Oligo-Miocene deltaic sequences, the project aims to establish technical and commercial viability for geologic CO2 storage within a major industrial corridor. This study provides a rigorous hydrogeological assessment to support Class VI permitting by quantifying the high degree of containment security. Utilizing a compositional reservoir simulator, we developed a suite of 27 distinct simulation cases to evaluate vertical plume dynamics near both planned injection wells and proximal legacy infrastructure. To ensure numerical accuracy near wellbores, we implemented a refined mesh strategy, determining that a 5.6 ft × 5.6 ft grid refinement offered the optimal balance between computational efficiency and descriptive precision. The modeling framework utilized a systematic sensitivity-based approach to evaluate the mechanical redundancy of the subsurface system by performing a bounding analysis of wellbore interfaces against hypothetical high-permeability microannuli. By systematically isolating competing physical drivers, including permeability, porosity, gas hysteresis, thermal gradients, salinity, and solubility trapping (quantified via Henry’s law with dynamically adjusted coefficients), this work moves beyond binary assessments to establish a nuanced hierarchy of containment factors. The results confirm that primary trapping mechanisms (e.g., gas hysteresis and solubility), combined with the site's unique geomechanical stratigraphy, significantly restrict vertical mobility and reinforce the robust containment security of the reservoir. Baseline results demonstrate substantial vertical separation between the CO2 plume and the upper confining system, ensuring robust containment. Sensitivity analysis reveals that even under highly conservative bounding scenarios—assuming theoretical 10-Darcy pathways at specific wellbore locations—the 2,900-ft thick multi-layered confining zone remains a reliable barrier. In these hypothetical upper-bound cases, peak upward fluxes of CO2 and saltwater after 15 years of injection remain localized and dissipate rapidly within the lower sections of the confining interval, leaving the integrity of the seal uncompromised. Furthermore, the study identifies that while localized wellbore pathways define theoretical upper bounds of vertical migration, the Area of Review (AoR) is primarily sensitive to regional thermal gradients and hysteresis, which can influence the AoR by over 3,000 acres in pessimistic configurations. Also, primary trapping mechanisms, specifically gas hysteresis and solubility, work in tandem with the Gulf Coast’s unique geomechanical stratigraphy to significantly restrict vertical mobility. Ductile, smectite-rich mudstones facilitate natural borehole convergence and the self-healing of potential conduits, creating a natural geomechanical bridge that effectively mitigates migration potential at both current injection points and legacy-well locations. This comprehensive modeling effort demonstrates that the integration of high-resolution wellbore simulations and regional geomechanical observations confirms the long-term storage security of the studied site, providing a physics-based foundation for industrial-scale CCS deployments. This modeling framework establishes a baseline for future research into coupled geomechanical effects, such as time-dependent borehole convergence, to further refine long-term containment projections. Acknowledgements We thank the Gulf Coast Carbon Center (GCCC) at the Bureau of Economic Geology for foundational research support. We appreciate Alex Bump for technical guidance and David Hoffman for model mesh generation. This work used TACC’s Frontera cluster for simulations and CMG Ltd. software licenses provided to UT-Austin. This material is based upon work supported by the Department of Energy under Award Number DE-FE0032338. Disclaimer This material is based upon work supported by the U.S. Department of Energy’s Fossil Energy and Carbon Management Office under the CarbonSAFE program, award Number DE-FE0032338. The views expressed herein do not necessarily represent the views of the U.S. Department of Energy or the United States Government.

58 GEOSCIENCES↗

NASA NTRS2012

Investigation of Transonic Wake Dynamics for Mechanically Deployable Entry Systems

A numerical investigation of transonic flow around a mechanically deployable entry system being considered for a robotic mission to Venus has been performed, and preliminary results are reported. The flow around a conceptual representation of the vehicle geometry was simulated at discrete points along a ballistic trajectory using Detached Eddy Simulation (DES). The trajectory points selected span the low supersonic to transonic regimes with freestream Mach numbers from 1:5 to 0:8, and freestream Reynolds numbers (based on diameter) between 2:09 x 10(exp 6) and 2:93 x 10(exp 6). Additionally, the Mach 0:8 case was simulated at angles of attack between 0 and 5 . Static aerodynamic coefficients obtained from the data show qualitative agreement with data from 70deg sphere-cone wind tunnel tests performed for the Viking program. Finally, the effect of choices of models and numerical algorithms is addressed by comparing the DES results to those using a Reynolds Averaged Navier-Stokes (RANS) model, as well as to results using a more dissipative numerical scheme.

Stern, Eric↗

DOE OSTI2021

Electron dynamics in radio frequency magnetron sputtering argon discharges with a dielectric target

Abstract We demonstrate a self-consistent and complete description of electron dynamics in a typical electropositive radio frequency magnetron sputtering (RFMS) argon discharge with a dielectric target. The electron dynamics, including the electron power absorption dynamics in one radio frequency (RF) period, is studied via a fully kinetic 2d3v particle-in-cell/Monte Carlo collision (PIC/MCC) electrostatic simulation. The interplay between the fundamental plasma parameters is analyzed through their spatiotemporal dynamics. Due to the influence of magnetic trap on the electron transport, a spatially dependent charging that perturbs the electric potential is observed on the dielectric target surface, resulting in a spatially dependent ion energy distribution along the target surface. The E × B drift-to-discharge current ratio is in approximate agreement with Bohm diffusion. The electron power absorption can be primarily decoupled into the positive Ohmic power absorption in the bulk plasma region and the negative pressure-induced power absorption near the target surface. Ohmic power absorption is the dominant electron power absorption mechanism, mostly contributed by the azimuthal electron current. The power absorption due to electron inertial effects is negligible on time-average. Both the maximum power absorption and dissipation of electrons appear in the bulk plasma region during the second half of the RF period, implying a strong electron trapping in magnetron discharges. The contribution of secondary electrons is negligible under typical RFMS discharge conditions.

Physics↗

NASA NTRS1993

Accuracy considerations in the computational analysis of jet noise

The application of computational fluid dynamics methods to the analysis of problems in aerodynamic noise has resulted in the extension and adaptation of conventional CFD to the discipline now referred to as computational aeroacoustics (CAA). In the analysis of jet noise accurate resolution of a wide range of spatial and temporal scales in the flow field is essential if the acoustic far field is to be predicted. The numerical simulation of unsteady jet flow has been successfully demonstrated and many flow features have been computed with reasonable accuracy. Grid refinement and increased solution time are discussed as means of improving accuracy of Navier-Stokes solutions of unsteady jet flow. In addition various properties of different numerical procedures which influence accuracy are examined with particular emphasis on dispersion and dissipation characteristics. These properties are investigated by using selected schemes to solve model problems for the propagation of a shock wave and a sinusoidal disturbance. The results are compared for the different schemes.

Scott, James N.↗

NASA NTRS1988

On the prediction of equilibrium states in homogeneous turbulence

A comparison of several commonly used turbulence models (including the Kappa-epsilon and two second-order closures) is made for the test problem of homogeneous turbulent shear flow in a rotating frame. The time evolution of the turbulent kinetic energy and dissipation rate is calculated for a variety of models and comparisons are made with previously published experiments and numerical simulations. Particular emphasis is placed on examining the ability of each model to accurately predict equilibruim states for a range of the parameter Omega/S (the ratio of the rotation rate to the shear rate). It is found that none of the commonly used second-order closure models yield substantially improved predictions for the time evolution of the turbulent kinetic energy and dissipation rate over the somewhat defective results obtained from the simpler Kappa-epsilon model for the turbulent flow regime. There is also a problem with the equilibrium states predicted by the various models. For example, the Kappa-epsilon model erroneously yields equilibrium states that are independent of Omega/S while the Launder, Reece, and Rodi model predicts a flow relaminarization when Omega/S is greater than 0.39 - a result which is contrary to numerical simulations and linear spectral analysis which indicate flow instability for at least the range 0 less than or = Omega/S less than or = 0.5. The physical implications of the results obtained from the various turbulence models considered here are discussed in detail along with proposals to remedy the deficiencies based on a dynamical systems approach.

Speziale, Charles G.↗