Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamic 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

TEAMER: Supporting model output files for Environmental Compliance Framework for Floating Tidal Turbines, Cook Inlet, AK

Orbital Marine Power (Orbital) is seeking to deploy their floating tidal technology in US waters and has considered the possibility of deploying in temperate waters including the Pacific Northwest (PNW) and the Western Passage, Maine. It has become apparent that some of the most promising tidal sites in the US are located in high latitudes, within the State of Alaska. Within the state, the most likely sites for grid-scale tidal energy are in Cook Inlet, with the largest city (Anchorage) located on the shores, ready and able to absorb the electrical output of a commercial scale tidal development. Through the TEAMER project, Pacific Northwest National Laboratory (PNNL) developed a framework to select environmentally compliant sites partially based on hydrodynamic model output. This uploaded dataset summarizes the key hydrodynamic model output files discussed in the TEAMER project report. Please unzip the file to see actual data files and the ReadMe file. This project is part of the TEAMER RFTS 4 (request for technical support) program.

16 TIDAL AND WAVE POWER↗

Characterizing in-stream turbulent flow for tidal energy converter siting in Cook Inlet, Alaska

Cook Inlet in Alaska is the most promising location for tidal energy development in the U.S. due to its significant tidal range of approximately 10 meters and high volume flux. The inlet's unique geometry and flow characteristics make it the most energetic tidal stream in the nation, with GW-scale potential energy capacity. With the growing interest in tidal energy converter (TEC) deployment in this area, we implemented a regional-scale, 3D hydrodynamic modeling framework to predict tidal current and turbulence characteristics that can assist TEC designers and project managers. We validated the model results extensively using various datasets collected with bottom-mounted acoustic Doppler current profilers and velocimeters. The comparison between the model outputs and observational data highlighted the effectiveness of the 3D FVCOM model and the Mellor-Yamada Level 2.5 Turbulence Model in accurately assessing macro-scale kinetic energy, turbulence intensity, and the production and dissipation rates at a prospective TEC site. Using two months of model simulation data, we examined the channel cross-section for TEC deployment, focusing on undisturbed power density and macro-scale turbulent properties. Further, our findings indicate that understanding the turbulence characteristics and flow properties can enhance Stage I/II resource characterization by identifying optimal locations for TECs and their layouts within the channel. Furthermore, we demonstrated that TEC designers can utilize macro-scale turbulence data from 3D coastal models as boundary conditions for other turbulence models, allowing for a more detailed resolution of the turbulence structure at TEC siting locations. Ultimately, this work emphasizes the importance of estimating flow and turbulence conditions in energetic systems to understand turbulent sites better and improve resource characterization.

16 TIDAL AND WAVE POWER↗

Modelling time variable and total eclipses of the millisecond pulsar PSR 1744-24A

The 11 ms pulsar PSR 1744-24A in the globular clusters Terzan 5 is the first and only known example of a binary eclipsing pulsar showing time-variable eclipses with irregular and nonsymmetric eclipse patterns. Occasionally, the pulsar eclipse becomes total, and no radio emission is detected for an extended period of time. We present here numerical results of a hydrodynamical model aimed at explaining the remarkable eclipse properties of PSR 1744-24A. We show that the hydrodynamic time scale for the partial or total pulsar enshrouding is of order of the orbital period for a relatively large mass loss rate, in agreement with observations. The time-variable eclipse egress positions imply that the mass loss rate from the companion varies within one orbital period. Hydrodynamic modelling of the mass outflow of PSR 1744-24A is of crucial importance for the understanding of a class of millisecond pulsars whose radio emission might be occasionally or permanently quenched by mass outflows from their companions. The properties of PSR 1744-24A suggest a new way of studying millisecond pulsars involving radio, X-ray, gamma-ray, and optical observations.

Tavani, Marco↗

Evaluating the Impact of Land Use Change on Submerged Aquatic Vegetation Stressors in Mobile Bay

