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.
We present three-dimensional smoothed particle hydrodynamics models for the evolution of magnetically supported clouds by ambipolar diffusion.
Using the 2017 Hurricane Harvey flood event as a test case, this study set up a series of sensitivity analyses to highlight three challenges associated with large-scale flood inundation modeling, including (a) model parameterization, (b) errors in digital elevation models, and (c) effects of reservoir retention. Driven by radar-based hourly rainfall data, a series of hydrologic-hydraulic models including the VIC hydrologic model, RAPID routing model, and Flood2D-GPU hydrodynamic model are set up over Harris County, Texas, to simulate flood inundation and hazards. The results demonstrate the importance of hydrologic parameters in improving flood modeling. For a large flood event such as Hurricane Harvey, the effect of the initial water depths is insignificant. The Manning's n values may increase the peak water depth by ~1%, the flood extents by 65km 2 , and the high danger zone by ~6%. On the contrary, the bathymetry correction factors may reduce the flood extent by ~1.4% and the high-danger zone by ~4%. Reducing the reservoir storage capacity to 1% may increase the flood extent by ~4% and the high-danger zone by ~17%. This study may provide supporting information to guide and prioritize the development of future high-performance computing hydrodynamic large-scale flood simulations.
During the last 10 months of the grant, both laminar and anomalous plasma processes occurring during the refilling of the outer plasmasphere after magnetic storms are investigated. Theoretical investigations were based on two types of models: (1) two-stream hydrodynamic model in which plasma flows from the conjugate ionospheres are treated as separate fluids and the ion temperature anisotropies are treated self-consistently; and (2) large-scale particle-in-cell code.
Here, we report on methods and results for a model-based flood hazard assessment we have conducted for the Hudson River from New York City to Troy/Albany at the head of tide. Our recent work showed that neglecting freshwater flows leads to underestimation of peak water levels at up-river sites and neglecting stratification (typical with two-dimensional modeling) leads to underestimation all along the Hudson. As a result, we use a three-dimensional hydrodynamic model and merge streamflows and storm tides from tropical and extratropical cyclones (TCs, ETCs), as well as wet extratropical cyclone (WETC) floods (e.g. freshets, rain-on-snow events). We validate the modeled flood levels and quantify error with comparisons to 76 historical events. A Bayesian statistical method is developed for tropical cyclone streamflows using historical data and consisting in the evaluation of (1) the peak discharge and its pdf as a function of TC characteristics, and (2) the temporal trend of the hydrograph as a function of temporal evolution of the cyclone track, its intensity and the response characteristics of the specific basin. A k-nearest-neighbors method is employed to determine the hydrograph shape. Out of sample validation tests demonstrate the effectiveness of the method. Thus, the combined effects of storm surge and runoff produced by tropical cyclones hitting the New York area can be included in flood hazard assessment. Results for the upper Hudson (Albany) suggest a dominance of WETCs, for the lower Hudson (at New York Harbor) a case where ETCs are dominant for shorter return periods and TCs are more important for longer return periods (over 150 years), and for the middle-Hudson (Poughkeepsie) a mix of all three flood events types is important. However, a possible low-bias for TC flood levels is inferred from a lower importance in the assessment results, versus historical event top-20 lists, and this will be further evaluated as these preliminary methods and results are finalized. Future funded work will quantify the influences of sea level rise and flood adaptation plans (e.g. surge barriers). It would also be valuable to examine how streamflows from tropical cyclones and wet cool-season storms will change, as this factor will dominate at upriver locations.
Altimeter measurements are corrected for several geophysical parameters in order to access ocean signals of interest, like mesoscale or sub-mesoscale variability. The ocean tide is one of the most critical corrections due to the amplitude of the tidal elevations and to the aliasing phenomena of high-frequency signals into the lower-frequency band, but the internal-tide signatures at the ocean surface are not yet corrected globally. Internal tides can have a signature of several centimeters at the surface with wavelengths of about 50–250 km for the first mode and even smaller scales for higher-order modes. The goals of the upcoming Surface Water Ocean Topography (SWOT) mission and other high-resolution ocean measurements make the correction of these small-scale signals a challenge, as the correction of all tidal variability becomes mandatory to access accurate measurements of other oceanic signals. In this context, several scientific teams are working on the development of new internal-tide models, taking advantage of the very long altimeter time series now available, which represent an unprecedented and valuable global ocean database. The internal-tide models presented here focus on the coherent internal-tide signal and they are of three types: empirical models based upon analysis of existing altimeter missions, an assimilative model and a three-dimensional hydrodynamic model. A detailed comparison and validation of these internal-tide models is proposed using existing satellite altimeter databases. The analysis focuses on the four main tidal constituents: M2, K1, O1 and S2. The validation process is based on a statistical analysis of multi-mission altimetry including Jason-2 and Cryosphere Satellite-2 data. The results show a significant altimeter variance reduction when using internal-tide corrections in all ocean regions where internal tides are generating or propagating. A complementary spectral analysis also gives some estimation of the performance of each model as a function of wavelength and some insight into the residual non-stationary part of internal tides in the different regions of interest. This work led to the implementation of a new internal-tide correction (ZARON'one) in the next geophysical data records version-F (GDR-F) standards.
Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.
Abstract Interest in wave energy converters to provide autonomous power to various ocean-bound systems, such as autonomous underwater vehicles, sensor systems, and even aquaculture farms, has grown in recent years. The Monterey Bay Aquarium Research Institute has developed and deployed a small two-body point absorber wave energy device suitable to such needs. This paper provides a description of the system to support future open-source access to the device and further the general development of similar wave energy systems. Additionally, to support future control design and system modification efforts, a set of hydrodynamic models are presented and cross-compared. To test the viability of using a linear frequency-domain admittance model for controller tuning, the linear model is compared against four WEC-Sim models of increasing complexity. The linear frequency-domain model is found to be generally adequate for capturing system dynamics, as the model agreement is good and the degree of nonlinearity introduced in the WEC-Sim models is generally less than 2.5%.
This dataset includes modeled tidal current velocities, direction and depth at two locations in East and North Forelands (60.716, -151.434 and 61.024, -151.157) near Nikiski and Tyonek, respectively, in Cook Inlet, Alaska. Data from two grid cells were provided by the Pacific Northwest National Laboratory based on a tidal hydrodynamic model that characterized the tidal stream resources in Cook Inlet for a period from May 1 to September 1, 2005 (Wang and Yang 2020). The model grid size had a horizontal spatial resolution of 100 m at East Forelands and 200 m at Tyonek; mean sea level (MSL) depth was 47.9 m and 23.7 m at each respective site, and there were 10 depth bins that ranged in size with the tide from 4.3-5.2 m and 1.9-2.8 m, respectively (Wang and Yang 2020).
Evaluating extreme precipitation events (EPEs) is essential for building climate-resilient water management strategies, but it remains a major challenge in ungauged basins. Using the 26,070 km 2 Wadi al-Rummah basin in central Saudi Arabia as a case study, we developed an alternative, reliable, cost-effective satellite-based framework that combines empirically derived EPE thresholds, imagery-calibrated 2D hydrodynamic modeling, GRACE water-storage diagnostics, and bias-corrected CMIP6 projections to assess flood hazards and recharge potential under current and future climate scenarios in ungauged basins. The integrated approach and the resulting findings followed four key steps: (1) Identified a 22.5 mm EPE threshold, the 80th percentile of 3-day GPM/IMERG rainfall (2000–2024), aligned with flood-triggering events (Nov 2018: 23–28 mm; Apr, 2023: 42 mm); (2) Developed and calibrated a RiverFlow2D model using Sentinel-2 and PlanetScope imagery for the November 2018 flood, accurately reproducing flood depth and extent (RMSE ≤0.31 m; fuzzy-Dice ≥0.91), and estimating runoff (41 %), infiltration (25 %), and evaporation (34 %); (3) Independently validated the model with the April 2023 event (RMSE ≤0.35 m; fuzzy-Dice ≥0.86); (4) Conducted climate projections (2025–2100) from five bias-corrected NEX-GDDP CMIP6 models that revealed a 34 % increase in EPE intensity under SSP2-4.5 and 48 % under SSP5-8.5 scenarios, relative to 20th-century baselines. Our findings indicate that while intensifying extremes in the 21st century increase flood risk, the results highlight the potential for episodic recharge if effective retention strategies are employed, and offer a transferable model for climate-informed planning in data-scarce arid regions.
Bimetallic tubes are used when a component requires more than one performance requirement, for instance, strength/creep resistance and oxidation/corrosion resistance. Shear-assisted processing and extrusion (ShAPE) can be used to fabricate extrudates that are comparable or superior in performance relative to conventionally manufactured extrudates. For the first time, ShAPE has been successfully employed in producing bimetallic Al tubing consisting of 6061 and 7075 Al alloys. Both light and electron microscopy techniques were used to investigate the integrity of the tubes, especially the interface between the core and cladding of the bimetallic tubes. Void-free bimetallic tubes were produced using ShAPE. Quantification of tube integrity was carried out with tensile testing under as-extruded and post-aged conditions. All the bimetallic tubes in the as-extruded samples exhibited uniform elongation above 5% with good tensile strength. Key insights such as material flow during bimetallic tube extrusion were obtained from the characterization of remnant billet and simulation results from a smoothed particle hydrodynamics model.
A hydrodynamic model was developed of a spherically symmetric dusty gas flow in a cometary atmosphere assuming a single fluid, inviscid, perfect gas. The hydrodynamics for gas and dust, which involves the gas drag force (momentum transfer), heat exchange between gas and dust, photodissociation energy for H2O gas, and radiative heating and cooling terms for dust particles, are solved using the Gear method for stiff, coupled differential equations. Calculations were done with a dust size distribution for radii alpha = 0.01 micron to 10 cm with densities variable with the size. A nucleus size of 4.0 km radius with a density of 0.5 g/cu cm and a total dust-to-gas mass ratio chi = 1 was adopted. There are indications from in situ observations that dust particle fragment into smaller ones. Fragmentation of dust particles was incorporated into the model. This is done by adding source and sink terms in the continuity equations for the dust. Lifetimes for the decay of dust particles were assumed as a function of particle size. It is also assumed that dust particles always fragment only into the next smaller size.
Abstract Charging of particles having diameters of tens of microns has been extensively studied at atmospheric pressure in the context of, for example, electrostatic precipitators where the focus was on unipolar charging. The ambipolar charging of particles in atmospheric pressure plasmas, and of droplets in particular, has received less attention. The plasma activation of droplets is of interest for water purification, fertilizer production and materials synthesis, all of which depend on the transport of the droplets through the plasma, which in turn depends on their charging. In this paper, we report on the transport dynamics of water droplets, tens of microns in diameter, carried by the gas flow through an atmospheric pressure radiofrequency glow discharge sustained in helium. The droplets pass through the plasma with minimal evaporation and without reaching the Rayleigh limit. The droplet trajectory in the presence and absence of the plasma provides insights on the forces acting on the droplet. The measurements were analyzed using results from a three-dimensional fluid model and a two-dimensional plasma hydrodynamics model. We found that the transport dynamics as the droplet enters and leaves the plasma are due to differential charging of the droplet in the plasma gradients of the bounding sheaths to the plasma.
This project focused on developing an automated workflow to evaluate and optimize the iProTech Pitching Inertial Pump (PIP) wave energy converter (WEC) using open-source Python packages and the MATLAB/Simulink tool, WEC-Sim. The process involved parameterizing key design variables, running time-domain simulations, and performing sensitivity analyses to determine their impact on power output. The workflow, designed for the PIP device, is generalized and can be extended to optimize other WECs that can be simulated in WEC-Sim. This work establishes a foundation for future time-domain-based WEC design optimizations. Included in this submission are all figures from the final report and the model inputs required to generate them. This includes Python scripts with inputs that produce the meshes, boundary element method (BEM) models, hydrodynamic coefficients, and the WEC-Sim models used for time-domain analyses. Although data for every single run is not included to save space, all of it can be reproduced using the provided models. Detailed instructions for setting up the environment and running the codes are also included.
The physical processes occurring as a result of the transfer of energy and momentum from the primary solar flare energy release site in the corona to the underlying chromosphere and transition region during the course of the flare are investigated through a comparison of theoretical models and observational data. Static, dynamic and hydrodynamic models of the lower-temperature chromospheric flare are reviewed. The roles of thermal conduction, radiation, fast particles and mass motion in chromosphere-corona interactions are analyzed on the basis of Skylab UV, EUV and X-ray data, and empirical and synthetic models of the chromospheric and upper photospheric responses to flares are developed. The canonical model of chromospheric heating during flares as a result of primary energy release elsewhere is found to be justified in the chromosphere as a whole, although not entirely as the temperature minimum, and a simplified model of horizontal chromospheric flare structure based on results obtained is presented.
Shuttle Imaging Radar (SIR-B) imagery is Fourier filtered to remove the estimated system-transfer function, reduce speckle noise, and produce ocean scenes with a gray scale that is proportional to wave height. The SIR-B system response to speckled scenes of uniform surfaces yields an estimate of the stationary wavenumber response of the imaging radar, modeled by the 15 even terms of an eighth-order two-dimensional polynomial. Speckle can also be used to estimate the dynamic wavenumber response of the system due to surface motion during the aperture synthesis period, modeled with a single adaptive parameter describing an exponential correlation along track. A Fourier filter can then be devised to correct for the wavenumber response of the remote sensor and scene correlation, with subsequent subtraction of an estimate of the speckle noise component. A linearized velocity bunching model, combined with a surface tilt and hydrodynamic model, is incorporated in the Fourier filter to derive estimates of wave height from the radar intensities corresponding to individual picture elements.
The worldwide effort to design and commission floating offshore wind turbines (FOWT) is motivating the need for reliable numerical models that adequately represent their physical behavior under realistic sea states. However, properly representing the hydrodynamic quadratic damping for FOWT remains uncertain, because of its dependency on the choice of drag coefficients (dimensionless or not). It is hypothesized that the limited degree of freedom (DoF) drag coefficient formulation that uses only translational drag coefficients causes mischaracterization of the rotational DoF drag, leading to underestimation of FOWT global loads, such as tower base fore-aft shear. To address these hydrodynamic modeling uncertainties, different quadratic drag models implemented in the open-source mid-fidelity simulation tool, OpenFAST, were investigated and compared with the experimental data from the Offshore Code Comparison Collaboration, Continued, with Correlation (OC5) project. The tower base fore-aft shear and up-wave mooring line tension were compared under an irregular wave loading condition to demonstrate the effects of the different damping models. Two types of hydrodynamic quadratic drag formulations were considered: (1) member-based dimensionless drag coefficients applied only at the translational DoF (namely limited-DoF drag model) and (2) quadratic drag matrix model (in dimensional form). Based on the results, the former consistently underestimated the 95th percentile peak loads and spectral responses when compared to the OC5 experimental data. In contrast, the drag matrix models reduced errors in estimates of the tower base shear peak load by 7–10% compared to the limited-DoF drag model. The underestimation in the tower base fore-aft shear was thus inferred be related to mischaracterization of the rotational pitch drag and the heave motion/drag by the limited-DoF model.
Ultrafast modulation of semiconductor quantum well (QW) laser is of technological importance for information technology. Improvement by order(s) of magnitude in data transfer rate is possible as terahertz (THz) radiation is available for heating the laser at picosecond time scale. Optical gain modulation in the QW is achieved via temperature modulation of electron-hole plasma (EHP). Applications include free-space THz communication, optical switching, and pulse generation. The EHP in the semiconductor QW is described with a two-band model. Semiconductor Bloch equations with many-body effects are used to derive a hydrodynamical model for the active QW region. Because of ultrafast carrier-carrier scatterings in the order of 50 fs, EHP follows quasiequilibrium Fermi-Dirac distributions and THz field interacts incoherently with it. Carrier-longitudinal optical (LO) phonon scatterings and coherent laser-EHP interaction are treated microscopically in our physical model. A set of hydrodynamical equations for plasma density, temperature, and laser envelop amplitude are derived and Runge-Kutta method is adopted for numerical simulation. A typical 8 nm GaAs/Al(0.3)Ga(0.7) As single QW at 300 K is used. Additional information is contained in the original extended abstract.
This project explored how applied magnetic fields can improve inertial confinement fusion (ICF), specifically the indirect-drive approach that uses a hohlraum. This has been proposed for several decades as potentially beneficial, due to thermal insulation (reduced losses) from the imploded hotspot. We performed the most advanced radiation-magneto-hydrodynamic modeling to date of magnetized ICF designs in the ignition regime. We found that adding technologically feasible fields up to 60 – 70 Tesla could increase the fusion yield of current igniting designs for the National Ignition Facility (NIF) by up to 8x. Also, simulations show that in certain cases relatively small fields of 3 – 5 Tesla could double the yield, and be implemented at much lower cost. Early work on re-optimizing NIF designs with magnetic fields, namely by using a thicker ablator with more mass remaining, could increase the yield of a sub-ignition target by 18x and bring it into the ignition regime. A separate benefit of magnetization besides reduced thermal loss is reduced hydrodynamic instability. Modeling work under this project shows this could be significant, though early experiments at NIF and the Omega Laser proved inconclusive. We designed and proposed an improved NIF experiment on magnetized mix, based on a large-amplitude imposed perturbation. The project also supported basic physics research into magnetized laser-plasma interactions, namely cross-beam energy transfer, both with experiments at Omega and theory / modeling.