Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 1,207 records · Page 67
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Quantum dynamics in strongly driven random dipolar magnets

The random dipolar magnet LiHo x Y 1– x F 4 enters a strongly frustrated regime for small Ho 3+ concentrations with x < 0.05. In this regime, the magnetic moments of the Ho 3+ ions experience small quantum corrections to the common Ising approximation of LiHo x Y 1– x F 4 , which lead to a $\mathbb{Z}_2$-symmetry breaking and small, degeneracy breaking energy shifts between different eigenstates. Here we show that destructive interference between two almost degenerate excitation pathways burns spectral holes in the magnetic susceptibility of strongly driven magnetic moments in LiHo x Y 1– x F 4 . Such spectral holes in the susceptibility, microscopically described in terms of Fano resonances, can already occur in setups of only two or three frustrated moments, for which the driven level scheme has the paradigmatic Λ shape. For larger clusters of magnetic moments, the corresponding level schemes separate into almost isolated many-body Λ schemes, in the sense that either the transition matrix elements between them are negligibly small or the energy difference of the transitions is strongly off-resonant to the drive. This enables the observation of Fano resonances, caused by many-body quantum corrections to the common Ising approximation also in the thermodynamic limit. Finally, we discuss its dependence on the driving strength and frequency as well as the crucial role that is played by lattice dissipation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS1992

The evolution equation for the flame surface density in turbulent premixed combustion

One central ingredient in flamelet models for turbulent premixed combustion is the flame surface density. This quantity conveys most of the effects of the turbulence on the rate of energy release and is obtained via a modeled transport equation, called the Sigma-equation. Past theoretical work has produced a rigorous approach that leads to an exact, but unclosed, formulation for the turbulent Sigma-equation. In this exact Sigma-equation, it appears that the dynamical properties of the flame surface density are determined by a single parameter, namely the turbulent flame stretch. Unfortunately, the flame surface density and the turbulent flame stretch are not available from experiments and, in the absence of experimental data, little is known on the validity of the closure assumptions used in current flamelet models. Direct Numerical Simulation (DNS) is the obvious, complementary approach to get basic information on these fundamental quantities. Three-dimensional DNS of premixed flames in isotropic turbulent flow is used to estimate the different terms appearing in the Sigma-equation. A new methodology is proposed to provide the source and sink terms for the flame surface density, resolved both temporally and spatially throughout the turbulent flame brush. Using this methodology, the effects of the Lewis number on the rate of production of flame surface area are described in great detail and meaningful comparisons with flamelet models can be performed. The analysis reveals in particular the tendency of the models to overpredict flame surface dissipation as well as their inability to reproduce variations due to thermo-diffusive phenomena. Thanks to the detailed information produced by a DNS-based analysis, this type of comparison not only underscores the shortcomings of current models but also suggests ways to improve them.

Trouve, A.↗
NASA NTRS1999

Emergence of postural patterns as a function of vision and translation frequency

Emergence of postural patterns as a function of vision and translation frequency. We examined the frequency characteristics of human postural coordination and the role of visual information in this coordination. Eight healthy adults maintained balance in stance during sinusoidal support surface translations (12 cm peak to peak) in the anterior-posterior direction at six different frequencies. Changes in kinematic and dynamic measures revealed that both sensory and biomechanical constraints limit postural coordination patterns as a function of translation frequency. At slow frequencies (0.1 and 0.25 Hz), subjects ride the platform (with the eyes open or closed). For fast frequencies (1.0 and 1.25 Hz) with the eyes open, subjects fix their head and upper trunk in space. With the eyes closed, large-amplitude, slow-sway motion of the head and trunk occurred for fast frequencies above 0.5 Hz. Visual information stabilized posture by reducing the variability of the head's position in space and the position of the center of mass (CoM) within the support surface defined by the feet for all but the slowest translation frequencies. When subjects rode the platform, there was little oscillatory joint motion, with muscle activity limited mostly to the ankles. To support the head fixed in space and slow-sway postural patterns, subjects produced stable interjoint hip and ankle joint coordination patterns. This increase in joint motion of the lower body dissipated the energy input by fast translation frequencies and facilitated the control of upper body motion. CoM amplitude decreased with increasing translation frequency, whereas the center of pressure amplitude increased with increasing translation frequency. Our results suggest that visual information was important to maintaining a fixed position of the head and trunk in space, whereas proprioceptive information was sufficient to produce stable coordinative patterns between the support surface and legs. The CNS organizes postural patterns in this balance task as a function of available sensory information, biomechanical constraints, and translation frequency.

