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 1,225 records · Page 68

Insights into Slip-Rate Time Functions, Rupture Parameter Correlations, and Ground Motions from Validated Multicycle Earthquake Ruptures

Earthquake strong-motion predictions using kinematic source modeling require knowledge of the slip-rate functions (SRFs) along the rupture and their distinct characteristics in asperities, background (off-asperity) areas and near the surface. Here, in this study we analyzed SRFs from well-validated, self-consistent, and fully dynamic rupture models from earthquake cycles obeying a rate-and-state friction law, from our companion study (Galvez et al., 2021). The shapes of SRFs in asperities are well described by the regularized Yoffe function (RYF), which has only two parameters: rise time T r and smoothing time T s , which control the generation of long- and short-period ground motions, respectively. In background areas, we demonstrate that, in addition to the primary rupture, multiple secondary ruptures may also nucleate from rupture heterogeneities related to asperities, resulting in SRFs with multiple peaks. Because it is impossible to fit a multiple-peak SRF by the single-peak RYF, we describe SRFs in background areas in an effective way by fitting their amplitude spectra with the RYF spectra. Such spectrally effective RYFs capture salient aspects of seismic-wave generation and can be used in rupture generators for strong motion prediction. We found that small T s values correlate with small characteristic weakening distances, large peak slip rates (PSRs), and large rupture velocities. T r values are larger in background areas and smaller in asperities. Within the shallow aseismic zone, T s values approximately quadruple whereas T r values approximately double. Because of this dominant T s increase, PSR values decrease in the near-surface zone. These features indicate that the generation of strong motions by the near-surface portions of the rupture is negligible in the studied scenarios.

Geosciences↗

Prediction Skill of the 2012 U.S. Great Plains Flash Drought in Subseasonal Experiment (SubX) Models

Flash droughts refer to droughts that develop much more rapidly than normal (i.e., on the order of weeks to a few months). Such droughts can have devastating impacts on agriculture, water resources, and ecosystems. The ability to predict flash droughts in advance would greatly enhance our preparation for them and potentially mitigate their impacts. We investigated the prediction skill of U.S. flash droughts at subseasonal lead times in global forecast systems participating in the Subseasonal Experiment (SubX) project. An additional comprehensive set of hindcasts with NASA?s GEOSv2.1, a model with relatively high prediction skill, was performed to investigate the separate contributions of atmospheric and land initial conditions to flash drought prediction skill. Here we focus on results for the 2012 Great Plains flash drought, noting that the findings based on this event are generally applicable to other U.S. flash droughts. The prediction skill of the SubX models is quite variable. While the skill is limited to less than 2 weeks in most models, it is considerably higher (3-4 weeks or more) for certain models and initialization dates. The enhanced prediction skill is found to originate from two robust sources: 1) accurate soil moisture initialization, and 2) the satisfactory representation of quasi-stationary cross-North Pacific Rossby wave trains that lead to the rapid intensification of flash droughts. Our results corroborate earlier findings that accurate soil moisture initialization is important for skillful subseasonal forecasts and highlight the need for additional research on the sources and predictability of drought-inducing quasi-stationary Rossby waves.

DeAngelis, Anthony M.↗

Optimal Growth in Hypersonic Boundary Layers

The linear form of the parabolized linear stability equations is used in a variational approach to extend the previous body of results for the optimal, nonmodal disturbance growth in boundary-layer flows. This paper investigates the optimal growth characteristics in the hypersonic Mach number regime without any high-enthalpy effects. The influence of wall cooling is studied, with particular emphasis on the role of the initial disturbance location and the value of the spanwise wave number that leads to the maximum energy growth up to a specified location. Unlike previous predictions that used a basic state obtained from a self-similar solution to the boundary-layer equations, mean flow solutions based on the full Navier-Stokes equations are used in select cases to help account for the viscous- inviscid interaction near the leading edge of the plate and for the weak shock wave emanating from that region. Using the full Navier-Stokes mean flow is shown to result in further reduction with Mach number in the magnitude of optimal growth relative to the predictions based on the self-similar approximation to the base flow.

