Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “unstable waves”

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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Data-driven Detection and Early Prediction of Thermoacoustic Instability in a Multi-nozzle Combustor

Thermoacoustic instability (TAI) is a critical issue in modern lean-burn gas-turbine combustors, which is induced by a strong coupling between the resonant combustor acoustics and fluctuations in the heat release rate. This instability may lead to high-amplitude pressure waves that generate undesirable noise levels as well as fatigue stresses in mechanical structures of the combustor. The intense pressure fluctuations due to TAI may also cause large flow perturbations and possibly flow reversal that may lead to flame oscillations, flame liftoff, and even flame blow-out. Hence, there is a strong need for exercising control actions in a timely fashion to mitigate the TAI phenomena. Anomaly detection is an essential prerequisite to the design of a good controller and such a detector must be able to reliably predict a forthcoming TAI. To detect and predict the onset of a TAI from an ensemble of pressure time series, this paper investigates three data-driven methods: Fast Fourier transform (FFT), symbolic time series analysis (STSA), and hidden Markov modeling (HMM). The main focus of the paper is to make a comparative evaluation of these three anomaly detection methods for classification of the current regime of operation into stable and unstable categories as well as for real-time identification of precursors to impending instabilities with short-length time series of measured variables (e.g., pressure oscillations). Finally, the results, generated on experimental data from a multi-nozzle combustor apparatus, have been compared to evaluate the performance of FFT, STSA, and HMM methods for TAI analysis.

42 ENGINEERING↗
DOE OSTI2023

Suppressed superexchange interactions in the cuprates by bond-stretching oxygen phonons

We study a multiorbital Hubbard–Su-Schrieffer-Heeger model for the one-dimensional (1D) corner-shared cuprates in the adiabatic and nonadiabatic limits using exact diagonalization and determinant quantum Monte Carlo. At half filling and in the adiabatic limit, lattice dimerization can be achieved only over a narrow range of couplings slightly below a critical coupling g c . Beyond this value, the sign of the effective hopping changes, and the lattice becomes unstable. Strong lattice fluctuations replace the dimerization state in the nonadiabatic case. We also examine the model's temperature-dependent uniform and dynamical magnetic susceptibilities and compare them to the results of an effective spin-1/2 Heisenberg model. In doing so, we demonstrate that lattice fluctuations induced by the e-ph interaction suppress the effective superexchange interaction when g < g c . Furthermore, our results elucidate the effect of bond-stretching phonons in the parent cuprate compounds in general and are particularly relevant to 1D cuprates, where strong e-ph interactions have recently been inferred.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Competing orders at higher-order Van Hove points

Van Hove points are special points in the energy dispersion, where the density of states exhibits analytic singularities. When a Van Hove point is close to the Fermi level, tendencies towards density wave orders, Pomeranchuk orders, and superconductivity can all be enhanced, often in more than one channel, leading to a competition between different orders and unconventional ground states. In this study we consider the effects from higher-order Van Hove points, around which the dispersion is flatter than near a conventional Van Hove point, and the density of states has a power-law divergence. We argue that such points are present in intercalated graphene and other materials. We use an effective low-energy model for electrons near higher-order Van Hove points and analyze the competition between different ordering tendencies using an unbiased renormalization-group approach. For purely repulsive interactions, we find that two key competitors are ferromagnetism and chiral superconductivity. For a small attractive interaction, we find an unconventional spin Pomeranchuk order, in wich the spin oder parameter winds around the Fermi surface. The supermetal state, predicted for a single higher-order Van Hove point, is an unstable fixed point in our case.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Excitation of fishbone-like mode in tokamaks due to bounce resonances of trapped thermal ions

When the drift kinetic effect of thermal ions is taken into account, a high-frequency fishbone-like mode (FLM) is found to be driven unstable by trapped thermal ions (TTIs) in tokamak plasmas, according to self-consistent magneto-hydrodynamic (MHD)-kinetic hybrid modeling utilizing the MARS-K code [Liu et al., Phys. Plasmas 15, 112503 (2008)] as well as an analytic theory. It is found that, similar to energetic particles, TTIs can also stabilize the internal kink mode, whereas the FLM is excited when the effective beta of TTIs exceeds a threshold value. The real frequency of the FLM is comparable to the bounce frequency of TTIs. The mode structure of the FLM can be significantly different from the conventional step-like function for the associated plasma radial displacement. This drift kinetic induced modification of the mode structure near the q = 1 surface is captured by non-perturbative MHD-kinetic hybrid computations with MARS-K. Furthermore, the FLM can only be triggered by TTIs at sufficiently high thermal temperatures. Both the FLM and the internal kink can be stabilized by sufficiently fast plasma toroidal rotation and parallel sound wave damping. These two conditions of high thermal temperature and (fast) flow stabilization, though making it challenging to observe the TTI-driven FLM in present day experiments, are favorable for the mode excitation in future reactor scale devices.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Deviations from the Isobaric Multiplet Mass Equation due to Threshold States

