SEARCH · Engineering Papers

Results for “Fluid Structure Interaction”

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.

SourceResearch connectionFrom yearThrough year

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

At least 1,387 records · Page 77

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Nanoparticle localization within chiral liquid crystal defect lines and nanoparticle interactions

Self-assembly of colloidal particles into predefined structures is a promising way to design inexpensive manmade materials with advanced macroscopic properties. Doping of nematic liquid crystals (LCs) with nanoparticles has a series of advantages in addressing these grand scientific and engineering challenges. It also provides a very rich soft matter platform for the discovery of unique condensed matter phases. The LC host naturally allows the realization of diverse anisotropic interparticle interactions, enriched by the spontaneous alignment of anisotropic particles due to the boundary conditions of the LC director. Here we demonstrate theoretically and experimentally that the ability of LC media to host topological defect lines can be used as a tool to probe the behavior of individual nanoparticles as well as effective interactions between them. LC defect lines irreversibly trap nanoparticles enabling controlled particle movement along the defect line with the use of a laser tweezer. Minimization of Landau–de Gennes free energy reveals a sensitivity of the ensuing effective nanoparticle interaction to the shape of the particle, surface anchoring strength, and temperature, which determine not only the strength of the interaction but also its repulsive or attractive character. Theoretical results are supported qualitatively by experimental observations. Furthermore, this work may pave the way toward designing controlled linear assemblies as well as one-dimensional crystals of nanoparticles such as gold nanorods or quantum dots with tunable interparticle spacing.

36 MATERIALS SCIENCE↗

DOE OSTI2021

Gate-tunable plasmons in mixed-dimensional van der Waals heterostructures

Abstract Surface plasmons, collective electromagnetic excitations coupled to conduction electron oscillations, enable the manipulation of light–matter interactions at the nanoscale. Plasmon dispersion of metallic structures depends sensitively on their dimensionality and has been intensively studied for fundamental physics as well as applied technologies. Here, we report possible evidence for gate-tunable hybrid plasmons from the dimensionally mixed coupling between one-dimensional (1D) carbon nanotubes and two-dimensional (2D) graphene. In contrast to the carrier density-independent 1D Luttinger liquid plasmons in bare metallic carbon nanotubes, plasmon wavelengths in the 1D-2D heterostructure are modulated by 75% via electrostatic gating while retaining the high figures of merit of 1D plasmons. We propose a theoretical model to describe the electromagnetic interaction between plasmons in nanotubes and graphene, suggesting plasmon hybridization as a possible origin for the observed large plasmon modulation. The mixed-dimensional plasmonic heterostructures may enable diverse designs of tunable plasmonic nanodevices.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1995

Large HAWT wake measurement and analysis

From the theoretical fluid dynamics point of view, the wake region of a large horizontal-axis wind turbine has been defined and described, and numerical models of wake behavior have been developed. Wind tunnel studies of single turbine wakes and turbine array wakes have been used to verify the theory and further refine the numerical models. However, the effects of scaling, rotor solidity, and topography on wake behavior are questions that remain unanswered. In the wind tunnel studies, turbines were represented by anything from scaled models to tea strainers or wire mesh disks whose solidity was equivalent to that of a typical wind turbine. The scale factor compensation for the difference in Reynolds number between the scale model and an actual turbine is complex, and not typically accounted for. Though it is wise to study the simpler case of wakes in flat topography, which can be easily duplicated in the wind tunnel, current indications are that wind turbine farm development is actually occurring in somewhat more complex terrain. Empirical wake studies using large horizontal-axis wind turbines have not been thoroughly composited, and, therefore, the results have not been applied to the well-developed theory of wake structure. The measurement programs have made use of both in situ sensor systems, such as instrumented towers, and remote sensors, such as kites and tethered, balloonborne anemometers. We present a concise overview of the work that has been performed, including our own, which is based on the philosophy that the MOD-2 turbines are probably their own best detector of both the momentum deficit and the induced turbulence effect downwind. Only the momentum deficit aspects of the wake/machine interactions have been addressed. Both turbine power output deficits and wind energy deficits as measured by the onsite meteorological towers have been analyzed from a composite data set. The analysis has also evidenced certain topographic influences on the operation of spatially diverse wind turbines.