NASA Discipline Neuroscience↗
DOE OSTI2022

Simulation of 15% and 50% Thermal Power Dispatch to an Industrial Facility Using a Flexible Generic Full-Scope Pressurized Water Reactor Plant Simulator

Nuclear power plants in the United States are increasingly challenged to compete in wholesale electricity markets due to the low electricity costs and increasingly dynamic grid conditions from competing generation sources. An alternative market for nuclear power is industrial facilities that can use the thermal and/or electrical power generated by a nuclear power plant to offset the economic losses incurred by electricity market challenges. A generic pressurized water reactor (PWR) simulator was used to show the results of a basic design for a generic thermal power extraction system and tests were run using a set of procedures to show what happens when a nuclear power plant transitions from full electrical power dispatch to 15% and 50% thermal power dispatch. This type of operation leads to losses in turbine performance efficiency due to the deviation from the design operating point, but because the thermal power is also used by the industry load without conversion losses, the combined thermal efficiency of the PWR increases. For the 15% case, the thermal efficiency increased from 32% to 41.9%, while for the 50% case, the efficiency increased up to 60.1%. In addition, for 50% thermal power dispatch, the power dissipated by the condenser decreased from approximately 2000 to approximately 1300 MW (thermal), indicating a substantially diminished impact on the environment in terms of releasing heat into the cooling water reservoir.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2025

Distinctive features of structural evolution and thermodynamic response in wide-bandgap semiconductors driven by intense electronic excitation

Radiation-tolerant material selection requires balancing lattice rigidity, defect dynamics, and electronic stability, as shown by covalent SiC outperforming ionic Ga 2 O 3 and GaN under extreme environments. Responding to intense electronic excitation, irradiation-driven phase segregation (β → δ/κ in Ga 2 O 3 ) and core–shell track (disordered structure in GaN), accompanied by elemental redistribution, contrastingly, exceptional radiation tolerance manifested by comparatively minimal lattice distortion (0.17 % strain variation) was demonstrated in SiC. These differential responses are primarily attributed to two fundamental mechanisms: (i) thermodynamic driving forces governing defect migration and phase separation, and (ii) the synergistic effects of robust covalent bonding composition coupled with efficient defect recombination processes. Here, the stronger electron–phonon (e-ph) coupling in Ga 2 O 3 (4.34 × 1018 W m −3 K −1 ) and GaN (3.55 × 10 18 W m −3 K −1 ) enhances lattice energy deposition, triggering thermal spikes (ΔT ≫ T m ) and structural transition behaviors, whereas weaker e-ph coupling in SiC (3.69 × 10 18 W m −3 K −1 ), relatively high thermodynamic parameters and efficient energy dissipation suppress thermal spikes to maintaining lattice integrity. The photoresponse degradation driven by enhanced radiative recombination is dominant in N-doped SiC, while V-doped systems achieve defect-mediated photoconduction optimization characterized by abrupt current transitions, matching fluorescence yield evolutions, and directly connecting defect engineering to optoelectronic performance.

Intense electronic excitation↗
DOE OSTI2024

Vortex bursting and associated twist dynamics on helical vortex tubes and vortex rings

