Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “modulational instability”

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.

431 records · Page 24
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Mesosphere Dynamics with Gravity Wave Forcing: Planetary Waves - 2

We present results from a non-linear, 3D, time dependent numerical spectral model (NSM) which extends from the ground up into the thermosphere and incorporates Hines' Doppler Spread Parameterization for small-scale gravity waves (GW). Our focal point is the mesosphere where wave interactions are playing a dominant role. We discuss planetary waves in the present paper and diurnal and semi-diurnal tides in the companion paper. Without external time dependent energy or momentum sources, planetary waves (PWs) are generated in the model for zonal wavenumbers 1 to 4, which have amplitudes in the mesosphere above 50 km as large as 30 m/s and periods between 2 and 50 days. The waves are generated primarily during solstice conditions, which indicates that the baroclinic instability (associated with the GW driven reversal in the latitudinal temperature gradient) is playing an important role. Results from a numerical experiment show that GWs are also involved directly in generating the PWs. For the zonal wavenumber m = 1, the predominant wave periods in summer are around 4 days and in winter between 6 and 10 days. For m = 2, the periods are in summer and close to 2.5 and 3.5 days respectively For m = 3, 4 the predominant wave periods are in both seasons close to two days. The latter waves have the characteristics of Rossby gravity waves with meridional winds at equatorial latitudes. A common feature of the PWs (m = 1 to 4) generated in summer and winter is that their vertical wavelengths throughout the mesosphere are large which indicates that the waves are not propagating freely but are generated throughout the region. Another common feature is that the PWs propagate preferentially westward in summer and eastward in winter, being launched from the westward and eastward zonal winds that prevail respectively in summer and winter altitudes below 80 km. During spring and fall, for m = 1 and 2 eastward propagating long period PWs are generated that are launched from the smaller eastward zonal winds that prevail in these seasons. PWs generated in the model produce large amplitude modulations of the diurnal tides at altitudes above 80 km and contribute to their seasonal variations.

Mayr, H. G.↗
DOE OSTI2022

Albedo susceptibility of northeastern Pacific stratocumulus: the role of covarying meteorological conditions

Abstract. Quantification of the radiative adjustment of marine low clouds to aerosol perturbations, regionally and globally, remains the largest source of uncertainty in assessing current and future climate. One of the important steps towards quantifying the role of aerosol in modifying cloud radiative properties is to quantify the susceptibility of cloud albedo and liquid water path (LWP) to perturbations in cloud droplet number concentration (Nd). We use 10 years of spaceborne observations from the polar-orbiting Aqua satellite to quantify the albedo susceptibility of marine low clouds to Nd perturbations over the northeast (NE) Pacific stratocumulus (Sc) region. Mutual information analysis reveals a dominating control of cloud state (e.g., LWP and Nd) on low-cloud albedo susceptibility, relative to the meteorological states that drive these cloud states. Through a LWP–Nd space decomposition of albedo susceptibilities, we show clear separation among susceptibility regimes (brightening or darkening), consistent with previously established mechanisms through which aerosol modulates cloud properties. These regimes include (i) thin non-precipitating clouds (LWP < 55 g m−2) that exhibit brightening (occurring 37 % of the time), corresponding to the Twomey effect; (ii) thicker non-precipitating clouds, corresponding to entrainment-driven negative LWP adjustments that manifest as a darkening regime (36 % of the time); and (iii) another brightening regime (22 % of the time) consisting of mostly precipitating clouds, corresponding to precipitation-suppression LWP positive adjustments. Overall, we find an annual-mean regional low-cloud brightening potential of 20.8±2.68 W m−2 ln(Nd)−1, despite an overall negative LWP adjustment for non-precipitating marine stratocumulus, owing to the high occurrence of the Twomey–brightening regime. Over the NE Pacific, clear seasonal covariabilities among meteorological factors related to the large-scale circulation are found to play an important role in grouping conditions favorable for each susceptibility regime. When considering the covarying meteorological conditions, our results indicate that for the northeastern Pacific stratocumulus, clouds that exhibit the strongest brightening potential occur most frequently within shallow marine boundary layers over a cool ocean surface with a stable atmosphere and a dry free troposphere above. Clouds that exhibit a darkening potential associated with negative LWP adjustments occur most frequently within deep marine boundary layers in which the atmospheric instability and the ocean surface are not strong and warm enough to produce frequent precipitation. Cloud brightening associated with warm-rain suppression is found to preferably occur either under unstable atmospheric conditions or humid free-tropospheric conditions that co-occur with a warm ocean surface.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework [1, 2, 3]. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick [4, 5]. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, no slip boundary conditions, and the hybrid outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, thrust, ect.). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the Acoustic Velocity Potential Equation (AVPE). The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob [2, 3]. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, Sean R.↗
NASA NTRS2023

