Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “wave prediction”

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 991 records · Page 55

Subseasonal Forecasting and MJO Teleconnections in Machine Learning Weather Prediction Models

Abstract In recent years, machine‐learning (ML) models trained on reanalysis data have rivaled physics‐based forecast models in terms of performance skill for global weather forecasting. With increased rollout stability, the question of how these models perform for subseasonal to seasonal (S2S, week 3–8) forecasting has emerged. In this study we run a large set of subseasonal hindcasts over 2004–2023 to evaluate two ML weather forecast models at the S2S time scale, SFNO‐HENS (Nvidia, fully ML) and NeuralGCM (Google Research, hybrid). Corresponding hindcasts from the European Centre for Medium‐Range Weather Forecasts (ECMWF) are used as a baseline for comparison to a physics‐based model. Because our focus is on predicting moisture transport over the Western United States between October and March, we evaluate the models' prediction skill for the Madden‐Julian Oscillation (MJO) and its associated teleconnections in the North Pacific. We find that both ML models are competitive with the ECWMF model, with comparable skill in predicting the North Pacific large‐scale circulation and the MJO at week 3 and beyond. Even though overall the mid‐latitude subseasonal prediction skill remains low, the ML models exhibit interesting behavior such as a realistic propagation of the MJO across the Maritime Continent and realistic teleconnections. A SFNO‐HENS sensitivity experiment with altered initial conditions in the tropics demonstrates the stability of the model, and it illustrates the capability of ML models to represent important physical processes of the atmosphere at the S2S time scale. Plain Language Summary Predicting weather patterns and precipitation a few weeks in advance (subseasonal time scale) is of great interest for stakeholders such as water managers in the Southwest United States (US), where arid conditions prevail. Subseasonal forecasts from traditional weather forecast models exhibit low skill in the region, limiting their applicability. Here we examine whether the recent breakthrough in weather forecasting made with machine learning/artificial intelligence models can translate to improved subseasonal forecasts. Recently‐developed machine learning models exhibit comparable skill to a state‐of‐the‐art physics‐based model for predicting weather patterns in the North Pacific/North America region, and associated moisture transport. The same applies to their skill in predicting the tropical pattern, the Madden‐Julian Oscillation, and its important remote perturbations over the midlatitude East Pacific and Southwest US. Additionally, a perturbation experiment carried out with one of the machine learning models illustrates their ability to not only predict the evolution of atmospheric fields, but also to learn and represent physical processes such as tropics‐extratropics Rossby wave propagation. Key Points Two machine learning weather forecast models exhibit state‐of‐the‐art prediction skill at the subseasonal time scale in the Pacific sector The models equal ECWMF in terms of Madden‐Julian oscillation (MJO) prediction skill, and they accurately predict the MJO propagation and associated teleconnections The two machine‐learning models represent key physical processes for subseasonal prediction, despite being trained for weather forecasting

Peings, Yannick↗

The evolution of methods for noise prediction of high speed rotors and propellers in the time domain

Linear wave equation models which have been used over the years at NASA Langley for describing noise emissions from high speed rotating blades are summarized. The noise sources are assumed to lie on a moving surface, and analysis of the situation has been based on the Ffowcs Williams-Hawkings (FW-H) equation. Although the equation accounts for two surface and one volume source, the NASA analyses have considered only the surface terms. Several variations on the FW-H model are delineated for various types of applications, noting the computational benefits of removing the frequency dependence of the calculations. Formulations are also provided for compact and noncompact sources, and features of Long's subsonic integral equation and Farassat's high speed integral equation are discussed. The selection of subsonic or high speed models is dependent on the Mach number of the blade surface where the source is located.

Farassat, F.↗

Simulations of ion cyclotron anisotropy instabilities in the terrestrial magnetosheath

Enhanced transverse magnetic fluctuations observed below the proton cyclotron frequency in the terrestrial magnetosheath have been identified as due to the proton cyclotron and helium cyclotron instabilities driven by the T-perpendicular greater than T-parallel condition of the sheath ions. One-dimensional hybrid computer simulations are used here to examine the nonlinear properties of these two growing modes at relatively weak fluctuation energies and for wave vectors parallel to the background magnetic field. Second-order theory predicts fluctuating magnetic field energies at saturation of the proton cyclotron anisotropy instability in semiquantitative agreement with the simulation results. Introduction of the helium component enhances the wave-particle exchange rate for proton anisotropy reduction by that instability, thereby reducing the saturation energy of that mode. The simulations demonstrate that wave-particle interactions by the proton cyclotron and helium cyclotron instabilities lead to the anticorrelation observed by Anderson and Fuselier (1993).