Navier Stokes equations↗

Prediction of jet mean flow structure in support of HSCT noise suppression concepts

The paper describes the application of techniques based on computational fluid dynamics to the simulation of jet flowfields. A solution code for the Reynolds-averaged Navier-Stokes equations is supplemented by conventional two-equation turbulence models based on the Boussinesq approximation. The axisymmetric SCIPVIS code is enhanced with the PARCH and CRAFT codes to examine plug-jet flowfields and imperfectly expanded axisymmetric free round jets. The sensitivity of shock/boundary layer interactions is observed in simulations of the plug case, and the adaptive gridding in the disk region and turbulence levels generated at the triple point are identified as areas in the Mach case that require improvement. Jet-wave structure in the region beyond the first several shock cells can be predicted, and turbulence modeling can be undertaken with respect to improving compressibility, length scale, vorticity, and energy budget. The mean flow structure of imperfectly expanded jets can be studied to develop related noise suppression concepts for the High-Speed Civilian Transport (HSCT).

Sinha, N.↗

Osmotic forces and gap junctions in spreading depression: a computational model

In a computational model of spreading depression (SD), ionic movement through a neuronal syncytium of cells connected by gap junctions is described electrodiffusively. Simulations predict that SD will not occur unless cells are allowed to expand in response to osmotic pressure gradients and K+ is allowed to move through gap junctions. SD waves of [K+]out approximately 25 to approximately 60 mM moving at approximately 2 to approximately 18 mm/min are predicted over the range of parametric values reported in gray matter, with extracellular space decreasing up to approximately 50%. Predicted waveform shape is qualitatively similar to laboratory reports. The delayed-rectifier, NMDA, BK, and Na+ currents are predicted to facilitate SD, while SK and A-type K+ currents and glial activity impede SD. These predictions are consonant with recent findings that gap junction poisons block SD and support the theories that cytosolic diffusion via gap junctions and osmotic forces are important mechanisms underlying SD.

Spreading Cortical Depression/physiology↗

Concerning the Perturbations Generated by Flux Transfer Events

The interaction of the solar wind with the Earth's magnetosphere is often highly unsteady. Bursts of magnetic reconnection at multiple locations on the dayside equatorial magnetopause generate flux transfer events or FTEs: twisted ropes of interconnected magnetosheath and magnetospheric magnetic field lines. Once formed, the events move antisunward, displacing and perturbing the ambient media. This talk explores the perturbations predicted by both global numerical simulations and analytical models. Results from global hybrid code simulation confirm the predictions of analytical models indicating that the events generate standing forward slow mode waves as their speeds relative to the magnetosheath flow approach the Alfven velocity. Geometric considerations lead to the conclusion that events generated by component reconnection on the dayside magnetopause move poleward and exhibit strong signatures during intervals of southward IMF orientation, but move towards the flanks and exhibit weak signatures during intervals of northward IMF orientation. Changing event orientations and magnetosheath/magnetospheric magnetic field orientation can enhance the amplitudes of events reaching the flanks. Although the orientations of events on the flanks inferred from multispacecraft timing techniques are consistent with the predictions of the component reconnection model, occurrence patterns versus latitude and IMF orientation require an explanation in terms of both the component and antiparallel reconnection models.

Sibeck, David G.↗

Tuning Neural Network Models for Improved Prediction of Boundary Layer Transition