Heat Transfer and Interfacial Flow Physics of Microgravity Flow Boiling in Single-Side-Heated Rectangular Channel with Subcooled Inlet Conditions – Experiments Onboard the International Space Station

This study is the culmination of a long-term collaborative effort between researchers from the Purdue University Boiling and Two-Phase Flow Laboratory (PU-BTPFL) and the NASA Glenn Research Center to investigate gravitational effects on flow boiling and flow condensation. The science and design concepts for this large-scale effort were initiated in 2011 and included several studies detailing various aspects of two-phase fluid physics in both Earth gravity and microgravity, culminating in construction of the large-scale experimental facility named “Flow Boiling and Condensation Experiment (FBCE)”. The experiment was launched to the International Space Station (ISS) in August 2021. Following the successful installation of FBCE, equipped with the Flow Boiling Module (FBM), onboard the ISS and completion of several safety checks, flow boiling experiments were performed for five months from February 2022 until July 2022. This resulted in a large flow boiling database covering broad ranges of operating parameters and heating configurations spanning several research objectives. This study investigates microgravity flow boiling of n-perfluorohexane with subcooled inlet in a single-side-heated rectangular channel of dimensions 114.6-mm heated length, 2.5-mm heated width, and 5.0-mm height. Key operating parameters investigated are mass velocity (199.90 – 3200.13 kg/m 2 s), inlet subcooling (0.10 – 45.76°C), and inlet pressure (113.30 – 164.29 kPa). Images and image sequences acquired via high-speed-video are presented to elucidate the interfacial flow physics. To analyze and explain the effects of various parameters in microgravity, heat transfer results are presented as flow boiling curves, streamwise profiles of wall temperature and heat transfer coefficient, and parametric trends of local and averaged heat transfer coefficient. Mass velocity and inlet subcooling significantly influenced most of the aforementioned aspects of flow boiling, whereas effects of inlet pressure were comparatively insignificant. Although the data and observed flow physics might be different, the parametric effects and trends in microgravity are similar to vertical upflow in Earth gravity. Some cases, especially low mass velocities, high heat fluxes, and large degrees of inlet subcooling, experienced temporally anomalous flow behaviors caused by two-phase flow instabilities manifesting as flow reversals and resulted in deviations in overall trends. Severe thermodynamic non-equilibrium is observed throughout the channel. Overall, FBCE’s ISS experiments were successful for subcooled inlet with single-sided heating of rectangular channel, and the collected data well established the various effects on flow boiling physics in highly controlled long-term microgravity conditions.

flow boiling↗
DOE OSTI2022

Extension of the flux fit method for estimating power deposition profiles

The flux fit method is used to self-consistently estimate the power deposition profile and heat transport profiles from temperature measurements originating from perturbative experiments with a modulated source. Furthermore, this letter improves on this method by addressing the limitations and assumptions. The most crucial improvement is the additional freedom in the source deposition profile. Allowing for a variable central deposition location and height and including a skewness parameter produces deposition profiles more consistent with the measurement data, but still wider than equilibrium ray tracing in two different DIII-D discharges. Moreover, we show that the quality of the estimated deposition profile is key to the accuracy of diffusivity and connectivity estimates, but inversely, the estimated transport parameters hardly affect the quality of the power deposition estimate. Using this method, we show that the power deposition profile estimate is broadened with respect to ray-tracing by about 1.7-1.8 times in two DIII-D discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS2023