Miller, A. H.↗

DOE OSTI2015

Computational and experimental study of laminar flames

During the past three years, our research has centered on an investigation of the effects of complex chemistry and detailed transport on the structure and extinction of hydrocarbon flames in coflowing axisymmetric configurations. We have pursued both computational and experimental aspects of the research in parallel on both steady-state and time-dependent systems. The computational work has focused on the application of accurate and efficient numerical methods for the solution of the steady-state and time-dependent boundary value problems describing the various reacting systems. Detailed experimental measurements were performed on axisymmetric coflow flames using two-dimensional imaging techniques. Previously, spontaneous Raman scattering, chemiluminescence, and laser-induced fluorescence were used to measure the temperature, major and minor species profiles. Particle image velocimetry (PIV) has been used to investigate velocity distributions and for calibration of time-varying flames. Laser-induced incandescence (LII) with an extinction calibration was used to determine soot volume fractions, while soot surface temperatures were measured with three-color optical pyrometry using a color digital camera. A blackbody calibration of the camera allows for determination of soot volume fraction as well, which can be compared with the LII measurements. More recently, we have concentrated on a detailed characterization of soot using a variety of techniques including time-resolved LII (TiRe-LII) for soot primary particles sizes, multi-angle light scattering (MALS) for soot radius of gyration, and spectrally-resolved line of sight attenuation (spec-LOSA). Combining the information from all of these soot measurements can be used to determine the soot optical properties, which are observed to vary significantly depending on spatial location and fuel dilution. Our goal has been to obtain a more fundamental understanding of the important fluid dynamic and chemical interactions in these flames so that this information can be used effectively in combustion modeling.

02 PETROLEUM↗

NASA NTRS2000

Physics of Hard Spheres Experiment: Significant and Quantitative Findings Made