Gary, S. P.↗

A NASTRAN DMAP alter for linear buckling analysis under dynamic loading

A modification to the NASTRAN solution sequence for transient analysis with direct time integration (COSMIC NASTRAN rigid format 9) was developed and incorporated into a DMAP alter. This DMAP alter calculates the buckling stability of a dynamically loaded structure, and is used to predict the onset of structural buckling under stress-wave loading conditions. The modified solution sequence incorporates the linear buckling analysis capability (rigid format 5) of NASTRAN into the existing Transient solution rigid format in such a way as to provide a time dependent eigensolution which is used to assess the buckling stability of the structure as it responds to the impulsive load. As a demonstration of the validity of this modified solution procedure, the dynamic buckling of a prismatic bar subjected to an impulsive longitudinal compression is analyzed and compared to the known theoretical solution. In addition, a dynamic buckling analysis is performed for the analytically less tractable problem of the localized dynamic buckling of an initially flawed composite laminate under transverse impact loading. The addition of this DMAP alter to the transient solution sequence in NASTRAN facilitates the computational prediction of both the time at which the onset of dynamic buckling occurs in an impulsively loaded structure, and the dynamic buckling mode shapes of that structure.

Aiello, Robert A.↗

A study of the limitations of linear theory methods as applied to sonic boom calculations

Current sonic boom minimization theories have been reviewed to emphasize the capabilities and flexibilities of the methods. Flexibility is important because it is necessary for the designer to meet optimized area constraints while reducing the impact on vehicle aerodynamic performance. Preliminary comparisons of sonic booms predicted for two Mach 3 concepts illustrate the benefits of shaping. Finally, for very simple bodies of revolution, sonic boom predictions were made using two methods - a modified linear theory method and a nonlinear method - for signature shapes which were both farfield N-waves and midfield waves. Preliminary analysis on these simple bodies verified that current modified linear theory prediction methods become inadequate for predicting midfield signatures for Mach numbers above 3. The importance of impulse is sonic boom disturbance and the importance of three-dimensional effects which could not be simulated with the bodies of revolution will determine the validity of current modified linear theory methods in predicting midfield signatures at lower Mach numbers.

Darden, Christine M.↗

Pulsar average wave forms and hollow-cone beam models

Pulsar wave forms have been analyzed from observations conducted over a wide radio-frequency range to assess the wave-form morphologies and to measure wave-form widths. The results of the analysis compare favorably with the predictions of a model with a hollow-cone beam of fixed dimensions and with random orientation of both the observer and the cone axis with respect to the pulsar spin axis. A class of three-component wave forms is included in the model by adding a central pencil beam to the hollow-cone hypothesis. The consequences of a number of discrepancies between observations and quantitative predictions of the model are discussed.

Backer, D. C.↗

Soliton formation and topology manipulation of coupled spins via ultrafast re-magnetization

The major goal of the project was to explore the properties of magnetic order driven far out of equilibrium by optical excitations. These include fundamental questions related to the interplay of magnetic, structural and electronic degrees of freedom in materials that are optically excited. These problems break down in short term quenching of the magnetization as a results of the energy of the optical pulse studying how the energy of the optical pulse and the angular momentum of the magnetization flow between the different degrees of freedom. This is followed by the longer time re-emergence of the magnetism as the system cools. Within this goal we explored the formation of solitons, both with and without non-trivial topology, via rapid re-magnetization processes after optical-driven ultrafast demagnetization. Guided by theory, we predict that turbulence and modulational instabilities will drive the formation of solitons, including dispersive shock waves, magnon droplets, and skyrmions. In the case of topological defect formation, i.e. skyrmion generation, related processes have long been predicted based on general principles of symmetry-breaking phase transitions; the density of topological defects in a long-range-ordered phase can be controlled by varying the quench rate through the second-order phase transition, i.e. the Kibble-Zurek (KZ) mechanism. We, more broadly, had the goal to study a range of other emergent magnetic behaviors after optical excitation. While great attention has been given to ultrafast demagnetization, less is known of the subsequent spin dynamics and coupling to the lattice. What is increasingly appreciated is that ultrafast demagnetization leads to spin currents that can carry angular momentum from the rapidly demagnetized sample. These spin currents appear in many ways but fundamentally controls demagnetization, drives interactions between regions of the material, and control the spin resulting structure. They can be probed many ways including THz emission which, in turn, gives insight into the demagnetization processes. Simultaneously if there is magneto-elastic coupling ultrafast demagnetization can drive structural excitations that are not expected from thermal energy added to the lattice. We explored a range of complex phenomena that arise from photo-excitation of magnetic systems with the coupling of electronic, magnetic and structural degrees of freedom. We have made a number of fundamental discoveries that are reflected in our publications list with additional work still being prepared for publication. Highlights of this work include (i) ultra-efficient, nonlinear THz surface acoustics (ii) spin-current-mediated rapid magnon localization and coalescence, (iii) spin-wave soliton formation in ferromagnetic FePt nanoparticles, (iv) dynamic phonon coupling in elemental antiferromagnetic Cr, (v) THz emission from Co/Pt bilayers and FeRh/Pt bilayers, (vi) theoretical investigation of spin hydrodynamics, solitons and shock waves, and (vii) ultrafast perturbation of magnetic domains by optical pumping.