Effects of Heating Configuration and Operating Parameters on Heat Transfer and Interfacial Physics of Microgravity Flow Boiling With Subcooled Inlet Conditions –Experiments Onboard the International Space Station

This study is part of the Flow Boiling and Condensation Experiment (FBCE), a collaborative effort between the Purdue University Boiling and Two-Phase Flow Laboratory (PU-BTPFL) and the NASA Glenn Research Center. The FBCE fitted with the Flow Boiling Module (FBM) was launched to the International Space Station (ISS) in August 2021 and experiments were successfully performed from February to July 2022 to amass a large microgravity-flow-boiling database. This study is focused on heat transfer and flow visualization of microgravity flow boiling of n-Perfluorohexane in a rectangular channel of 5.0 mm height, 2.5 mm width (heated), and 114.6 mm length, with subcooled inlet conditions. High-speed-video photography is utilized to present flow patterns and temporal interfacial behavior. Heat transfer results are presented in the form of flow boiling curves and both parametric curves and streamwise profiles of wall temperature and heat transfer coefficient. Firstly, the parametric effects of mass velocity (199.4 – 3200.0 kg/m 2 s), inlet subcooling (0.2 – 46.0°C), and inlet pressure (124.2 – 176.7 kPa), on the aforementioned aspects are assessed for double-sided heating to establish them for a microgravity environment. Of these three parameters, mass velocity and inlet subcooling mostly determine the microgravity flow boiling behavior, while inlet pressure plays an insignificant role. Flow patterns for double-sided heating are more complex than those for single-sided heating due to interaction between the two vapor layers. Vapor interaction is minimized at high subcoolings and high mass velocities due to strong condensation offered by the subcooled bulk liquid layer separating them. Despite the different flow patterns, both single- and double-sided heating generally result in similar parametric trends and local heat transfer coefficients for similar operating conditions. Flow instabilities manifest as temporal flow anomalies and temperature oscillations, and their severity increases with increasing boiling number. Secondly, the effects of heating configuration are analyzed by comparing and contrasting several aspects of single- and double-sided heating data. The heat fluxes at which onset of nucleate boiling degradation (ONBD) and critical heat flux (CHF) occur are distinctly different for single- and double-sided heating. There exists a threshold inlet subcooling demarcating the dominance of flow acceleration and condensation effects in vapor removal from the near-wall region and replenishment of fresh liquid for boiling. Above the threshold, condensation from the near-wall region is dominant and single-sided heating yields higher heat fluxes, and below it, acceleration is dominant and double-sided yields higher heat fluxes. At mass velocity in the range of 200 – 2400 kg/m 2 s, the threshold inlet subcooling lies in the approximate range of 20 – 30°C (corresponding inlet quality of roughly -0.40 – -0.20).

Microgravity↗
NASA NTRS2007

Inter-annual and Long-term Temperature Variations in the Mesopause Region at High Latitudes Generated by the Stratospheric QBO

The Numerical Spectral Model (NSM) simulates the Quasi-biennial Oscillation (QBO) that dominates the zonal circulation of the lower stratosphere at low latitudes. In the model, the QBO is generated with parameterized small-scale gravity waves (GW), which are partially augmented in 3D with planetary waves owing to baroclinic instability. Due to GW filtering, the QBO extends into the upper mesosphere, evident in UARS zonal wind and TIMED temperature measurements. While the QBO zonal winds are confined to equatorial latitudes, even in simulations with latitude-independent wave source, the associated temperature variations extend to high latitudes. The meridional circulation redistributes some of the QBO energy to focus it partially onto the Polar Regions. The resulting QBO temperature variations away from the equator tend to increase at higher altitudes to produce inter-annual variations that can exceed 5 K in the polar mesopause region -- and our 3D model simulations show that the effect is variable from year to year and can produce large differences between the two hemispheres, presumably due to interactions involving the seasonal variations. Modeling studies with the NSM have shown that long-term variations can also be generated by the QBO interacting with the seasonal cycles through OW node-filtering. A 30-month QBO, optimally synchronized by the 6-month Semi-Annual Oscillation (SAO), thus produces a 5-year or semi-decadal (SD) oscillation -- and observational evidence for that has been provided by a recent analysis of stratospheric NCEP data. In a simulation with the 2D version of the NSM, this SD oscillation extends into the upper mesosphere, and we present results to show that the related temperature variations could contribute significantly to the long-term variations of the polar mesopause region. Quasi-decadal variations could furthermore arise from the modeled solar cycle modulations of the QBO and 12-month annual oscillation. Our numerical results are discussed in the context of the observed low summer temperatures reproduced by the model, to demonstrate that the above interannual and long-term variations could contribute significantly to the climatology of Polar Mesospheric Clouds (PMC) investigated by the Aeronomy of Ice in the Mesosphere (AIM) mission.

