Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

New Developments and Capabilities Within WEC-Sim: Preprint

WEC-Sim is an open-source software for simulating wave energy converters, which has been actively developed and applied since its initial release in 2014 to simulate a wide variety of device archetypes. WEC-Sim is developed jointly by the National Renewable Energy Laboratory (NREL) and Sandia National Laboratories (Sandia) within the MATLAB/SIMULINK environment. A general wave-to-wire model begins with a deployment site resource characterization, which is used to complete the hydrodynamic simulation of wave energy converters (WEC), with the power generation profile imported to a grid simulator to understand the influence on the local electrical network. While modeling the entire wave-to-wire is difficult and encompasses multiple time scales and physics, WEC-Sim is focused on the hydrodynamics simulation to predict, analyze, and optimize WEC dynamics and power performance. WEC-Sim simulations are performed in the time domain based on the radiation and diffraction method using hydrodynamics coefficients derived from boundary element method (BEM)-based frequency-domain potential flow solvers (e.g., WAMIT, NEMOH, Capytaine, or ANSYS-AQWA). With this level of modeling fidelity, WEC-Sim can handle floating body hydrodynamics, mechanical and electrical power generation methods, advanced control implementation, mooring systems, and other unique applications such as desalination. Additional WEC-Sim functionalities include pre-built Simulink blocks and MATLAB scripts that can simulate a wide range of floating systems and the corresponding auxiliary subsystems. The developers of WEC-Sim continue to release new versions of the software, at least annually, with our latest release in September 2022. These releases include bug fixes, updates to software documentation, as well as new features to expand WEC-Sim's capabilities to model a wide range of WEC concepts. This publication will highlight the new features added to WEC-Sim between versions 4.1.0 to 5.0.1 which spans over a two year period from June 2020 to September 2022. New features to be described will include topics such as continuous integration checks, revised Morison Element and nonlinear hydro implementations, run directly from Simulink (required for hardware-in-the-loop execution), BEMIO updates to import Capytaine BEM hydrodynamics, addition of cable blocks, and new wave visualization features.

TIDAL AND WAVE POWER↗
NASA NTRS1986

Three-dimensional, time-dependent, MHD model of a solar flare-generated interplanetary shock wave

A three-dimensional time-dependent MHD model of the propagation of an interplanetary shock wave into an ambient three-dimensional heliospheric solar wind is initialized with a peak velocity of 1000 km/s at the center of a right circular cone of 18 deg included angle at 18 solar radii. Differences from a previous 2-1/2 simulation (Wu et al., 1983; Gislason et al., 1984; Dryer et al., 1984) include diminuation of the solar peak velocity and concentration of the peak density at each radius. The IMF magnitude starts with high-latitude peaks, and helical-like IMF rotation is noted due to a large-amplitude nonlinear Alfven wave in the shocked plasma.

Dryer, M.↗
NASA NTRS2013

Analyzing Tropical Waves Using the Parallel Ensemble Empirical Model Decomposition Method: Preliminary Results from Hurricane Sandy

In this study, we discuss the performance of the parallel ensemble empirical mode decomposition (EMD) in the analysis of tropical waves that are associated with tropical cyclone (TC) formation. To efficiently analyze high-resolution, global, multiple-dimensional data sets, we first implement multilevel parallelism into the ensemble EMD (EEMD) and obtain a parallel speedup of 720 using 200 eight-core processors. We then apply the parallel EEMD (PEEMD) to extract the intrinsic mode functions (IMFs) from preselected data sets that represent (1) idealized tropical waves and (2) large-scale environmental flows associated with Hurricane Sandy (2012). Results indicate that the PEEMD is efficient and effective in revealing the major wave characteristics of the data, such as wavelengths and periods, by sifting out the dominant (wave) components. This approach has a potential for hurricane climate study by examining the statistical relationship between tropical waves and TC formation.

PEEMD↗
NASA NTRS1984

