Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “turbulent correlation length”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Comparison of Corcos-based and experimentally-derived coherence factors for BFFs estimation

In this paper, high-spatial-resolution unsteady Pressure Sensitive Paint (uPSP) data are utilized to compare two methods for panel Buffet Forcing Functions (BFF) estimation for the Space Launch System (SLS). Such methods are based on discrete pressure measurements within a panel but employ coherence factors to account for partially-correlated fluctuating pressures across the whole panel. In one method, coherence factors are derived based on the Corcos model, whereas the second method utilizes experimentally-derived coherence factors. To simulate discrete measurements using uPSP data, suitable subsets of the data are extracted. When full uPSP resolution is retained, uPSP data provide a benchmark to assess discrete-measurements-based methods. The analysis focuses on the peak SLS buffet environment located downstream of the Forward Attachment Hardware (FAH) between the core stage and solid rocket boosters. Trends of Corcos-based and experimentally-derived coherence factors are in reasonable agreement with the benchmark. However, at certain frequencies, experimentally-derived coherence factors are sensitive to the separation distance between pressure measurements utilized to compute coherence lengths. Such sensitivity originates from deviation of the experimental-based coherence function from an exponential decay assumption. On the other hand, the present implementation of the Corcos model fails to capture certain non-turbulent boundary layer related environments, such as a subharmonic of FAH vortex-shedding. For all methods presented in this paper, at near transonic conditions, increased pressure coherence and spatial nonuniformity lead to BFF overestimation and sensitivity to the pressure measurement location within the panel.

transonic buffet↗
NASA NTRS2019

Optimization of a Remote Sensing Energy Balance Method over Different Canopy Applied at Global Scale

Parameterization methods which calculate turbulent heat and water fluxes with thermal remote sensing data were evaluated in the revised remote sensing surface energy balance system (SEBS) model (Chen et al., 2013). The model calculates sensible heat (H) based on the Monin-Obukhov similarity theory (MOST) and determines latent heat (LE) as the residual of energy balance. We examined the uncertainties of H and LE in the SEBS model due to five key parameters at the local station point scale. Observations at 27 flux towers located in seven land cover types (needle-leaf forest, broad leaf forest, shrub, savanna, grassland, cropland, and sparsely vegetated land) and an artificial intelligence particle swarm optimization (PSO) algorithm was combined to calibrate the five parameters (leaf drag coefficient, leaf heat transfer coefficients, roughness length for soil, and two parameters for ground heat calculation) in the SEBS model. The root-mean-square error at the site scale was reduced by 9 W/sq.m for H, and 92 W/sq.m for LE, and their correlation coefficients were increased by 0.07 (H) and 0.11 (LE) after using the calibrated parameters. The updated model validation was further conducted globally for the remotely sensed evapotranspiration (ET) calculations. Overestimation of SEBS global ET was significantly improved by using the optimized values of the parameters. The results suggested PSO was able to consistently locate the global optimum of the SEBS model, and appears to be capable of solving the ET model optimization problem.

Chen, Xuelong↗
NASA NTRS1995

Turbulence in the solar wind: Kinetic effects

Although a casual look at the fluctuating magnetic and velocity fields in the solar wind may be reminiscent of a chaotic and disordered flow, there is, nonetheless. considerable organization and structure in the temporal and spatial evolution of those fluctuations. Much of that evolution is controlled by processes operating on rather large scales for example, in the inner heliosphere, the fluctuations in magnetic and velocity are highly correlated in the sense of outward propagating Alfven waves. This correlation can be destroyed both in time and distance by the velocity gradients present between fast and slow streams and by other nonlinear processes which stir the medium, producing a turbulent cascade of energy from large to small scales. Many aspects of this turbulent evolution can be described using fluid models; however, at some scale the fluid approximation breaks down and a more detailed paradigm is necessary. The breakdown is evident in the power spectrum of magnetic fluctuations at scales approaching the wavelength of ion cyclotron waves. At those scales, as evident in Mariner 10 and other magnetometer data, the spectrum bends over and the fluctuations damp, possibly heating the ambient plasma. Some evidence for heating of the solar wind is present in the Voyager data. Fluid models can be modified to some extent to incorporate aspects of a kinetic treatment. This is done by modifying the dissipation terms in the fluid equations and by including extra terms, such as the Hall term. As the scale lengths of phenomena shrink further and approach the spatial and temporal scales characteristic of electron phenomena, the fluid description must be abandoned altogether and a fully kinetic treatment is required. One example is the generation of Langmuir solitons produced by the electron beams that generate type 3 solar radio bursts.