The interaction of opposite-signed twist waves on vortex tubes can lead to vortex bursting, a process where the core expands into a double ring-like structure with strong swirling flows. Previous works have studied vortex bursting on rectilinear vortices by axially perturbing the initial core size to generate the twist waves, and observed largely axisymmetric bursting dynamics. In this work, we numerically study bursting on vortical structures with curved centrelines, analysing the interaction between the centreline dynamics, twist wave generation and propagation, and vortex bursting. We focus on axially perturbed helical vortex tubes with small radius-to-pitch ratios up to$0.0625$, as well as vortex rings with a large radius-to-core size ratio$10$, both at a circulation-based Reynolds number$5000$. The results show that though the initial twist wave propagation speeds are relatively unaffected by the curvature and torsion of the centreline, the bursting process is altered significantly compared with rectilinear vortices. The self-induced rotation of the centreline of the helical tube induces a non-axisymmetric distortion of the bursting structure, which rapidly breaks up the vortex core into small-scale helical structures. A similar destabilization of the bursting structure also occurs on vortex rings. The enstrophy increase and accelerated energy decay associated with bursting are predominantly determined by the twist wave strength, rather than the curvature and torsion of the centreline. Combined, our findings imply that bursting could play an important role in transferring and dissipating energy of vortical structures in wakes, and turbulent flows in general.

Mechanics↗
DOE OSTI2023

Simulation of radioactive plume transport in the atmosphere including dynamics of particle aggregation and breakup

Accurate prediction of the atmospheric transport of debris particles relies heavily on our knowledge of the size distribution of the particles within a debris cloud. Assuming a fixed particle size during simulations is not always viable since the size distribution of the debris can change during transport. Various microphysical processes, such as aggregation and breakup, influence debris particles and dictate any changes to the size distribution. To track those changes that can occur, a population balance model can be adopted and instituted within a model framework. Nonetheless, many of the models that simulate the transport of radioactive debris following a device-driven fission incident have historically neglected to consider these processes. As such, this work describes our effort to develop a modeling framework capable of simulating the transport and deposition of a radioactive plume generated from a fission incident with a dynamic population balance including particle aggregation and breakup. Herein the impact of aggregation and breakup, individually and collectively, on the particle size distribution is explored using the developed framework. When simulating aggregation, for example, six mechanisms, including Brownian coagulation, the convective enhancement to Brownian coagulation, van der Waals-viscous force correction for Brownian coagulation, gravitational collection, turbulent inertial motion, and turbulent shear, are considered. Brownian coagulation and its corrections have, as one would expect, a large impact on relatively small aggregates. Aggregates with a diameter that is less than or equal to 1.0 μm, for instance, comprise 50.6 vol % of all aggregates in the absence of aggregation and 31.2 vol % when Brownian coagulation and its corrections are considered. Gravitational collection and, to a much lesser extent, turbulent shear and turbulent inertial motion are, conversely, of great importance to relatively large aggregates (i.e., diameter greater than 3.0 μm). Additionally, the individual effects of atmospheric and particle parameters, such as wind speed and particle density, are examined. Of the parameters examined, turbulent energy dissipation and aggregate fractal dimension (i.e., aggregate shape with lower values representing more irregular particles) were of substantial importance since both terms directly impact aggregate stability and, by extension, the breakup rate. Large-scale transport and deposition simulations in a dry atmosphere are also presented and discussed as a proof of concept.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
NASA NTRS2019

Consequences of Tidal Dissipation in a Putative Venusian Ocean

The solar tide in an ancient Venusian ocean is simulated using a dedicated numerical tidal model. Simulations with varying ocean depth and rotational periods ranging from minus 243 to 64 sidereal Earth days are used to calculate the tidal dissipation rates and associated tidal torque. The results show that the tidal dissipation could have varied by more than 5 orders of magnitude, from 0.001 to 780 gigawatts (GW), depending on rotational period and ocean depth. The associated tidal torque is about 2 orders of magnitude below the present day Venusian atmospheric torque, and could change the Venusian daylength by up to 72 days per million years depending on rotation rate. Consequently, an ocean tide on ancient Venus could have had significant effects on the rotational history of the planet. These calculations have implications for the rotational periods of similarly close-in exoplanetary worlds and the location of the inner edge of the liquid water habitable zone.