Recent studies have completed the A=16 isospin quintets for states with J π =0 + and 2 + . The dependence of their masses as a function of isospin projection shows evidence for deviations from quadratic behavior indicating isospin violation beyond the expectation from two-body forces. The deviation is most pronounced for the 2 + states. Predictions from the shell model embedded in the continuum (SMEC) allow us to explain that this isospin violation is associated with a modification of the nuclear structure due to the open-quantum-system nature of the proton-rich members of the quintet. In particular, the 0 + and 2 + states in 16 Ne and the 2 + state in 16 F are threshold resonances located just above a proton-decay threshold where s-wave coupling to the continuum is expected. The measured deviations of these threshold states from the quadratic behavior of the remaining members of the multiplets makes it possible to obtain information on the magnitude and the energy dependence of the continuum-coupling energy correction. Finally, continuum coupling is also indicated for the ground state of 8 C, but this time through p-wave coupling

coulomb energies & analogue states↗
DOE OSTI2022

Design of a high energy density experiment to measure the suppression of hydrodynamic instability in an applied magnetic field

A transverse magnetic field can suppress the hydrodynamic instability growth of an unstable plasma interface. This effect, of interest to inertial confinement fusion and astrophysics, has mostly been studied via simulation. Here, we present the design of an experiment at the National Ignition Facility (NIF) to demonstrate this effect in a laboratory. Simulations indicate that the timescale for the diffusion of the magnetic field across the mixing region should be at least comparable to the timescale of the instability growth in order to have a measurable suppression effect. This motivates the use of lower density target materials than usual high energy density (HED) hydrodynamics experiments to permit faster hydrodynamics and higher plasma conductivities (through higher temperature), for a given laser drive and magnetic field. We discuss a target design for creating a Rayleigh–Taylor unstable HED plasma interface that uses 320 mg/cc iodine-doped carbon foam as the heavy material, 20 mg/cc carbon foam as the light material, and a 6 lm amplitude, a 120 lm wavelength ripple machined at the interface, which shows a measurable suppression effect with a nominal NIF drive and a 30 T magnetic field (the present facility limit). Here, models indicate lower density foams that may display even larger suppression effects, as the Hohlraum drive also radiatively preheats the foam to permit even higher temperatures and, hence, higher conductivities.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Neutron-producing gas puff Z-pinch experiments on a fast, low-impedance, 0.5 MA linear transformer driver

A study on the neutron production from single and double gas puff Z-pinches on the CESZAR linear transformer driver with ~0.45 MA current and 170 ns rise time is presented. Total neutron yield measurements made with a LaBr activation detector are compared for three configurations, using a double nozzle setup. When a single, hollow, deuterium gas shell was used, reliable implosions could only be attained at higher load mass than the optimal value to match implosion time with the driver rise time, with neutron yields of ~10 6 per pulse. The use of a double gas puff configuration with a deuterium center jet allowed a reduction in the shell density and operation closer to machine-matched conditions, recording up to (4.1 ± 0.3) × 10 7 neutrons/pulse when either Kr or D 2 was used in the shell. For a comparable mass and implosion time, using a higher atomic-number gas in the outer shell results in more unstable plasma surface and smaller plasma radius at the location of instability bubbles, which, however, do not seem to consistently correlate with a higher neutron yield. Comparing implosion dynamics with models and neutron yields with literature scaling suggests that the machine current is not well coupled to the plasma during the final stages of compression. Optimizing current and energy coupling to the pinched plasma is critical to improving performance, particularly in low-impedance drivers.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Feasibility and performance of the staged Z-pinch: A one-dimensional study with FLASH and MACH2