Boundary layer transition can strongly impact flight vehicle performance as it influences surface skin friction and aerodynamic heating, making accurate transition prediction a key to designing next generation aircraft. Artificial neural networks (ANNs) have shown promise toward predicting laminar-turbulent transition based on linear stability correlations. The computational efficiency of ANNs and the substantially reduced user involvement in relation to direct computations based on the linear stability theory (LST) makes them an attractive methodology for integrating the LST based correlations in computational fluid dynamics codes. Tollmien-Schlichting (TS) waves correspond to the dominant transition mechanism in 2D or weakly 3D subsonic boundary layers, such as those encountered in general aviation applications. Improvements to neural network model accuracy in predicting the amplification rates of TS instability waves have been investigated by leveraging recent machine learning developments in conjunction with surrogate optimization techniques and via suitable augmentation of the data used to train the networks. The optimized models trained on the modified dataset reduced the average transition location errors on different airfoils at several flow conditions by 51% of the original manually-tuned network’s errors on the same flow cases. The actual transition locations were derived from the Langley Stability and Transition Analysis Code (LASTRAC).

Machine Learning↗

Nonlinear theory of type I irregularities in the equatorial electrojet

By allowing for the effect of wave electric fields on electron orbits, a nonlinear dispersion relation for type I irregularities is obtained. This relation predicts (1) isotropy of the Doppler shift with elevation, (2) a limiting phase velocity equal to the ion acoustic speed, and (3) a saturation amplitude which is maximum for horizontally propagating waves and decreases with elevation. With the theory presented here, the in situ electric field measurements of Pfaff et al. (1980) in the electrojet environment could provide a quantitative check for the theory of orbit diffusion by stochastic electric fields.

Sudan, R. N.↗

Acceleration of auroral electrons in parallel electric fields

Rocket observations of auroral electrons are compared with the predictions of a number of theoretical acceleration mechanisms that involve an electric field parallel to the earth's magnetic field. The theoretical models are discussed in terms of required plasma sources, the location of the acceleration region, and properties of necessary wave-particle scattering mechanisms. We have been unable to find any steady state scatter-free electric field configuration that predicts electron flux distributions in agreement with the observations. The addition of a fluctuating electric field or wave-particle scattering several thousand kilometers above the rocket can modify the theoretical flux distributions so that they agree with measurements. The presence of very narrow energy peaks in the flux contours implies a characteristic temperature of several tens of electron volts or less for the source of field-aligned auroral electrons and a temperature of several hundred electron volts or less for the relatively isotropic 'monoenergetic' auroral electrons. The temperature of the field-aligned electrons is more representative of the magnetosheath or possibly the ionosphere as a source region than of the plasma sheet.

Kaufmann, R. L.↗

Robust ab initio predictions for dimensionless ratios of 𝐸⁢2 and radius observables. II. Estimation of 𝐸⁢2 transition strengths by calibration to the charge radius

Converged results for 𝐸⁢2 observables are notoriously challenging to obtain in ab initio no-core configuration interaction approaches. Matrix elements of the 𝐸⁢2 operator are sensitive to the large-distance tails of the nuclear wave function, which converge slowly in an oscillator basis expansion. Similar convergence challenges beset ab initio prediction of the nuclear charge radius. However, we exploit systematic correlations between the calculated 𝐸⁢2 and radius observables to yield meaningful predictions for relations among these observables. In particular, we examine ab initio predictions for dimensionless ratios of the form 𝐵⁡(𝐸⁢2)/(𝑒 2 ⁢𝑟 4 ) for nuclei throughout the 𝑝 shell. Finally, meaningful predictions for 𝐸⁢2 transition strengths may then be made by calibrating to the ground-state charge radius if experimentally known.

ab initio calculations↗

Modeling of shock-induced force on an isolated particle in water and air

The prediction of force on an isolated particle, while a shock is passing over it, is an important problem in many natural and industrial applications. Although the flow monotonically changes from the pre-shock to the post-shock state, the particle's force has been observed to behave nonmonotonically with a sharp peak when the shock is located halfway across the particle. This nonmonotonic behavior is due to the unsteady nature of the compression and rarefaction waves that radiate as the shock diffracts around the particle and, therefore, cannot be predicted by a quasi-steady model. An accurate force model must account for the unsteady nature of the flow and the sharp discontinues in the flow properties across the shock. In this work, we test four different inviscid models and observe that the compressible Maxey–Riley–Gatignol (C-MRG) model is the most accurate based on comparison with results from particle-resolved inviscid simulations at two different Mach numbers for both water and air as the medium. The C-MRG model is first demonstrated to predict the force on a stationary particle accurately and then extended to capture the force on a moving particle. Numerical complexities regarding the implementation of the C-MRG model are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Planetary Transmission Diagnostics

