Search NASA⌕ Search

SEARCH · Search NASA

Results for “wave modeling”

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 145 records · Page 8

Improving tide-estuary representation in MPAS-Ocean (Final Report)

The overarching goal of this research is to produce a skill-assessed tide-estuary representation, including BGC (biogeochemistry) components of the ecosystem, into the MPAS-O global climate modeling system. Specifically, in the original proposal we proposed one-way nesting, into global MPAS-O, of a high-resolution model for TES (tidal and estuary systems) including the coupling of the BGC components. We have completed the proposed tasks during the past 4 years. Specifically, after extensive consultation with the MPAS-O group during the first year, we have decided to go beyond the original proposal and to two-way nest SCHISM into MPAS-O. As a result, SCHISM-CoSiNE (with CoSiNE being the BGC model) is now fully integrated into MPAS-O as a subroutine (more precisely, an analysis member) and can be invoked as MPAS-O is running with two-way exchange of information between the TES and global models (with proper spatial interpolation done on the fly). As originally proposed, we have also developed the new finite-volume (FV) model MPAS-OI (aka ELCIRC-sub) for TES. We have conducted several tests using the coupled MPAS-SCHISM- CoSiNE model and also MPAS-OI for Columbia River, San Francisco Bay and Chesapeake Bay. Four journal publications have been produced so far as the result of the project and 1 is under review (cf. Publication section). With these achievements, we will be ready to conduct fully coupled simulations of TES and global ocean anywhere, as soon as some missing capabilities in MPAS-O (e.g. tides, improved ability to handle the boundary conditions, and streamlining of the global BGC model) are worked out. The significance of the outcome of this project cannot be overstated. It marks the first time, as far as we know, that a global climate model is tightly coupled to a bona fide nearshore and terrestrial model (SCHISM system); in fact SCHISM is the first bona fide compound flooding model that seamlessly covers the coastal ocean, estuaries, rivers and watersheds (Ye et al. 2020; Zhang et al. 2020). The coupled model will reduce the errors and uncertainties of the MPAS-O global climate model system. It also complements and extends on-going E3SM project, as the latter seeks in the long-term to include components required to simulate impacts of sea-level change and storm surge on coastal regions, including wave models and focusing resolution in coastal and storm-track regions, as an important component of the Cryosphere-ocean system. Note that a spectral wave model (the unstructured-grid version of Wave Watch III, aka Wind Wave Model) is already incorporated inside the SCHISM system and is therefore directly accessible in the coupled system.

58 GEOSCIENCES↗

The Atmospheric Effects of HSCT Emissions Simulated by a 3-Wave Interactive Model

An interactive model which couples a semi-spectral dynamical model, a radiative transfer code and a two-dimensional chemistry transport model (2-D CTM), is used to assess the atmospheric effects of the High-Speed Civil Transport (HSCT) engine emissions. The residual mean meridional circulation, the zonal-mean temperature and the eddy diffusion coefficients are calculated using zonal means and three longest zonal waves of dynamical variables integrated in the semi-spectral dynamical model. They are used in the 2-D CTM to simulate the distribution of trace gases in the atmosphere. The simulated ozone is sent to the radiative transfer code to calculate the heating rates, which drive the dynamics. This radiative coupling connects the dynamical and photochemical processes and creates feedback when the atmosphere is perturbed. It is found that in most areas the ozone depletion caused by HSCT emissions calculated using the 3-wave model has the features similar to, but with significantly larger magnitude than that calculated by the AER 2-D CTM with prescribed transport parameters and temperature. The difference is mostly due to the differences in the circulation in the two models. The radiative feedback effects are investigated by comparing the ozone depletion calculated with the baseline dynamics and with the dynamics perturbed by the HSCT emissions. The feedback through changes in the residual mean meridional circulation and the eddy diffusion coefficients has moderate effects on the simulated ozone depletion. It reduced the ozone depletion by 20-30% in northern mid and high-latitudes. However, the feedback through changes in the zonal- mean temperature is negligible.