Direct examination of atomic interactions is difficult. One powerful approach to visualizing atomic interactions is to study near-index-matched colloidal dispersions of microscopic plastic spheres, which can be probed by visible light. Such spheres interact through hydrodynamic and Brownian forces, but they feel no direct force before an infinite repulsion at contact. Through the microgravity flight of the Physics of Hard Spheres Experiment (PHaSE), researchers have sought a more complete understanding of the entropically driven disorder-order transition in hard-sphere colloidal dispersions. The experiment was conceived by Professors Paul M. Chaikin and William B. Russel of Princeton University. Microgravity was required because, on Earth, index-matched colloidal dispersions often cannot be density matched, resulting in significant settling over the crystallization period. This settling makes them a poor model of the equilibrium atomic system, where the effect of gravity is truly negligible. For this purpose, a customized light-scattering instrument was designed, built, and flown by the NASA Glenn Research Center at Lewis Field on the space shuttle (shuttle missions STS 83 and STS 94). This instrument performed both static and dynamic light scattering, with sample oscillation for determining rheological properties. Scattered light from a 532- nm laser was recorded either by a 10-bit charge-coupled discharge (CCD) camera from a concentric screen covering angles of 0 to 60 or by sensitive avalanche photodiode detectors, which convert the photons into binary data from which two correlators compute autocorrelation functions. The sample cell was driven by a direct-current servomotor to allow sinusoidal oscillation for the measurement of rheological properties. Significant microgravity research findings include the observation of beautiful dendritic crystals, the crystallization of a "glassy phase" sample in microgravity that did not crystallize for over 1 year in 1g (Earth's gravity), and the emergence of face-centered-cubic (FCC) crystals late in the coarsening process (as small crystallites lost particles to the slow ripening of large crystallites). Significant quantitative findings from the microgravity experiments have been developed describing complex interactions among crystallites during the growth process, as concentration fields overlap in the surrounding disordered phase. Time-resolved Bragg scattering under microgravity captures one effect of these interactions quite conclusively for the sample at a volume fraction of 0.528. From the earliest time until the sample is almost fully crystalline, the size and overall crystallinity grow monotonically, but the number of crystallites per unit volume (number density) falls. Apparently nucleation is slower than the loss of crystallites because of the transfer of particles from small to large crystals. Thus, coarsening occurs simultaneously with growth, rather than following the completion of nucleation and growth as is generally assumed. In the same sample, an interesting signature appears in the apparent number density of crystallites and the volume fraction within the crystallites shortly before full crystallinity is reached. A brief upturn in both indicates the creation of more domains of the size of the average crystallite simultaneous with the compression of the crystallites. Only the emergence of dendritic arms offers a reasonable explanation. The arms would be "seen" by the light scattering as separate domains whose smaller radii of curvature would compress the interior phase. In fiscal year 1999, numerous papers, a doctoral dissertation, and the PHaSE final report were produced. Although this flight project has been completed, plans are in place for a follow-on colloid experiment by Chaikin and Russel that employs a light microscope within Glenn's Fluids and Combustion Facility on the International Space Station. PHaSE is providing us with a deeper understanding of the nure of phase transitions. The knowledge derived has added to the understanding of condensed matter. In addition, the burgeoning study of the dynamics of colloidal self-assembly may lead to the development of a range of photonic materials that control the desirable properties of light. Thus, applications of ordered colloidal structures include not only ultrastructure ceramics, but also photonic crystals and photothermal nanosecond light-switching devices. Industries dealing with semiconductors, electro-optics, ceramics, and composites stand to benefit from such advancements.

Doherty, Michael P.↗

DOE OSTI2016

Review: magnetically assisted resistance spot welding

Currently, the use of advanced high strength steels (AHSSs) is the most cost effective means of reducing vehicle body weight and maintaining structural integrity at the same time. However, AHSSs present a big challenge to the traditional resistance spot welding (RSW) widely applied in automotive industries because the rapid heating and cooling procedures during RSW produce hardened weld microstructures, which lower the ductility and fatigue properties of welded joints and raise the probability of interfacial failure under external loads. Changing process parameters or post-weld heat treatment may reduce the weld brittleness, but those traditional quality control methods also increase energy consumption and prolong cycle time. In recent years, a magnetically assisted RSW (MA-RSW) method was proposed, in which an externally applied magnetic field would interact with the conduction current to produce a Lorentz force that would affect weld nugget formation. This paper is a review of an experimental MA-RSW platform, the mode of the external magnetic field and the mechanism that controls nugget shape, weld microstructures and joint performance. In conclusion, the advantages of the MA-RSW method in improving the weldability of AHSSs are given, a recent application of the MA-RSW process to light metals is described and the outlook for the MA-RSW process is presented.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Thermodynamic and Kinetic Modulation of Methylammonium Lead Bromide Crystallization Revealed by In Situ Monitoring

Hybrid organic–inorganic perovskite (HOIP) crystals are promising optoelectronic materials, but little is known about either the thermodynamic and kinetic controls on crystal growth or the underlying growth mechanism(s). Herein, we use fluid cell atomic force microscopy (AFM) and solution nuclear magnetic resonance (NMR) spectroscopy to investigate growth of the model HOIP crystal CH 3 NH 3 PbBr 3 (MAPbBr 3 ) and to determine how formic acid (HCOOH) modulates the thermodynamics and kinetics of growth. The results show that growth of MAPbBr 3 in dimethylformamide (DMF) proceeds through the classical pathway by the spreading of atomic crystal steps generated at screw dislocations on the {100} surface. Temperature dependent step velocity measurements demonstrate that with increasing concentration, HCOOH decreases both the solubility of MAPbBr 3 and the kinetic coefficient (b) of step movement. 1 H-NMR measurements indicate that HCOOH increases the lifetime of the methylammonium (MA + ) ions and promotes the association of MAPbBr 3 , thus tuning the solubility of the perovskite. HCOOH also alters the molecular tumbling motion and bulk diffusion of the MA + ions, possibly via H-bonding. Further, our findings establish a direct correlation between the mesoscale crystal growth kinetics and the molecular-scale interactions between organic additives and constituent ions, providing unprecedented insights for developing predictive syntheses of HOIP crystals with defined size, crystal habit and shape, and defect distribution.

36 MATERIALS SCIENCE↗

NASA NTRS1999

Novel Approaches to Cellular Transplantation from the US Space Program

Research in the treatment of type I diabetes is entering a new era that takes advantage of our knowledge in an ever increasing variety of scientific disciplines. Some may originate from very diverse sources, one of which is the Space Program at National Aeronautics and Space Administration (NASA). The Space Program contributes to diabetes-related research in several treatment modalities. As an ongoing effort for medical monitoring of personnel involved in space exploration activities NASA and the extramural scientific community investigate strategies for noninvasive estimation of blood glucose levels. Part of the effort in the space protein crystal growth program is high-resolution structural analysis insulin as a means to better understand the interaction with its receptor and with host immune components and as a basis for rational design of a "better" insulin molecule. The Space Program is also developing laser technology for potential early cataract detection as well as a noninvasive analyses for addressing preclinical diabetic retinopathy. Finally, NASA developed an exciting cell culture system that affords some unique advantages in the propagation and maintenance of mammalian cells in vitro. The cell culture system was originally designed to maintain cell suspensions with a minimum of hydrodynamic and mechanical sheer while awaiting launch into microgravity. Currently the commercially available NASA bioreactor (Synthecon, Inc., Houston, TX) is used as a research tool in basic and applied cell biology. In recent years there is continued strong interest in cellular transplantation as treatment for type I diabetes. The advantages are the potential for successful long-term amelioration and a minimum risk for morbidity in the event of rejection of the transplanted cells. The pathway to successful application of this strategy is accompanied by several substantial hurdles: (1) isolation and propagation of a suitable uniform donor cell population; (2) management of host immune rejection; (3) protection from the autoimmune component of the disease; and (4) anatomic placement of the engrafted cells that permits timely response to blood sugar levels as well as effective release and deployment of insulin. Bioreactor technology may provide some critical advances for surmounting some of these scientific hurdles. The NASA bioreactor is a horizontally rotating cylinder that is completely filled with culture medium. Gaseous exchange is maintained by a concentric cylinder of permeable silicon. In slow rotation (15-25 rpm) particles of small mass such as cells and tissue aggregates remain suspended in the rotating body of fluid. This novel approach suspends cells without stirring thus, allowing objects of different size mass to colocate and interact in a very low shear environment. In fact, analysis of the forces acting on individual cells in the rotating bioreactor reveals that the cells are continuously falling through the fluid medium. The conditions in the bioreactor permit assembly of cells into aggregates. three-dimensional tissue growth, synthesis of intercellular matrix,10 differentiation, and some sinusoid formations that may serve as a surrogate vasculature for larger tissue segments.

Pellis, Neal R.↗

DOE OSTI2019

Oxidation state and coordination environment of Pb in U-bearing minerals

This study examines the oxidation and coordination of lead (Pb) in the uranium (U) minerals uraninite (UO 2 ), coffinite (USiO 4 ), and brannerite (UTi 2 O 6 ), and minerals that have U as a minor constituent, zircon (ZrSiO 4 ) and titanite (CaTiSiO 5 ). The characterization of Pb was conducted through a combination of micro-synchrotron X-ray fluorescence mapping (µSXRF) and micro-X-ray absorption spectroscopy (µXAS), collected at the Pb L II and L III edge on characterized grains (electron microscopy and X-ray diffraction) from a range of geological settings and ages. For all minerals considered in this study, Pb L II and L III X-ray absorption near-edge structure (XANES) and extended X-ray absorption fine structure (EXAFS) data indicate that Pb exists as Pb 2+ . In the case of uraninite, DFT structure optimization calculations suggest that the substitution of Pb 2+ for U 4+ is facilitated by the stereochemically active electron lone pair indicative of Pb 2+ , which effectively takes the place of an oxygen vacancy (coupled substitution Pb 2+ + = U 4+ + O 2- ). In detail, however, the EXAFS and XANES data, supported by first principle calculations of the XANES associated with Pb 2+ replacing U 4+ in the crystal lattice of uraninite, indicate that Pb 2+ is present in a less distorted coordination than suggested by DFT, which may reflect a yet to be characterized link between interstitial oxygen defects and Pb 2+ incorporation, rather than a simple oxygen vacancy. For coffinite, the XANES data demonstrate that Pb exists mainly as Pb 2+ in the form of galena, indicating limited stability of the isomorphous substitution in coffinite. In contrast, Pb is homogeneously distributed in the metamict mineral brannerite, which displays a distinct XANES spectrum suggestive of Pb 2+ bonded to oxygens in a highly distorted coordination environment. Pb L II XANES indicate that radiogenic Pb exists predominantly as Pb 2+ within the crystal lattice of zircon and as distinct mineral PbO nano-scale inclusions or domains. Additionally, Pb L II XANES measurements also indicate that Pb 2+ exists within titanite although the coordination environment is unclear. Overall, these data suggest that the previous indications of the presence of Pb 4+ in radioactive minerals need to be reconsidered. The nature of radiogenic Pb 2+ incorporation in U-bearing minerals controls the mobilization of Pb in hydrothermal fluids and the utility of a mineral as a geochronometer. The Pb valence and coordination data presented in this study provide fundamental constraints imperative for our understanding of the controls governing the mobilization of radiogenic Pb during water-rock interaction and metamorphism.

58 GEOSCIENCES↗

DOE OSTI2021

Activated penetrant dynamics in glass forming liquids: size effects, decoupling, slaving, collective elasticity and correlation with matrix compressibility

In this work, we employ the microscopic self-consistent cooperative hopping theory of penetrant activated dynamics in glass forming viscous liquids and colloidal suspensions to address new questions over a wide range of high matrix packing fractions and penetrant-to-matrix particle size ratios. The focus is on the mean activated relaxation time of smaller tracers in a hard sphere fluid of larger particle matrices. This quantity also determines the penetrant diffusion constant and connects directly with the structural relaxation time probed in an incoherent dynamic structure factor measurement. The timescale of the non-activated fast dissipative process is also studied and is predicted to follow power laws with the contact value of the penetrant–matrix pair correlation function and the penetrant–matrix size ratio. For long time penetrant relaxation, in the relatively lower packing fraction metastable regime the local cage barriers are dominant and matrix collective elasticity effects unimportant. As packing fraction and/or penetrant size grows, much higher barriers emerge and the collective elasticity associated with the correlated matrix dynamic displacement that facilitates penetrant hopping becomes important. This results in a non-monotonic variation with packing fraction of the degree of decoupling between the matrix and penetrant alpha relaxation times. The conditions required for penetrant hopping to become slaved to the matrix alpha process are determined, which depend mainly on the penetrant to matrix particle size ratio. By analyzing the absolute and relative importance of the cage and elastic barriers we establish a mechanistic understanding of the origin of the predicted exponential growth of the penetrant hopping time with size ratio predicted at very high packing fractions. A dynamics-thermodynamics power law connection between the penetrant activation barrier and the matrix dimensionless compressibility is established as a prediction of theory, with different scaling exponents depending on whether matrix collective elasticity effects are important. Quantitative comparisons with simulations of the penetrant relaxation time, diffusion constant, and transient localization length of tracers in dense colloidal suspensions and cold viscous liquids reveal good agreements. Multiple new predictions are made that are testable via future experiments and simulations. Extension of the theoretical approach to more complex systems of high experimental interest (nonspherical molecules, semiflexible polymers, crosslinked networks) interacting via variable hard or soft repulsions and/or short range attractions is possible, including under external deformation.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1994

Towards Understanding the Fluid Dynamic Phenomenon of Interest to Rocket Base Heating: A Review

The significance of the base heating problem for rockets during ascent is due to the complex interaction between the rocket nozzle plumes and the external-flow which can change the flow characteristics in the base region dramatically. At lower altitudes the external-flow merges with the plume-flow, without the formation of a large separated flow region, and the cooler external-flow promotes convective cooling of the base wall. Under these conditions the majority of the base heating is due to radiative heating from the shock heated plume gases. At higher altitudes, however, the process of base heating is not so straightforward. The plume and the base flow expands dramatically and separated flow regions occur in the base area. Hot exhaust gases from the rocket nozzle will be entrained into the separated flow regions and produce a convective component to the base wall heating. Further, if the rocket exhaust-gas contains soot, the soot can increase the emission from the gas and dramatically increase the wall absorption coefficient for radiative heating if it is deposited on the walls . In addition, if the rocket exhaust gas is fuel rich, the fuel can bum in the separated flow regions and further increase the base heating. The base burning phenomenon, and the increased base heating caused by it at higher altitudes, have been observed for the Space Shuttle and Saturn Rocket. Under these conditions, the total heating is significantly higher than the heating without separated flow in the base region, and the increase in heating is directly attributable to the fluid dynamic complexity of the base region. Realistic simulation of the base heating requires that the calculated flow environment reproduce the fluid dynamic flow features accurately. Thus, it will be necessary to introduce into the CFD codes the capability for the flow to respond to the complex vehicle geometry, the effect of turbulence, the ability to accurately reproduce the plume shock/shear layer structures and capture such phenomenon as shock induced base separation and base-burning phenomenon. A survey of experimental, theoretical and computational work that details the fluid dynamics of the base flow environment will be presented in the proposed paper. CFD simulations of rocket base flows using standard CFD codes such as OVERFLOW or GASP will be explored to capture these phenomenon accurately. Merits and limitations of these codes for base flow environment predictions will be explored.

Venkatapathy, E.↗

NASA NTRS1999

Double Diffusive Convection in Materials Processing

A great number of crystals grown in space are plagued by convective motions which contribute to structural flaws. The character of these instabilities is not well understood but is associated with density variations in the presence of residual gravity (g-jitter). As a specific example, past HgCdTe crystal growth space experiments by Lehoczky and co-workers indicate radial compositional asymmetry in the grown crystals. In the case of HgCdTe the rejected component into the melt upon solidification is HgTe which is denser than the melt. The space grown crystals indicate the presence of three dimensional flow with the heavier HgTe-rich material clearly aligned with the residual gravity (0.55-1.55 micro g) vector. This flow stems from double-diffusive convection, namely, thermal and solutal buoyancy driven flow in the melt. The study of double-diffusive convection is multi-faceted and rather vast. In our investigation, we seek to focus on one specific aspect of this discipline that is of direct relevance to materials processing especially crystal growth, namely, the side ways heating regime. This problem has been widely studied, both experimentally and numerically, in the context of solar ponds wherein the system is characterized by a linear salt (solutal) gradient with an imposed lateral temperature gradient. The induced flow instabilities arise from the wide disparity between the fluid thermal diffusivity and the solute diffusivity. The extension of the analysis to practical crystal growth applications has however not been rigorously made and understood. One subtle but important difference in crystal growth systems is the fact that die system solute gradient is non-linear (typically exponential). Besides, the crystal growth problem has the added complexities of solidification, both lateral and longitudinal thermal gradients and segregation phenomena in systems where binary and ternary compounds are being grown. This paper treats the side ways heating problem alone in a model fluid system. Results from detailed numerical calculations, mainly two dimensional are provided. The interactions between a non-linear solute gradient and an imposed transverse thermal gradient are investigated. The buoyancy effects are treated in the traditional Boussinesq approximation and also in a more complete density formulation to address recent concerns of the first approach especially in simulations of the system response in a reduced gravity environment. Detailed flow, temperature and solute field plots along with heat and mass transfer results are presented in the paper. Implications to practical crystal growth systems as discerned from the modeling results are also explored and reported.

Ramachandra, Narayanan↗

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

Shebalin, John V↗

DOE OSTI2019

Dynamics of potential vorticity staircase evolution and step mergers in a reduced model of beta-plane turbulence

A two-field model of potential vorticity (PV) staircase structure and dynamics relevant to both beta-plane and drift-wave plasma turbulence is studied numerically and analytically. The model evolves averaged PV whose flux is both driven by and regulates a potential enstrophy field, ε. The model employs a closure using a mixing length model. Its link to bistability, vital to staircase generation, is analyzed and verified by integrating the equations numerically. Long-time staircase evolution consistently manifests a pattern of metastable quasiperiodic configurations, lasting for hundreds of time units, yet interspersed with abrupt ( Δt << 1 ) mergers of adjacent steps in the staircase. The mergers occur at the staircase lattice defects where the pattern has not completely relaxed to a strictly periodic solution that can be obtained analytically. Other types of stationary solutions are solitons and kinks in the PV gradient and ε-profiles. The waiting time between mergers increases strongly as the number of steps in the staircase decreases. This is because of an exponential decrease in interstep coupling strength with growing spacing. The long-time staircase dynamics is shown numerically be determined by local interaction with adjacent steps. In conclusion, mergers reveal themselves through the explosive growth of the turbulent PV flux, which, however, abruptly drops to its global constant value once the merger is completed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2012

Derivation of nonlinear Schroedinger equation for electrostatic and electromagnetic waves in fully relativistic two-fluid plasmas by the reductive perturbation method

The reductive perturbation method is used to derive a generic form of nonlinear Schroedinger equation (NLSE) that describes the nonlinear evolution of electrostatic (ES)/electromagnetic (EM) waves in fully relativistic two-fluid plasmas. The matrix eigenvector analysis shows that there are two mutually exclusive modes of waves, each mode involving only either one of two electric potentials, A and {phi}. The general result is applied to the electromagnetic mode in electron-ion plasmas with relativistically high electron temperature (T{sub e} Much-Greater-Than m{sub e}c{sup 2}). In the limit of high frequency (ck Much-Greater-Than {omega}{sub e}), the NLSE predicts bump type electromagnetic soliton structures having width scaling as {approx}kT{sub e}{sup 5/2}. It is shown that, in electron-positron pair plasmas with high temperature, dip type electromagnetic solitons can exist. The NLSE is also applied to electrostatic (Langmuir) wave and it is shown that dip type solitons can exist if k{lambda}{sub D} Much-Less-Than 1, where {lambda}{sub D} is the electron's Debye length. For the k{lambda}{sub D} Much-Greater-Than 1, however, the solution is of bump type soliton with width scaling as {approx}1/(k{sup 5}T{sub e}). It is also shown that dip type solitons can exist in cold plasmas having relativistically high streaming speed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2013

Magnetohydrodynamics and Heat Transfer in a Free-Surface, Flowing Liquid Metal Experiment

This thesis examines magnetohydrodynamic (MHD) effects and heat transfer in free-surface, liquid metal flow. The interaction of the flow with the applied magnetic field introduces an anisotropy to the fluid motion, resulting in starkly different behavior depending on the orientation and strength of the applied magnetic field. The transition from a hydrodynamic state to an MHD state was experimentally found to be characterized by the interaction parameter, N = &sigma L B 2 /&rho v0 , crossing above 0.5. The experiment consisted of a channel 1 m in length, up to 16 cm in width and 1 to 2.5 cm in depth situated within an electromagnet capable of producing a uniform magnetic field up to 2.7 kG. The fluid velocity was measured with an array of 25 potential probes that resolved the spanwise velocity profile, while waves on the surface of the metal were monitored via a diagnostic that tracked the motion of a laser beam reflected off of the free surface. Lastly, the temperature response of the fluid was recorded with an array of 32 thermocouples embedded in the bottom of the channel and a mid-wavelength infrared camera used to image the free surface. Three distinct sets of experiments were performed to investigate the effects of the applied magnetic field on the dynamics within the flow, studying: 1) changes to vortices in the wake of a cylinder, 2) behavior of surface fluctuations due to the turbulent flow, and 3) vertical heat transfer resulting from heat deposition on the free surface. The first set of experiments examined the wake of a cylinder inserted into the flow with its axis parallel to the magnetic field. Measurements with the potential probes indicated that vortices in the wake became laminarized as the field was increased due to the alignment of the injected vorticity with the applied field. The second set of experiments showed that surface fluctuations were suppressed at high magnetic fields, but the precise nature of the damping depended on the orientation of the applied field. Lastly, when heat was injected at the free surface via a resistive heater, the temperature response of the fluid was found to vary considerably depending on the strength and direction of the applied field. In all experiments, it was observed that the angle between the vorticity of structures in the flow and the applied magnetic field dictated the macroscopic behavior of the fluid. This work provides the first experimental results of several phenomena associated with free-surface, liquid metal flows in an applied magnetic field. Further, the results have direct implications for proposed flowing liquid metal divertors for fusion reactors, and offer insight into the nature of turbulent MHD flows.