Dynamical evolution and stability↗
DOE OSTI2023

Spark Channel Dynamics of Electrostatic Discharges

When two differently-charged objects are brought in close proximity to each other, the resulting high electric fields can cause electron avalanche breakdown of the air gap separating the objects, a process known as electrostatic discharge (ESD). If enough initial charge is stored on the objects, the electrical breakdown can proceed to ionize the air to such a degree that a highly conductive filament of plasma forms in the gap, known as a spark channel. The spark electrically bridges the air gap, resulting in a rapid pulse of current that neutralizes the charge difference. The current pulse produces significant heating of the gas in the spark, resulting in dissociation, ionization, thermal radiation, and hydrodynamic expansion. ESD presents a hazard to electrically-sensitive devices, with consequences such as economic losses (e.g. damaged electronics) or unsafe response (e.g. unintended ignition of flammable gas mixtures, initiation of detonators, etc.). For this thesis, the ESD spark is taken to occur between two conducting electrodes, with the spark channel being axisymmetric in a cylindrical coordinate system centered on the channel. An RLC-type circuit is used for the discharge model of the ESD event. The spark is treated as a time-dependent resistance that is in series with a capacitance, an inductance, and (optionally) a load resistance representing a “victim” component under threat from the ESD event. The primary motivation of this work is to use a numerical hydrodynamic model to understand the energy dissipation and transport processes in the spark. The model consists of the compressible Euler equations of mass, momentum, and energy conservation together with an Eddington/P1 approximation for thermal radiation transport. To close the hydrodynamic system, an equation of state (EOS) was fitted from tabular data for air that accounts for the dissociation and ionization of air species. The hydrodynamic equations are solved using a conservative Lagrangian finite volume method. These partial differential equations are coupled to the circuit equations by calculation of the spark resistance via numerical integration of the electrical conductivity of the channel. Computational results are compared against experimental measurements of discharge current and radial density of the spark channel.

42 ENGINEERING↗
DOE OSTI2021

Pressure Response Optimization of an Eddy Current-Driven Flyer Plate Valve for the ITER Shattered Pellet Injection System

One technique for mitigating disruptions in a tokamak is shattered pellet injection (SPI). SPI is a process in which a large solid pellet consisting of deuterium, neon, or argon is desublimated in a pipe gun barsrel and launched downstream. Pellets are shattered just before entering the plasma by an impact with an angled tube. Injection of these materials into the plasma radiates stored thermal energy, limits current decay rates, suppresses the generation of runaway electrons, and dissipates runaway electrons if necessary. A critical element of the SPI system is a fast-acting valve that releases high-pressure gas to dislodge and accelerate pellets directly, or indirectly via a mechanical punch. A prototype valve sized for the ITER SPI system has been designed and fabricated. A pulsed high-voltage power supply energizes the valve’s internal magnetic coil, which induces eddy currents in the adjacent flyer plate resulting in a repulsive force between the flyer plate and the coil. The flyer plate action lifts a valve seat, allowing high-pressure gas to flow from the valve plenum to the downstream (breech) location of the pellet or mechanical punch. All of the valve’s internal components are designed to operate in ITER-level static background magnetic fields. Here, a study was conducted to optimize the downstream pressure response for a range of valve sizes and operating pressures. In particular, the study analyzes the breech pressure response associated with varying plenum pressures as well as varying breech volumes. A computational fluid dynamics simulation was built in STAR-CCM+ and validated against data from laboratory experiments. The resulting simulation outputs, in the form of downstream responses for a variety of initial plenum pressures and breech volumes, will be used as a complement to experimental data to ensure the pressure pulse is suitable for pellet survivability. These data, combined with studies on pellet shear strength and shock response, will be applied to optimization of overall operating parameters of the ITER SPI system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Space-time generalization of mutual information