Shia, R.-L.↗

From Safer Space Travel to Improved Fisheries Monitoring: Wave Energy Modeling Tool Helps Design Tomorrow's Tech

What does wave energy have to do with space exploration? It turns out, plenty. The National Aeronautics and Space Administration (NASA) is working toward launching the Artemis I mission, an uncrewed lunar test flight, in 2021. It will be the maiden flight of the Space Launch System heavy-lift launch vehicle, as well as the Orion crew module (CM). Orion will orbit the moon for multiple days before returning to Earth and landing in the Pacific Ocean.

Tannen S Vanzwieten↗

Effects of Surface Turbulence Flux Parameterizations on the MJO: The Role of Ocean Surface Waves

This study investigates the sensitivity of the Madden–Julian oscillation (MJO) to changes to the bulk flux parameterization and the role of ocean surface waves in air–sea coupling using a fully coupled ocean–atmosphere–wave model. The atmospheric and ocean model components of the Energy Exascale Earth System Model (E3SM) are coupled to a spectral wave model, WAVEWATCH III (WW3). Two experiments with wind speed–dependent bulk algorithms (NCAR and COARE3.0a) and one experiment with wave-state-dependent flux (COR3.0a-WAV) were conducted. We modify COARE3.0a to include surface roughness calculated within WW3 and also account for the buffering effect of waves on the relative difference between air-side and ocean-side momentum flux. Differences in surface fluxes, primarily caused by discrepancies in drag coefficients, result in significant differences in MJO’s properties. While COARE3.0a has better convection–circulation coupling than NCAR, it exhibits anomalous MJO convection east of the date line. The wave-state-dependent flux (COR3.0-WAV) improves the MJO representation over the default COARE3.0 algorithm. Strong easterlies over the Pacific Ocean in COARE3.0a enhance the latent heat flux (LHFLX). This is responsible for the anomalous MJO propagation after the date line. In COR3.0a-WAV, waves reduce the anomalous easterlies, leading to a decrease in LHFLX and MJO dissipation after the date line. These findings highlight the role of surface fluxes in MJO simulation fidelity. Most importantly, we show that the proper treatment of wave-induced effects in bulk flux parameterization improves the simulation of coupled climate variability.

54 ENVIRONMENTAL SCIENCES↗

Model of Wave Driven Flow Oscillation for Solar Cycle

At low latitudes in the Earth's atmosphere, the observed zonal flow velocities are dominated by the semi-annual and quasi-biennial oscillations with periods of 6 months and 20 to 32 months respectively. These terrestrial oscillations, the SAO and QBO respectively, are driven by wave-mean flow interactions due to upward propagating planetary-scale waves (periods of days) and small-scale gravity waves (periods of hours). We are proposing (see also Mayr et al., GRL, 2001) that such a mechanism may drive long period oscillations (reversing flows) in stellar and planetary interiors, and we apply it to the Sun. The reversing flows would occur below the convective envelope where waves can propagate. We apply a simplified, one dimensional, analytical flow model that incorporates a gravity wave parameterization due to Hines (1997). Based on this analysis, our estimates show that relatively small wave amplitudes less than 10 m/s can produce zonal flow amplitudes of 20 m/s, which should be sufficient to generate the observed variations in the magnetic field. To produce the 22-year period of oscillation, a low buoyancy frequency must be chosen, and this places the proposed flow in a region that is close to (and below) the base of the convective envelope. Enhanced turbulence associated with this low stability should help to generate the dynamo currents. With larger stability at deeper levels in the solar interior, the model can readily produce also oscillations with much longer periods. To provide an understanding of the fluid dynamics involved, we present numerical results from a 2D model for the terrestrial atmosphere that exemplify the non-linear nature of the wave interaction for which a mechanical analog is the escapement mechanism of the clock.