Strong turbulence of plasma waves

This paper reviews recent work related to modulational instability and wave envelope self-focusing in dynamical and statistical systems. After introductory remarks pertinent to nonlinear optics realizations of these effects, the author summarizes the status of the subject in plasma physics, where it has come to be called 'strong Langmuir turbulence'. The paper treats the historical development of pertinent concepts, analytical theory, numerical simulations, laboratory experiments, and spacecraft observations. The role of self-similar self-focusing Langmuir envelope wave packets is emphasized, both in the Zakharov equation model for the wave dynamics and in a statistical theory based on this dynamical model.

Goldman, M. V.↗
NASA NTRS1993

Nonhydrostatic effects in numerical modeling of mesoscale convective systems and baroclinic waves

The present investigation is concerned with the role of convection upon mesoscale modeling results, particularly when the grid resolution becomes small enough that there is not a clear scale separation between the explicitly resolved circulations and the parameterized clouds. In those situations, the vertical accelerations in explicitly resolved circulations become strong enough that the hydrostatic assumption may no longer be valid. These concerns arise from interests in improving mesoscale modeling per se and in improving the subgrid-scale parameterizations in global models. The hydrostatic and the nonhydrostatic options of the Colorado State University Regional Atmospheric Modeling System were used to simulate dry gravity currents in two dimensions, using several different horizontal grid sizes. With horizontal grid intervals of 10 km or less, nonhydrostatic simulations produce wider and colder heads and weaker but wider forced updrafts than do the hydrostatic simulations. Comparing the hydrostatic and nonhydrostatic models show that the difference between the vertical mass fluxes is much less than the difference between the vertical velocities. When the grid is fine enough to resolve the head of the gravity current, horizontal convergence at the gust front extends upwards almost to the head of the cold air. Vertical mass flux in the forced updraft at the front varies with horizontal grid size mainly as a function of the height of the simulated head. For coarser grids, which do not resolve the head, vertical mass flux at all heights decreases with increasing horizontal grid size. A comparison on nonhydrostatic simulations with horizontal grid intervals of 1 km and 2 km illustrates how decreasing the grid size does not necessarily increase the intensity of the resolved circulation. The smaller grid enables the simulated gravity current to entrain a bubble of warm air behind the head, which results in a weaker circulation with a shorter head and weaker updraft.

Cohen, Charles↗
NASA NTRS1988

The GISS global climate-middle atmosphere model. II - Model variability due to interactions between planetary waves, the mean circulation and gravity wave drag

The variability which arises in the GISS Global Climate-Middle Atmosphere Model on two time scales is reviewed: interannual standard deviations, derived from the five-year control run, and intraseasonal variability as exemplified by statospheric warnings. The model's extratropical variability for both mean fields and eddy statistics appears reasonable when compared with observations, while the tropical wind variability near the stratopause may be excessive possibly, due to inertial oscillations. Both wave 1 and wave 2 warmings develop, with connections to tropospheric forcing. Variability on both time scales results from a complex set of interactions among planetary waves, the mean circulation, and gravity wave drag. Specific examples of these interactions are presented, which imply that variability in gravity wave forcing and drag may be an important component of the variability of the middle atmosphere.

Rind, D.↗
DOE OSTI2020

Ground state wave functions for single-band Hubbard models from the Gutzwiller conjugate gradient minimisation theory

The Gutzwiller conjugate gradient minimisation (GCGM) theory is an ab initio quantum many-body theory for computing the ground-state properties of infinite systems. Previous applications of GCGM provides satisfying accuracy of ground-state energy of Hubbard models. In the current work, we address the problem of whether the obtained wave function is a good approximation for the true ground state by comparing the correlation functions with the benchmark data. Additionally, our results confirms the accuracy of the reproduced ground state of the regular Hubbard model, but with some exception for the frustrated Hubbard model.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2019

High-Resolution Regional Wave Hindcast for the U.S. Alaska Coast