This report presents a methodology for detecting and diagnosing gear faults in the planetary stage of a helicopter transmission. This diagnostic technique is based on the constrained adaptive lifting algorithm. The lifting scheme, developed by Wim Sweldens of Bell Labs, is a time domain, prediction-error realization of the wavelet transform that allows for greater flexibility in the construction of wavelet bases. Classic lifting analyzes a given signal using wavelets derived from a single fundamental basis function. A number of researchers have proposed techniques for adding adaptivity to the lifting scheme, allowing the transform to choose from a set of fundamental bases the basis that best fits the signal. This characteristic is desirable for gear diagnostics as it allows the technique to tailor itself to a specific transmission by selecting a set of wavelets that best represent vibration signals obtained while the gearbox is operating under healthy-state conditions. However, constraints on certain basis characteristics are necessary to enhance the detection of local wave-form changes caused by certain types of gear damage. The proposed methodology analyzes individual tooth-mesh waveforms from a healthy-state gearbox vibration signal that was generated using the vibration separation (synchronous signal-averaging) algorithm. Each waveform is separated into analysis domains using zeros of its slope and curvature. The bases selected in each analysis domain are chosen to minimize the prediction error, and constrained to have the same-sign local slope and curvature as the original signal. The resulting set of bases is used to analyze future-state vibration signals and the lifting prediction error is inspected. The constraints allow the transform to effectively adapt to global amplitude changes, yielding small prediction errors. However, local wave-form changes associated with certain types of gear damage are poorly adapted, causing a significant change in the prediction error. The constrained adaptive lifting diagnostic algorithm is validated using data collected from the University of Maryland Transmission Test Rig and the results are discussed.

Lewicki, David G.↗

Relativistic calculation of atomic M-shell ionization by protons

Relativistic plane-wave Born-approximation calculations of cross sections for M-shell ionization of Ho-67, Au-79, U-92 by protons with incident energies from 0.05 to 1 MeV are reported. Dirac-Hartree-Slater wave functions were employed and binding-energy change and Coulomb deflection were taken into account. Associated X-ray production cross sections were also computed. Results are compared with previous theoretical predictions and with experimental data. Definite improvement in the theory has been attained by the use of realistic wave functions and consistent inclusion of the effects of relativity.

Chen, M. H.↗

A review of computational methods for studying oscillating water columns – the Navier-Stokes based equation approach

This review evaluates the state-of-the-practice numerical tools used to predict the performance of Oscillating Water Column (OWC). The OWC is a widely studied Wave Energy Converter that provides a reliable form of renewable form of electricity that can potentially meet global energy needs. However, the fluid-flow phenomena affecting its hydrodynamic performance are not fully understood. While there has been the successful full-scale deployment of OWCs, various computational methods are being explored to optimize this technology. Potential flow theory is commonly used to evaluate the efficiency of OWCs; however, this assumption tends to over-predict the hydrodynamic performance. Recently, numerical studies using a diverse set of commercial, open-source, or in-house Computational Fluid Dynamics (CFD) using Reynolds Averaged Navier Stokes (RANS) and Large-Eddy Simulation codes show a better comparison to available experimental results but are computationally expensive. ANSYS Fluent was found to be the most widely used CFD code applied to the study of the OWCs, with a high degree of accuracy in terms of experimental validation of numerical results.

16 TIDAL AND WAVE POWER↗

Public Reference Data for Megawatt-Scale Hydrogen Electrolysis - Simulated Wave