The mutual information characterizes correlations between spatially separated regions of a system. Yet, in experiments we often measure dynamical correlations, which involve probing operators that are also separated in time. Here, we introduce a space-time generalization of mutual information which, by construction, satisfies several natural properties of the mutual information and at the same time characterizes correlations across subsystems that are separated in time. In particular, this quantity, that we call the space-time mutual information, bounds all dynamical correlations. We construct this quantity based on the idea of the quantum hypothesis testing. As a by-product, our definition provides a transparent interpretation in terms of an experimentally accessible setup. We draw connections with other notions in quantum information theory, such as quantum channel discrimination. Finally, we study the behavior of the space-time mutual information in several settings and contrast its long-time behavior in many-body localizing and thermalizing systems.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Clogging, diode and collective effects of skyrmions in funnel geometries

Abstract Using a particle-based model, we examine the collective dynamics of skyrmions interacting with a funnel potential under dc driving as the skyrmion density and relative strength of the Magnus and damping terms are varied. For driving in the easy direction, we find that increasing the skyrmion density reduces the average skyrmion velocity due to jamming of skyrmions near the funnel opening, while the Magnus force causes skyrmions to accumulate on one side of the funnel array. For driving in the hard direction, there is a critical skyrmion density below which the skyrmions become trapped. Above this critical value, a clogging effect appears with multiple depinning and repinning states where the skyrmions can rearrange into different clogged configurations, while at higher drives, the velocity-force curves become continuous. When skyrmions pile up near the funnel opening, the effective size of the opening is reduced and the passage of other skyrmions is blocked by the repulsive skyrmion–skyrmion interactions. We observe a strong diode effect in which the critical depinning force is higher and the velocity response is smaller for hard direction driving. As the ratio of Magnus force to dissipative term is varied, the skyrmion velocity varies in a non-linear and non-monotonic way due to the pile up of skyrmions on one side of the funnels. At high Magnus forces, the clogging effect for hard direction driving is diminished.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1998

Pulsing Inertial Oscillation, Supercell Storms, and Surface Mesonetwork Data

The pulsing inertial oscillation (PIO) model is a nonlinear, time-dependent, translating vortex solution of the inviscid, compressible fluid dynamic equations in the middle troposphere. The translation of this vortex during a pulse is strikingly similar to that of a supercell storm -- a rotating thunderstorm that can generate tornadoes and hail. Two studies were performed to test the hypothesis that some supercell storms are manifestations of a PIO pulse. The first study applied the model to an intense interior draft whose buoyancy was bounded by a temperature excess of +/- 12 K. The peak updraft speed achieved was 41.5 m/ s and the peak Rossby number was 92.9. The study also pointed to an advanced concept for attaining higher values. The second study applied the PIO model to a supercell storm as a whole and succeeded in replicating its bulk properties, such as mesocyclonic circulation, net mass and moisture influxes, and time track. This study also identified a critical feature of the PIO model that could be tested against storm data: The average vertical draft is downward before the turn in the storm track and upward afterwards. In the conventional theory, the average vertical draft is upward from storm inception until dissipation. These differing draft predictions were compared with the best available data, which are surface mesonetwork data. These data were found to support the PIO model. However, surface data alone are not conclusive, and further measurements are warranted.

Costen, R. C.↗
DOE OSTI2026

Novel application of the modal strain energy technique for state-of-the-art damping predictions

