Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics simulations”

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 1,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
NASA NTRS1990

Monte Carlo simulation of nonadiabatic expansion in cometary atmospheres - Halley

Monte Carlo methods developed for the characterization of velocity-dependent collision processes and ballistic transports in planetary exospheres form the basis of the present computer simulation of icy comet atmospheres, which iteratively undertakes the simultaneous determination of velocity distribution for five neutral species (water, together with suprathermal OH, H2, O, and H) in a flow regime varying from the hydrodynamic to the ballistic. Experimental data from the neutral mass spectrometer carried by Giotto for its March, 1986 encounter with Halley are compared with a model atmosphere.

Hodges, R. Richard, Jr.↗
NASA NTRS1994

Synthesized spectra of turbulent clouds

A comparison has been made of the velocity structure of non-star-forming dense molecular clouds and the output of a 512(exp 3) hydrodynamic code for compressible gas displaying features of turbulence. A remarkable similarity exists between the interstellar line profiles and the velocity profiles found by integrating, with an optically thin assumption, through subsections of the simulation. This similarity is in terms of both the structure of individual spectra, showing skewness, multiple peaks, and wings, and also the variations of the spectra from point to point. It is possible to see the origin of the various velocity substructures and how they correlate with regions of high density and velocity divergence, or regions of high vorticity, by deprojecting individual spectra of the simulation. It is found that the interstellar medium spectra are best represented by the epoch of the simulation in which the energy spectrum is essentially Kolmogorov-like and initially developed shocks have been largely converted into vorticity structures. Further, it is found that the vorticity peaks ascribed to intermittency contribute to the non-Gaussian features and wings of the line profile, but at an intensity level that is weak compared to the contribution of the bulk of the flow.

Falgarone, E.↗
NASA NTRS1995

Unsteady numerical simulations of the stability and dynamics of flames

In this report we describe the research performed at the Naval Research Laboratory in support of the NASA Microgravity Science and Applications Program over the past three years (from Feb. 1992) with emphasis on the work performed since the last microgravity combustion workshop. The primary objective of our research is to develop an understanding of the differences in the structure, stability, dynamics and extinction of flames in earth gravity and in microgravity environments. Numerical simulations, in which the various physical and chemical processes can be independently controlled, can significantly advance our understanding of these differences. Therefore, our approach is to use detailed time-dependent, multi-dimensional, multispecies numerical models to perform carefully designed computational experiments. The basic issues we have addressed, a general description of the numerical approach, and a summary of the results are described in this report. More detailed discussions are available in the papers published which are referenced herein. Some of the basic issues we have addressed recently are (1) the relative importance of wall losses and gravity on the extinguishment of downward-propagating flames; (2) the role of hydrodynamic instabilities in the formation of cellular flames; (3) effects of gravity on burner-stabilized flames, and (4) effects of radiative losses and chemical-kinetics on flames near flammability limits. We have also expanded our efforts to include hydrocarbon flames in addition to hydrogen flames and to perform simulations in support of other on-going efforts in the microgravity combustion sciences program. Modeling hydrocarbon flames typically involves a larger number of species and a much larger number of reactions when compared to hydrogen. In addition, more complex radiation models may also be needed. In order to efficiently compute such complex flames recent developments in parallel computing have been utilized to develop a state-of-the-art parallel flame code. This is discussed below in some detail after a brief discussion of the numerical models.

Kailasanath, K.↗
NASA NTRS1993

Hydrodynamic models of the Cartwheel ring galaxy

A series of increasingly sophisticated models of the Cartwheel ring galaxy is studied in order to test the collisional model for the galaxy formation and examine the star formation processes in this unique environment, using new data acquired in the last decade. The simulations provided some possible answers to a number of questions about the Cartwheel. First, an explanation for the wide spacing between inner and outer rings is suggested by the simple epicyclic kinematics within the dark matter-dominated potential implied by H I rotation curve. These models and the kinematic model of Struck-Marcell and Lotan (1990) also predict that the outer ring should be relatively weak, while the second inner ring should be stronger, with a dense orbit-crossing region of significant width bounded by sharp, caustic edges. The collisional model is given support by the agreement between the observations and the morphological and kinematic properties of the numerical simulations presented.

Struck-Marcell, Curtis↗
NASA NTRS1995

Solar wind/local interstellar medium interaction including charge exchange with neural hydrogen