Goldstein, M. L.↗
NASA NTRS1991

Heat transfer measurements from a smooth NACA 0012 airfoil

Local convective heat transfer coefficients were measured from a smooth NACA 0012 airfoil having a chord length of 0.533 m. Flight data were taken for the smooth airfoil at Reynolds numbers based on chord in the range 1.24 to 2.50 million and at various angles of attack up to 4 deg. During these flight tests, the freestream velocity turbulence intensity was found to be very low. Wind tunnel data were acquired in the Reynolds number range 1.20 to 4.52 million and at angles of attack from -4 to +8 deg. The turbulence intensity in the IRT was 0.5-0.7 percent with the cloud-generating sprays off. A direct comparison between the results obtained in flight and in the IRT showed that the higher level of turbulence intensity in the IRT had little effect on the heat transfer for the lower Reynolds numbers but caused a moderate increase in heat transfer at the higher Reynolds numbers. Turning on the cloud-generating spray nozzle atomizing air in the IRT did not alter the heat transfer. The present data were compared with leading-edge cylinder and flat plate heat transfer correlations that are often used to estimate airfoil heat transfer in computer codes.

Poinsatte, Philip E.↗
NASA NTRS2006

Shuttle Damage/Repair from the Perspective of Hypersonic Boundary Layer Transition - Experimental Results

An overview is provided of the experimental wind tunnel program conducted at the NASA Langley Research Center Aerothermodynamics Laboratory in support of an agency-wide effort to prepare the Shuttle Orbiter for Return-to-Flight. The effect of an isolated protuberance and an isolated rectangular cavity on hypersonic boundary layer transition onset on the windward surface of the Shuttle Orbiter has been experimentally characterized. These experimental studies were initiated to provide a protuberance and cavity effects database for developing hypersonic transition criteria to support on-orbit disposition of thermal protection system damage or repair. In addition, a synergistic experimental investigation was undertaken to assess the impact of an isolated mass-flow entrainment source (simulating pyrolysis/outgassing from a proposed tile repair material) on boundary layer transition. A brief review of the relevant literature regarding hypersonic boundary layer transition induced from cavities and localized mass addition from ablation is presented. Boundary layer transition results were obtained using 0.0075-scale Orbiter models with simulated tile damage (rectangular cavities) of varying length, width, and depth and simulated tile damage or repair (protuberances) of varying height. Cavity and mass addition effects were assessed at a fixed location (x/L = 0.3) along the model centerline in a region of near zero pressure gradient. Cavity length-to-depth ratio was systematically varied from 2.5 to 17.7 and length-to-width ratio of 1 to 8.5. Cavity depth-to-local boundary layer thickness ranged from 0.5 to 4.8. Protuberances were located at several sites along the centerline and port/starboard attachment lines along the chine and wing leading edge. Protuberance height-to-boundary layer thickness was varied from approximately 0.2 to 1.1. Global heat transfer images and heating distributions of the Orbiter windward surface using phosphor thermography were used to infer the state of the boundary layer (laminar, transitional, or turbulent). Test parametrics include angles-of-attack of 30 deg and 40 deg, sideslip angle of 0 deg, freestream Reynolds numbers from 0.02x106 to 7.3x106 per foot, edge-to-wall temperature ratio from 0.4 to 0.8, and normal shock density ratios of approximately 5.3, 6.0, and 12 in Mach 6 air, Mach 10 air, and Mach 6 CF4, respectively. Testing to simulate the effects of ablation from a proposed tile repair concept indicated that transition was not a concern. The experimental protuberance and cavity databases highlighted in this report were used to formulate boundary layer transition correlations that were an integral part of an analytical process to disposition observed Orbiter TPS damage during STS- 114.

Horvath, Thomas J.↗
NASA NTRS1989

Computational fluid dynamics analysis of space shuttle main propulsion feed line 17-inch disconnect valves

