Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “dynamical orthogonality”

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 289 records · Page 16
All sourcesNASA NTRSDOE OSTI
NASA NTRS1988

Stability analysis of multiple rigid robot manipulators holding a common rigid object

The authors consider several possible control structures for multiple-arm systems by regarding either joint torques, tip force, or a generalized acceleration of the control input. They emphasize the first case, since a class of relatively model-independent control laws can be generated for both motion and internal force control. The recently developed move/squeeze orthogonal subspace decomposition coupled with the energy Lyapunov function formulation provides a basic analytical framework within which motion and force control are considered as independent problems. Simulation results of two three-link planar arms are included to demonstrate good transient performance for both motion and force that can be attained with the full dynamics control paradigm.

Wen, John T.↗
NASA NTRS2001

System Identification and POD Method Applied to Unsteady Aerodynamics

The representation of unsteady aerodynamic flow fields in terms of global aerodynamic modes has proven to be a useful method for reducing the size of the aerodynamic model over those representations that use local variables at discrete grid points in the flow field. Eigenmodes and Proper Orthogonal Decomposition (POD) modes have been used for this purpose with good effect. This suggests that system identification models may also be used to represent the aerodynamic flow field. Implicit in the use of a systems identification technique is the notion that a relative small state space model can be useful in describing a dynamical system. The POD model is first used to show that indeed a reduced order model can be obtained from a much larger numerical aerodynamical model (the vortex lattice method is used for illustrative purposes) and the results from the POD and the system identification methods are then compared. For the example considered, the two methods are shown to give comparable results in terms of accuracy and reduced model size. The advantages and limitations of each approach are briefly discussed. Both appear promising and complementary in their characteristics.

Tang, Deman↗
NASA NTRS2008

Mechanisms of Interannual Variations of the Meridional Overturning Circulation of the North Atlantic Ocean

The authors investigate the nature of the interannual variability of the meridional overturning circulation (MOC) of the North Atlantic Ocean using an Estimating the Circulation and Climate of the Ocean (ECCO) assimilation product for the period of 1993-2003. The time series of the first empirical orthogonal function of the MOC is found to be correlated with the North Atlantic Oscillation (NAO) index, while the associated circulation anomalies correspond to cells extending over the full ocean depth. Model sensitivity experiments suggest that the wind is responsible for most of this interannual variability, at least south of 40(deg)N. A dynamical decomposition of the meridional streamfunction allows a further look into the mechanisms. In particular, the contributions associated with 1) the Ekman flow and its depth-independent compensation, 2) the vertical shear flow, and 3) the barotropic gyre flowing over zonally varying topography are examined. Ekman processes are found to dominate the shorter time scales (1.5-3 yr), while for longer time scales (3-10 yr) the MOC variations associated with vertical shear flow are of greater importance. The latter is primarily caused by heaving of the pycnocline in the western subtropics associated with the stronger wind forcing. Finally, how these changes in the MOC affect the meridional heat transport (MHT) is examined. It is found that overall, Ekman processes explain a larger part of interannual variability (3-10 yr) for MHT (57%) than for the MOC (33%).

meridional heat transport↗
DOE OSTI2023

Tutorial for Non-Orthogonal Configuration Interaction Fragment Studies Using GronOR

A brief tutorial is presented for the GronOR non-orthogonal configuration interaction application illustrating the download and build process, setting up calculations with OpenMolcas and GronOR, and a description of the file structure and input options. An overview of the computational workflow includes the generation of fragment wave functions and integrals with an adapted version of OpenMolcas, running GronOR calculations, and accounting for dynamic correlation corrections. Three end-to-end use cases are provided, representing a molecular dimer, an ion pair, and overlapping molecular fragments.

Straatsma, T.P.↗
DOE OSTI2022

Do ambient shear and thermal stratification impact wind turbine tip-vortex breakdown?