Structural damping, which measures the energy dissipation of a vibrating structure, is a key modeling input for lightweight structures but is notoriously hard to predict. This work utilizes vibration-based measurements of centimeter-scale coupons and the modal strain energy approach to predict structural damping of a lightly damped structure. The approach was originally validated with panels shorter than a meter in length. This work extends the validation to a 2.75 m beam made of unidirectional and biaxial glass fiber laminates bonded by adhesive. The comparison between three-dimensional finite element model predictions and full-scale experimental measurements of damping show an average error of 5.2% for the first five modes. Additionally, the modal strain energy approach is newly applied with a one-dimensional geometrically exact beam theory model and a two-dimensional sectional analysis solver. This beam approach accurately predicts the damping behavior of the first bending modes but loses accuracy for higher order modes that are dominated by three-dimensional effects. This novel approach provides faster simulations while allowing arbitrary beam cross sections. The paper also investigates traditional and high-force dynamic mechanical analysis to measure structural damping of coupons. Both alternatives show significant errors in attempted validation against the theoretical thermoelastic damping of aluminum coupons.

17 WIND ENERGY↗
DOE OSTI2020

Initial Atmospheric Transport of Particles

A new version of the Initial Atmospheric Transport (IAT) model has been developed. This report fully characterizes the new model. IAT is designed to model heat release into the atmosphere that contains contaminated particles. The heat release forms a buoyant puff or plume. Buoyancy is quickly dissipated in 10's to 100's of seconds through drag and entrainment of ambient air. The final location of contaminant particles in the atmosphere after dissipation of abnormal heat is an important input to long term transportation codes such as the National Oceanic and Atmospheric Administration's (NOAA) Hybrid Single Particle Lagrangian Integrated Trajectory (HYSPLIT) code. Ultimately the settling of contaminants after long term transport provides a basis for estimating health consequences. In comparison to the old version, the new model includes explicit simulation of particles, newly formulated equations of motion for a rising puff, and a new plume model. The plume model releases many instances of a single puff along a single trajectory in a way that balances energy released over a longer period than would be appropriate for a single puff. The puff model has additions for vorticity, turbulence, virtual mass, drag, and entrainment. The particle model leads to direct interfacing between IAT and the long-term transport code HYSPLIT. Previously, particle locations were given a pattern based on IAT's puff trajectory with no particle simulation. Now, particles are explicitly simulated, and 3-D spatial distributions of particles are output based on final positions. Ballistic particles are thrown outward at high velocity from explosive events and tend to hit the ground before the end of the simulation. Non-ballistic particles tend to get pulled into the rising puff s ring of vorticity with some escaping to ambient conditions as buoyant energy dissipates. A new model verification and validation effort was made for the new plume model. The Wallops Island Solid Rocket Propellant (SRP) fire plume was extracted to a 3-D transient dataset. After calibration, the model performs better than previous single puff attempts. Even so, error levels in the process used to estimate location and ambient weather conditions approach 30%. This leads to the conclusion that unplanned photogrammetry evaluations of the IAT model are useful but not sufficiently accurate to quantify model parameters. Like previous attempts, the vertical motion entrainment parameters come out much too high in comparison to values for other studies in the literature. It is suspected that a new model of shear entrainment would alleviate this problem. Unfortunately, IAT already has more calibration parameters than the data comparisons thus far can resolve with statistical certainty. This underscores the need to verify and validate the IAT model using detailed CFD studies that allow ideal statistical comparison with the IAT model for hundreds of cases. Even so, validation by comparison to carefully planned experiments (i.e. with acceptable error levels) is also necessary because the dynamics of entrainment are at a sub-grid level for CFD. Though better validation and verification are most important, many new enhancements are still envisioned for future IAT development. These include particle rain out due to condensation and agglomeration of water, direct connection to the North American Mesoscale (NAM) weather dataset, dynamic weather conditions, including ambient particle entrainment into the puff, and detrainment at the puff to ambient mixing interface. In general, including detailed aerosol dynamics is the greatest coming challenge which will require an entire new sub-model in IAT. This is very important if contaminant particles change in such a way that they increase or decrease in their potential harm to the environment.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1994

Field-Aligned Poynting Flux Observations in the High-Latitude Ionosphere