Mayr, Hans G.↗

Deep Learning for Full Waveform Inversion of Elastic Active-Source Seismic Data to Estimate P-Wave Velocity Models

Seismic imaging methods are critical for Global Security and Energy & Homeland Security missions and activities that rely on subsurface characterization, but traditional methods remain computationally expensive and require significant labor hours and expertise to execute. Within the past few years, machine learning (ML), namely deep learning (DL), has been used to develop data-driven end-to-end full waveform inversion (FWI) methods to estimate 2D P-wave velocity (Vp) models in a fraction of the time as conventional FWI. These methods, however, are trained on simplistic acoustic wave seismic data and Vp models that are not realistic nor representative of real-world observations, leaving a large gap between the state-of-the-art and deployable, feasible, and practical DL FWI methods. Here, we generate a synthetic active-source, 3D, elastic wave seismic data set and a variety of Vp models with realistic geologic structure for training DL FWI methods. We evaluate six different methods that have performed well for acoustic DL FWI or medical imaging tasks using our more realistic dataset. We find that these six trained models do not match the performance of published acoustic end-to-end DL FWI methods, indicating more training data may be needed, physics may need to be incorporated to achieve good accuracy at the sacrifice of the end-to-end advantage, and/or novel methods need to be developed to enable end-to-end DL FWI methods to perform well for real-world seismic data.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Lower stratosphere waves during 1986-1989 southern springs

Global NMC height analyses and Nimbus 7 TOMS Grid-T Version 6 data are used to investigate the horizontal structure of ozone and 50-mb height waves for the 1986-1989 springs of the Southern Hemisphere. Empirical orthogonal functions (EOFs) are obtained for the ozone and height data independently to identify the meridional structure of observed long waves as well as their propagation characteristics. The leading EOFs for both the ozone and 50-mb height waves exhibit wavenumber 1 and 2 patterns. The correlation coefficients between ozone and 50-mb height principal components are high. They indicate high/low ozone values associated with high/low values of 50-mb heights. This distribution is consistent with the advection of ozone by planetary waves south of the ozone maximum, as shown by a simple linear wave model. The observed total ozone amplitudes for waves 1 and 2 for most years are consistent with the total ozone predictions by the wave model.

Nogues-Paegle, Julia↗

Recent improvements to and validation of the one dimensional NASA wave rotor model

A numerical model has been developed at the NASA Lewis Research Center which can predict both the unsteady flow quantities within a wave rotor passage and the steady averaged flows in the ports. The model is based on the assumptions of one-dimensional, unsteady, perfect gas flow. The model assesses not only the dominant wave behavior, but the loss effects of finite passage opening time, leakage from the passage ends, viscosity, and heat transfer to and from the passages. The model operates in the rotor reference frame; however, until recently no account was made for the often significant effect of the rotor circumferential velocity component. The present model accounts for this by modifying the passage boundary conditions, allowing the internal computational scheme to remain the rotor reference frame, while quantities such as inlet duct stagnation properties may be specified in the fixed or absolute reference frame. Accurate modeling of this effect is critical to successful wave rotor analysis and design, particularly in off-design predictions where the flows in the inlet ducts are mismatched with the rotor passages and significant turning may take place (i.e., work is done on the gas). The relative simplicity of the model makes it useful for design and optimization, as well as analysis, of wave rotor cycles for many applications. This report, building on several earlier papers, describes the most recent modifications to the model. These include accounting for the relative/absolute transition at the passage boundaries and refinements to the viscous source term correlation which resulted from this accounting. Comparison of model predictions with measured data is then presented and discussed.

Paxson, Daniel E.↗

A finite element model for wave propagation in an inhomogeneous material including experimental validation