Modern wind turbines experience uneven inflow conditions across the rotor, due to the ambient flow’s shear and thermal stratification. Such conditions alter the shape and length of turbine wakes and thus impact the loads and power generation of downstream turbines. To this end, understanding the spatial evolution of the individual wakes under different atmospheric conditions is key to controlling and optimising turbine arrays. With this numerical study we aim to obtain a better understanding of the fundamental physics governing the near-wake dynamics of wind turbines under shear and thermal stability, by examining their tip-vortex breakup mechanisms. Our approach considers scale-resolving simulations of a single turbine wake under a linear shear profile as well as the application of harmonic tip perturbations to trigger flow instabilities. For the subsequent analysis we use the proper orthogonal decomposition (POD) method to extract coherent structures from the flow, and we also calculate mean kinetic energy fluxes to quantify each coherent structure’s contribution to wake recovery. The wake’s helical spiral is found to hinder wake recovery for all studied ambient flow conditions, whereas the mutual inductance instability has positive MKE flux leading to an enhanced wake recovery. Finally, the ambient shear has the largest impact on the local MKE flux with respect to downstream location by changing the shape of the curve and location of extrema, whereas thermal stratification has only a minimal impact on the magnitude of the near-wake local MKE flux distribution.

17 WIND ENERGY↗
NASA NTRS1994

Correlation of finite element free vibration predictions using random vibration test data

Finite element analysis is regularly used during the engineering cycle of mechanical systems to predict the response to static, thermal, and dynamic loads. The finite element model (FEM) used to represent the system is often correlated with physical test results to determine the validity of analytical results provided. Results from dynamic testing provide one means for performing this correlation. One of the most common methods of measuring accuracy is by classical modal testing, whereby vibratory mode shapes are compared to mode shapes provided by finite element analysis. The degree of correlation between the test and analytical mode shapes can be shown mathematically using the cross orthogonality check. A great deal of time and effort can be exhausted in generating the set of test acquired mode shapes needed for the cross orthogonality check. In most situations response data from vibration tests are digitally processed to generate the mode shapes from a combination of modal parameters, forcing functions, and recorded response data. An alternate method is proposed in which the same correlation of analytical and test acquired mode shapes can be achieved without conducting the modal survey. Instead a procedure is detailed in which a minimum of test information, specifically the acceleration response data from a random vibration test, is used to generate a set of equivalent local accelerations to be applied to the reduced analytical model at discrete points corresponding to the test measurement locations. The static solution of the analytical model then produces a set of deformations that once normalized can be used to represent the test acquired mode shapes in the cross orthogonality relation. The method proposed has been shown to provide accurate results for both a simple analytical model as well as a complex space flight structure.

Chambers, Jeffrey A.↗
NASA NTRS2020

Analysis of Sample Acquisition Dynamics Using Discrete Element Method

The analysis presented in this paper is conducted in the framework of the Ocean Worlds Autonomy Testbed for Exploration Research and Simulation (OceanWATERS) project, currently under development at NASA Ames Research Center. OceanWATERS aims at designing a simulation environment which allows for testing autonomy of scientific lander missions to the icy moons of our solar system. Mainly focused on reproducing the end effector interaction with the inherent terrain, this paper introduces a novel discrete element method (DEM)-based approach to determine forces and torques acting on the lander’s scoop during the sample acquisition process. An accurate force feedback from the terrain on the scoop is required by fault-detection and autonomous decision-making algorithms to identify when the requested torque on the robotic arm’s joints exceeds the maximum available torque. Knowledge of the terrain force feedback significantly helps evaluating the arm’s links structural properties and properly selecting actuators for the joints. Models available in literature constitute a partial representation of the dynamics of the interaction. As an example, Balovnev derived an analytical expression of the vertical and horizontal force acting on a bucket while collecting a sample as a function of its geometry and velocity, soil parameters and reached depth. Although the model represents an adequate approximation of the two force components, it ignores the direction orthogonal to the scoop motion and neglects the torque. This work relies on DEM analysis to compensate for analytical models’ deficiencies and inaccuracies, i. e. provide force and torque 3D vectors, defined in the moving reference (body) frame attached to the scoop, at each instant of the sample collection process. Results from the first presented analysis relate to the specific OceanWATERS sampling strategy, which consists of collecting the sample through five consecutive passes with increasing depth, each pass following the same circularlinear- circular trajectory. Data is collected given a specific scoop design interacting with two types of bulk materials, which may characterize the surface of icy planetary bodies: snow and ice. Although specifically concerned with the OceanWATERS design, this first analysis provides the expected force trends for similar sampling strategies and allows to deduce phenomenological information about the general scooping process. In order to further instruct the community on the use of DEM tools as a solution to the sampling collection problem, two more analyses have been carried out, mainly focused on reducing the DEM computation time, which increases with a decrease in particle size. After running a set of identical simulations, where the only changing parameter is the size of the spherical particle, it is observed that the resulting force trajectories, starting from a given particle size, converge to the true trend. It is deducible that a further decrease in size yields negligible improvements in the accuracy, while it sensibly increases computation time. A final analysis aims at discussing limitations of approximating bulk material particles having a complex shape, e. g. ice fragments, with spheres, by comparing force trends resulting in the two cases for the same simulation scenario.