Regional long-term wave hindcast with state-of-art third-generation models is essential for characterization of regional wave resources. This report summarizes modeling efforts for the simulation of the wave climate along southern Alaska, which used an unstructured, nested-grid modeling approach that incorporates a global-regional nested grid using WAVEWATCH III ® and the high-resolution SWAN (Simulating Waves Nearshore) model. Wave resource and wave bulk parameters were simulated for a 32-year period from 1979 to 2010 and were subsequently validated with wave buoy data in the model domain. Error statistics of model skills were calculated at all buoy stations. Overall, model results, both resource and bulk parameters, match observations well. The model hindcast was able to reproduce the seasonal variation of the sea state with large waves that occur in the winter and early spring months when wind forcing is strong and the calm sea state during the summer when wind forcing is weak. The hindcast also considers the effect of sea ice on wave growth, which affects the southern Bering Sea. The hindcast results, including resource and bulk parameters, were archived at every model grid point for the entire model domain at three-hour intervals. In addition, spectral results were saved at National Data Buoy Center buoy stations and selected virtual buoy locations at hourly intervals. The nested-grid modeling framework employed in this study provides a powerful and efficient modeling approach for accurately simulating wave climate at regional and long-term temporal scales with sufficiently fine resolutions in the nearshore region.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1996

A study of dusty nonradiative shock waves

A dynamic theoretical model of grain evolution through steady-state nonradiative shock waves in the Cygnus Loop and Puppis A supernova remnants was constructed. Data from the Rosat position sensitive proportional counter (PSPC) and the Infrared Astronomical Satellite (IRAS) were used. A theoretical model of the grain evolution was developed by using the steady-state nonradiative shocks. Predictions were made concerning the X-ray spectrum emanating from a post-shock plasma were made. The fitting of the model predictions to the spectral shape of the Cygnus Loop X-ray emissions is discussed.

Vancura, O.↗
NASA NTRS2005

Ring Current Dynamic in the Presence of EMIC Waves

The effect of EMIC waves, generated by a positive ion temperature anisotropy on Earth s RC ions is one of the best known examples of wave-particle interaction in the magnetosphere and the most controversial mechanism of RC losses. Under certain conditions, relativistic electrons with energy 21 MeV can be removed from the outer RB by EMIC wave scattering during a magnetic storm much faster than by any other loss mechanisms. That is why the calculation of EMIC waves is a very critical part of the NASA LWS program. The new RC model that we have developed and present for the first time has several new features that we have combine together in a one single model: (a) several lower frequency cold plasma wave modes are taken into account; (b) wave tracing of these wave has been incorporated in the energy EMIC wave equation; (c) no assumptions regarding wave shape spectra have been made; (d) no assumptions regarding the shape of particle distribution have been made to calculate the growth rate; (e) pitch- angle, energy, and mix diffusions are taken into account together for the first time; (f) the exact loss-cone RC analytical solution has been found and coupled with bounce-averaged numerical solution of kinetic equation; (g) LHW are included as an additional factor that contributes to saturation process of EMIC waves; and (h) the hot ions were included in the real part of dielectric permittivity tensor. We compare our theoretical results with the different EMIC waves models as well as RC experimental data.

Khazanov, G. V.↗
NASA NTRS2009

Self-Consistent Model of Magnetospheric Electric Field, Ring Current, Plasmasphere, and Electromagnetic Ion Cyclotron Waves: Initial Results

Further development of our self-consistent model of interacting ring current (RC) ions and electromagnetic ion cyclotron (EMIC) waves is presented. This model incorporates large scale magnetosphere-ionosphere coupling and treats self-consistently not only EMIC waves and RC ions, but also the magnetospheric electric field, RC, and plasmasphere. Initial simulations indicate that the region beyond geostationary orbit should be included in the simulation of the magnetosphere-ionosphere coupling. Additionally, a self-consistent description, based on first principles, of the ionospheric conductance is required. These initial simulations further show that in order to model the EMIC wave distribution and wave spectral properties accurately, the plasmasphere should also be simulated self-consistently, since its fine structure requires as much care as that of the RC. Finally, an effect of the finite time needed to reestablish a new potential pattern throughout the ionosphere and to communicate between the ionosphere and the equatorial magnetosphere cannot be ignored.