A finite element model was developed to solve for the acoustic pressure field in a nonhomogeneous region. The derivations from the governing equations assumed that the material properties could vary with position resulting in a nonhomogeneous variable property two-dimensional wave equation. This eliminated the necessity of finding the boundary conditions between the different materials. For a two media region consisting of part air (in the duct) and part bulk absorber (in the wall), a model was used to describe the bulk absorber properties in two directions. An experiment to verify the numerical theory was conducted in a rectangular duct with no flow and absorbing material mounted on one wall. Changes in the sound field, consisting of planar waves was measured on the wall opposite the absorbing material. As a function of distance along the duct, fairly good agreement was found in the standing wave pattern upstream of the absorber and in the decay of pressure level opposite the absorber.

Baumeister, Kenneth J.↗

A finite element model for wave propagation in an inhomogeneous material including experimental validation

A finite element model was developed to solve for the acoustic pressure field in a nonhomogeneous region. The derivations from the governing equations assumed that the material properties could vary with position resulting in a nonhomogeneous variable property two-dimensional wave equation. This eliminated the necessity of finding the boundary conditions between the different materials. For a two media region consisting of part air (in the duct) and part bulk absorber (in the wall), a model was used to describe the bulk absorber properties in two directions. An experiment to verify the numerical theory was conducted in a rectangular duct with no flow and absorbing material mounted on one wall. Changes in the sound field, consisting of planar waves, was measured on the wall opposite the absorbing material. As a function of distance along the duct, fairly good agreement was found in the standing wave pattern upstream of the absorber and in the decay of pressure level opposite the absorber.

Baumeister, Kenneth J.↗

Simulating Long–Term Dynamics of Radiation Belt Electrons Using DREAM3D Model

Here, we compared the performance of DREAM3D simulations in reproducing the long-term radiation belt dynamics observed by Van Allen Probes over the entire year of 2017 with various boundary conditions (BCs) and model inputs. Specifically, we investigated the effects of three different outer boundary conditions, two different low-energy boundary conditions for seed electrons, four different radial diffusion (RD) coefficients (DLL), four hiss wave models, and two chorus wave models from the literature. Using the outer boundary condition driven by GOES data, our benchmark simulation generally well reproduces the observed radiation belt dynamics inside L* = 6, with a better model performance at lower μ than higher μ, where μ is the first adiabatic invariant. By varying the boundary conditions and inputs, we find that: (a) The data-driven outer boundary condition is critical to the model performance, while adding in the data-driven seed population doesn't further improve the performance. (b) The model shows comparable performance with DLL from Brautigam and Albert (2000, https://doi.org/10.1029/1999ja900344), Ozeke et al. (2014, https://doi.org/10.1002/2013ja019204), and Liu et al. (2016, https://doi.org/10.1002/2015gl067398), while with DLL from Ali et al. (2016, https://doi.org/10.1002/2016ja023002) the model shows less RD compared to data. (c) The model performance is similar with data-based hiss models, but the results show faster loss is still needed inside the plasmasphere. (d) The model performs similarly with the two different chorus models, but better capturing the electron enhancement at higher μ using the Wang et al. (2019, https://doi.org/10.1029/2018ja026183) model due to its stronger wave power, since local heating for higher energy electrons is under-reproduced in the current model.

79 ASTRONOMY AND ASTROPHYSICS↗

A non-Gaussian statistical model for surface elevation of nonlinear random wave fields

Probability density function of the surface elevation of a nonlinear random wave field is obtained. The wave model is based on the Stokes expansion carried to the third order for both deep water waves and waves in finite depth. The amplitude and phase of the first-order component of the Stokes wave are assumed to be Rayleigh and uniformly distributed and slowly varying, respectively. The probability density function for the deep water case was found to depend on two parameters: the root-mean-square surface elevation and the significant slope. For water of finite depth, an additional parameter, the nondimensional depth, is also required. An important difference between the present result and the Gram-Charlier representation is that the present probability density functions are always nonnegative. It is also found that the 'constant' term in the Stokes expansion, usually neglected in deterministic studies, plays an important role in determining the details of the density function. The results compare well with laboratory and field experiment data.