We present results from a hydrodynamic model of the interaction of the solar wind with the local interstellar medium (LISM), self-consistently taking into account the effects of charge exchange between the plasma component and the interstellar neutrals. The simulation is fully time dependent, and is carried out in two or three dimensions, depending on whether the helio-latitudinal dependence of the solar wind speed and number density (both giving rise to three dimensional effects) are included. As a first approximation it is assumed that the neutral component of the flow can be described by a single, isotropic fluid. Clearly, this is not the actual situation, since charge exchange with the supersonic solar wind plasma in the region of the nose results in a 'second' neutral fluid propagating in the opposite direction as that of the LISM neutrals.

Pauls, H. Louis↗
DOE OSTI2020

General Correlations for Laminar Flow Friction Loss and Heat Transfer in Plain Rectangular Plate-Fin Cores

In this paper, new generalized correlations for predicting the average fanning friction factor f and average Nusselt number Nu for laminar flow in plain plate-fin compact cores of rectangular cross section are presented. These are based on extended experimental data, as well as three-dimensional computational simulations, obtained for a broad range of fin density and geometrical attributes. The results indicate that while the fully developed forced convection scales only with the interfin channel cross-sectional ratio α (fin spacing by fin height), the entrance region hydrodynamic and thermal performance is additionally a function of the fin-core length L, flow Reynolds number Re, and fluid Prandtl number Pr. The developing flow and convection is further shown to scale as: (fRe)~(L/d h Re) 1/2 , and Nu ~(L/d h Re) 1/2 Pr 1/3 Φ(α), where f, Re, and Nu are all based on the hydraulic diameter dh of the interfin flow channel. Generalized correlations for both (fRe) and Nu are developed by the corresponding scaling of the forced convection behavior and asymptotic matching of the entrance or developing flow (short fin-core flow length) and the fully developed flow (large fin-core flow length) region performance. Finally, the predictions from these correlations are found to be within ±15% of all available experimental data for air, water, and glycol (0.71 ≤ Pr ≤ 10), and fin cores with 0 < α ≤ 1.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Direct Measurement of Ice-Ablator Interface Motion for Instability Mitigation in Indirect Drive ICF Implosions

In indirect drive inertial confinement fusion (ICF) implosions hydrodynamic instability growth at the imploding capsule ablator-DT fuel interface can reduce fuel compressibility and inject ablator into the hot spot hence reducing hot spot pressure and temperature. As a mitigation strategy, a gentle acceleration of this interface is predicted by simulations and theory to significantly reduce this instability growth in the early stage of the implosion. We have performed high-contrast, time-resolved x-ray refraction enhanced radiography (RER) to accurately measure the level of acceleration as a function of the initial laser drive time history for indirect-drive implosions on the National Ignition Facility. Here, we demonstrate a transition from no acceleration to 20 ± 1.8μm ns -2 acceleration by tweaking the drive that should reduce the initial instabilities by an order of magnitude at high modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

The dynamics of suspensions of prolate spheroidal particles—Effects of volume fraction

Here, the effect of volume fraction on the dynamics of a gravity driven suspension of prolate spheroidal solid particles in a fully periodic domain is examined by fully resolved numerical simulations, for relatively modest Reynolds numbers (around 20, depending on the volume fraction). Three systems are examined, for volume fractions of 2.2% (20 particles), 5.5% (50 particles), and 9.9% (90 particles). The results show a transition from flow at low volume fractions dominated by the hydrodynamic interactions between the particles and the fluid, although modified by collisions, to flow dominated by collisions at higher volume fractions. The distribution of the particles with respect to each other, as measured by the probability distribution of nearest distances, is nearly random at low volume fractions but at high volume fractions, there is more clustering than for a random distribution. At lower volume fractions most of the particles fall broadside on whereas at the highest volume fraction, their orientation is essentially random.

42 ENGINEERING↗
DOE OSTI2023

As a Matter of Dynamical Range – Scale Dependent Energy Dynamics in MHD Turbulence