Z-pinch platforms constitute a promising pathway to fusion energy research. Here, we present a one-dimensional numerical study of the staged Z-pinch (SZP) concept using the FLASH and MACH2 codes. We discuss the verification of the codes using two analytical benchmarks that include Z-pinch-relevant physics, building confidence on the codes' ability to model such experiments. Then, FLASH is used to simulate two different SZP configurations: a xenon gas-puff liner (SZP1*) and a silver solid liner (SZP2). The SZP2 results are compared against previously published MACH2 results, and a new code-to-code comparison on SZP1* is presented. Using an ideal equation of state and analytical transport coefficients, FLASH yields a fuel convergence ratio (CR) of approximately 39 and a mass-averaged fuel ion temperature slightly below 1 keV for the SZP2 scheme, significantly lower than the full-physics MACH2 prediction. For the new SZP1* configuration, full-physics FLASH simulations furnish large and inherently unstable CRs (>300) but achieve fuel ion temperatures of many kilo-electron volts. While MACH2 also predicts high temperatures, the fuel stagnates at a smaller CR. The integrated code-to-code comparison reveals how magnetic insulation, heat conduction, and radiation transport affect platform performance and the feasibility of the SZP concept.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Resistive wall tearing mode disruptions in DIII-D and ITER tokamaks

Disruptions are a serious problem in tokamaks, in which thermal and magnetic energy confinement is lost. This paper uses data from the DIII-D experiment, theory, and simulations to demonstrate that resistive wall tearing modes (RWTMs) produce the thermal quench (TQ) in a typical locked mode shot. Analysis of the linear RWTM dispersion relation shows the parameter dependence of the growth rate, particularly on the resistive wall time. Linear simulations of the locked mode equilibrium show that it is unstable with a resistive wall and stable with an ideally conducting wall. Nonlinear simulations demonstrate that the RWTM grows to sufficient amplitude to cause a complete thermal quench. The RWTM growth time is proportional to the thermal quench time. Here, the nonlinearly saturated RWTM magnetic perturbation amplitude agrees with experimental measurements. The onset condition is that the q = 2 rational surface is sufficiently close to the resistive wall. Collectively, this identifies the RWTM as the cause of the TQ. In ITER, RWTMs will produce long TQ times compared to present-day experiments. ITER disruptions may be significantly more benign than previously predicted.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Competing mechanisms at vibrated interfaces of density-contrast fluids

Fluid-fluid interfacial instability and subsequent fluid mixing are ubiquitous in nature and engineering. The hydrodynamic instability of fluid interfaces has long centered on the pressure gradient-driven long-wavelength Rayleigh-Taylor instability and the resonance-induced short-wavelength Faraday instability. However, neither instability alone can explain the dynamics when both mechanisms are present. We identify a previously unseen multi-modal instability emerging from their coexistence. When the denser fluid is polydimethylsiloxane, the mixed region at a high density contrast (Atwood number = 0.9) spans a vibration amplitude range approximately twice the gravitational acceleration. Using Floquet stability analysis, we show how vibrations govern transitions between the RT and Faraday instabilities, leading to contention between these instabilities rather than resonant enhancement. Here, the initial transient growth is represented by the exponential modal growth of the most unstable Floquet exponent, along with its accompanying periodic behavior. Direct numerical simulations validate these findings and track interface breakup into the multiscale and nonlinear regimes. Specifically, we show that growing RT modes nonlinearly suppresses Faraday responses even when the initial growth rate of the Faraday instability is 3.63 times that of RT, so a bidirectional competition hinders their sustained coexistence.

Direct numerical simulations↗
DOE OSTI2023

Effect of outer-cylinder rotation on the radially heated Taylor–Couette flow

A Taylor–Couette setup with radial heating is considered where a Boussinesq fluid is sheared in the annular region between two concentric, independently rotating cylinders maintained at different temperatures. Linear stability analysis is performed to determine the Taylor number for the onset of instability. Two radius ratios corresponding to wide and thin gaps with several rotation rate ratios are considered. The rotation of the outer cylinder is found to have a general stabilizing effect on the stability threshold as compared to pure inner-cylinder rotation, with a few exceptions. The radial heating sets up an axial flow which breaks the reflection symmetry of isothermal Taylor–Couette flow in the axial coordinate. This symmetry breaking separates linear stability thresholds, and we find the fastest growing modes with both positive and negative azimuthal numbers for different parameters. Another important finding of the current study is the discovery of unstable modes in the Rayleigh-stable regime. Furthermore, closed disconnected neutral curves (CDNCs) are observed for both wide and thin gaps which can separate from or merge into open neutral–stability curves. Alternatively, CDNCs can also morph into open neutral stability curves as the rotation rate ratio is changed. CDNCs are observed to be sensitive to changes in control parameters, and their appearance/disappearance is shown to induce discontinuous jumps in the critical Taylor number. Finally, for both wide and thin gaps, the fastest-growing modes found in the pure corotation case are shown to have their origins in the instability islands at smaller values of rotation rate ratios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Effects of spin-orbit coupling and very large supercells on the description of acceptors in CdTe