Catanoso, Damiana↗
DOE OSTI2025

Gauge-fixing quantum density operators at scale

We provide a theory, algorithms, and simulations of nonequilibrium quantum systems using a one-dimensional (1D) completely positive (CP), matrix-product (MP) density-operator (𝜌) representation. By generalizing the matrix product state's orthogonality center, to additionally store positive classical mixture correlations, the MP⁢𝜌 factorization naturally emerges. In this setting, we analytically and numerically examine the virtual gauge freedoms associated with the representation of quantum density operators. Based on this perspective, we simplify algorithms in certain limits to speed up the integration of the canonical-form master-equation dynamics. This enables us to quickly evolve under the dynamics of two-body quantum channels without resorting to optimization-based methods. In addition to this technical advance, we also scale up numerical examples and discuss implications for accurately modeling hardware architectures and predicting their performance in the near term. This includes an example of the quantum to classical transition of informationally leaky, i.e., decohering, qubits. In this setting, because of loss from environmental interactions, nonlocal complex coherence correlations are converted into global incoherent classical statistical mixture correlations. Lastly, the representation of both global and local correlations is discussed. We expect this work to have applications in additional nonequilibrium settings, beyond qubit engineering.

Gangapuram, Amit Jamadagni [Oak Ridge National Lab↗
NASA NTRSDate not supplied

Modal Test and Model Correlation of NASA Plum Brook Station Mechanical Vibration Facility Head Expander –Lessons Learned from the Perspective of an Early-Career Engineer

In preparation for the Sierra Nevada Corporation’s (SNC) Dream Chaser spacecraft vibration test campaign at the Mechanical Vibration Facility (MVF) at NASA Plum Brook Station (PBS) in Sandusky, Ohio, a test-verified model of MVF is needed in order to be able to perform accurate pretest analysis used for determining response limits and abort levels. MVF was designed to vibration test MPCV Orion and was used to perform the system level vibration test of the European Service Module Structural Test Article (E-STA) in 2016. MVF is comprised of an 18 ft diameter annulus table that is driven with sixteen hydraulic vertical actuator assemblies and four hydraulic horizontal actuator assemblies, which allow it to perform single axis vibration testing in the vertical axis and in each of the two orthogonal horizontal axes without the need for reconfiguring the test article. A head expander for the MVF Table has been designed and built that fills in the center opening providing a continuous flat mounting surface with a maximum diameter of 16.25feet that expands the vibration testing capabilities of MVF. The MVF Table with this head expander will be used during the SNC Dream Chaser spacecraft vibration test campaign. Therefore, a critical element in a test-verified model of the MVF will be a test correlated finite element model (FEM) of the head expander. To obtain this, engineers from the Structural Dynamics Lab (SDL) at NASA Glenn Research Center (GRC) in Cleveland, Ohio performed a modal pretest analysis, conducted a modal test in July 2019, and most recently correlated the head expander finite element model to the modal test data up to 300 Hz. From the initial test preparations to the final delivery of a correlated finite element model, all efforts mentioned were led by the same early-career engineers at NASA GRC. From the viewpoint of an early-career engineer, lessons learned about modal pretest analysis, modal testing, and finite element model correlation of the MVF Table expander head will be presented and discussed. This will include the importance of understanding the limitations of using uncorrelated finite element models in the modal pretest analysis and planning, the importance of orthogonality metrics in judging adequacy and accuracy of test mode shapes, and the importance of having the FEM match the as built hardware in the model correlation effort.

Emma L Pierson↗
DOE OSTI2022

ORNL_AISD_NiNb

This dataset describes the nickel-niobium solid solution binary alloy, where the two constituent elements nickel (Ni) and niobium (Nb) are randomly placed on an underlying crystal lattice. This dataset for nickel-niobium (Ni-Nb) alloys available includes the formation energy and bulk modulus for each crystal structure. Each atomic sample has a disordered phase which is obtained starting from an initial regular crystal structure of type body-centered cubic (BCC), face-centered cubic (FCC), or hexagonal compact packed (HCP). The geometry optimization ensures that all the alloy samples reached the equilibrium with negative formation energy. We perform geometry optimizations using the LAMMPS simulation package [1], a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. We utilized the embedded atom model (EAM) potential for Ni and Nb developed in a previous study [2]. The potential could describe behaviors of the liquid and solid phases of Ni-Nb alloy. The structural factors and angular distributions of three atoms are well-matched with X-ray and ab initio-based molecular dynamics data. We prepared the three different crystals with different initial lattice parameters (3.52 Ã… for FCC, 3.32 Ã… for BCC, and 3.5 Ã… for HCP). We performed energy minimization in two steps. Firstly, we minimized the structures with an isotropic unit cell to minimize the side effects from our arbitrary lattice parameters for all other compositions. Then, we applied geometry optimization with a triclinic (non-orthogonal) unit cell to fully minimize the stress components to calculate the elastic constants. In this procedure, we chose 10,000 as the maximum number of allowable steps aimed at obtaining fully relaxed atomic geometries. The dataset consists of three sets of crystal structures. The first set contains 46,086 irregular crystal structures, each of them with 54 atoms, obtained through optimization starting from a regular BCC crystal structure. The second set contains 24,543 irregular crystal structures, each of them with 32 atoms, obtained through optimization starting from a regular FCC crystal structure. The third set contains 39,303 irregular crystal structures, each of them with 48 atoms, obtained through optimization starting from a regular HCP crystal structure. The atomic configurations within each set span the possible compositional range. The three sets have been unified in a global dataset, which is extremely heterogeneous in terms of crystal structures, lattice volumes, and atomic configurations. Organization of files inside the dataset: the dataset contains three subdirectories called • BCC_opt • FCC_opt • HCP_opt based on the type of initial regular structure used to start the geometry optimization. Inside each of these folders, every atomic structure is identified by a string “A_B_Câ€, where A denotes the number of Nb in the system, B denotes index of structure with a given Nb number, and C denotes the total number of structures generated with a given Nb number. For each optimized crystal structure identified by the unique string of characters “A_B_Câ€, three files are provided: • A_B_C_opt.xyz: The optimized geometries in xyz format • A_B_C_opt.cfg: The optimized geometries in cfg format. It includes cell information and atomic energy, and forces calculated from LAMMPS. • A_B_C.elastic: Raw data of 21 elastic constants from LAMMPS output. • A_B_C.bulk: Calculated upper and lower bounds of bulk modulus and averaged one based on Voigt-Reuss-Hill approach from *.elastic. References: [1] A. P. Thompson, H. M. Aktulga, R. Berger, D. S. Bolintineanu, W. M. Brown, P. S. Crozier, P. J. in 't Veld, A. Kohlmeyer, S. G. Moore, T. D. Nguyen, R. Shan, M. J. Stevens, J. Tranchida, C. Trott, and S. J. Plimpton. LAMMPS - a flexible simulation tool for particle-based materials modeling at the atomic, meso, and continuum scales. Comp. Phys. Comm., 271:108171, 2022. [2] Y Zhang, R Ashcraft, MI Mendelev, CZ Wang, and KF Kelton. Experimental and molecular dynamics simulation study of structure of liquid and amorphous ni62nb38 alloy. The Journal of chemical physics, 145(20):204505, 2016.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Properties of giant dipole resonances within an extended pairing model with a focus on spectral statistics

In this paper, we report on a study of the spectral features associated with dipole resonances in medium mass nuclei (E ≤ 12 MeV), as revealed in the framework of the spd-interacting boson model. The effect of pairing correlations on the theory follows from solutions obtained through an application of the Bethe Ansatz Equation. In general, calculated spectra around the critical point of the vibrational to γ-soft transitions appears to approach that of a Gaussian Orthogonal Ensemble, while near the rotational and vibrational limits of the theory the spectra show more regular behavior. Specifically, the results reveal that the statistical features of the spectra are sensitive to the vector boson pairing strength, c p , in the transition region; that is, when c p is zero, or when the system approaches one of its dynamical symmetries limits the spectrum display regular features, while for stronger c p values, or when near to the critical phase transition region, the spectral feature show more chaotic behavior. Overall, our results indicate that the statistical features are governed by the interplay between dipole resonant energies, pairing correlations, and interactions between and among the single and vector bosons modes of the theory. As part of this work we also found out that chaoticity occurs when results were fit to a Berry-Robnik distribution. Throughout our analyses, we used experimentally known information about both positive and negative parity states. Our findings suggest that dipole resonances appear to be best-described by Poisson statistics for A ≈ 32-138 nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1996

Dynamics and Control of Articulated Anisotropic Timoshenko Beams

The paper illustrates the use of continuum models in control design for stabilizing flexible structures. A 6-DOF anisotropic Timoshenko beam with discrete nodes where lumped masses or actuators are located provides a sufficiently rich model to be of interest for mathematical theory as well as practical application. We develop concepts and tools to help answer engineering questions without having to resort to ad hoc heuristic ("physical") arguments or faith. In this sense the paper is more mathematically oriented than engineering papers and vice versa at the same time. For instance we make precise time-domain solutions using the theory of semigroups of operators rather than formal "inverse Laplace transforms." We show that the modes arise as eigenvalues of the generator of the semigroup, which are then related to the eigenvalues of the stiffness operator. With the feedback control, the modes are no longer orthogonal and the question naturally arises as to whether there is still a modal expansion. Here we prove that the eigenfunctions yield a biorthogonal Riesz basis and indicate the corresponding expansion. We prove mathematically that the number of eigenvalues is nonfinite, based on the theory of zeros of entire functions. We make precise the notion of asymptotic modes and indicate how to calculate them. Although limited by space, we do consider the root locus problem and show for instance that the damping at first increases as the control gain increases but starts to decrease at a critical value, and goes to zero as the gain increases without bound. The undamped oscillatory modes remain oscillatory and the rigid-body modes go over into deadbeat modes. The Timoshenko model dynamics are translated into a canonical wave equation in a Hilbert space. The solution is shown to require the use of an "energy" norm which is no more than the total energy: potential plus kinetic. We show that, under an appropriate extension of the notion of controllability, rate feedback with a collocated sensor can stabilize the structure in the sense that all modes are damped and the energy decays to zero. An example, non-numeric, is worked out in some detail illustrating the concepts and theory developed.

Balakrishnan, A. V.↗
DOE OSTI2021

Improvement of Numerical Methods in Pronghorn

During the fiscal year of 2021, the finite-volume method (FVM) was deployed in Pronghorn to provide improved efficiency, stability, and accuracy for coarse-mesh, thermal-hydraulics problems. While the main goals of fiscal year 2021 were met, several issues emerged from the early deployment of the finite-volume method in Pronghorn. These issues were: The compressible and incompressible formulations that were implemented are inadequate for many nuclear reactor flow problems; Omission of terms accounting for the porosity and Darcy-Forchheimer body force discontinuities in the Rhie-Chow interpolation lead to oscillations in pressure and velocity at these discontinuities; The FVM lacks a correction for non-orthogonal grids for computing accurate pressure gradient leading to loss of accuracy in regions with skewed elements; The FVM currently uses a monolithic solver. Monolithic solvers have issues dealing with the saddle-point nature of the discretized fluid equations. The result are bad convergence if direct factorization is not used and large memory consumption when direct factorization is used. These four issues are addressed in this report. In particular, we report the completion of the following task: Implementation of a weakly compressible formulation in the MOOSE Navier-Stokes module; Implementation of Moukalled’s method for including body forces in the RCI. Additionally, we identified the need to smooth the porosity using Moukalled’s face-cell smoothing operator; Implementation of a non-orthogonal correction for the Green-Gauss gradient computation; Preliminary implementation of a SIMPLE segregated solver.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1982

Tropical cyclone track and genesis forecasting using satellite microwave sounder data

Although many dynamical and statistical prediction schemes are available to forecasters, tropical cyclone track errors are still large. One primary difficulty is that tropical cyclones exist over the data-sparse tropical oceans. Satellite sounders, however, routinely provide numerous data over these areas. Mean layer temperatures from the Scanning Microwave Spectrometer on board the Nimbus 6 satellite are decomposed using empirical orthogonal functions, and the expansion coefficients are related to deviations from the persistence forecast location, to speed change, to direction change and to intensity change. The significance of the regression equations is tested by a null hypothesis of zero correlation coefficient. It appears that significant information about tropical cyclone motion exists in the satellite-estimated mean layer temperatures, especially at upper levels. A physical interpretation of the statistical results is offered, and a one-storm-out independent test is used to test the stability of the equations. Finally, some further work is suggested.

Kidder, S. Q.↗
NASA NTRS1974

An ultrahigh-accuracy body pointing system for the Large Space Telescope

The Large Space Telescope (LST) program is aimed at placing a three-meter diffraction-limited telescope in a 270-nm orbit to perform astronomical observations that are not possible with earth-based telescopes. A complex simulation model is described which was developed to determine whether the stringent pointing stability requirements could be met. The model (programmed on a hybrid computer) included detailed dynamic representation of control moment gyros (CMGs) and reaction wheels (RWs), including their noise characteristics; dynamic sensor representation (including noise); shockmounts for the CMG actuators; detailed representation of an image motion compensation system; and a detailed flexible body vehicle model. Stability and performance studies based on the simulation model showed that the body pointing system will meet LST requirements in the presence of CMG vibrational disturbances and sensor noise. The recommended system consists of three orthogonally mounted RWs for primary short-term control, and a cluster of CMG actuators for continuous RW desaturation and vehicle maneuvering.

Rybak, S. C.↗
NASA NTRS1977

Vibration analysis of rotor blades with an attached concentrated mass

The effect of an attached concentrated mass on the dynamics of helicopter rotor blades is determined. The point transmission matrix method was used to define, through three completely automated computer programs, the natural vibrational characteristics (natural frequencies and mode shapes) of rotor blades. The problems of coupled flapwise bending, chordwise bending, and torsional vibration of a twisted nonuniform blade and its special subcase pure torsional vibration are discussed. The orthogonality relations that exist between the natural modes of rotor blades with an attached concentrated mass are derived. The effect of pitch, rotation, and point mass parameters on the collective, cyclic, scissor, and pure torsional modes of a seesaw rotor blade is determined.

Murthy, V. R.↗
NASA NTRS2014

Statistical Mechanics of Turbulent Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V.↗
NASA NTRS2014

MHD Turbulence and Magnetic Dynamos

Incompressible magnetohydrodynamic (MHD) turbulence and magnetic dynamos, which occur in magnetofluids with large fluid and magnetic Reynolds numbers, will be discussed. When Reynolds numbers are large and energy decays slowly, the distribution of energy with respect to length scale becomes quasi-stationary and MHD turbulence can be described statistically. In the limit of infinite Reynolds numbers, viscosity and resistivity become zero and if these values are used in the MHD equations ab initio, a model system called ideal MHD turbulence results. This model system is typically confined in simple geometries with some form of homogeneous boundary conditions, allowing for velocity and magnetic field to be represented by orthogonal function expansions. One advantage to this is that the coefficients of the expansions form a set of nonlinearly interacting variables whose behavior can be described by equilibrium statistical mechanics, i.e., by a canonical ensemble theory based on the global invariants (energy, cross helicity and magnetic helicity) of ideal MHD turbulence. Another advantage is that truncated expansions provide a finite dynamical system whose time evolution can be numerically simulated to test the predictions of the associated statistical mechanics. If ensemble predictions are the same as time averages, then the system is said to be ergodic; if not, the system is nonergodic. Although it had been implicitly assumed in the early days of ideal MHD statistical theory development that these finite dynamical systems were ergodic, numerical simulations provided sufficient evidence that they were, in fact, nonergodic. Specifically, while canonical ensemble theory predicted that expansion coefficients would be (i) zero-mean random variables with (ii) energy that decreased with length scale, it was found that although (ii) was correct, (i) was not and the expected ergodicity was broken. The exact cause of this broken ergodicity was explained, after much investigation, by greatly extending the statistical theory of ideal MHD turbulence. The mathematical details of broken ergodicity, in fact, give a quantitative explanation of how coherent structure, dynamic alignment and force-free states appear in turbulent magnetofluids. The relevance of these ideal results to real MHD turbulence occurs because broken ergodicity is most manifest in the ideal case at the largest length scales and it is in these largest scales that a real magnetofluid has the least dissipation, i.e., most closely approaches the behavior of an ideal magnetofluid. Furthermore, the effects grow stronger when cross and magnetic helicities grow large with respect to energy, and this is exactly what occurs with time in a real magnetofluid, where it is called selective decay. The relevance of these results found in ideal MHD turbulence theory to the real world is that they provide at least a qualitative explanation of why confined turbulent magnetofluids, such as the liquid iron that fills the Earth's outer core, produce stationary, large-scale magnetic fields, i.e., the geomagnetic field. These results should also apply to other planets as well as to plasma confinement devices on Earth and in space, and the effects should be manifest if Reynolds numbers are high enough and there is enough time for stationarity to occur, at least approximately. In the presentation, details will be given for both theoretical and numerical results, and references will be provided.

Shebalin, John V↗