Abstract Magnetized turbulence is ubiquitous in many astrophysical and terrestrial plasmas but no universal theory exists. Even the detailed energy dynamics in magnetohydrodynamic (MHD) turbulence are still not well understood. We present a suite of subsonic, super-Alfvénic, high plasma beta MHD turbulence simulations that only vary in their dynamical range, i.e., in their separation between the large-scale forcing and dissipation scales, and their dissipation mechanism (implicit large eddy simulation, ILES, and direct numerical simulation (DNS)). Using an energy transfer analysis framework we calculate the effective numerical viscosities and resistivities, and demonstrate that all ILES calculations of MHD turbulence are resolved and correspond to an equivalent visco-resistive MHD turbulence calculation. Increasing the number of grid points used in an ILES corresponds to lowering the dissipation coefficients, i.e., larger (kinetic and magnetic) Reynolds numbers for a constant forcing scale. Independently, we use this same framework to demonstrate that—contrary to hydrodynamic turbulence—the cross-scale energy fluxes are not constant in MHD turbulence. This applies both to different mediators (such as cascade processes or magnetic tension) for a given dynamical range as well as to a dependence on the dynamical range itself, which determines the physical properties of the flow. We do not observe any indication of convergence even at the highest resolution (largest Reynolds numbers) simulation at 2048 3 cells, calling into question whether an asymptotic regime in MHD turbulence exists, and, if so, what it looks like.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Unraveling the 2021 Central Tennessee flood event using a hierarchical multi-model inundation modeling framework

Flood prediction systems need hierarchical atmospheric, hydrologic, and hydraulic models to predict rainfall, runoff, streamflow, and floodplain inundation. The accuracy of such systems depends on the error propagation through the modeling chain, sensitivity to input data, and choice of models. In this study, we used multiple precipitation forcings (hindcast and forecast) to drive hydrologic and hydrodynamic models to analyze the impacts of various drivers on the estimates of flood inundation depth and extent. We implement this framework to unravel the August 2021 extreme flooding event that occurred in Central Tennessee, USA. We used two radar-based quantitative precipitation estimates (STAGE4 and MRMS) as well as quantitative precipitation forecasts (QPF) from the National Weather Service Weather Prediction Center (WPC) to drive a series of models in the hierarchical framework, including the Variable Infiltration Capacity (VIC) land surface model, the Routing Application for Parallel Computation of Discharge (RAPID) river routing model, and the AutoRoute and TRITON inundation models. An evaluation with observed high-water marks demonstrates that the framework can reasonably simulate flood inundation. Despite the complex error propagation mechanism of the modeling chain, we show that inundation estimates are most sensitive to rainfall estimates. Most notably, QPF significantly underestimates flood magnitudes and inundations leading to unanticipated severe flooding for all stakeholders involved in the event. Finally, we discuss the implications of the hydrodynamic modeling framework for real-time flood forecasting.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2026

Improving coastal water level estimation by merging nadir-only satellite altimetry data into a hydrodynamic model

Providing robust real time flood warnings is of paramount importance to coastal communities. Although state-of-the-art hydrodynamic models are capable of robustly predicting Coastal Water Levels (CWL), unresolved drivers affecting level fluctuations are often not represented by the model governing equations. This work evaluates a novel method to improve the performance of the ADvanced CIRCulation (ADCIRC) hydrodynamic model by assimilating observations from four nadir-only satellite altimetry missions against a set of National Oceanic and Atmospheric Administration (NOAA) gauge stations located across the entire U.S. East Coast. Two different types of simulations were performed – Open Loop (OL) and Data Assimilation (DA). Five different simulations were performed where four different satellite altimetry observations were assimilated individually and combined with two different scenarios – with and without considering the data quality flags. Results indicate that, despite their limited spatial coverage, merging nadir-only observations into ADCIRC from the newly launched Surface Water and Ocean Topography (SWOT)’s nadir altimeter can improve the model performance at 76% of the gauge locations, whereas Sentinel-6 improves 73% of the total locations, Jason-3 74%, and SARAL 21%. Furthermore, combining observations from SWOT-nadir, Jason-3, and Sentinel-6 can improve the ADCIRC performance at more than 80% of the gauge locations for 107-day simulation. Nadir-only satellite altimetry observations can be useful for improving the model performance even if flagged as “poor quality” near the coast. When the flagged data are disregarded, SWOT can improve ADCIRC at 78%, Sentinel-6 at 73%, Jason-3 at 53%, and SARAL at 21% of the gauge locations. The ability to improve the model simulations largely depends on the availability of a satellite overpass nearby. Therefore, model performance can be further enhanced if satellite observations are available during a storm surge event, stressing the importance of frequent satellite overpasses.

Aafnan Bhuiyan, Soelem↗
DOE OSTI2023

Prediction of Breakdown in Air and Solid Dielectrics: A Complete Plasma Model from Discharge Initiation to Flashover (Final Technical Report)