Mayr, Hans G.↗
DOE OSTI2020

Theory of edge localized mode suppression by static resonant magnetic perturbations in the DIII-D tokamak

According to a recent paper [Hu et al ., Phys. Plasmas 26, 120702 (2019)], mode penetration at the top of the pedestal is a necessary and sufficient condition for the suppression of edge localized modes (ELMs) by a resonant magnetic perturbation (RMPs) in an H-mode tokamak discharge. This paper employs asymptotic matching theory to model a particular DIII-D discharge in which ELMs were suppressed by an externally generated, static, n ¼ 2, RMP whose amplitude was modulated at a frequency of 1 Hz. It is demonstrated that the response of the plasma to the applied RMP, in the immediate vicinities of the rational (i.e., resonant) surfaces, is governed by nonlinear, rather than by linear, physics. This is the case because the magnetic island widths associated with driven reconnection exceed the linear layer widths, even in cases where driven reconnection is strongly suppressed by plasma rotation. The natural frequency at a given rational surface (i.e., the helical frequency at which the locally resonant component of the RMP would need to propagate in order to maximize driven reconnection) is found to be offset from the local E B frame in the ion diamagnetic direction. The size of the offset is mostly determined by neoclassical poloidal rotation. Finally, the predictions of a fully nonlinear plasma response model are found to be broadly consistent with the DIII-D experimental data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

PPPL Report on Reduced Modeling of Fusion Alpha Transport in ARC Burning Plasmas