36 MATERIALS SCIENCE↗

Dynamics in Layer Models of Solid Flame Propagation

Self-propagating high-temperature synthesis (SHS) is a process in which combustion waves, e.g., "solid flames", which are considered here, are employed to synthesize desired materials. Like many other systems, SHS is a pattern forming system. The problem of describing experimentally observed patterns and of predicting new, as yet unobserved, patterns continues to attract the attention of scientists and mathematicians due to the fundamental significance of the phenomena in combustion in particular, and in nonlinear science in general. Here, we analyze the dynamics of solid flame propagation in a 2D region by considering the region to be composed of parallel, identical layers aligned along the direction of propagation and having thermal contact. Each layer is then described by wave propagation in 1D, with the transverse Laplacian replaced by a term describing heat exchange between neighboring layers. This configuration is the simplest model of a 2D system because it accounts, in a simple way, for the principal feature of the problem, i.e., heat exchange between neighbors in the transverse direction. For simplicity, we describe the situation for two layers. Because the layers are identical, uniformly propagating waves in each layer must be identical, independent of the heat exchange rate alpha. When the Zeldovich number Z exceeds a critical value Z(sub c), which depends on alpha, uniformly propagating waves become unstable. The stability diagram for the two coupled layers reproduces that for the full 2D problem after appropriate identification of parameters in the two problems. Depending on parameter values, we determine three different steady-state dynamical behaviors (though additional behaviors are also expected to occur). The three behaviors are: (i) waves in each layer which pulsate in phase as they propagate, so that together they form a single pulsating propagating wave; (ii) the waves in each layer are no longer identical, and antiphase pulsations occur, with the two waves alternately advancing and receding as they propagate. This wave is the analog of the spinning wave on the surface of a circular cylinder; (iii) finally, there is a region of bistability between the in phase and antiphase waves. with each having its own domain of attraction, so that which of the two behaviors occur depends on the condition of initiation of the wave. The results of our computations indicate a qualitative similarity in the behavior of combustion waves in the layer model and in the full 2D model. Specifically, due to the similarity between the small alpha wave behavior in the layer model and the large diameter behavior in the model of waves on the surface of a cylinder, we are able to predict the behavior of the mean velocity for the waves on the cylinder, where computations of the full problem can be rather difficult. We also compute solutions for three or more layers. The results of our computations prompt us to predict that, while planar uniformly propagating waves are unstable, the wave will be quasiplanar, i.e., the resulting spinning waves have very low amplitude hot spots, and travel with the velocity close to that of the uniformly propagating wave. Such waves may be difficult to distinguish from uniformly propagating waves in experiments. We also find that for both the layer model and full 2D problem, steady-state time-dependent waves, e.g., pulsating and spinning wave, have a conserved quantity H which characterizes the excess energy in the wave, just as in the case of uniformly propagating waves. The quantity H, which is generated by dissipation, does not vary in time and is proportional to the diffusivity and caloricity of the system, and inversely proportional to the mean wave velocity.

Aldushin, A. P.↗

Acceleration of energetic charged particles by interplanetary and interstellar shock waves

Pioneer 11 observations of energetic proton acceleration in interplanetary shocks are in good agreement with the predictions of the Vector-V x Vector-B electric field model of shock acceleration. The Vector-V x Vector-B model predicts that relativistic particles are accelerated by a factor of 2-15 in kinetic energy during a single encounter with a quasi-perpendicular (not less than 70 deg) supernova shock wave, and that the post-encounter pitch angle distributions are very anisotropic.

Pesses, M. E.↗

Simulations of SAR wave spectra using high spectral resolution estimates from the SCR and ROWS instruments

A numerical model for predicting the synthetic aperture radar (SAR) image of a moving ocean surface is described, and results are presented for two SIR-B data sets collected off the coast of Chile. Wave height spectra measured by the NASA radar ocean wave spectrometer (ROWS) and surface contour radar (SCR) were used as inputs to this model, and results are compared with actual SIR-B image spectra from orbits 91 and 106.