Alabama coastal systems have been subjected to increasing pressure from a variety of activities including urban and rural development, shoreline modifications, industrial activities, and dredging of shipping and navigation channels. The impacts on coastal ecosystems are often observed through the use of indicator species. One such indicator species for aquatic ecosystem health is submerged aquatic vegetation (SAV). Watershed and hydrodynamic modeling has been performed to evaluate the impact of land use change in Mobile and Baldwin counties on SAV stressors and controlling factors (temperature, salinity, and sediment) in Mobile Bay. Watershed modeling using the Loading Simulation Package in C++ (LSPC) was performed for all watersheds contiguous to Mobile Bay for land use scenarios in 1948, 1992, 2001, and 2030. Landsat-derived National Land Cover Data (NLCD) were used in the 1992 and 2001 simulations after having been reclassified to a common classification scheme. The Prescott Spatial Growth Model was used to project the 2030 land use scenario based on current trends. The LSPC model simulations provided output on changes in flow, temperature, and sediment for 22 discharge points into the Bay. Theses results were inputted in the Environmental Fluid Dynamics Computer Code (EFDC) hydrodynamic model to generate data on changes in temperature, salinity, and sediment on a grid with four vertical profiles throughout Mobile Bay. The changes in the aquatic ecosystem were used to perform an ecological analysis to evaluate the impact on SAV habitat suitability. This is the key product benefiting the Mobile Bay coastal environmental managers that integrates the influences of temperature, salinity, and sediment due to land use driven flow changes with the restoration potential of SAVs.

Al-Hamdan, Mohammad↗

Multi-wavelength Observations and Modelling of a Solar Flare

Aims: The aim of this work is to investigate the dynamic behavior of a C-class solar flare through the evolution of temperature, emission measure, energy loss and velocity. In particular, the variation of these properties with time are studied using multi-wavelength observations in combination with a recently developed 0-D hydrodynamic model. Methods: The temperature and emission measure evolution were studied using several instruments covering a wide range of temperatures - the Reuven Ramaty High Energy Solar Spectroscopic Imager (RHESSI, >5 MK), GOES-12 (5- 30 MK), the Transition Region and Coronal Explorer (TRACE 171 A, 1 MK) and the Coronal Diagnostic Spectrometer (CDS, 0.03-8 MK). The temperature and emission measure were analysed through the systematic cooling of flare plasma through the response functions of these instruments. These parameters were then investigated using the Enthalpy Based Thermal Evolution of Loops model (EBTEL). The Doppler shifts at both flare footpoints were analysed using five emission lines seen by CDS. Results: The flare began with clear evidence for pre-flare heating. Upflows of approx.90 km/s and low level emission, both observed in Fe XIX before the main impulsive phase were explained by pre-flare gentle chromospheric evaporation. During the main impulsive phase, the flare plasma was heated to a temperature of >13 MK in approximately 10 minutes. Explosive chromospheric evaporation was observed, driving upflows of approx.80 km/s in Fe XIX and simultaneous downflows of approx.20 km/s in He I and O v. At the peak of the Rare, conduction modelled by EBTEL was found to be the dominant loss mechanism, working efficiently to both lower the temperatures and drive gentle chromospheric evaporation. As the temperature fell below approx.8 MK, radiation became the dominant loss mechanism. During the final stages of the decay phase, downflowing plasma was observed at the footpoints in He I, O v and Mg x at velocities of up to approx.40 km/s, suggesting loop draining occurred. Conclusions. This is the first extensive study of the evolution of flare plasma using both spectroscopic and broad-band instruments in conjunction with a comprehensive hydrodynamic model. The flare began with pre-flare heating and then evolved following the predictions of the standard flare model. Detailed analysis of the plasma heating mechanisms was carried out and the heating function most consistent with observations was found to be Gaussian in shape. The simulations suggested that both direct heating and heating by a non-thermal beam played significant roles in this event.

Raftery, Claire L.↗

Development of the Packed Bed Reactor ISS Flight Experiment