The U.S. Department of Energy and the National Laboratory of the Rockies (NLR) demonstrate hydrogen electrolysis, hydrogen compression and storage, and variable hydrogen fuel cell power production using megawatt-scale equipment at NLR’s Flatirons Campus as part of the Advanced Research on Integrated Energy Systems (ARIES) initiative. This dataset represents part of that effort and is intended for academic, national laboratory, industrial, and other stakeholders to plan, design, and validate models of megawatt-scale hydrogen technologies and diverse energy infrastructure nationwide. These data provide a baseline for how existing hydrogen electrolysis technologies perform when coupled with various energy technologies. Future datasets will demonstrate how existing hydrogen fuel cell technologies can provide controllable, dispatchable, and variable power output for artificial intelligence (AI) data centers and other variable loads. This dataset entry describes hydrogen production using a single, simulated wave energy conversion device. The electrolyzer is a 1.25-MW proton exchange membrane type MC250 system manufactured by Nel Hydrogen. While the unit supports up to 2.5 MW of electrolysis, NLR only has a single 1.25-MW electrolysis stack. For the wave energy, NLR used a wave energy converter model from PacWave. These devices can be equipped with accumulators and pressure relief values to smooth the power output by storing and releasing hydraulic energy. Using a peak power output of 10 MW, the model created two 25-minute profiles: one with and one without the accumulators and pressure relief valves. To down select the profile data from the native resolution of 20 Hz to 1 Hz, NLR took the mean of every 20 data points. NLR experimented with two simulated wave energy power plants: one that peaks at 10 MW, and one that peaks at 5 MW. These profiles were scaled for the physical 1.25 MW electrolyzer by multiplying the original profiles by one eighth and one quarter, respectively. The first profile matches the capacity rating of eight of the 1.25 MW electrolyzers, while the second matches four electrolyzers. Finally, NLR experimented with two settings for the electrolyzer power supply minimum and maximum current ramp rates (gain and slew): 200 and 400 amperes per second. The simulated profiles were translated from power (kilowatts) to current (amperes) using a curve fit with calibration data and sent to the electrolyzer power supply at 1-Hz frequency. These datasets report relevant hydrogen balance-of-plant and system data, all captured at 1 Hz, including hydrogen mass production measured with an Emerson Coriolis flow meter. Each .zip file represents a single wave electrolysis experiment and is formatted as follows: {technology}-{accumulator?}_{number of 1.25 MW electrolyzers connected}-{electrolyzer ramp rate in amperes/second} For instance, “wavePacWave-Noacc_4-400.zip” represents the 25 minute-long experiment using the PacWave’s wave energy converter model, equipped with no accumulator, connected to four 1.25-MW electrolyzers with their power supplies set to a maximum current ramp rate (gain and slew) of 400 A/s. Each .zip folder contains the following files: A .csv file containing raw data. An .xlsx file explaining all the fields in the raw data. A .png plot showing the time series of hydrogen production in kilograms per hour, electrolysis power consumption, and input wave power. An experiment, labeled “characterization_200.zip”, demonstrates the MC250 electrolyzer steady-state response with 30 minute load steps for a total duration of 5 hours. Finally, a .csv file is provided with all wave profiles combined into one dataset labeled "combined_wave_experiments.csv". NLR also built an AI/machine-learning predictive model based on these datasets. The model ingests the electrolyzer current command in amperes, as well as various pressures and temperatures across the system, and predicts hydrogen output in kilograms per hour. The complete model can be found at https://huggingface.co/NatLabRockies/ptmelt-hydrogen-electrolysis.

08 HYDROGEN↗

Consideration of some factors affecting low-frequency fuselage noise transmission for propeller aircraft