Huang, N. E.↗

Observations of ULF oscillations in the ion fluxes at small pitch angles with ATS 6

The ultra-low-frequency modulation of ion flux densities at small pitch angles observed by ATS 6 is examined, with particular attention given to a detailed analysis of a representative event. ULF modulation events with maximum modulation at small pitch angles were identified 14 times during the first eight months of operation of the NOAA low-energy particle detector on ATS 6. For the event of October 23, 1974, maximum flux modulation, with a maximum/minimum intensity ratio of 3.7, was observed in the 100 to 150 keV detector at an angle of 32 deg to the ambient field. Spectral analysis of magnetic field data reveals a right elliptically polarized magnetic perturbation with a 96-sec period and a 5-gamma rms amplitude, propagating in the dipole meridian at an angle of about 15 deg to the ambient field and the dipole axis. Proton flux modulation is found to lag the field by up to 180 deg for the lowest-energy channel. Observations are compared with the drift wave, MHD slow wave, and bounce resonant interaction associated with transverse wave models, and it is found that none of the wave models can adequately account for all of the correlated particle and field oscillations.

Su, S.-Y.↗

Spiral structure of M51: Streaming motions across the spiral arms

The atomic hydrogen (HI) and the H alpha emission line in the grand-design spiral galaxy M51 have been observed with the Westerbork Synthesis Radio Telescope and the Taurus Fabry-Perot imaging spectrometer, respectively. Across the inner spiral arms significant tangential and radial velocity gradients are detected in the H alpha emission after subtraction of the axi-symmetric component of the velocity field. The shift is positive on the inside and negative on the outside of the northern arm. Across the southern arm this situation is reversed. The direction of the shifts is such that the material is moving inward and faster compared to circular rotation in both arms, consistent with the velocity perturbations predicted by spiral density wave models for gas downstream of a spiral shock. The observed shifts amount to 20 to 30 km (s-1), corresponding to streaming motions of 60 to 90 km (s-1) in the plane of the disk (inclination angle 20 degrees). Comparable velocity gradients have also been observed by Vogel et al. in the CO emission from the inner northern arm of M51. The streaming motions in M51 are about 2 to 3 times as large as the ones found in HI by Rots in M81, and successfully modelled by Visser with a self-consistent density wave model. Researchers have not been able to detect conclusively streaming motions in the HI emission from the arms, perhaps due to the relatively poor angular resolution (approx. 15 seconds) of the HI observations.

Tilanus, R. P. J.↗

Accurate Time-Dependent Traveling-Wave Tube Model Developed for Computational Bit-Error-Rate Testing

The phenomenal growth of the satellite communications industry has created a large demand for traveling-wave tubes (TWT's) operating with unprecedented specifications requiring the design and production of many novel devices in record time. To achieve this, the TWT industry heavily relies on computational modeling. However, the TWT industry's computational modeling capabilities need to be improved because there are often discrepancies between measured TWT data and that predicted by conventional two-dimensional helical TWT interaction codes. This limits the analysis and design of novel devices or TWT's with parameters differing from what is conventionally manufactured. In addition, the inaccuracy of current computational tools limits achievable TWT performance because optimized designs require highly accurate models. To address these concerns, a fully three-dimensional, time-dependent, helical TWT interaction model was developed using the electromagnetic particle-in-cell code MAFIA (Solution of MAxwell's equations by the Finite-Integration-Algorithm). The model includes a short section of helical slow-wave circuit with excitation fed by radiofrequency input/output couplers, and an electron beam contained by periodic permanent magnet focusing. A cutaway view of several turns of the three-dimensional helical slow-wave circuit with input/output couplers is shown. This has been shown to be more accurate than conventionally used two-dimensional models. The growth of the communications industry has also imposed a demand for increased data rates for the transmission of large volumes of data. To achieve increased data rates, complex modulation and multiple access techniques are employed requiring minimum distortion of the signal as it is passed through the TWT. Thus, intersymbol interference (ISI) becomes a major consideration, as well as suspected causes such as reflections within the TWT. To experimentally investigate effects of the physical TWT on ISI would be prohibitively expensive, as it would require manufacturing numerous amplifiers, in addition to acquiring the required digital hardware. As an alternative, the time-domain TWT interaction model developed here provides the capability to establish a computational test bench where ISI or bit error rate can be simulated as a function of TWT operating parameters and component geometries. Intermodulation products, harmonic generation, and backward waves can also be monitored with the model for similar correlations. The advancements in computational capabilities and corresponding potential improvements in TWT performance may prove to be the enabling technologies for realizing unprecedented data rates for near real time transmission of the increasingly larger volumes of data demanded by planned commercial and Government satellite communications applications. This work is in support of the Cross Enterprise Technology Development Program in Headquarters' Advanced Technology & Mission Studies Division and the Air Force Office of Scientific Research Small Business Technology Transfer programs.