Packed bed reactors are compact, require minimum power and maintenance to operate, and are highly reliable. These features make this technology a leading candidate as a potential unit operation in support of long duration human space exploration. On earth, this type of reactor accounts for approximately 80% of all the reactors used in the chemical process industry today. Development of this technology for space exploration is truly crosscutting with many other potential applications (e.g., in-situ chemical processing of planetary materials and transport of nutrients through soil). NASA is developing an ISS experiment to address this technology with particular focus on water reclamation and air revitalization. Earlier research and development efforts funded by NASA have resulted in two hydrodynamic models which require validation with appropriate instrumentation in an extended microgravity environment. The first model developed by Motil et al., (2003) is based on a modified Ergun equation. This model was demonstrated at moderate gas and liquid flow rates, but extension to the lower flow rates expected in many advanced life support systems must be validated. The other model, developed by Guo et al., (2004) is based on Darcy s (1856) law for two-phase flow. This model has been validated for a narrow range of flow parameters indirectly (without full instrumentation) and included test points where the flow was not fully developed. The flight experiment presented will be designed with removable test sections to test the hydrodynamic models. The experiment will provide flexibility to test additional beds with different types of packing in the future. One initial test bed is based on the VRA (Volatile Removal Assembly), a packed bed reactor currently on ISS whose behavior in micro-gravity is not fully understood. Improving the performance of this system through an accurate model will increase our ability to purify water in the space environment.

Patton, Martin O.↗

Advances in constraining intrinsic alignment models with hydrodynamic simulations

We use galaxies from the illustristng, massiveblack-ii, and illustris-1 hydrodynamic simulations to investigate the behaviour of large scale galaxy intrinsic alignments. Our analysis spans four redshift slices over the approximate range of contemporary lensing surveys z = 0–1. We construct comparable weighted samples from the three simulations, which we then analyse using an alignment model that includes both linear and quadratic alignment contributions. Our data vector includes galaxy–galaxy, galaxy–shape, and shape–shape projected correlations, with the joint covariance matrix estimated analytically. In all of the simulations, we report non-zero IAs at the level of several σ. For a fixed lower mass threshold, we find a relatively strong redshift dependence in all three simulations, with the linear IA amplitude increasing by a factor of ~2 between redshifts z = 0 and z = 1. We report no significant evidence for non-zero values of the tidal torquing amplitude, A 2 , in TNG, above statistical uncertainties, although MBII favours a moderately negative A 2 ~ –2. Examining the properties of the TATT model as a function of colour, luminosity and galaxy type (satellite or central), our findings are consistent with the most recent measurements on real data. We also outline a novel method for constraining the TATT model parameters directly from the pixelized tidal field, alongside a proof-of-concept exercise using TNG. This technique is shown to be promising, although comparison with previous results obtained via other methods is non-trivial.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling the hydrodynamic impact on the tool influence function during hemispherical subaperture optical polishing

To fabricate high-precision and accurate optics relative to the optical design surface, a high level of deterministic control of material removal (i.e., the tool influence function, TIF) during subaperture tool polishing is required. In this study, a detailed analysis of the pressure distribution, which is a key component of the TIF, has been performed using finite element analysis to couple together solid mechanics and fluid dynamics. Modeling experimental parameters of recently published work reveals that, when considering tool deformation, which in turn influences the fluid film thickness between the tool and workpiece, the effective pressure profile has a flat-top distribution. This flat-top pressure profile differs from the parabolic pressure distributions predicted by Hertzian mechanics. Furthermore, the shear contribution is shown here to be a key contributor to material removal, inducing the removal at the periphery of the contact edge and even outside the generally accepted contact area. Finally, the simulated fluid velocities provide evidence of mixed-mode contact polishing, supporting recent experimental findings that also suggest that onset of hydroplaning contributions lead to material removal drop-off.

36 MATERIALS SCIENCE↗

Review of TEAMER Awards for WEC-Sim Support: Preprint