A major challenge for most of the technologies of today, and the technologies of the future, is the compact and reliable design of electrical insulation systems for their electrical components. Although controlled low-temperature plasmas (LTPs) or non-thermal plasmas in air at atmospheric pressure are of interest for a wide range of technologies, the appearance of such discharges is considered highly undesirable in insulation systems because they are precursors of a breakdown or flashover, where complete insulation failure takes place. For non-thermal plasmas, often only one mechanism, or the transition from one discharge mechanism to another, has typically been studied to date. On the other hand, in the case of thermal plasma applications, the conditions following the occurrence of an electric arc or flashover have always been considered. Here, there is a big technical gap: there is no paper describing a complete plasma model simulating all discharge mechanisms and their transitions from discharge initiation to electrical arc and flashover. To address this gap, a complete plasma model for both thermal and non-thermal plasmas should be developed using a hydrodynamic approach. The project aimed to build the model. The work has the potential to address industrial applications and is of great importance to higher voltage systems.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Numerical Modeling of a Two-Stage Ocean Current Turbine

The Equinox Ocean Turbines (EQOT) current energy converter has a unique design with power generation in two small-diameter turbines attached to the tips of a large-diameter passive rotor. This configuration offers some key advantages for capturing ocean currents. With no centrally placed generator, almost no reaction torque is required at the nacelle of the main large-diameter rotor, and the small-diameter tip turbine generators operate at a higher speed and lower torque. The physics that determine the performance and loads on the turbine are also unique. The interactions of the flow field between the two stages and the general architecture of the system cannot be captured with traditional mid-fidelity modeling tools. For design iterations and large sets of load cases, it is important to have mid-fidelity models that can capture the important phenomenon with enough accuracy to identify global trends. This work uses a limited set of high-fidelity computational fluid dynamics (CFD) simulations to help inform the selection of and construction of a custom mid-fidelity model. Mid-fidelity modeling approaches were verified by comparing key turbine performance quantities to those found with the CFD model. Hydrodynamic interactions of the two-stage rotor were identified through high-fidelity CFD modeling. This highlighted the impact of the main rotor tip vortex and wake on the secondary rotor apparent inflow. This results in a relative flow rotation and sharp deficit, that change the optimal secondary rotor rotation speed and adds unsteadiness to the blade loading respectively. Multiple mid-fidelity approaches were evaluated for their ability to capture these effects. A simple approximation of the combined-stage performance based on single-stage BEM provides a reasonable rough prediction, especially near the peak TSR values, with some larger discrepancy at higher TSRs. Predicting the combined-stage performance based on single-stage CFD data improves this prediction across the TSR range. Although the combined-stage modeling in OLAF was not successful in this stage of the project, it showed promise as a mid-fidelity method, assuming the parameters can be tuned to account for the significant differences in time and length scales between the main and secondary rotors. This may be addressed through code changes in future work. A significant finding from the OLAF work was the agreement between the vortex core radius values found independently via a parameter space search and via CFD. The technique of using single-stage secondary rotor BEM, with a custom inflow taken from single-stage main rotor CFD or OLAF, provides an efficient method to capture one-way coupled flow interactions. This method provided generally good predictions of the impact of the flow rotation on the secondary rotor but struggled to accurately predict the peaks of the unsteady load progression. Future work could include some superposition of a tuned main rotor trailing edge viscous wake into the custom inflow to better predict this interaction.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Numerical Model of IProTech PIP WEC Device

iProTech PIP wave energy converter (WEC) is a slack moored, single hull device with no moving parts in the water, joints or bearings. This submission includes data of the simulation, reports, and code for the iProTech PIP (WEC) project. The organization of the data included in the provided archive is detailed below and in the data description of the archive. The data teamer-iprotech-nrel folder includes and explains matlab and python code developed to hydrodynamically model the PIP WEC device in WEC-Sim. The subfolders cover the following steps: 1) report: explanatory information on device geometry 2) pip_mesher: python code to generate mesh panels from device profile data 3) wec-sim_models: matlab code to run WEC-Sim The data uploaded is a snapshot as of 11/02/2121 of code residing in a Github repository administered by David Ogden of NREL.

16 TIDAL AND WAVE POWER↗
DOE OSTI2021

Emergence and dynamics of unconfined self-organised vortices in active magnetic roller liquids