A steady incompressible three-dimensional (3-D) viscous flow analysis was conducted for the Space Shuttle Main Propulsion External Tank (ET)/Orbiter (ORB) propellant feed line quick separable 17-inch disconnect flapper valves for liquid oxygen (LO2) and liquid hydrogen (LH2). The main objectives of the analysis were to predict and correlate the hydrodynamic stability of the flappers and pressure drop with available water test data. Computational Fluid Dynamics (CFD) computer codes were procured at no cost from the public domain, and were modified and extended to carry out the disconnect flow analysis. The grid generator codes SVTGD3D and INGRID were obtained. NASA Ames Research Center supplied the flow solution code INS3D, and the color graphics code PLOT3D. A driver routine was developed to automate the grid generation process. Components such as pipes, elbows, and flappers can be generated with simple commands, and flapper angles can be varied easily. The flow solver INS3D code was modified to treat interior flappers, and other interfacing routines were developed, which include a turbulence model, a force/moment routine, a time-step routine, and initial and boundary conditions. In particular, an under-relaxation scheme was implemented to enhance the solution stability. Major physical assumptions and simplifications made in the analysis include the neglect of linkages, slightly reduced flapper diameter, and smooth solid surfaces. A grid size of 54 x 21 x 25 was employed for both the LO2 and LH2 units. Mixing length theory applied to turbulent shear flow in pipes formed the basis for the simple turbulence model. Results of the analysis are presented for LO2 and LH2 disconnects.

Kandula, Max↗
DOE OSTI2021

UNH TDP - Concurrent Measurements of Inflow, Power Performance, and Loads for a Grid-Synchronized Vertical Axis Cross-Flow Turbine Operating in a Tidal Estuary

This data was collected between October 12 and December 15 of 2021 at the University of New Hampshire (UNH) and Atlantic Marine Energy Center (AMEC) turbine deployment platform (TDP). This data set includes over 29 days of grid connected turbine operation during this 65 day time frame. The priority for this measurement campaign was to collect data while the turbine was electrically connected to the grid by means of a rectifier and inverter. The Fall_2021_UNH_Measurement_Timeline.png highlights when each instrument was functioning and the Fall_2021_UNH_Test_Log.jpg indicates the four main regions for analysis available from this measurement campaign. The TDP is a floating structure moored on the Portsmouth facing side of Memorial Bridge pier #2, which spans the Piscataqua River between Portsmouth, NH and Kittery, ME. The Piscataqua River connects the Great Bay Estuary to the Gulf of Maine and the river currents are dominated by tidal forcing with water velocities exceeding 2.5 m/s during spring ebb tides at this site which were previously characterized by Kaelin Chancey (Assessment Of The Localized Flow And Tidal Energy Conversion System At An Estuarine Bridge - UNH MS Thesis 2019). The turbine under test was a modified New Energy Corporation (Calgary, CA) model EVG-025 4-blade H-Darrius type vertical axis cross flow turbine that rotates in the clockwise direction with a rotor diameter of 3.2m and blade length of 1.7m. The hydro-foil profile was a NACA 0021 with a 10 inch chord length and a blade preset pitch angle of +4deg with a positive angle corresponding with the toe in direction. The standard EVG-025 has a rotor diameter of 3.4m and its rated power output is 25kW at 3 m/s. The rotor diameter was reduced to accommodate the size of the existing TDP moon-pool. This project was pursued to quantify device performance for cross flow turbines operating in a marine environment. Accurate physical models, to characterize cross flow turbine performance, require real operational data sets due to the complexity of blade fluid interactions. This data can help support model development which will help predict turbine performance when analyzing perspective project locations in the future. Instrumentation was deployed to measure; water speed/direction, electrical power output, turbine shaft speed, turbine thrust force, and platform motion. Concurrent measurements of these parameters allow for correlations (cause and affect) to be inferred, allowing for characterization of device performance over a range of operating conditions. Water currents were measured using Acoustic Doppler Current Profilers (ADCP's) and Acoustic Doppler Velocimeters (ADV's) directly upstream and downstream of the turbine for inflow, wake and turbulence measurements. Electrical power output was measured using the Voltsys rectifier and the Shark power meter. Shaft speed was calculated based on the Voltsys measurements of the permanent magnet three phase generator AC generation frequency, coupled directly to the cross flow turbine under test (i.e., no gear box). Platform motions were captured using a Yost IMU (inertial measurement unit). Turbine thrust loading was measured using a reaction arm about the turbine deployment platform spanning beam, where two bi-directional load cells were connected to the system via a pinned connection. This submission includes zipped folders for each instrument containing quality controlled (QC'd) data in daily .csv files for the relevant duration specific to each instrument, along with separate .csv file that contains the units for each variable. Some instrument daily files are quite large and can pose a challenge for a visual spreadsheet editor to open. A processing software like MATLAB or Python is recommended. Note the degree of QC varied between each instrument due to time constraints. Particular time and attention was given to perform quality control tests on the acoustic based instruments that are particularly suscep...

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Computational fluid dynamics modeling and analysis of silica nanoparticle synthesis in a flame spray pyrolysis reactor