Testing Expertise and Access for Marine Energy Research (TEAMER) is a U.S. Department of Energy Water Power Technologies Office sponsored program, overseen by the Pacific Ocean Energy Trust, which aims to advance the state of marine energy technologies. The program connects technology developers with experts at U.S. facilities, including numerical modeling and analysis facilities, like WEC-Sim. The WEC-Sim facility is supported by the WEC-Sim development team at Sandia National Laboratories and the National Renewable Energy Laboratory. WEC-Sim (Wave Energy Converter SIMulator) is an open-source software for simulating wave energy converters. WEC-Sim can model the multi-body dynamics of devices comprised of bodies, joints, power take-off systems, and mooring systems. Since TEAMER's first round of support in 2020, there have been eighteen TEAMER awards focused on numerical model development in WEC-Sim. TEAMER awards for WEC-Sim support have modeled a wide range of wave energy converter archetypes, including point absorbers, attenuators, oscillating water columns, and many other novel architectures. A wide variety of studies have been conducted, leading to important insights for TEAMER partners and software improvements for WECSim. This article highlights several successful WEC-Sim TEAMER awards. The awards described herein include TEAMER recipients Ocean Motion Technologies, AquaHarmonics, iProTech, East Carolina University, Virginia Tech, Maiden Wave Energy, and the University of Massachusetts Dartmouth. The awards of these seven partners contain a wide range of investigations and cover the creation of baseline hydrodynamic models, PTO modeling, geometry optimization in both boundary element methods and WECSim, and model tuning and validation.

industry support↗

On the X-ray spectrum of Kepler's supernova remnant

We have devised a method to do nonequilibrium ionization calculations on the results of two-dimensional hydrodynamical simulations, based on the algorithm of Hughes & Helfand (1985). We have calculated the ionization structure and X-ray emission for a two-dimensional numerical hydrodynamical simulation for the remnant of Kepler's supernova (SN); the hydrodynamical model was presented in a previous paper. In this model, the progenitor of Kepler's SN is assumed to have been a massive runaway star ejected from the Galactic plane. In its red supergiant stage, its dense stellar wind was distorted and compressed into a bow shock by the ram pressure of the tenuous interstellar medium. The subsequent interaction of the supernova ejecta with this asymmetric circumstellar matter produced a strongly asymmetric supernova remnant (SNR). In this paper, we present calculated X-ray spectra for this hydrodynamical model. A comparison with observations implies a moderate overabundance of Fe in Kepler's SNR (only 50% larger than its cosmic value), in contrast to a large (6 to 15) Fe overabundance derived previously. However, we confirm earlier conclusions that Si and S abundances are 2 to 3 times solar. These modest enhancements of Si, S, and Fe may be attributed either to heavy-element enriched SN ejecta or to the initial chemical abundances of the SN progenitor, which originated in the metal-rich inner Galaxy. The comparison of our models with the observed spectra confirm theoretical predictions that moderate electron heating occurs at strong collisionless shock fronts, with the implied electron/mean temperature ratio of approximately 0.5.

Borkowski, Kazimierz J.↗

San Juan Islands Tidal Energy Characterization

The San Juan Islands are an archipelago with multiple tidal channels that produce strong tidal currents that could be harnessed for electricity generation. At present, there is limited electrical generation on the islands, with power primarily provided by a subsea cable connection with the mainland. As part of the C-MIST program, NOAA's CO-OPS program collected ADCP data at 46 stations in the San Juan Islands and adjacent waterways, from April to August of 2017. PNNL has recently developed a 3-D tidal hydrodynamic model for tidal energy resource characterization and assessment in the Salish Sea (Yang et al. 2021). Several tidal channels in the San Juan Islands were identified as top hotspots for potential tidal energy extraction. However, due to the complex geometry and inter-connected waterways, tidal currents around San Juan Islands also exhibit strong spatial and temporal asymmetry. Therefore, it is important to analyze the NOAA ADCP data to define metrics relevant to tidal energy, and further refine PNNL's tidal hydrodynamic model to accurately simulate currents in small tidal channels and account for sharp bathymetry gradients. The outcome of this technical assistance will allow developers to identify the promising opportunities for tidal energy in the San Juan Islands.

16 TIDAL AND WAVE POWER↗

Vortex-sheet modeling of hydrodynamic instabilities produced by an oblique shock interacting with a perturbed interface in the HED regime