Lyzenga, D.↗

Pressure Fluctuation Spectral Peaks Due to ‘Guided Waves’ in Regions Upstream of the Jet Exit

Pressure fluctuation spectral peaks due to ‘guided’ (or ‘trapped’) waves in compressible jets in regions upstream of the nozzle exit are investigated experimentally. These spectral peaks are known to disappear in the far acoustic field on the sideline and in the downstream direction, where noise measurements are done typically. The data in this paper show that they are prevalent in the far upstream regions of the ‘forward arc’, at shallow angles relative to the jet axis. This is noted with a round as well as a rectangular jet. The result is significant since such spectral behavior is not predicted by existing noise prediction codes. The experimental results not only expand our understanding of the trapped wave phenomenon but may also aid in refining noise prediction codes for basic as well as application-oriented propulsion systems.

Jets↗

Twenty-two-year cycle of the upper limiting rigidity of Daly waves

The method of calculating energy losses along regular particle trajectories is applied to obtain the predicted cosmic ray anisotropies from 200 to 500 GV. The tilt angle of the interplanetary neutral sheet varies to simulate a 22-year cycle magnetic cycle. The calculated values of solar diurnal and semidiurnal, and sidereal diurnal intensity waves are compared with observations.

Erdos, G.↗

Session on validation of coupled models

The session on validation of coupled models is reviewed. The current use of a mesoscale model with a grid size of 20-km during STORM-FEST in 1992 has proven to be extremely valuable. The availability of forecast products at a much higher temporal and spatial resolution was very helpful for mesoscale forecasting, mission planning, and the guidance of research aircraft. Recent numerical simulation of ocean cyclones and mesoscsle convective systems using nonhydrostatic cloud/mesoscale models with a grid size as small as 2-km have demonstrated the potential of these models for predicting mesoscale convective systems, squall lines, hurricane rainbands, mesoscale gravity waves, and mesoscale frontal structures embedded within an extratropical cyclone. Although mesoscale/cloud scale models have demonstrated strong potential for use in operational forecasting, very limited quantitative evaluation (and verification) of these models were performed. As a result, the accuracy, the systematic biases, and the useful forecasts limits were not properly defined for these models. Also, no serious attempts were made to use these models for operational prediction of mesoscale convective systems.

Kuo, Bill↗

Progress on Computational Modeling of Water-Based NSTF (FY2021)

This report documents the FY21 progress and achievements made in the computational analyses of the water-based NSTF. Both system-level and high-fidelity Computational Fluid Dynamics (CFD) analyses were performed to gain a complete understanding of the complex flow and heat transfer phenomena in natural convection systems. The progress on the waterbased NSTF experimental testing is summarized in a companion report (ANL-ART-230). As a continuation of progress from previous years, in FY21 the RELAP5 model of the NSTF was first updated to include a representative heated cavity that incorporated convection and conduction means of heat transfer. Two different cavity models were developed, one that models the radiation heat exchange among all surfaces in one enclosure (Cavity Model 1), and another that employs multiple radiation enclosures with each considering the surfaces at one axial level (Cavity Model 2). This updated RELAP5 deck was benchmarked with single-phase test data to tune form loss coefficients of the elbows and tees, along with heat transfer coefficients within the cavity and off the external insulation panels. This tuned model was then found to accurately predict fluid temperatures within 2.5% and both system and riser flow rates within 6% during single-phase, steady-state simulations. For transient simulations, there is a slight over prediction of the flow rate, fluid temperatures, and heater temperatures during the heat-up period, which indicate that the transient heat loss may be underestimated. Following, the capabilities of RELAP5-3D were examined, which feature a conduction enclosure model that is not present in the current version of RELAP5-MOD3.3 used for all previous NSTF simulations. As an initial comparison of the two codes, single-phase simulations were run that used similar input decks. The results indicate that the single-phase flow predicted by the two codes was very similar, however RELAP5-3D appears to predict a higher system flow during the two-phase transient when compared to the flow rates predicted by RELAP5-MOD3.3. Furthermore, RELAP5-3D run did not capture the density wave oscillations or geysering phenomena, which are large oscillation flow instabilities observed in the experimental facility. To further the prediction capabilities for capturing two-phase instabilities, Options 55, 57, 58, 61, which are implemented by default in RELAP5-MOD3.3 but not in RELAP5-3D, were investigated. When enabling Option 61, RELAP5-3D agrees well with MOD3.3 in predicting the total system flow, while other options have negligible effect on the predicted total system flow by RELAP5-3D. However, none of these options enables RELAP5-3D to produce the same flow instabilities. For future steps, efforts are currently underway to investigate RELAP5-3D’s different behavior in predicting low-pressure two-phase instabilities. STAR-CCM+ was used to perform CFD analysis of the water-based test facility and provide simulation results for comparisons to experimental data. Efforts first focused on describing the characteristics of flow within the heated cavity. A study was performed on whether possible leakage flow between the two sides of the cavity could substantially impact and explain the results. An additional second assumption was examined which defined adiabatic side walls, so physical insulation structures were added to the models to allow conduction across this outer insulation. The results suggest that potential porosity of the insulation could be a major factor in the temperature distribution within the cavity and should be considered in any studies where it is desirable to know the cavity wall temperatures to higher accuracy.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A Gapped Phase in Semimetallic T d -WTe 2 Induced by Lithium Intercalation