Gamayunov, K. V.↗
NASA NTRS1996

A Simulated Spectrum of Convectively Generated Gravity Waves: Propagation from the Tropopause to the Mesopause and Effects on the Middle Atmosphere

This work evaluates the interaction of a simulated spectrum of convectively generated gravity waves with realistic middle atmosphere mean winds. The wave spectrum is derived from the nonlinear convection model described by Alexander et al. that simulated a two-dimensional midlatitude squall line. This spectrum becomes input to a linear ray tracing model for evaluation of wave propagation as a function of height through climatological background wind and buoyancy frequency profiles. The energy defined by the spectrum as a function of wavenumber and frequency is distributed spatially and temporally into wave packets for the purpose of estimating wave amplitudes at the lower boundary of the ray tracing model. A wavelet analysis provides an estimate of these wave packet widths in space and time. Without this redistribution of energies into wave packets the Fourier analysis alone inaccurately assumes the energy is evenly distributed throughout the storm model domain. The growth with height of wave amplitudes is derived from wave action flux conservation coupled to a convective instability saturation condition. Mean flow accelerations and wave energy dissipation profiles are derived from this analysis and compared to parameterized estimates of gravity wave forcing, providing a measure of the importance of the storm source to global gravity wave forcing. The results suggest that a single large convective storm system like the simulated squall line could provide a significant fraction of the zonal mean gravity wave forcing at some levels, particularly in the mesosphere. The vertical distributions of mean flow acceleration and energy dissipation do not much resemble the parameterized profiles in form because of the peculiarities of the spectral properties of the waves from the storm source. The ray tracing model developed herein provides a tool for examining the role of convectively generated waves in middle atmosphere physics.

Alexander, Joan↗
NASA NTRS1996

A Simulated Spectrum of Convectively Generated Gravity Waves: Propagation from the Tropopause to the Mesopause and Effects on the Middle Atmosphere

This work evaluates the interaction of a simulated spectrum of convectively generated gravity waves with realistic middle atmosphere mean winds. The wave spectrum is derived from the nonlinear convection model described by Alexander et al. [1995] that simulated a two-dimensional midlatitude squall line. This spectrum becomes input to a linear ray tracing model for evaluation of wave propagation as a function of height through climatological background wind and buoyancy frequency profiles. The energy defined by the spectrum as a function of wavenumber and frequency is distributed spatially and temporally into wave packets for the purpose of estimating wave amplitudes at the lower boundary of the ray tracing model. A wavelet analysis provides an estimate of these wave packet widths in space and time. Without this redistribution of energies into wave packets the Fourier analysis alone inaccurately assumes the energy is evenly distributed throughout the storm model domain. The growth with height of wave amplitudes is derived from wave action flux conservation coupled to a convective instability saturation condition. Mean flow accelerations and wave energy dissipation profiles are derived from this analysis and compared to parameterized estimates of gravity wave forcing, providing a measure of the importance of the storm source to global gravity wave forcing. The results suggest that a single large convective storm system like the simulated squall line could provide a significant fraction of the zonal mean gravity wave forcing at some levels, particularly in the mesosphere. The vertical distributions of mean flow acceleration and energy dissipation do not much resemble the parameterized profiles in form because of the peculiarities of the spectral properties of the waves from the storm source. The ray tracing model developed herein provides a tool for examining the role of convectively generated waves in middle atmosphere physics.

Alexander, M. Joan↗
NASA NTRS1983

A theoretical model study of observed correlations between whistler mode waves and energetic electron precipitation events in the magnetosphere

