SEARCH · Engineering Papers
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.
The attached zip file includes a SolidWorks pack-and-go assembly of NREL's HERO WEC (hydraulic and electric reverse osmosis wave energy converter) V1.0. This model does not include all aspects of the design (i.e. RO (reverse osmosis) system, electrical enclosure, hose, cable) it only includes the WEC and PTO (power take-off) design.
Protection systems are a critical feature for the proper operation of a power grid. Protection systems should be able to quickly detect, identify, locate, and isolate the faults within a time frame that prevents damage to grid components and reduces customer outages. As more inverter-based photovoltaics (PV) are integrated into the grid, phasor domain based protection schemes become more inaccurate due to low levels of fault current from inverters or respond too slowly to accurately identify fault conditions. This points to the need to identify and develop novel protection approaches for distribution systems with high penetration of inverter based PV. The main objective of the project is to develop new time domain-based protection schemes to enable more accurate and faster detection, identification and location of faults in distribution systems and/or microgrids with high penetration of inverter-based PV. Our approach aims to exploit the traveling waves that are generated during a fault or other disturbances in the grid.
Current research focused on the examination of the effects of ozone heating and Newtonian cooling on the linear stability of transient planetary waves. A finite difference numerical model was developed in which the basic state wind, temperature, and ozone distributions were constructed using ground-based and satellite data. Preliminary results indicate that ozone dynamics interaction may play a more important role in the stability and maintenance of planetary waves than previously thought, particularly in summer when the mean solar zenith angle is smallest and thus ozone heating largest. Examination of the finite amplitude interactions among radiation, ozone, and dynamics was carried out. Self-consistent, coupled equations governing the weakly nonlinear interactions are currently being analyzed to provide a better understanding of wave-mean flow interactions and ozone transport in a continuously stratified model of the troposphere-stratosphere coupled system.
Marine energy resources such as river, tidal, and ocean currents can provide abundant, predictable, and clean power to both densely populated and remote communities globally. However, marine turbine technologies are still in an early stage of development. Recently, features have been added to the National Renewable Energy Laboratory's open-source wind turbine modeling tool OpenFAST to support the simulation of axial-flow marine turbines. These new features enable marine energy developers to predict the performance and safety of their turbines and allow researchers and students to advance the field through the exploration of novel designs. Developing and demonstrating the new capabilities of OpenFAST for marine turbines requires a baseline design that can be used as an example and test case. There are limited open-source turbine models suitable for this purpose. One existing baseline design is the Reference Model 1 (RM1), a fixed-bottom, dual-rotor marine turbine that was published in 2014 as part of the Department of Energy's Water Power Reference Model Project. To test and demonstrate OpenFAST's ability to model floating marine turbines, several modifications were made to the RM1. A floating platform and mooring system were designed to support a single RM1 rotor in a wide range of possible current and wave conditions. The system and model are intended to be a stable reference case and starting point for research efforts. The platform features all cylindrical members and familiar design elements from the offshore wind industry; the mooring system uses simple catenary chains. While not optimized for the lowest platform cost, the design is understandable and robust and can serve as a baseline for floating marine turbine technologies. This is a useful starting point for studies to explore new technologies or design methodologies. Additionally, the floating RM1 design enables cross-collaboration among industry, academia, and national lab researchers and serves as an educational platform for newcomers to marine energy. This short paper presents the specifications of the turbine, including its floating support structure and mooring system. The rationale behind the selection of certain design parameters is described. To enable adoption of this baseline case for research and educational efforts, an OpenFAST model of the full system has been made publicly available at https://github.com/OpenFAST/r-test/tree/main/glue-codes/openfast/MHK_RM1_Floating.
One of the strongest plasma wave signals observed during the Voyager 2 encounter with Neptune is a narrowband emission between 3.0 and 4.3 kHz that was detected over a period of roughly 2 hours around closest approach. The emission occurs below the electron cyclotron frequency and the low-frequency cutoff of the radio continuum radiation. Of the naturally occurring signals in the earth's auroral zone and in Jupiter's magnetosphere this emission most resembles trapped Z mode waves found near the left-hand cutoff frequency. Using this identification, a plasma density profile is obtained that is independent of the plasma temperature. These densities greatly exceed those measured by the plasma science instrument on Voyager but are lower than estimates based on other models of Neptunian plasma wave phenomenology. If this wave mode is not a natural emission, it might arise from an unusual interaction of the spacecraft with the cold, dense ambient plasma.
A statistical study of diffuse energetic ion events and their related waves upstream of the Earth's bow shock was performed using data from the Active Magnetospheric Particle Tracer Explorers/Ion Release Module (AMPTE/IRM) satellite over two 5-month periods in 1984 and 1985. The data set was used to test the assumption in the self-consistent model of the upstream wave and particle populations by Lee (1982) that the particle acceleration through hydromagnetic waves and the wave generation are directly coupled. The comparison between the observed wave power and the wave power predicted on the observed energetic particle energy density and solar wind parameters results in a high correlation coefficient of about 0.89. The intensity of diffuse ions falls off approximately exponentially with the distance upstream from the bow shock parallel to the magnetic field with e-folding distances which vary from approximately 3.3 R(sub E) to approximately 11.7 R(sub E) over the energy range from 10 keV/e to 67.3 keV/e for both protons and alpha particles. After normalizing the upstream particle densities to zero bow shock distance by using these exponential variations, a good correlation (0.7) of the density of the diffuse ions with the solar wind density was found. This supports the suggestion that the solar wind is the source of the diffuse ions. Furthermore, the spectral slope of the diffuse ions correlates well with the solar wind velocity component in the direction of the interplanetary magnetic field (0.68 and 0.66 for protons and alpha particles) which concurs with the notion that the solar wind plays an important role in the acceleration of the upstream particles.
We studied the shock-induced collapse of butane gas bubbles in the homogeneous explosive nitromethane (NM) to investigate the effects of hot spot formation on the detonation process. A butane bubble was injected into a sample of NM, and a shock wave from a flat plate impactor compressed the bubble, creating a localized hot spot. We measured shock and detonation wave speeds with optical velocimetry, and we used a high-speed camera to image the shock propagation and bubble collapse processes. A multiband optical fiber pyrometer measured the time-resolved thermal radiance, and we used the results and emissivity values extracted from spectral fits to estimate temperatures. We measured the characteristics of the shock-to-detonation transition in NM with and without a bubble. All experiments were performed at shock pressures near 8 GPa, where neat NM can detonate. A single bubble in this system was shown to sensitize NM, leading to a reduced run-to-detonation time. We used hydrodynamic modeling to predict shock wave propagation, the extent of chemical reaction, and subsequent temperature rise from the collapsing bubble. We used a temperature-dependent Arrhenius burn model for simulations, and it yielded much better results than reactive burn models that depend only on pressure and density.
The Goddard Institute for Space Studies (GISS) climate/middle atmosphere model has been used to study the impacts of increasing greenhouse gases, polar ozone depletion, volcanic eruptions, and solar cycle variability. We focus on the projection of the induced responses onto Northern Hemisphere winter surface climate. Changes in the model's surface climate take place largely through enhancement of existing variability patterns, with greenhouse gases, polar ozone depletion and volcanic eruptions primarily affecting the Arctic Oscillation (AO) pattern. Perturbations descend from the stratosphere to the surface in the model by altering the propagation of planetary waves coming up from the surface, in accord with observational evidence. Models lacking realistic stratospheric dynamics fail to capture these wave flux changes. The results support the conclusion that the stratosphere plays a crucial role in recent AO trends. We show that in our climate model, while ozone depletion has a significant effect, greenhouse gas forcing is the only one capable of causing the large, sustained increase in the AO observed over recent decades. This suggests that the AO trend, and a concurrent strengthening of the stratospheric vortex over the Arctic, are very likely anthropogenic in origin.
Steady-state gasdynamical studies, previously limited to tightly wound normal spiral galaxies, are extended to models of barred spirals with a 5% to 10% perturbing potential. The models show that a strong wave manifestation is an important constituent of the bar structure in many barred spirals and that a density-wave shock wave can form a bar structure as pronounced as the narrow bars often evident in optical photographs of barred spirals. The dark narrow dust lanes often observed along the leading edges of bar structures are identified as tracers of shocks, and it is found that strong shocks along a bar structure during even a small part of a galaxy's lifetime might easily deplete a large enough proportion of the gas to cause a lack of gas in the inner annuli encompassing the bar by the time of the present epoch. It is emphasized that even moderate-amplitude barlike perturbations in the disk can drive large noncircular gas motions, typically 50 to 150 km/s.
We analyse Starobinsky inflation in the presence of the Brout Englert Higgs (BEH) boson with a non-minimal coupling to the Ricci scalar, $R$. The latter induces a coupling of the massive scaleron associated with the $R^2$ term to the BEH boson and this leads to a radiative correction to the BEH mass that must be fine tuned to keep the scalar light. For the case of $R^2$ driven inflation this requires a high level of fine tuning of order 1 part in $10^{8}$; for the case of Higgs inflation it is very much greater. We consider a scale invariant extension of the $R^2$/Higgs model and find that for $R^2$ driven inflation but not for Higgs inflation the required fine tuning is significantly reduced to one part in $10^{3-4}$. We consider the vacuum stability of the fine tuned model and its reheating and dilaton abundance after inflation. We also discuss possible gravitational wave signals associated with the model and the constraint on the mass of scalar or fermion dark matter candidates if they are produced by the gravitational couplings of the scalaron.
This work deals with an instability arising from the shear-induced migration of particles in dense suspensions coupled with a dependence of viscosity on particle concentration. The analysis summarized here treats the inertialess (Re = O) linear stability of homogeneous simple shear flows for a Stokesian suspension model of the type proposed by Leighton and Acrivos (1987). Depending on the importance of shear-induced migration relative to concentration-driven diffusion, this model admits short-wave instability arising from wave-vector stretching by the base flow and evolving into particle-depleted shear bands. Moreover, this instability in the time-dependent problem corresponds to loss of ellipticity in the associated static problem (Re = O, Pe = O). While the isotropic version of the Leighton-Acrivos model is found to be stable with their experimentally determined parameters for simple shear, it is known that the stable model does not give a good quantitative description of particle clustering in the core of pipe flow (Nott and Brady 1994). This leads to the conjecture that an appropriate variant on the above model could explain such clustering as a two-phase bifurcation in the base flow.
We studied the collapse of individual helium gas bubbles in the homogeneous explosive nitromethane (NM) to investigate effects of hot-spot formation on the detonation process. A bubble was injected into a NM sample, and a shock wave from an explosive detonator compressed the bubble, creating a localized hot spot. We measured shock and detonation wave speeds with optical velocimetry, and we used a high-speed camera to image the shock propagation and the pre- and post-bubble collapse processes. An infrared camera image showed the residual radiance temperature distribution after the bubble collapse, and an optical fiber pyrometer measured the time-resolved thermal radiance. We measured the optical spectra of light emitted from detonating NM without a bubble and from a collapsing bubble in shocked, undetonated NM. We estimated temperatures of the detonation fronts and of the hot spots formed by bubble collapse. To study the incipient detonation process, we performed all bubble collapse experiments at pressures below the threshold for creating a sustained detonation. Where the bubble collapsed, we observed an opaque, thermally emissive region believed to be chemical reaction products. Chemical reactions in NM can be produced with lower shock pressures (~1 GPa) when a helium bubble is present than without a bubble (~10 GPa). We used hydrodynamic modeling to predict shock wave propagation, extent of chemical reaction, and subsequent temperature rise from the collapsing bubble. Simulations using a temperature-dependent Arrhenius burn model gave much better results than reactive burn models that depend only on pressure and density.
Not Available
The following information and metadata applies to both the Phase I (Hydrodynamics) and Phase II (Full System Power Take-Off) zip folders which contain testing data from the OSU (Oregon State University) O.H. Hinsdale Wave Research Laboratory, from both OSU and the University of Hawaii at Manoa (UH). See zip folders provided further below in the downloads section. For experimental data of the full system, including PTO, see Phase II dataset. There are two main directories in each Phases's zip folder: "OSU_data" and "UH_data". The "OSU_data" directory contains data collected from their DAQ (data acquisition system), which includes all wave gauge observations, as well as body motions derived from their Qualisys motion tracking system. The organization of the directory follows OSU's convention. Detailed information on the instrument setup can be found under "OSU_data/docs/setup/instm_locations". The experiments conducted are documented in the "OSU_data/docs/daq_logs", which provides the trial number to the corresponding data located under "OSU_data/data" in several formats (e.g., ".mat" and ".txt"). Inside the trial directory, data is provided for each of the instruments defined in "OSU_data/docs/setup/instm_locations". The "UH_data" directory contains data collected from their DAQ. The data is stored in a ".tdms" file format. There are free plug-ins for Microsoft Excel and MathWorks MATLAB to read the ".tdms" format. Below are a few links providing methods to read in the data, but a Google search should identify alternatives sources if these no longer exist (valid as of January 2024): Excel: http://www.ni.com/example/27944/en/ MATLAB: https://www.mathworks.com/matlabcentral/fileexchange/30023-tdms-reader The Excel plugin is recommend to get a quick overview of the data. The UH data is organized by directory name, in which the sub-directories for each experiment contains a directory whose name defines the wave height and period for the experimental data within. For example, a directory name "H02_T0275" corresponds to an experiment with wave height 0.1m and a period of 2.75s. For random wave data, the gamma value is also included in the directory name. For example, a directory name "H02_T0225_G18" corresponds to an experiment with a significant wave height of 0.2m, a peak period of 2.25s, and a gamma value of 1.8, with each spectra being a TMA spectrum. For the free decay experiments, the directory name is defined by the initial angular displacement. For example, a directory name "ang05_run01" corresponds to an experiment with an initial angular displacement of 5 degrees. There is a dataset in the UH data for each corresponding experiment defined in the OSU DAQ logs. The ".tdms" data is output from the DAQ at fixed intervals. Therefore, if multiple files are contained within the folder, the data will need to be stitched together. Within the UH dataset, there are two input channels from the OSU DAQ providing a random square wave signal for time synchronization ("ENV-WHT-0010") and a high/low signal ("ENV-WHT-0012") to identify when the wave maker is active (+5V). The UH data is logged as a collection of channel outputs. Channels not in use for the OSU testing (either Phase I or Phase II) are marked "nan" below. If the sensor is disconnected, it will record noise throughout the experiment. Below are the channel definitions in terms of what they measure: GPS Time = time CYL-POS-0001 = position between flap and fixed reference CYL-LCA-0001 = force between flap and hydraulic cylinder REC-LPT-0001 = nan REC-HPT-0001 = nan REC-HPT-0002 = nan REC-HPT-0003 = nan HHT-HPT-0001 = pressure at exhaust ("head" only) REC-FQC-0001 = nan REC-FQC-0002 = nan HHT-FQC-0001 = flow at exhaust ("head" only) ENV-WHT-0001 = nan ENV-WHT-0002 = nan ENV-WHT-0003 = nan ENV-WHT-0010 = random signal from OSU DAQ ENV-WHT-0012 = high/low signal from OSU DAQ Also included is a calibration curve to convert the string pot data to flap pi...
We present observations and modeling of spatial eigen-functions of resonating waves within fault zone waveguide, using data recorded on a dense seismic array across the San Jacinto Fault Zone (SJFZ) in southern California. The array consists of 5-Hz geophones that cross the SJFZ with ~10-30 m spacing at the Blackburn Saddle near the Hemet Stepover. Wavefield snapshots after the S wave arrival are consistent for more than 50 near-fault events, suggesting that this pattern is controlled by the fault zone structure rather than source properties. Data from example event with high signal to noise ratio shows three main frequency peaks at ~1.3 Hz, ~2.0 Hz, and ~2.8 Hz in the amplitude spectra of resonance waves averaged over stations near the fault. The data are modeled with analytical expressions for eigen-functions of resonance waves in a low-velocity layer (fault zone) between two quarter-spaces. Using a grid search-based method, we investigate the possible width of the waveguide, location within the array, and shear wave velocities of the media that fit well the resonance signal at ~1.3 Hz. The results indicate a ~300 m wide damaged fault zone layer with ~65% S wave velocity reduction compared to the host rock. The SW edge of the low-velocity zone is near the mapped fault surface trace, indicating that the damage zone is asymmetrically located at the regionally faster NE crustal block. The imaging resolution of the fault zone structure can be improved by modeling fault zone resonance modes and trapped waves together.
The frequency and intensity of summer heat waves in East Asia have increased sharply in recent decades, significantly impacting public health and the economy. The Arctic-Siberian Plain (ASP) teleconnection pattern has been identified as a key driver, with ASP warming amplifying atmospheric circulation patterns conducive to extreme temperatures. This study evaluates the ability of Coupled Model Inter-comparison Project phase 6 models to simulate the ASP pattern across interannual variability (IAV) and intra-seasonal variability (ISV) timescales using the Common Basis Function method. The multi-model mean shows statistically significant pattern correlations with ERA5 reanalysis, with correlation coefficients of 0.90 and 0.99 for IAV and ISV, respectively. While the ASP pattern is generally well captured, models exhibit substantial inter-model diversity in the intensity and position of anticyclonic anomalies over the ASP and East Asia. Models with ASP pattern variability similar to reanalysis better reproduce extreme East Asian temperatures, whereas those over- or underestimating ASP variability exhibit lower skill. These performance differences are related to differences in simulating key variables associated with the development of the ASP pattern. Our findings highlight the role of the ASP pattern in modulating extreme heat events, as models with improved ASP simulations align more closely with observed temperature extremes. Refining ASP representations in models could enhance seasonal heat wave predictions, improving climate adaptation strategies.
The interactions between stationary planetary waves 1 and 2 are investigated using a truncated quasi-geostrophic model in a midlatitude beta channel. The model includes eight interacting horizontal modes with three zonal wave numbers and two meridional scales in a quasi-geostrophic beta channel 40 deg wide; all wave-wave interactions in the model are a result of dissipation. The procedures for solving the model equations are examined. The theory of dissipation-induced interactions is discussed in terms of interactions between weak waves under severe truncation. The behavior of small amplitude waves 1 and 2 in basic states is discussed; it is observed that wave 1 is amplified in a wave 2 basic state, and a wave 2 disturbance in a wave 1 basic state is affected by the relative phases of the waves. The dynamics of these wave interactions are analyzed using potential enstrophy budgets of the waves. The relevance of the waves to the atmosphere is studied. It is noted that the wave behavior is a combination of the linear results with an amplification of wave 1 and a strong dependence on the relative phases of waves 1 and 2 in the stratosphere.
Electromagnetic Ion Cyclotron (EMIC) wave scattering has been proved to be responsible for the fast loss of both radiation belt (RB) electrons and ring current (RC) protons. However, its role in the concurrent dropout of these two co–located populations remains to be quantified. In this work, we study the effect of EMIC wave scattering on both populations during the 27 February 2014 storm by employing the global physics–based RAM–SCB model. Throughout this storm event, MeV RB electrons and 100s keV RC protons experienced simultaneous dropout following the occurrence of intense EMIC waves. By implementing data–driven initial and boundary conditions, we perform simulations for both populations through the interplay with EMIC waves and compare them against Van Allen Probes observations. Notably, the results indicate that by including EMIC wave scattering loss, especially by the He–band EMIC waves, the model aligns closely with data for both populations. Additionally, we investigate the simulated pitch angle distributions (PADs) for both populations. Including EMIC wave scattering in our model predicts a 90° peaked PAD for electrons with stronger losses at lower pitch angles, while protons exhibit an isotropic PAD with enhanced losses at pitch angles above 40°. Furthermore, our model predicts considerable precipitation of both particle populations, predominantly confined to the afternoon to midnight sector (12 hr < MLT < 24 hr) during the storm's main phase, corresponding closely with the presence of EMIC waves.