We have used data from Dynamics Explorer 2 to investigate the rate of conversion of electromagnetic energy into both thermal and bulk flow particle kinetic energy in the high-latitude ionosphere. The flux tube integrated conversion rate E(dot)J can be determined from spacecraft measurements of the electric and magnetic field vectors by deriving the field-aligned Poynting flux, S(parallel) = S(dot)B(sub 0), where B(sub 0) is in the direction of the geomagnetic field. Determination of the Poynting flux from satellite observations is critically dependent upon the establishment of accurate values of the fields and is especially sensitive to errors in the baseline (unperturbed) geomagnetic field. We discuss our treatment of the data in some detail, particularly in regard to systematically correcting the measured magnetic field to account for attitude changes and model deficiencies. S(parallel) can be used to identify the relative strengths of the magnetosphere and thermospheric winds as energy drivers and we present observations demonstrating the dominance of each of these. Dominance of the magnetospheric driver is indicated by S(parallel) directed into the ionosphere. Electromagnetic energy is delivered to and dissipated within the region. Dominance of the neutral wind requires that the conductivity weighted neutral wind speed in the direction of the ion drift be larger than the ion drift, resulting in observations of an upward directed Poynting flux. Electromagnetic energy is generated within the ionospheric region in this case. We also present observations of a case where the neutral atmosphere motion may be reaching a state of sustained bulk flow velocity as evidenced by very small Poynting flux in the presence of large electric fields.

Slavin, J. A.↗
NASA NTRS2007

Shell Models of RMHD Turbulence and the Heating of Solar Coronal Loops

A simplified nonlinear numerical model for the development of incompressible magnetohydrodynamics in the presence of a strong magnetic field B|| and stratification, nicknamed 'Shell-Atm,' is presented. In planes orthogonal to the mean field, the nonlinear incompressible dynamics is replaced by two-dimensional shell models for the complex variables u and b, allowing one to reach large Reynolds numbers while at the same time carrying out sufficiently long integrations to obtain good statistics at moderate computational cost. The shell models of different planes are coupled by Alfve'n waves propagating along B||. The model may be applied to open or closed magnetic field configurations where the axial field dominates and the plasma pressure is low; here we apply it to the specific case of a magnetic loop of the solar corona heated by means of turbulence driven by photospheric motions, and we use statistics for its analysis. The Alfven waves interact nonlinearly and form turbulent spectra in the directions perpendicular and, through propagation, also parallel to the mean field. A heating function is obtained and shown to be intermittent; the average heating is consistent with values required for sustaining a hot corona and is proportional to the aspect ratio of the loop to the -1.5 power, and haracteristic properties of heating events are distributed as power laws. Crosscorrelations show a delay of dissipation compared with energy content.

corona↗
DOE OSTI2020

Turbulent properties from Doppler lidar

Permanent and mobile facilities deployed by the ARM program cover a wide range of climatic regions, with a suite of active remote-sensing of cloud and aerosol profiles complemented by comprehensive in-situ measurements of the chemical, physical and optical properties of aerosol measured at the surface. To be able to fully exploit these datasets requires an understanding of the transport from the surface and within the atmosphere, especially within the planetary boundary layer. This requires a method of identifying when the surface and atmosphere are dynamically-linked. Determining the state of the boundary layer, and defining a well-mixed boundary-layer, is therefore suitably described in terms of mixing; here, we will assume that mixing essentially means turbulent mixing. This necessitates accurate measurements of the winds and turbulent mixing in the atmosphere. Doppler lidar is an excellent tool for observing turbulent properties, and ARM has deployed commercially available Doppler lidars at many of the fixed and mobile facilities. Doppler lidars have been shown to provide accurate retrievals of wind and turbulent properties within the boundary layer , however, robust characterization and calculation of the Doppler lidar measurement uncertainties is required before deriving higher level data products.

54 ENVIRONMENTAL SCIENCES↗