Flame Spray Pyrolysis (FSP) is a method for large-scale production of nanoparticles and nanoscale powders employed in a wide range of industrial applications. Particle size and morphology are complex functions of the physicochemical phenomena occurring in the FSP reactor. An extensive study of FSP-related phenomena can be utilized to develop effective strategies for achieving desired particle size/morphology and scaling up the overall yield of an FSP system. In this work, a computational fluid dynamics (CFD) model of an FSP reactor is developed to simulate the coupling of key phenomena involved in the particle synthesis process: liquid spray breakup and evaporation, mixing, combustion, and particle formation/growth of silica nanoparticles. Herein, the particle sizes and their distributions from the CFD simulations are validated against experimental data. Subsequently, the simulations are utilized to investigate the impact of process parameters on the resultant flame dynamics and particle growth. Firstly, the CFD results show that the particle sizes are strongly correlated with the precursor concentration in the solvent. At lower precursor concentrations, the spread of the distribution is relatively insensitive to the value of the concentration. At higher concentrations, the spread is higher as the collision probability between particles is higher. Secondly, increasing the pilot flow rate increases the length of the pilot flames impacting the local ignition location of the spray flame. Lastly, it is shown that the dispersion gas flow rate strongly influences the spray flame shape. This shape can be used for control of particle growth as it helps determine the regions of high temperature and the residence time of the particles in the high temperature region enabling the design and process optimization of the FSP reactor.

42 ENGINEERING↗
NASA NTRS1997

Atmospheric Angular Momentum Fluctuations During 1979-1988 Simulated by Global Circulation Models

Changes in major global dynamical phenomena in the Earth's atmosphere are manifested in the time series of atmospheric angular momentum (AAM), as determined directly from meteorological observations and indirectly from geodetic observations of small fluctuations in the rotation of the solid Earth which are proportional to length of day. AAM fluctuations are intimately linked with energetic processes throughout the whole atmosphere and also with the stresses at the Earth's surface produced largely by turbulent momentum transport in the oceanic and continental boundary layers and by the action of normal pressure forces on orographic features. A stringent test of any numerical global circulation model (GCM) is therefore provided by a quantitative assessment of its ability to represent AAM fluctuations on all relevant timescales, ranging from months to several years. From monthly data provided by the Atmospheric Model Intercomparison Project (AMIP) of the World Climate Research Programme, we have investigated seasonal and interannual fluctuations and the decadal mean in the axial component of AAM in 23 AMIP GCMs over the period 1979-1988. The decadal means are generally well simulated, with the model median value (1.58 x 10(exp 26) kg sq m/s) being only 3.5% larger than the observed mean and with 10 of the models being within 5% of the observed. The seasonal cycle is well reproduced, with the median amplitude of the models seasonal standard deviations being only 2.4% larger than observed. Half the seasonal amplitudes lie within 15% of the observed, and the median correlation found between the observed and model seasonal cycles is 0.95. The dominant seasonal error is an underestimation of AAM during northern hemisphere winter associated with errors in the position of subtropical jets. Less robust are the modeled interannual variations, although the median correlation of 0.61 between model simulations and observed AAM is statistically significant. The two El Nino-Southem Oscillation events that occurred during the AMIP decade 1979-1988 have the expected positive AAM anomalies, although the AAM signature of the 1982-1983 event tends to be underestimated and that of the 1986-1987 event overestimated.

Hide, R.↗
NASA NTRS1997

Atmospheric Angular Momentum Fluctuations During 1979-1988 Simulated by Global Circulation Models