Predicting accurate shallow donor and acceptor levels in semiconductors has been quite challenging using periodic boundary conditions as implemented in current density functional theory codes. The reason is that the wave functions associated with the shallow centers are quite extended and are not fully contained in the supercell (typically with a few hundred atoms) so impurities in the periodically repeated image cells interact with each other. Errors of ~0.1 eV are expected, and these are of the same order of magnitude as the ionization energies themselves. In the case of acceptors in CdTe, this problem is exacerbated by the strong spin-orbit coupling that split the Te-related states at the top of the valence band, making the calculations at least 8 times more expensive. CdTe is an important solar-cell material with record high efficiency of 22%. One of the main limiting factors to increasing the efficiency towards the theoretical limit of ~30% is the often reported very low hole concentration. Pushing to the limit of computational capability by using very large supercells with spin-orbit coupling we report the results of hybrid functional calculations of group-V acceptors in CdTe. We show that extrapolation to the dilute limit leads to an interpretation of the experimental data that is qualitatively different from previous DFT and hybrid functional calculation reports. We find that the group-V impurities indeed behave as shallow acceptors and that the corresponding compensating AX-centers are unstable and do not limit p-type doping. We address the differences between our results and previous theoretical predictions and show that our calculated ionization energies predict hole concentrations that are in excellent agreement with recent temperature-dependent Hall measurements on high-quality single-crystal samples.

14 SOLAR ENERGY↗
DOE OSTI2021

Adaptive Critical Balance and Firehose Instability in an Expanding, Turbulent, Collisionless Plasma

Abstract Using a hybrid-kinetic particle-in-cell simulation, we study the evolution of an expanding, collisionless, magnetized plasma in which strong Alfvénic turbulence is persistently driven. Temperature anisotropy generated adiabatically by the plasma expansion (and consequent decrease in the mean magnetic-field strength) gradually reduces the effective elasticity of the field lines, causing reductions in the linear frequency and residual energy of the Alfvénic fluctuations. In response, these fluctuations modify their interactions and spatial anisotropy to maintain a scale-by-scale “critical balance” between their characteristic linear and nonlinear frequencies. Eventually the plasma becomes unstable to kinetic firehose instabilities, which excite rapidly growing magnetic fluctuations at ion-Larmor scales. The consequent pitch-angle scattering of particles maintains the temperature anisotropy near marginal stability, even as the turbulent plasma continues to expand. The resulting evolution of parallel and perpendicular temperatures does not satisfy double-adiabatic conservation laws, but is described accurately by a simple model that includes anomalous scattering. Our results have implications for understanding the complex interplay between macro- and microscale physics in various hot, dilute, astrophysical plasmas, and offer predictions concerning power spectra, residual energy, ion-Larmor-scale spectral breaks, and non-Maxwellian features in ion distribution functions that may be tested by measurements taken in high-beta regions of the solar wind.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Cooling and instabilities in colliding flows

ABSTRACT Collisional self-interactions occurring in protostellar jets give rise to strong shocks, the structure of which can be affected by radiative cooling within the flow. To study such colliding flows, we use the AstroBEAR AMR code to conduct hydrodynamic simulations in both one and three dimensions with a power-law cooling function. The characteristic length and time-scales for cooling are temperature dependent and thus may vary as shocked gas cools. When the cooling length decreases sufficiently and rapidly, the system becomes unstable to the radiative shock instability, which produces oscillations in the position of the shock front; these oscillations can be seen in both the one- and three-dimensional cases. Our simulations show no evidence of the density clumping characteristic of a thermal instability, even when the cooling function meets the expected criteria. In the three-dimensional case, the nonlinear thin shell instability (NTSI) is found to dominate when the cooling length is sufficiently small. When the flows are subjected to the radiative shock instability, oscillations in the size of the cooling region allow NTSI to occur at larger cooling lengths, though larger cooling lengths delay the onset of NTSI by increasing the oscillation period.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Modification of wind turbine wakes by large-scale, convective atmospheric boundary layer structures