fusion↗

DOE OSTI2015

Study of the influence of heat sources on the out-of-pile calibration curve of calorimetric cells used for nuclear energy deposition quantification

At present the Jules Horowitz Reactor is under construction in Cadarache research center of CEA 'French Alternative Energies and Atomic Energy Commission' center located in the south-east of France. This new Material Testing Reactor (MTR) will be operational in late 2019 and will allow the generation of a new experimental potential (up to 20 irradiation devices simultaneously) and new harsh conditions such as higher neutron fluxes (5.10{sup 14} n.cm{sup -2}.s{sup -1} for E≥1 MeV), faster material ageing and higher nuclear heating (up to 20 W/g for nominal capacity of 100 MW). In nuclear research field, the control and the measurement of the nuclear heating (energy deposition rate per mass unit induced by the interactions of radiations with matter) is crucial to carry out accurate studies on ageing of materials and on the behavior of nuclear fuels under irradiation. Several experiments need to know precisely this key parameter in order to establish dedicated thermal conditions. The measurement of the nuclear heating inside MTRs is realized by three kinds of sensors: single-cell calorimeter, differential calorimeter and gamma thermometer. One scientific objective of the IN-CORE program, between CEA and Aix-Marseille University in 2009, is to improve the nuclear heating measurement. In this context a new multi-sensor device, called CARMEN, was made. This device contains in particular a differential calorimeter which was designed to measure the nuclear heating in the periphery of OSIRIS reactor (a MTR located at Saclay, France) up to 2 W/g and tested during two irradiation campaigns. Results obtained during these campaigns showed that temperatures reached inside the calorimeter are higher than ones obtained during the preliminary out-of-pile calibration experiments. For instance for 1.74 W/g, the in-pile temperature of the calorimeter rod is equal to 305 deg. C against 225 deg. C in laboratory conditions by simulating the nuclear heating by Joule Effect inside the calorimeter cell head. This discrepancy is higher than in previous experiments because the calorimeter owns a high sensitivity. Consequently, a new prototype was created and instrumented by other heat sources in order to impose an energy deposition on the calorimetric cell structure (in particular in the base) and to improve the calibration step in out-of-pile conditions. In this paper, on the first part a detailed description of the new calorimetric sensor will be given. On the second part, the experimental response of the sensor obtained for several internal heating conditions will be shown. The influence of these conditions on the calibration curve will be discussed. Then the response of this prototype will be also presented for different external cooling fluid conditions (in particular flow temperature). In this part, the comparison between the in-pile and out-of-pile experimental results will be performed. On the last part, these out-of-pile experiments will be completed by 2D axisymmetrical thermal simulations with the CEA code CAST3M using Finite Elements Method. After a comparison between experimental and numerical works, improvements of the sensor prototype will be studied (new heat sources). (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