In this work, we consider hydrodynamic instabilities produced by the interaction of an oblique shock with a perturbed material interface under high-energy-density (HED) conditions. During this interaction, a baroclinic torque is generated along the interface due to the misalignment between the density and pressure gradients, thus leading to perturbation growth. Our objective is to understand the competition between the impulsive acceleration due to the normal component of the shock velocity, which drives the Richtmyer–Meshkov instability, and the shear flow across the interface due to the tangential component of the shock velocity, which drives the Kelvin–Helmholtz instability, as well as its relation to perturbation growth. Since the vorticity resulting from the shock-interface interaction is confined to the interface, we describe the perturbation growth using a two-dimensional vortex-sheet model. We demonstrate the ability of the vortex-sheet model to reproduce roll-up dynamics for non-zero Atwood numbers by comparing to past laser-driven HED experiments. We determine the dependence of the interface dynamics on the tilt angle and propose a time scaling for the perturbation growth at early time. Eventually, this scaling will serve as a platform for the design of future experiments. This study is the first attempt to incorporate into a vortex-sheet model the time-dependent interface decompression and the deceleration (as well as the corresponding Rayleigh–Taylor instability) arising from laser turn-off.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Tidal energy resource characterization measurements at Cook Inlet’s East Foreland: Velocity and turbulence

To characterize tidal current and turbulence at a top tidal energy site off the East Foreland in Cook Inlet, Alaska, United States, three moorings were deployed for two months between July and August 2021, and a transect survey was conducted over the course of two tidal cycles at the end of the deployment period. Measurements of velocity and turbulence were then analyzed to better understand the site's hydrodynamics and power potential. Analysis reveals that swift, north-flowing flood currents peak at 4~m/s, while south-flowing ebb currents reach just over 3~m/s. Turbulence intensity ranges from 23\% at the seafloor to 8\% near the surface, and the presence of the foreland creates more intense turbulence near-shore during ebb tide than flood. Power availability at the site could be as high as 720~MW, or 13~kW/m$^2$, though the energy available to a marine energy device will be smaller than this estimate because of water-to-wire efficiency and wake losses. The results from this measurement campaign will inform the validation of a high-resolution tidal hydrodynamic model, as well as early tidal energy projects that are beginning to move beyond the prototyping and demonstration stages to full-scale deployments.

McVey, James R.↗

Spectral Evidence for Heating at Large Column Mass in Umbral Solar Flare Kernels. I. IRIS Near-UV Spectra of the X1 Solar Flare of 2014 October 25

The GOES X1 flare SOL2014-10-25T17:08:00 was a three-ribbon solar flare observed with the Interface Region Imaging Spectrograph (IRIS)in the near-UV(NUV)and far-UV. One of the flare ribbons crossed a sunspot umbra, producing a dramatic, ∼1000% increase in the NUV continuum radiation. We comprehensively analyze the UV spectral data of the umbral flare brightenings, which provide new challenges for radiative−hydrodynamic modeling of the chromospheric velocity field and the white-light continuum radiation. The emission line profiles in the umbral flare brightenings exhibit redshifts and profile asymmetries, but these are significantly smaller than in another, well-studied X-class solar flare. We present a ratio of the NUV continuum intensity to the FeIIλ2814.45intensity. This continuum-to-line ratio is a new spectral diagnostic of significant heating at high column mass (log m/[gcm−2]>−2) during solar flares because the continuum and emission line radiation originate from relatively similar temperatures but moderately different optical depths. The full spectral readout of these IRIS data also allow for a comprehensive survey of the flaring NUV landscape: in addition to many lines of Fe II and Cr II, we identify a new solar flare emission line, HeIλ2829.91(as previously identified in laboratory and early-type stellar spectra). The Fermi/GBM hard X-ray data provide inputs to radiative−hydrodynamic models (which will be presented in Paper II) in order to better understand the large continuum-to-line ratios, the origin of the white-light continuum radiation, and the role of electron beam heating in the low atmosphere.

Adam F Kowalski↗

Propagation of positive discharges in an air bubble having an embedded water droplet