In this study, we consider the impact of large-scale, convective structures in an unstable atmospheric boundary layer on wind turbine wakes. Simulation data from a high-fidelity large-eddy simulation (LES) of the AWAKEN wind farm site matching unstable atmospheric conditions were analyzed, and both turbine performance and wake behavior were affected based on their location relative to the convective structures. Turbines located in updraft regions of the flow experienced lower inflow velocity and generated less power, but their wakes were observed to recover faster and saw greater turbulent kinetic energy mixing higher in the boundary layer. The opposite effect was found for turbines in the downdraft regions of the convective structures. A simplified model of this wake behavior was also developed based on a two-dimensional k–ε Reynolds-Averaged Navier–Stokes formulation. This simplified model included the effects of vertical transport, but could be efficiently solved as a parabolic system, and was found to capture similar wake modifications observed in the high-fidelity LES computations.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Updated systematics of intermediate-energy single-nucleon removal cross sections

The body of experimental measurements of intermediate-energy reactions that remove a single nucleon from a secondary beam of neutron- or proton-rich nuclei continues to grow. These data have been analyzed consistently using an approximate, eikonal-model treatment of the reaction dynamics combined with appropriate shell-model descriptions of the projectile initial state, the bound final states spectrum of the reaction residue, and single-particle removal strengths computed from their wave-function overlaps. Here, the systematics of the ratio R s of the measured inclusive cross section to all bound final states and the calculated cross section to bound shell-model states—in different regions of the nuclear chart and involving both very weakly bound and strongly bound valence nucleons—is important in relating the empirically deduced orbital occupancies to those from the best available shell-model predictions. Importantly, several new higher-energy measurements, for which the sudden-approximation aspect of the dynamical description is placed on an even stronger footing, now supplement the previously analyzed measurements. These additional data sets are discussed. Their R s values are shown to conform to and reinforce the earlier-observed systematics, with no indication that the approximately linear reduction in R s with increasing nucleon separation energy is a consequence of a breakdown of the sudden approximation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2025

Microscopic optical potentials from a Green's function approach

Optical potentials are a standard tool in the study of nuclear reactions, as they describe the interaction between a target nucleus and a projectile. The use of phenomenological optical potentials built using experimental data on stable isotopes is widespread. Although successful in their dedicated domain, it is unclear whether these phenomenological potentials can provide reliable predictions for unstable isotopes. To address this problem, optical potentials based on microscopic nuclear structure input calculations prove to be crucial and are an important current line of research. In this work we present an explicit implementation of the Feshbach formalism for the systematic derivation of optical potentials using input from nuclear structure models. Numerical tools for the derivation of Green's functions associated with nonlocal potentials are presented. In conclusion, the new optical potential, based on the valence shell model, is applied to the calculations of 𝑛 + 24 Mg elastic scattering and yields a close agreement with the experimental data.

Direct reactions↗
DOE OSTI2022

The role of plasma instabilities in relativistic radiation-mediated shocks: stability analysis and particle-in-cell simulations

ABSTRACT Relativistic radiation-mediated shocks are likely formed in prodigious cosmic explosions. The structure and emission of such shocks are regulated by copious production of electron–positron pairs inside the shock-transition layer. It has been pointed out recently that substantial abundance of positrons inside the shock leads to a velocity separation of the different plasma constituents, which is expected to induce a rapid growth of plasma instabilities. In this paper, we study the hierarchy of plasma microinstabilities growing in an electron-ion plasma loaded with pairs and subject to a radiation force. Linear stability analysis indicates that such a system is unstable to the growth of various plasma modes which ultimately become dominated by a current filamentation instability driven by the relative drift between the ions and the pairs. These results are validated by particle-in-cell simulations that further probe the non-linear regime of the instabilities, and the pair-ion coupling in the microturbulent electromagnetic field. Based on this analysis, we derive a reduced-transport equation for the particles via pitch-angle scattering in the microturbulence and demonstrate that it can couple the different species and lead to non-adiabatic compression via a Joule-like heating. The heating of the pairs and, conceivably, the formation of non-thermal distributions, arising from the microturbulence, can affect the observed shock-breakout signal in ways unaccounted for by current single-fluid models.

79 ASTRONOMY AND ASTROPHYSICS↗