Changes in major global dynamical phenomena in the Earth's atmosphere are manifested in the time series of atmospheric angular momentum (AAM), as determined directly from meteorological observations and indirectly from geodetic observations of small fluctuations in the rotation of the solid Earth which are proportional to length of day. AAM fluctuations are intimately linked with energetic processes throughout the whole atmosphere and also with the stresses at the Earth's surface produced largely by turbulent momentum transport in the oceanic and continental boundary layers and by the action of normal pressure forces on orographic features. A stringent test of any numerical global circulation model (GCM) is therefore provided by a quantitative assessment of its ability to represent AAM fluctuations on all relevant timescales, ranging from months to several years. From monthly data provided by the Atmospheric Model Intercomparison Project (AMIP) of the World Climate Research Programme, we have investigated seasonal and interannual fluctuations and the decadal mean in the axial component of AAM in 23 AMIP GCMs over the period 1979-1 988. The decadal means are generally well simulated, with the model median value (1.58 x 10(exp 26) kg sq m/s) being only 3.5% larger than the observed mean and with 10 of the models being within 5% of the observed. The seasonal cycle is well reproduced, with the median amplitude of the models' seasonal standard deviations being only 2.4% larger than observed. Half the seasonal amplitudes lie within 15% of the observed, and the median correlation found between the observed and model seasonal cycles is 0.95. The dominant seasonal error is an under- estimation of AAM during northern hemisphere winter associated with errors in the position of subtropical jets. Less robust are the modeled interannual variations, although the median correlation of 0.61 between model simulations and observed AAM is statistically significant. The two El Nino-Southern Oscillation events that occurred during the AMIP decade 1979-1988 have the expected positive AAM anomalies, although the AAM signature of the 1982-1983 event tends to be underestimated and that of the 1986-1987 event overestimated.

Hide, R.↗
NASA NTRS1952

Experimental Investigation of Average Heat-Transfer and Friction Coefficients for Air Flowing in Circular Tubes Having Square-Thread-Type Roughness

An investigation of forced-convection heat transfer and associated pressure drops was conducted with air flowing through electrically heated Inconel tubes having various degrees of square-thread-type roughness, an inside diameter of 1/2 inch, and a length of 24 inches. were obtained for tubes having conventional roughness ratios (height of thread/radius of tube) of 0 (smooth tube), 0.016, 0.025, and 0.037 over ranges of bulk Reynolds numbers up to 350,000, average inside-tube-wall temperatures up to 1950deg R, and heat-flux densities up to 115,000 Btu per hour per square foot. Data The experimental data showed that both heat transfer and friction increased with increase in surface roughness, becoming more pronounced with increase in Reynolds number; for a given roughness, both heat transfer and friction were also influenced by the tube wall-to-bulk temperature ratio. Good correlation of the heat-transfer data for all the tubes investigated was obtained by use of a modification of the conventional Nusselt correlation parameters wherein the mass velocity in the Reynolds number was replaced by the product of air density evaluated at the average film temperature and the so-called friction velocity; in addition, the physical properties of air were evaluated at the average film temperature. The isothermal friction data for the rough tubes, when plotted in the conventional manner, resulted in curves similar to those obtained by other investigators; that is, the curve for a given roughness breaks away from the Blasius line (representing turbulent flow in smooth tubes) at some value of Reynolds number, which decreases with increase in surface roughness, and then becomes a horizontal line (friction coefficient independent of Reynolds number). A comparison of the friction data for the rough tubes used herein indicated that the conventional roughness ratio is not an adequate measure of relative roughness for tubes having a square-thread-type element. The present data, as well as those of other investigators, were used to isolate the influence of ratios of thread height to width, thread spacing to width, and the conventional roughness ratio on the friction coefficient. A fair correlation of the friction data was obtained for each tube with heat addition when the friction coefficient and Reynolds number were defined on the basis of film properties; however, the data for each tube retained the curve characteristic of that particular roughness. The friction data for all the rough tubes could be represented by a single line for the complete turbulence region by incorporating a roughness parameter in the film correlation. No correlation was obtained for the region of incomplete turbulence.

Sams, E. W.↗
DOE OSTI2024

On the subject of large-scale pool fires and turbulent boundary layer interactions