The Weyl semimetal WTe 2 has shown several correlated electronic behaviors, such as the quantum spin Hall effect, superconductivity, ferroelectricity, and a possible exciton insulator state, all of which can be tuned by various physical and chemical approaches. A new electronic phase in WTe 2 induced by lithium intercalation is discovered. The new phase exhibits an increasing resistivity with decreasing temperature and its carrier density is almost two orders of magnitude lower than the carrier density of the semimetallic T d phase, probed by in situ Hall measurements as a function of lithium intercalation. The theoretical calculations predict the new lithiated phase to be a potential charge density wave (CDW) phase with a bandgap of ≈0.14 eV, in good agreement with the in situ transport data. The new phase is structurally distinct from the initial T d phase, characterized by polarization-angle-dependent Raman spectroscopy, and large lattice distortions close to 6% are predicted in the new phase. This finding of a new gapped phase in a 2D semimetal demonstrates electrochemical intercalation as a powerful tuning knob for modulating electron density and phase stability in 2D materials.

36 MATERIALS SCIENCE↗

Why Are There so Few Reports of High-Energy Electron Drift Resonances? Role of Radial Phase Space Density Gradients

Models of monochromatic Pc5 (2–7 mHz) ultralow frequency (ULF) wave interactions with high energy (greater than ~1 MeV) electrons predict drift resonant interactions that can cause rapid radial transport and acceleration. There are few reports of electron drift resonance at energies greater than ~1 MeV, in contrast to lower energies; moreover, all previous reports occur in the aftermath of interplanetary shocks. These two facts are difficult to reconcile with theory and numerical simulations predicting that greater than ~1 MeV drift resonances should occur more often and in a wider variety of driving conditions. In this study, we show that a combination of observational sampling biases and nominal radial phase space density gradients is one explanation for this discrepancy between theory and observations. In particular, we examine electron dynamics in two case studies with very similar satellite coverage, solar wind conditions, and Pc5 wave properties, yet with different radial phase space density profiles. Using global wave and particle observations, we show that the events have vastly different particle responses despite having similar wave properties. Placing these results in context with past studies, we further show here that nominal radial PSD gradients near geostationary orbit can mask the expected drift resonance particle response and explain (1) the small number of past greater than ~1 MeV drift resonance reports and (2) the restriction of these reports to interplanetary shock events. We argue that future observational studies characterizing radial transport via drift resonance should examine global particle dynamics, including observations of the radial phase space density profile.

79 ASTRONOMY AND ASTROPHYSICS↗

Do black holes remember what they are made of?

We study the ringdown signal of black holes formed in prompt-collapse binary neutron star mergers. We analyze data from 47 numerical relativity simulations. We show that the ($l$ = 2, $m$ = 2) and ($l$ = 2, $m$ = 1) multipoles of the gravitational wave signal are well fitted by decaying damped exponentials, as predicted by black-hole perturbation theory. We show that the ratio of the amplitude in the two modes depends on the progenitor binary mass ratio q and reduced tidal parameter $\overline{Λ}$. Unfortunately, the numerical uncertainty in our data is too large to fully quantify this dependency. If confirmed, these results will enable novel tests of general relativity in the presence of matter with next-generation gravitational-wave observatories.

79 ASTRONOMY AND ASTROPHYSICS↗

Confirmation of the predicted Hartree-Fock limit in BeH sub 2.

We have recently estimated the total energy of the BeH 2 molecule as -15.918 ± 0.015 au. One of the steps in the analysis was to estimate the difference between the SCF energies computed with a minimal basis set of Slater orbits and with the best possible single-determinant wave function, the so-called Hartree-Fock limit and was predicted to be -15.773 ± 0.007 au for BeH 2 .

Molecular energy↗