Low temperature plasmas generated inside gas bubbles immersed in water is an effective method of rapidly transferring plasma generated reactive species to the water for applications in biomedicine, agriculture and environment. Reactive species are generally produced in the gas phase plasma and then solvate into the liquid. The large surface-to-volume ratio (SVR) of the bubble accelerates this process. In generating bubbles in water, aerosols and droplets are also contained within the bubble. These droplets also have a large SVR and so can be rapidly plasma activated. However, the presence of the droplets can also impact the propagation of the plasma in the bubble. In this paper, results are discussed from computational and experimental investigations of the formation and evolution of discharges in an air bubble immersed in water with an embedded water droplet. The computations were performed with a two-dimensional plasma hydrodynamics model. Experiments were performed with a quasi-2D bubble apparatus. In bubbles having a droplet, a plasma filament typically bridges from the powered electrode to the droplet, and then from the droplet to the bubble surface. A surface-hugging streamer also occurs on the inner bubble surface and on the surface of the droplet. Both surface streamers result in part from surface charge accumulation and can dominate the formation of reactive species that transport into the droplet. Increasing droplet conductivity suppresses propagation of the surface discharge and leads to a lower density of aqueous reactive species. Increasing conductivity of the surrounding water does not change the overall structure of the discharge but does slightly elevate the discharge intensity. Finally, the size and shape of the embedded droplet can significantly affect the formation and propagation of the streamer.

2D fluid modelling↗

Shock physics in warm dense matter: A quantum hydrodynamics perspective

Abstract Warm dense matter (WDM), an exotic, highly compressed state of matter between solid and plasma phases, is of high current interest, in particular for astrophysics and inertial confinement fusion. For the latter, in particular the propagation of compression shocks is crucial. The main unknown in the shock propagation in WDM is the behaviour of the electrons since they are governed by correlations, quantum and spin effects that need to be accounted for simultaneously. Here we describe the shock dynamics of the warm dense electron gas using a quantum hydrodynamic model. From the numerical hydrodynamic simulations, we observe that the quantum Bohm pressure induces shear force that weakens the formation and strength of the shock. In addition, the Bohm pressure induces an electron density response that takes the form of oscillations. This is confirmed by the theoretical analysis of the early stage of the shock formation. Our theoretical and numerical analyses allow us to identify characteristic dimensionless shock propagation parameters at which the effect of the Bohm force is important.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

TRITON: A Multi-GPU open source 2D hydrodynamic flood model

A new open source multi-GPU 2D flood model called TRITON is presented in this work. The model solves the 2D shallow water equations with source terms using a time-explicit first order upwind scheme based on an Augmented Roe's solver that incorporates a careful estimation of bed strengths and a local implicit formulation of friction terms. Here, the scheme is demonstrated to be first order accurate, robust and able to solve for flows under various conditions. TRITON is implemented such that the model effectively utilizes heterogeneous architectures, from single to multiple CPUs and GPUs. Different test cases are shown to illustrate the capabilities and performance of the model, showing promising runtimes for large spatial and temporal scales when leveraging the computer power of GPUs. Under this hardware configuration, communication and input/output subroutines may impact the scalability. The code is developed under an open source license and can be freely downloaded in https://code.ornl.gov/hydro/triton.

2D flood model↗

Evaluating precipitation, streamflow, and inundation forecasting skills during extreme weather events: A case study for an urban watershed

Integrated forecasting systems for precipitation, streamflow, and floodplain inundation are of critical importance in mitigating the impacts of destructive floods caused by extreme weather events. However, the skills of streamflow and floodplain inundation forecasts derived from various Quantitative Precipitation Forecasts (QPF) require a greater level of understanding. In this paper, a set of QPF developed by the National Weather Service (NWS) were used to drive a flood modeling system obtained utilizing offline coupling of a physics-based distributed hydrological model, the Distributed Hydrology Soil and Vegetation Model (DHSVM), and a hydrodynamic model, the Two-dimensional Runoff Inundation Toolkit for Operational Needs (TRITON). This flood modeling system was used to produce forecasts of streamflow and floodplain inundation maps during three major flood events in the Brays Bayou Watershed (Houston, Texas, USA) for a range of QPF durations (6–72 h). Then, to investigate the effects of increasing QPF durations on the forecasts, the forecasting skills of precipitation, streamflow, and floodplain inundation were quantified. The results show that: 1) QPF skills for more intense and sustained events such as hurricanes and tropical storms are higher than for shorter, less intense events; 2) while QPF and streamflow forecasting skills decrease as QPF durations increase, inundation forecasts under longer QPF durations (24 or 72 h) show higher skills; 3) extending the maximum QPF duration in operational hydrologic modeling from 24 h (under normal circumstances) to 72 h (for extreme events) may increase the skills of long lead time forecasts for large-scale events like Hurricane Harvey.

54 ENVIRONMENTAL SCIENCES↗