The role to which a realistic inflow turbulent boundary layer (TBL) influences transient and mean large-scale pool fire quantities of interest (QoIs) is numerically investigated. High-fidelity, low-Mach large-eddy simulations that activate low-dissipation, unstructured numerics are conducted using an unsteady flamelet combustion modeling approach with mutiphysics coupling to soot and participating media radiation transport. Three inlet profile configurations are exercised for a large-scale, high-aspect rectangular pool that is oriented perpendicular to the flow direction: a time-varying, TBL inflow profile obtained from a periodic precursor simulation, the time-mean of the transient TBL, and a steady power-law inflow profile that replicates the mean TBL crosswind velocity of 10.0 m/s at a vertical height of 10 m. Results include both qualitative transient flame evolution and quantitative flame shape with ground-level temperature and convective/radiative heat flux profiles. While transient fire events, which are driven by burst-sweep TBL coupling, such as blow-off and reattachment are vastly different in the TBL case (contributing to increased root mean square QoI fluctuation prediction and disparate flame lengths), mean surface QoI magnitudes are similar. Quadrant analysis demonstrates that the TBL configuration modifies burst-sweep phenomena at windward pool locations, while leeward recovery is found. In conclusion, positive fluctuations of convective heat flux correlate with fast moving fluid away from the pool surface due to intermittent combustion events.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS2001

Initial Circulation and Peak Vorticity Behavior of Vortices Shed from Airfoil Vortex Generators

An extensive parametric study of vortices shed from airfoil vortex generators has been conducted to determine the dependence of initial vortex circulation and peak vorticity on elements of the airfoil geometry and impinging flow conditions. These elements include the airfoil angle of attack, chord length, span, aspect ratio, local boundary layer thickness, and free stream Mach number. In addition, the influence of airfoil-to-airfoil spacing on the circulation and peak vorticity has been examined for pairs of co-rotating and counter-rotating vortices. The vortex generators were symmetric airfoils having a NACA-0012 cross-sectional profile. These airfoils were mounted either in isolation, or in pairs, on the surface of a straight pipe. The turbulent boundary layer thickness to pipe radius ratio was about 17 percent. The circulation and peak vorticity data were derived from cross-plane velocity measurements acquired with a seven-hole probe at one chord-length downstream of the airfoil trailing edge location. The circulation is observed to be proportional to the free-stream Mach number, the angle-of-attack, and the span-to-boundary layer thickness ratio. With these parameters held constant, the circulation is observed to fall off in monotonic fashion with increasing airfoil aspect ratio. The peak vorticity is also observed to be proportional to the free-stream Mach number, the airfoil angle-of-attack, and the span-to-boundary layer thickness ratio. Unlike circulation, however, the peak vorticity is observed to increase with increasing aspect ratio, reaching a peak value at an aspect ratio of about 2.0 before falling off again at higher values of aspect ratio. Co-rotating vortices shed from closely spaced pairs of airfoils have values of circulation and peak vorticity under those values found for vortices shed from isolated airfoils of the same geometry. Conversely, counter-rotating vortices show enhanced values of circulation and peak vorticity when compared to values obtained in isolation. The circulation may be accurately modeled with an expression based on Prandtl's relationship between finite airfoil circulation and airfoil geometry. A correlation for the peak vorticity has been derived from a conservation relationship equating the moment at the airfoil tip to the rate of angular momentum production of the shed vortex, modeled as a Lamb (ideal viscous) vortex. This technique provides excellent qualitative agreement to the observed behavior of peak vorticity for low aspect ratio airfoils typically used as vortex generators.

Wendt, Bruce J.↗
DOE OSTI2025

Turbulence statistical analysis of the L-H transition and RMPs in KSTAR

Here, we investigate the turbulence statistics associated with low-to-high confinement (L-H) transitions and externally applied resonant magnetic perturbations (RMPs) in KSTAR. Time-series fluctuations of electron density n e , electron temperature T e , and the time derivative of the poloidal magnetic field dB θ /dt (Mirnov coils) are analysed using information-geometric measures (information rate Γ and information length $\mathcal{L}$ = ∫ Γ dt), together with kurtosis κ and variance σ 2 . In low-density upper single-null plasmas (n e ~ 1.2 x 10 19 m -3 ), a ~80 kHz magnetic mode coupling n e , T e , dB θ /dt and emerges prior to the L-H transition and persists into the edge-localised modes H-mode. Edge-localised RMPs (ERMPs) suppress this coherent mode but enhance intermittency, producing frequent bursts that abruptly reshape the time-dependent probability density functions (PDFs) and generate large spikes in Γ (with smaller changes in κ), signalling ERMP-driven departures from quasi-stationarity. The impact of ERMPs on background fluctuation levels depends on density, radial location, and the fluctuating variable itself ($\tilde{n}$, $\tilde{T}$, $\dot{B}$ θ ), whereas $\mathcal{L}$ provides a robust, regime-agnostic measure of cumulative statistical reorganisation and spatial decorrelation. In particular, at low density we observe weaker coupling between $\tilde{n}$ and $\tilde{T}$, along with a tendency toward decreased radial correlation-most clearly for $\tilde{T}$-under ERMPs. Overall, information geometry cleanly captures intermittent events, quantifies non-equilibrium PDF evolution, and offers a compact, cross-diagnostic metric for assessing resonant magnetic perturbation effects on edge transport and correlation across densities, radial locations, and confinement states.