Actively driven colloids demonstrate complex out-of-equilibrium dynamics often rivaling self-organized patterns and collective behavior observed in living systems. Recent studies revealed the emergence of steady macroscopic states with multiple interacting vortices in an unconfined environment that emerge from the coupling between microscale particle rotation and translation. Yet, insights into the microscopic behavior during the vortex emergence, growth, and formation of a multi-vortical state remain lacking. Here, we investigate in experiments and simulations how the microscale magnetic roller behavior leads to the emergence of seed vortices, their aggregation or annihilation, and the formation of stable large-scale vortical structures. We reveal that the coupling of roller-induced hydrodynamic flows guides the local self-densifications and self-organization of the micro-rollers into seed vortices. Here, the resulting multi-vortical state is sensitive to the external magnetic field amplitude and allows tuning the rollers' number density in a vortex and its characteristic size.

36 MATERIALS SCIENCE↗
NASA NTRS2005

On the Origin of the Saturnian Satellite System: Did Iapetus Form In-Situ?

Current models of planet and satellite formation are marred by our lack of understanding regarding the turbulent state of accretion disks. According to the Rayleigh criterion, Keplerian disks are hydrodynamically stable. Indeed, it has been argued that a carefully designed Taylor Couette experiment shows stability in the case of positive radial gradients in specific angular momentum even for high Reynolds numbers [1], in agreement with numerical simulations which consistently show turbulence decay [2]. Other possible sources of turbulence may fail due to low ionization, may decay as the optical depth decreases due to dust coagulation, may involve unrealistic boundary conditions, or result in limited transport. The difficulty stems not only from the degree of turbulence, but also from the kind of turbulence, and whether it may be characterized by an alpha parameter.

Mosqueira, I.↗
DOE OSTI2023

An innovative and alternative approach toward gear fabrication

In this paper, we present the first demonstration of single-step gear manufacturing via an innovative in-situ friction stir forging approach. This novel process is a natural extension to the friction stir processing technique and relies on shear and normal stresses to form the complex shapes. Friction stir welding pin tool with shoulder is plunged at high revolutions per minute (rpm) into the sample to heat it to the desired temperature and enable plastic flow. Then, the tool is forged at high plunge rate to form the desired structurally-sound and defect free complex shapes. To better visualize the process the experiments are complemented by a 3D thermomechanically coupled, smoothed particle hydrodynamics (SPH) model to understand the material flow pattern both radially and through thickness in the gear teeth, plastic strain distribution, and temperature profile. We have demonstrated application of I-FSF to inherently poor formable AZ31 Mg, AA 7075-T6, and AA5083–10 vol% TiB 2 composite and fabricated 14-tooth spur gears from these materials with minimal farther machining. In conclusion, microstructure of the forged gears exhibited a directional material flow pattern due to the normal forging action and refinement of grains to 6–13 μm through dynamically recrystallization across the radial, through-thickness direction, which is in good agreement with the corresponding SPH simulation results.

42 ENGINEERING↗
DOE OSTI2021

Influence of Wetting on Viscous Fingering Via 2D Lattice Boltzmann Simulations

We present simulations of two-phase flow using the Rothman and Keller colour gradient Lattice Boltzmann method to study viscous fingering when a “red fluid” invades a porous model initially filled with a “blue” fluid with different viscosity. We conducted eleven suites of 81 numerical experiments totalling 891 simulations, where each suite had a different random realization of the porous model and spanned viscosity ratios in the range $$M\in [0.01,100]$$ M ϵ [ 0.01 , 100 ] and wetting angles in the range $$\theta _w\in [180^\circ ,0^\circ ]$$ θ w ϵ [ 180 ° , 0 ° ] to allow us to study the effect of these parameters on the fluid-displacement morphology and saturation at breakthrough (sweep). Although sweep often increased with wettability, this was not always so and the sweep phase space landscape, defined as the difference in saturation at a given wetting angle relative to saturation for the non-wetting case, had hills, ridges and valleys. At low viscosity ratios, flow at breakthrough is localized through narrow fingers that span the model. After breakthrough, the flow field continues to evolve and the saturation continues to increase albeit at a reduced rate, and eventually exceeds 90% for both non-wetting and wetting cases. The existence of a complicated sweep phase space at breakthrough, and continued post-breakthrough evolution suggests the hydrodynamics and sweep is a complicated function of wetting angle, viscosity ratio and time, which has major potential implications to Enhanced Oil Recovery by water flooding, and hence, on estimates of global oil reserves. Validation of these results via experiments is required to ensure they translate to field studies.

Engineering↗