Possible reasons for disagreement between measured and predicted trends of sidewall noise transmission at low frequency are investigated using simplified analysis methods. An analytical model combining incident plane acoustic waves with an infinite flat panel is used to study the effects of sound incidence angle, plate structural properties, frequency, absorption, and the difference between noise reduction and transmission loss. Analysis shows that these factors have significant effects on noise transmission but they do not account for the differences between measured and predicted trends at low frequencies. An analytical model combining an infinite flat plate with a normally incident acoustic wave having exponentially decaying magnitude along one coordinate is used to study the effect of a localized source distribution such as is associated with propeller noise. Results show that the localization brings the predicted low-frequency trend of noise transmission into better agreement with measured propeller results. This effect is independent of low-frequency stiffness effects that have been previously reported to be associated with boundary conditions.

Mixson, J. S.↗

Full Facility Shock Frame Simulations of the Electric Arc Shock Tube

Radiative heating computations are performed for high speed lunar return experiments conducted in the Electric Arc Shock Tube (EAST) facility at NASA Ames Research Center. The nonequilibrium radiative transport equations are solved via NASA's in-house radiation code NEQAIR using flow field input from US3D flow solver. The post-shock flow properties for the 10 km/s Earth entry conditions are computed using the stagnation line of a blunt-body and a full facility CFD (Computational Fluid Dynamics) simulation of the EAST shock tube. The shocked gas in the blunt-body flow achieves a thermochemical equilibrium away from the shock front whereas EAST flow exhibits a nonequilibrium behavior due to strong viscous dissipation of the shock by boundary layer. The full-tube flow calculations capture the influence of the boundary layer on the shocked gas state and provide a realistic fluid dynamic input for the radiative predictions. The integrated radiance behind the shock is calculated in NEQAIR for wavelength regimes from Vacuum-UltraViolet (VUV) to InfraRed (IR), which are pertinent to the emission characteristics of high enthalpy shock waves in air. These radiance profiles are validated against corresponding EAST shots. The full-tube simulations successfully predict a sharp radiance peak at the shock front which gets smeared in the test data due to the spatial resolution in the measurements. The full facility based radiance behind the shock shows a slightly better match with the test data in the VUV and Red spectral regions, as compared to that from a blunt-body based predictions. The UV radiance is very similar for both geometries and under-predicts the test behavior. The IR test data matches better with the blunt-body based predictions where the full-tube simulations show a significant over-prediction.

Radiation↗

Sea bottom topography imaging with SAR

It is well known that under favorable meteorological and hydrodynamical conditions the bottom topography of shallow seas can be mapped with airborne or spaceborne imaging radar. This phenomenon was observed for the first time in 1969 by de Loor and co-workers in Q-band Side Looking Airborne Radar (SLAR) imagery of sandwaves in the North Sea. It is now generally accepted that the imaging mechanism consists of three steps: (1) interaction between (tidal) current and bottom topography causes spatial modulations in the surface current velocity; (2) modulations in the surface current velocity give rise to variations in the spectrum of wind-generated waves, as described by the action balance equation; and (3) variations in the wave spectrum show up as intensity modulations in radar imagery. In order to predict radar backscatter modulations caused by sandwaves, an imaging model, covering the three steps, was developed by the Dutch Sea Bottom Topography Group. This model and some model results will be shown. On 16 Aug. 1989 an experiment was performed with the polarimetric P-, L-, and C-band synthetic aperture radar (SAR) of NASA/JPL. One scene was recorded in SAR mode. On 12 Jul. 1991 another three scenes were recorded, of which one was in the ATI-mode (Along-Track Interferometer). These experiments took place in the test area of the Sea Bottom Topography Group, 30 km off the Dutch coast, where the bottom topography is dominated by sand waves. In-situ data were gathered by a ship in the test area and on 'Measuring Platform Noordwijk', 20 km from the center of the test area. The radar images made during the experiment were compared with digitized maps of the bottom. Furthermore, the profiles of radar backscatter modulation were compared with the results of the model. During the workshop some preliminary results of the ATI measurements will be shown.

Vanderkooij, M. W. A.↗