Kim, Eun-jin [Coventry Univ. (United Kingdom); Seo↗
DOE OSTI2024

Direct numerical simulations of three-component Rayleigh–Taylor mixing and an improved model for multicomponent reacting mixtures

We present direct numerical simulations of a three-layer Rayleigh–Taylor instability (RTI) problem with a configuration based on the experiments of Suchandra & Ranjan ( J. Fluid Mech. , vol. 974, 2023, A35) and Jacobs & Dalziel ( J. Fluid Mech. , vol. 542, 2005, pp. 251–279). The problem consists of a layer of light fluid between two layers of heavy fluid with an Atwood number of 0.3. These simulations are first validated through comparison with available experimental data. The validated simulations are then utilized to analyse statistics in this three-component flow. First, length scales are examined utilizing spectra and two-point spatial correlations of velocity and species concentration fluctuations. Next, joint probability density functions (p.d.f.s) of species concentration are compared against several model p.d.f.s representing generalizations of the bivariate beta distribution. Notably, the joint p.d.f.s do not appear to be accurately described by a Dirichlet distribution, indicating the marginal distributions do not conform to a beta distribution. Finally, similarity of the present configuration to three-component mixing found in inertial confinement fusion (ICF) applications is exploited to develop and validate an improved model for the impact of multicomponent mixing on thermonuclear (TN) reaction rates. A single time instant from the present simulations is chosen for a TN burn calculation under the hypothetical assumption of ICF materials and temperatures. Total TN output from this second calculation is then compared against the prediction of the improved model. The new model is found to accurately predict TN reaction rates in both premixed and non-premixed configurations.

42 ENGINEERING↗
DOE OSTI2022

Dependence of scrape-off layer profiles and turbulence on gas fuelling in high density H-mode regimes in TCV

A set of high density, highly shaped H-mode discharges has been performed in the TCV tokamak with the aim of assessing the effects of increasing divertor neutral recycling on the properties of upstream inter-ELM scrape-off layer (SOL) profiles and transport. An increase of divertor neutral pressure has been correlated with the evolution of separatrix properties and turbulence level. The latter has been quantified by means of the α t parameter introduced in (Eich 2020 Nucl. Fusion 60 056016), describing the contribution of resistive-interchange turbulence in the SOL relative to drift wave transport. The analysis reveals a general broadening of the upstream SOL profiles as α t increases, with the SOL power width measured by the vertical IR thermography system increasing significantly. In a similar way, the upstream density profile widens in the near SOL, whereas in the far SOL a density shoulder is observed to progressively form and increase in amplitude. This behaviour is associated with an enhancement of far SOL turbulent transport in the form of blob-filaments travelling radially faster across the far SOL and becoming bigger at higher α t . Finally, the detected filaments, evaluated from the fast reciprocating probe at the outer midplane, are determined to mostly belong to the resistive ballooning and resistive X-point regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2024

Azimuthal Size Scales of Solar Wind Periodic Density Structures

Periodic Density Structures (PDSs) are quasi-periodic variations of solar wind density ranging from a few minutes to a few hours. PDSs advect with the solar wind and have radial length scales (Lx) of tens to several thousand Mm, thus belonging to the class of “mesoscale structures”. Current interplanetary multi-spacecraft observations are not at spatial separations capable of directly measuring the 3D size scale of PDSs or other mesoscale structures. Instead, previous investigations estimated characteristic spatial scales in solar wind parameters using cross-correlation and/or coherence analysis applied to multi-spacecraft observations. For solar wind density and interplanetary magnetic field (IMF) intensity, the reported size scales perpendicular to the Sun-Earth line (Ly) ranged between ≈30 and ≈200 Earth Radii (RE). Here, we implemented a similar approach on the same parameters but focused on high density slow solar wind intervals with PDSs observed by the Wind and ARTEMIS-P1 spacecraft. Additionally, this is the first statistical study of the IMF intensity periodicities in relation to PDSs. We identified intervals in which the two spacecraft observed the same periodicity obtaining two PDS groups based on their radial length scale: Lx1 ≈86RE and Lx2 ≈35RE. Then, we classified the events based on the periodic variations’ coherence level. Reproducing the results with simulations of the PDSs’ transit, we inferred the Ly order of magnitudes for the two PDS groups: Ly1 ≈340RE and Ly2 ≈187RE. Knowing the PDSs’ size scales is fundamental to constrain models aimed at reproducing these structures and is critical to better understand the PDS-magnetosphere coupling.

Interplanetary physics↗
NASA NTRS1995

Numerical Simulations of High Enthalpy Pulse Facilities

Axisymmetric flows within shock tubes and expansion tubes are simulated including the effects of finite rate chemistry and both laminar and turbulent boundary layers. The simulations demonstrate the usefulness of computational fluid dynamics for characterizing the flows in high enthalpy pulse facilities. The modeling and numerical requirements necessary to simulate these flows accurately are also discussed. Although there is a large body of analysis which explains and quantifies the boundary layer growth between the shock and the interface in a shock tube, there is a need for more detailed solutions. Phenomena such as thermochemical nonequilibrium. or turbulent transition behind the shock are excluded in the assumptions of Mirels' analysis. Additionally there is inadequate capability to predict the influence of the boundary layer on the expanded gas behind the interface. Quantifying the gas in this region is particularly important in expansion tubes because it is the location of the test gas. Unsteady simulations of the viscous flow in shock tubes are computationally expensive because they must follow features such as a shock wave over the length of the facility and simultaneously resolve the small length scales within the boundary layer. As a result, efficient numerical algorithms are required. The numerical approach of the present work is to solve the axisymmetric gas dynamic equations using an finite-volume formulation where the inviscid fluxes are computed with a upwind TVD scheme. Multiple species equations are included in the formulation so that finite-rate chemistry can be modeled. The simulations cluster grid points at the shock and interface and translate this clustered grid with these features to minimize numerical errors. The solutions are advanced at a CFL number of less than one based on the inviscid gas dynamics. To avoid limitations on the time step due to the viscous terms, these terms are treated implicitly. This requires a block tri-diagonal matrix inversion along each line of cells normal to the wall. The cost of this inversion is more than offset by the larger allowable time step. The source terms representing the finite-rate chemical kinetics are also treated implicitly. An algebraic turbulence model for compressible flow is used. The flow in a low pressure shock tube is computed and the results are compared with Mirels'analysis. The driven gas is nitrogen at 70 Pa, and the incident shock speed is approximately 2.9 km/sec so that there is little dissociation. The simulations include a laminar boundary layer and are run until the limiting flow regime is achieved. At this limit, the shock and interface travel at the same velocity because the amount of driven gas between these two features remains the same: the mass flow across the shock is equal to the mass of gas being entrained at the interface by the boundary layer. Simulations with several grids are presented to establish the grid independence of the solution, Good agreement is achieved between Mirels' correlations and the computations. This is expected since the flow conditions are chosen to be consistent with the assumptions used in Mirels' analysis. This comparison adds credibility to the numerical approach and highlights some of the differences between the theory and the detailed simulations. In addition, simulations of the HYPULSE expansion tube are presented for two operating conditions and the computations are compared to experimental data. The operating gas for both cases is nitrogen. One test condition is at a total enthalpy of 15.2 MJ/Kg and a relatively low pressure of 2 kPa. This case is characterized by a laminar boundary layer and significant chemical nonequilibrium. in the acceleration gas. The second test condition is at a total enthalpy of 10.2 MJ/Kg and a pressure of 38 kPa and is characterized by a turbulent boundary layer. The simulations compare well with experiment and reveal that the nonuniformity in pressure observed during the test time is related to variations in the boundary layer displacement thickness.

Wilson, Gregory J.↗