Kory, Carol L.↗

Quantifying modeling uncertainties when combining multiple gravitational-wave detections from binary neutron star sources

Here, with the increasing sensitivity of gravitational-wave detectors, we expect to observe multiple binary neutron-star systems through gravitational waves in the near future. The combined analysis of these gravitational-wave signals offers the possibility to constrain the neutron-star radius and the equation of state of dense nuclear matter with unprecedented accuracy. However, it is crucial to ensure that uncertainties inherent in the gravitational-wave models will not lead to systematic biases when information from multiple detections is combined. To quantify waveform systematics, we perform an extensive simulation campaign of binary neutron-star sources and analyze them with a set of four different waveform models. For our analysis with 38 simulations, we find that statistical uncertainties in the neutron-star radius decrease to ± 250 m (2% at 90% credible interval) but that systematic differences between currently employed waveform models can be twice as large. Hence, it will be essential to ensure that systematic biases will not become dominant in inferences of the neutron-star equation of state when capitalizing on future developments.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Sunspot models with Alfven wave emission

Sunspot models have been computed on the assumption that the missing flux is transported by undissipated Alfven waves. In order to estimate the flux of these waves, we propose an extension of Opik's cellular model of convection to include the effects of a vertical magnetic field on horizontal gas flow. Horizontal motions, not vertical motions, are impeded more or less severely depending on the electrical conductivity, and this reduces the convective flux. These motions, however, shake the field lines; this is assumed to be a source of Alfven waves, compensating for the reduction in the convective flux. The free parameter D/H (ratio of cell diameter to cell depth) is adjusted such that the total sum of radiative, convective, and Alfven wave fluxes remains constant at all depths and equals the undisturbed solar flux entering the spot from below. Our treatment is speculative to some extent, but it leads to a unique spot model for a given field strength.

Mullan, D. J.↗

Time-Dependent Traveling Wave Tube Model for Intersymbol Interference Investigations

For the first time, a computational model has been used to provide a direct description of the effects of the traveling wave tube (TWT) on modulated digital signals. The TWT model comprehensively takes into account the effects of frequency dependent AM/AM and AM/PM conversion, gain and phase ripple; drive-induced oscillations; harmonic generation; intermodulation products; and backward waves. Thus, signal integrity can be investigated in the presence of these sources of potential distortion as a function of the physical geometry of the high power amplifier and the operational digital signal. This method promises superior predictive fidelity compared to methods using TWT models based on swept-amplitude and/or swept-frequency data. The fully three-dimensional (3D), time-dependent, TWT interaction model using the electromagnetic code MAFIA is presented. This model is used to investigate assumptions made in TWT black-box models used in communication system level simulations. In addition, digital signal performance, including intersymbol interference (ISI), is compared using direct data input into the MAFIA model and using the system level analysis tool, SPW.

Kory, Carol L.↗