A recently extended test particle computer model of the gyroresonance wave-particle interaction in the magnetosphere is applied to previously reported cases of observed correlations between whistler mode waves and ionospheric responses to particle precipitation. Three different ionospheric effects, namely, X-ray bursts, photoemissions, and D region perturbations, all correlated with VLF waves and believed to be caused by precipitated particles, are considered. The precipitation flux level, the pulse shape, and the associated time delays are computed for the parameters relevant to each case and are compared with values deduced from the data. The results demonstrate that the existing theoretical model can be useful for interpreting experimental results of this kind. Furthermore, the model results and observations, used together, provide a basis for additional diagnostics of the various parameters of the cold and energetic particle distributions in the magnetosphere. For example, when applied to the observed photoemission case (Helliwell et al., 1980) the model results imply that the trapped energetic particle distribution function at the time could be modeled as proportional to E exp -n/2 with n about 3.5 to 6, where E is the particle energy.

Chang, H. C.↗
NASA NTRS2002

A Self-Consistent Model of the Interacting Ring Current Ions and Electromagnetic Ion Cyclotron Waves, Initial Results: Waves and Precipitating Fluxes

Initial results from a newly developed model of the interacting ring current ions and ion cyclotron waves are presented. The model is based on the system of two kinetic equations: one equation describes the ring current ion dynamics, and another equation describes wave evolution. The system gives a self-consistent description of the ring current ions and ion cyclotron waves in a quasilinear approach. These equations for the ion phase space distribution function and for the wave power spectral density were solved on aglobal magnetospheric scale undernonsteady state conditions during the 2-5 May 1998 storm. The structure and dynamics of the ring current proton precipitating flux regions and the ion cyclotron wave-active zones during extreme geomagnetic disturbances on 4 May 1998 are presented and discussed in detail.

Khazanov, G. V.↗
NASA NTRS1986

Shock waves and shock tubes; Proceedings of the Fifteenth International Symposium, Berkeley, CA, July 28-August 2, 1985

A detailed survey is presented of shock tube experiments, theoretical developments, and applications being carried out worldwide. The discussions explore shock tube physics and the related chemical, physical and biological science and technology. Extensive attention is devoted to shock wave phenomena in dusty gases and other multiphase and heterogeneous systems, including chemically reactive mixtures. Consideration is given to techniques for measuring, visualizing and theoretically modeling flowfield, shock wave and rarefaction wave characteristics. Numerical modeling is explored in terms of the application of computational fluid dynamics techniques to describing flowfields in shock tubes. Shock interactions and propagation, in both solids, fluids, gases and mixed media are investigated, along with the behavior of shocks in condensed matter. Finally, chemical reactions that are initiated as the result of passage of a shock wave are discussed, together with methods of controlling the evolution of laminar separated flows at concave corners on advanced reentry vehicles.

Bershader, D.↗
NASA NTRS1991

A grid-independent approximate Riemann solver with applications to the Euler and Navier-Stokes equations

A new two-dimensional approximate Riemann solver has been developed that obtains fluxes on grid faces via wave decomposition. By utilizing information propagation in the velocity-difference directions rather than in the grid-normal directions, this flux function more appropriately interprets and hence more sharply resolves shock and shear waves when they lie oblique to the grid. The model uses five waves to describe the difference in states at a grid face. Two acoustic waves, one shear wave, and one entropy wave propagate in the direction defined by the local velocity difference vector, while the fifth wave is a shear wave that propagates at a right angle to the other four. Test cases presented include a shock reflecting off a wall, a pure shear wave, supersonic flow over an airfoil, and viscous separated airfoil flow. Results using the new model give significantly sharper shock and shear contours than a grid-aligned solver. Navier-Stokes computations over an aifoil show reduced pressure distortions in the separated region as a result of the grid-independent upwinding.

Rumsey, Christopher L.↗