We are reporting on the modeling of fusion alpha particle transport in the planned ARC fusion device being designed by the CFS (Commonwealth Fusion Systems: https://cfs.energy). The ARC tokamak is designed to operate in a burning-plasma regime characterized by a substantial population of fusion-born alpha particles. Alfvén eigenmode (AE) stability is assessed both analytically and numerically, incorporating alpha-particle drive, ion Landau, and radiative damping from thermal species and collisional damping from trapped electrons. Regions of unstable and near-threshold AE activity are mapped across ARC’s operational parameter space. Linear stability analysis with NOVA indicates multiple, often marginally unstable AEs, extending to toroidal mode numbers up to n= 30. The present report focuses on the ARC flat-top operating point prior to the sawtooth event. Alpha-particle transport on timescales exceeding the neoclassical slowing-down time is assessed using the NUBEAM module [1][2] of the TRANSP code [3], employing transport coefficients derived from the RBQ quasilinear modeling (cf. Appendix B). These global simulations identify favorable and unfavorable operating regimes with respect to alpha confinement, pressure redistribution, and overall alpha-heating efficiency. We also evaluate additional transport mechanisms—including neoclassical tearing mode (TM)–induced stochasticity, sawtooth-driven redistribution, and toroidal-field ripple using the kick model (cf. Appendix C) which makes use of the guiding-center code ORBIT, see Section 5. The kick model is integrated into TRANSP to enable self-consistent predictions of alpha-driven current formation and sustainment within the ARC scenario. Sensitivity scans are performed over the mode frequency, rational-surface alignment, island width, mode amplitude, and proximity of the limiter to the plasma. Our study provides an initial, physics-based guidance for machine design, operational planning, and equilibrium control, ensuring adequate alpha confinement and robust self-heating performance in ARC. Our simulations mostly targeted worst case scenarios, e.g. for TMs and sawteeth. Overall, we expect benign effects for the ARC scenario investigated in this work on fusion alpha confinement and losses in the presence of AEs, tearing modes and sawteeth. This report addresses three thrusts identified at the outset. The first thrust focuses on analytic estimates of the parametric dependencies of EP relaxation based on local AE stability simulations (Section 3). The second thrust involves global evaluations of AE stability using the NOVA, RBQ, and NUBEAM codes (Section 4). Finally, we investigate alpha-particle transport driven by low-frequency instabilities associated with sawteeth and tearing modes (Section 5).

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Electron Temperature Gradient Driven Transport Model for Tokamak Plasmas

A new model for electron temperature gradient (ETG) modes is developed as a component of the Multi-Mode anomalous transport module [T. Rafiq \textit{et al.,} Phys Plasmas \textbf{20}, 032506 (2013)] to predict a time dependent electron temperature profile in conventional and low aspect ratio tokamaks. This model is based on two-fluid equations that govern the dynamics of low-frequency short- and long-wavelength electromagnetic toroidal ETG driven drift modes. A low collisionality NSTX discharge is used to scan the plasma parameter dependence on the ETG real frequency, growth rate, and electron thermal diffusivity. Electron thermal transport is discovered in the deep core region where modes are more electromagnetic in nature. Several previously reported gyrokinetic trends are reproduced, including the dependencies of density gradients, magnetic shear, $\beta$ and gradient of $\beta$ $(\betap)$, collisionality, safety factor, and toroidicity, where $\beta$ is the ratio of plasma pressure to the magnetic pressure. The electron heat diffusivity associated with the ETG mode is discovered to be on a scale consistent with the experimental diffusivity determined by power balance analysis.

ETG↗
DOE OSTI2023

Theoretical investigation of braking of tearing mode rotation by resistive walls in ITER

The locking of the 2/1 tearing mode to the resistive wall in the ITER tokamak (15 MA inductive scenario 2) is investigated theoretically using a cylindrical asymptotic matching model. The model takes into account the fact that ITER plasmas will effectively be surrounded by two walls; the inner blanket module layer with a time constant of about 23 ms, and the outer vacuum vessel with a time constant of about 380 ms. The model also takes cognizance of the fact that neither the blanket module layer nor the vacuum vessel can be accurately described as “thin” walls (in the ordinarily accepted sense). The model incorporates changes in both the plasma poloidal and the toroidal angular velocity profiles, in response to the electromagnetic braking torque that develops at the rational surface, because it turns out that neoclassical poloidal flow-damping is not strong enough to completely suppress changes in the poloidal velocity. Finally, the model accurately calculates changes in the poloidal and toroidal plasma angular velocity profiles by evolving the full angular equations of motion, taking the electromagnetic braking torque, plasma inertia, plasma viscosity, and poloidal flow-damping into account. The time required for the 2/1 tearing mode to grow from a small amplitude to a sufficient one to lock to the walls is found to be about 3.5 s. As a result, the critical full radial island width at which wall locking is triggered is found to be about 9% of the plasma minor radius.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Resolving ECRH deposition broadening due to edge turbulence in DIII-D

Microwave heat pulse propagation experiments have demonstrated a correlation between millimeter-scale turbulence and deposition profile broadening of electron cyclotron (EC) waves on the DIII-D tokamak. In a set of discharges in DIII-D, a variation in edge density fluctuations on the mm-scale is associated with 40%–150% broader deposition profiles, expressed in terms of normalized minor radius, as compared with equilibrium ray tracing. The 1D power profile is determined from transport analysis of the electron temperature response to EC power modulation using perturbative analysis with a square wave power modulation at 20–70 Hz, producing a series of Fourier harmonics that are fit collectively to resolve transport. Fitting an integrated heat flux expressed in the Fourier basis of the modulation to diffusive, convective, and coupled transport terms in a linear model can resolve the broadened EC deposition width from the power perturbation to resolve a broadening in each case. The best fit degree of beam broadening observed scales approximately linearly with the Doppler backscattering measured fluctuation level in the steep gradient region. Quantifying the effect of edge fluctuation broadening on EC current drive power needs of future devices will require 3D full-wave codes that can be validated on the current generation of machines. These DIII-D experiments provide a quantitative measure of fluctuation effects and a dataset to benchmark full-wave simulations that can model and eventually predict nonlinear effects neglected by 1D equilibrium beam and ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Patterns and rates of soil movement and shallow failures across several small watersheds on the Seward Peninsula, Alaska

Abstract. Thawing permafrost can alter topography, ecosystems, and sediment and carbon fluxes, but predicting landscape evolution of permafrost-influenced watersheds in response to warming and/or hydrological changes remains an unsolved challenge. Sediment flux and slope instability in sloping saturated soils have been commonly predicted from topographic metrics (e.g., slope, drainage area). In addition to topographic factors, cohesion imparted by soil and vegetation and melting ground ice may also control spatial trends in slope stability, but the distribution of ground ice is poorly constrained and hard to predict. To address whether slope stability and surface displacements follow topographic-based predictions, we document recent drivers of permafrost sediment flux present on a landscape in western Alaska that include creep, solifluction, gullying, and catastrophic hillslope failures ranging in size from a few meters to tens of meters, and we find evidence of rapid and substantial landscape change on an annual timescale. We quantify the timing and rate of surface movements using a multi-pronged, multi-scalar dataset including aerial surveys, interannual GPS surveys, synthetic aperture radar interferometry (InSAR), and climate data. Despite clear visual evidence of downslope soil transport of solifluction lobes, we find that the interannual downslope surface displacement of these features does not outpace downslope displacement of soil in locations where lobes are absent (downslope movement means: 7 cm yr−1 for lobes over 2 years vs. 10 cm yr−1 in landscape positions without lobes over 1 year). Annual displacements do not appear related to slope, drainage area, or modeled total solar radiation but are likely related to soil thickness, and volumetric sediment fluxes are high compared to temperate landscapes of comparable bedrock lithology. Time series of InSAR displacements show accelerated movement in late summer, associated with intense rainfall and/or deep thaw. While mapped slope failures do cluster at slope–area thresholds, a simple slope stability model driven with hydraulic conductivities representative of throughflow in mineral and organic soil drastically overpredicts the occurrence of slope failures. This mismatch implies permafrost hillslopes have unaccounted-for cohesion and/or throughflow pathways, perhaps modulated by vegetation, which stabilize slopes against high rainfall. Our results highlight the breadth and complexity of soil transport processes in Arctic landscapes and demonstrate the utility of using a range of synergistic data collection methods to observe multiple scales of landscape change, which can aid in predicting periglacial landscape evolution.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS2015

Solid Rocket Motor Combustion Instability Modeling in COMSOL Multiphysics

Combustion instability modeling of Solid Rocket Motors (SRM) remains a topic of active research. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process, acoustics, and steady-state gas dynamics. Recent advances in defining the energy transport of disturbances within steady flow-fields have been applied by combustion stability modelers to improve the analysis framework. Employing this more accurate global energy balance requires a higher fidelity model of the SRM flow-field and acoustic mode shapes. The current industry standard analysis tool utilizes a one dimensional analysis of the time dependent fluid dynamics along with a quasi-three dimensional propellant grain regression model to determine the SRM ballistics. The code then couples with another application that calculates the eigenvalues of the one dimensional homogenous wave equation. The mean flow parameters and acoustic normal modes are coupled to evaluate the stability theory developed and popularized by Culick. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The current study employs the COMSOL Multiphysics finite element framework to model the steady flow-field parameters and acoustic normal modes of a generic SRM. This work builds upon previous efforts to verify the use of the acoustic velocity potential equation (AVPE) laid out by Campos. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, del squared psi - (lambda/c) squared psi - M x [M x del((del)(psi))] - 2((lambda)(M)/c + M x del(M) x (del)(psi) - 2(lambda)(psi)[M x del(1/c)] = 0. with M as the Mach vector, c as the speed of sound, and λ as the complex eigenvalue. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluates the gas flow inside of a SRM using St. Robert's law to model the solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are verified by comparing the pertinent ballistics parameters with the industry standard code outputs (i.e. pressure drop, axial velocity, exit velocity). These results are then used by the coefficient form of the mathematics module to determine the complex eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. The one way coupled analysis is perform four times utilizing geometries determined through traditional SRM modeling procedures. The results of the steady-state CFD and AVPE analyses are used to calculate the linear acoustic growth rate as is defined by Flandro and Jacob. In order to verify the process implemented within COMSOL we first employ the Culick theory and compare the results with the industry standard. After the process is verified, the Flandro/Jacob energy balance theory is employed and results displayed.

Fischbach, S. R.↗
DOE OSTI2020

Proton Bunch Self-Modulation in Plasma with Density Gradient

We study experimentally the effect of linear plasma density gradients on the self-modulation of a 400 GeV proton bunch. Results show that a positive or negative gradient increases or decreases the number of microbunches and the relative charge per microbunch observed after 10 m of plasma. The measured modulation frequency also increases or decreases. With the largest positive gradient we observe two frequencies in the modulation power spectrum. Results are consistent with changes in wakefields’ phase velocity due to plasma density gradients adding to the slow wakefields’ phase velocity during self-modulation growth predicted by linear theory.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Simulation and experimental study of proton bunch self-modulation in plasma with linear density gradients

We present numerical simulations and experimental results of the self-modulation of a long proton bunch in a plasma with linear density gradients along the beam path. Simulation results agree with the experimental results reported [F. Braunmller, T. Nechaeva et al. (AWAKE Collaboration), Phys. Rev. Lett. 125, 264801 (2020)]: with negative gradients, the charge of the modulated bunch is lower than with positive gradients. In addition, the bunch modulation frequency varies with gradient. Simulation results show that dephasing of the wakefields with respect to the relativistic protons along the plasma is the main cause for the loss of charge. The study of the modulation frequency reveals details about the evolution of the self-modulation process along the plasma. In particular for negative gradients, the modulation frequency across time-resolved images of the bunch indicates the position along the plasma where protons leave the wakefields. Simulations and experimental results are in excellent agreement.

43 PARTICLE ACCELERATORS↗
DOE OSTI2024

Updraft Width Modulates Ambient Atmospheric Controls on Convective Cloud Depth

Abstract The depth of convective clouds affects vertical transport of atmospheric constituents, influencing downstream weather and climate. Atmospheric controls on the maximum depth reached by moist convection are investigated with radar‐tracked convective cells tagged with sounding‐derived atmospheric parameters from a field campaign in central Argentina. Regression analyses show that narrow (<12‐km diameter) and wide (>16‐km diameter) cell depths respond to disparate factors, where cell areas are defined using composite reflectivity signatures. Undiluted lifted parcel indices including convective available potential energy (CAPE) and level of neutral buoyancy (LNB) are top predictors of wide cell maximum depth while mid‐tropospheric relative humidity is the top predictor of narrow cell maximum depth. Because narrow cells are more numerous than wide cells, the overall outcome of the full cell population does not strongly correlate with CAPE and LNB conditions. Tracked cells and atmospheric conditions in a simulation with 3‐km grid spacing covering the field campaign produce similar results to those observed. Narrow cells that are relatively deep have a cooler and moister mid‐troposphere with weaker free tropospheric subsidence, while relatively deep wide cells have much warmer and moister lower tropospheric conditions. These atmospheric differences are present 1 hr before cell initiation at both a fixed observing site and variable cell initiation locations. Simulated narrow cell maximum equivalent potential temperature decreases with height at a rate similar to the ambient vertical gradient, causing these cells to fall short of their LNB and supporting the view that entrainment‐driven dilution is a dominant control on their depth.

54 ENVIRONMENTAL SCIENCES↗