Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamic efficiency”

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 235 records · Page 13

CFD Study of Countercurrent Flow in Triply Periodic Minimal Surfaces with CO 2 BOL Solvent

Triply periodic minimal surfaces (TPMS) are promising structures for heat exchangers and ultrafiltration due to their high heat and mass transfer efficiency compared to traditional devices. Because of this, the TPMS structures are a potential candidate for building the multifunctional, intensified devices in CO 2 capture tasks. With the Lawrence Livermore National Laboratory’s (LLNL’s) additive manufacturing (AM) capability, these complicated TPMS structures can be printed in the packed column for experimental tests. To optimize the TPMS structure and enable fast prototyping, a computationally-efficient computational fluid dynamics (CFD) framework was established to explore the TPMS structure’s hydrodynamics and mass transfer performance. The established model was able to simulate the countercurrent flows by integrating the Pacific Northwest National Laboratory (PNNL) CO 2 binding organic liquids (CO 2 BOL) solvent into the TPMS structures. Three typical TPMS structures - namely, Gyroid, Schwarz-D, and Schwarz-P - were investigated using the periodic boundary set up in all three directions. Two different unit cell size (1 cm and 2 cm) geometries are tested for each type of TPMS. The solvent flow rate in this study covered a range of [0, 0.1] m/s with an initial liquid film thickness of [0.1, 1.8] mm. These simulations provided a preliminary understanding of the TPMS structure behaviors in countercurrent flow conditions. The study also related the TPMS geometry requirement to the given CO 2 BOL solvent physical property and covered reasonable packed column operation ranges.

36 MATERIALS SCIENCE↗

Three-dimensional Hydrodynamical Simulations of Mass Transfer in Binary Systems by a Free Wind

A large fraction of stars in binary systems are expected to undergo mass and angular momentum exchange at some point in their evolution, which can drastically alter the chemical and dynamical properties and fates of the systems. Interaction by stellar wind is an important process in wide binaries. However, the details of wind mass transfer are still not well understood. We perform three-dimensional hydrodynamical simulations of wind mass transfer in binary systems to explore mass-accretion efficiencies and geometries of mass outflows, for a range of mass ratios from 0.05 to 1.0. In particular, we focus on the case of a free wind, in which some physical mechanism accelerates the expelled wind material balancing the gravity of the mass-losing star with the wind velocity comparable to the orbital velocity of the system. We find that the mass-accretion efficiency and accreted specific angular momentum increase with the mass ratio of the system. For an adiabatic wind, we obtain that the accretion efficiency onto the secondary star varies from about 0.1% to 8% for mass ratios between 0.05 and 1.0.

79 ASTRONOMY AND ASTROPHYSICS↗

Observation of hydrodynamic processes of radiation-ablated plasma in a small hole

In the hohlraum used in laser indirect-drive inertial confinement fusion experiments, hydrodynamic processes of radiation-ablated high-Z plasma have a great effect on laser injection efficiency, radiation uniformity, and diagnosis of hohlraum radiation field from diagnostic windows (DW). To study plasma filling in the DWs, a laser-irradiated Ti disk was used to generate 2–5 keV narrow energy band X-ray as the intense backlighter source, and laser-produced X-ray in a hohlraum with low-Z foam tamper was used to heat a small hole surrounded by gold wall with 150 μm in diameter and 100 μm deep. The hydrodynamic movement of the gold plasma in the small hole was measured by an X-ray framing camera and the results are analyzed. Quantitative measurement of the plasma areal density distribution and evolution in the small hole can be used to assess the effect of plasma filling on the diagnosis from the DWs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis of eccentric annular incompressible seals. I - A new solution using fast Fourier transforms for determining hydrodynamic force

A new analysis procedure is presented which solves for the flow variables of an incompressible-flow annular pressure seal in which the rotor has a large static displacement from the centered position. The analysis begins with a set of governing equations based on a turbulent bulk-flow model and Moody's friction equation. No simplification of these bulk-flow equations is required for the solution procedure. Perturbation of the flow variables yields a set of zeroth and first-order equations. The zeroth-order equations (which model the large static displacement) are integrated by means of an efficient new method which employs fast Fourier transforms. Further integration of the zeroth-order pressures yields the hydrodynamic reactive force. Predictions for the hydrodynamic forces from this analysis procedure are in excellent agreement with available experimental results.

Nelson, C. C.↗

Blow-off momentum from melt and vapor in nuclear deflection scenarios

For Earth-impacting objects that are large in size or have short warning times nuclear explosives are an effective threat mitigation response. Nuclear-based deflection works by means of conservation of momentum: as material is heated by incoming photons and neutrons it is ejected from the body which imparts momentum to the remaining mass of the asteroid. Predicting the complete response of a particular object is difficult, since the ejecta size and velocity distributions rely heavily on the unknown, complicated internal structure of the body. However, lower bounds on the blow-off momentum can be estimated using the melted and vaporized surface material. In this paper, we model the response of a one-dimensional SiO2 surface to monoenergetic soft X-ray, hard X-ray and neutron sources using Arbitrary Lagrangian–Eulerian radiation/hydrodynamic simulations. Errors in the blow-off momentum due to our hydrodynamic mesh resolution are quantified and inform zone sizing that balances numerical discretization error with computational efficiency. We explore deposited energy densities ranging from 1.1 to 200 times the melt energy density for SiO2, and develop an approximate relation that gives the mesh resolution needed for a desired percent error in the blow-off momentum as a function of deposited energy density and melt depth. Using these mesh constraints, the response of our one-dimensional SiO2 surface to the energy sources is simulated, and lower bounds are placed on the melt/vapor blow-off momentum as a function of deposited energy density and source energy type.

79 ASTRONOMY AND ASTROPHYSICS↗

The Quantum Hydrodynamic Description of Tunneling

The quantum hydrodynamic approach is based on the de Broglie-Bohm formulation of quantum mechanics. The resulting fluid-like equations of motion describe the flow of probability and an accurate solution to these equations is equivalent to solving the time-dependent Schroedinger equation. Furthermore, the hydrodynamic approach provides new insight into the mechanisms as well as an alternative computational approach for treating tunneling phenomena. New concepts include well-defined 'quantum trajectories', 'quantum potential', and 'quantum force' all of which have classical analogues. The quantum potential and its associated force give rise to all quantum mechanical effects such as zero point energy, tunneling, and interference. A new numerical approach called the Iterative Finite Difference Method (IFDM) will be discussed. The IFDM is used to solve the set of non-linear coupled hydrodynamic equations. It is 2nd-order accurate in both space and time and exhibits exponential convergence with respect to the iteration count. The stability and computational efficiency of the IFDM is significantly improved by using a 'smart' Eulerian grid which has the same computational advantages as a Lagrangian or Arbitrary Lagrangian Eulerian (ALE) grid. The IFDM is also capable of treating anharmonic potentials. Example calculations using the IFDM will be presented which include: a one-dimensional Gaussian wave packet tunneling through an Eckart barrier, a one-dimensional bound-state Morse oscillator, and a two-dimensional (2D) model collinear reaction using an anharmonic potential energy surface. Approximate treatments of the quantum hydrodynamic equations will also be discussed which could allow scaling of the calculations to hundreds of degrees of freedom which is important for treating tunneling phenomena in condensed phase systems.

36 MATERIALS SCIENCE↗

Planet Formation by Gas-assisted Accretion of Small Solids

We compute the accretion efficiency of small solids, with radii 1 cm ≤ R s ≤ 10 m, on planets embedded in gaseous disks. Planets have masses 3 ≤ M p ≤ 20 Earth masses (M ⊕ ) and orbit within 10 au of a solar mass star. Disk thermodynamics is modeled via 3D radiation-hydrodynamics calculations that typically resolve the planetary envelopes. Both icy and rocky solids are considered, explicitly modeling their thermodynamic evolution. The maximum efficiencies of 1 ≤ R s ≤ 100 cm particles are generally ≲10%, whereas 10 m solids tend to accrete efficiently or be segregated beyond the planet’s orbit. A simplified approach is applied to compute the accretion efficiency of small cores, with masses M p ≤ 1 M ⊕ and without envelopes, for which efficiencies are approximately proportional to $M^{2/3}_{p}$. The mass flux of solids, estimated from unperturbed drag-induced drift velocities, provides typical accretion rates dM p /dt ≲ 10 -5 M ⊕ yr -1 . In representative disk models with an initial gas-to-dust mass ratio of 70–100 and total mass of 0.05–0.06 M ⊙ , the solids’ accretion falls below 10 -6 M ⊕ yr -1 after 1–1.5 Myr. The derived accretion rates, as functions of time and planet mass, are applied to formation calculations that compute dust opacity self-consistently with the delivery of solids to the envelope. Assuming dust-to-solid coagulation times of ≈0.3 Myr and disk lifetimes of ≈3.5 Myr, heavy-element inventories in the range 3–7 M ⊕ require that ≈90–150 M ⊕ of solids cross the planet’s orbit. The formation calculations encompass a variety of outcomes, from planets a few times M ⊕ , predominantly composed of heavy elements, to giant planets. The peak luminosities during the epoch of the solids’ accretion range from ≈10 -7 to ≈10 -6 L ⊙ .

79 ASTRONOMY AND ASTROPHYSICS↗

Origin and dynamics of emission line clouds in cooling flow environments

The author suggests that since clouds are born co-moving in a turbulent intra-cluster medium (ICM), the allowed parameter space can now be opened up to a more acceptable range. Large-scale motions can be driven in the central parts of cooling flows by a number of mechanisms including the motion of the central and other galaxies, and the dissipation of advected, focussed rotational and magnetic energy. In addition to the velocity width paradox, two other paradoxes (Heckman et al. 1989) can be solved if the ICM is turbulent. Firstly, the heating source for the emission line regions has always been puzzling - line luminosities are extremely high for a given (optical or radio) galaxy luminosity compared to those in non-cooling flow galaxies, therefore a mechanism peculiar to cooling flows must be at work. However most, if not all, previously suggested heating mechanisms either fail to provide enough ionization or give the wrong line ratios, or both. The kinetic energy in the turbulence provides a natural energy source if it can be efficiently converted to cloud heat. Researchers suggest that this can be done via magneto-hydrodynamic waves through plasma slip. Secondly, while the x ray observations indicate extended mass deposition, the optical line emission is more centrally concentrated. Since many of the turbulence-inducing mechanisms are strongest in the central regions of the ICM, so is the method of heating. In other words material is dropping out everywhere but only being lit up in the center.

Loewenstein, Michael↗

Simulation of 100-300 GHz solid-state harmonic sources

Accurate and efficient simulations of the large-signal time-dependent characteristics of second-harmonic Transferred Electron Oscillators (TEO's) and Heterostructure Barrier Varactor (HBV) frequency triplers have been obtained. This is accomplished by using a novel and efficient harmonic-balance circuit analysis technique which facilitates the integration of physics-based hydrodynamic device simulators. The integrated hydrodynamic device/harmonic-balance circuit simulators allow TEO and HBV circuits to be co-designed from both a device and a circuit point of view. Comparisons have been made with published experimental data for both TEO's and HBV's. For TEO's, excellent correlation has been obtained at 140 GHz and 188 GHz in second-harmonic operation. Excellent correlation has also been obtained for HBV frequency triplers operating near 200 GHz. For HBV's, both a lumped quasi-static equivalent circuit model and the hydrodynamic device simulator have been linked to the harmonic-balance circuit simulator. This comparison illustrates the importance of representing active devices with physics-based numerical device models rather than analytical device models.

NONLINEAR CIRCUITS↗

Sub-Scale Testing and Development of the J-2X Fuel Turbopump Inducer

In the early stages of the J-2X upper stage engine program, various inducer configurations proposed for use in the fuel turbopump (FTP) were tested in water. The primary objectives of this test effort were twofold. First, to obtain a more comprehensive data set than that which existed in the Pratt & Whitney Rocketdyne (PWR) historical archives from the original J-2S program, and second, to supplement that data set with information regarding the cavitation induced vibrations for both the historical J-2S configuration as well as those tested for the J-2X program. The J-2X FTP inducer, which actually consists of an inducer stage mechanically attached to a kicker stage, underwent 4 primary iterations utilizing sub-scaled test articles manufactured and tested in PWR's Engineering Development Laboratory (EDL). The kicker remained unchanged throughout the test series. The four inducer configurations tested retained many of the basic design features of the J-2S inducer, but also included variations on leading edge blade thickness and blade angle distribution, primarily aimed at improving suction performance at higher flow coefficients. From these data sets, the effects of the tested design variables on hydrodynamic performance and cavitation instabilities were discerned. A limited comparison of impact to the inducer efficiency was determined as well.

Sargent, Scott R.↗

Hydrogen and Electric Power Generation from Liquid Microjets: Design Principles for Optimizing Conversion Efficiency

Historically, Liquid water microjets have been successfully employed for both electrical power generation and gaseous hydrogen production, but the demonstrated efficiencies have been low. Here, we employ a combination of a modified Poisson-Boltzmann description, continuum hydrodynamic equations, and microjet electrokinetic experiments to gain detailed insight into the origin of the streaming currents produced in pure water. We identify the contributions to the streaming current from specific ion adsorption at the solid/liquid interface and from long-ranged electrostatic interactions, finding that the portion originating from the latter dominate at charged surfaces. The detailed understanding afforded by theory and the close agreement with experimental results elucidates design principles for optimizing hydrogen production and power generation. Changing the sign of the surface charge density through targeted use of surface coatings via silanization switches the primary charge carrier between hydronium and hydroxide and therefore switches the corresponding production of molecular hydrogen to oxygen at the target electrode. Moreover, hydrophobic surface coatings reduce dissipation due to fluid/solid friction, thereby increasing the conversion efficiency.

08 HYDROGEN↗

Effects of Contact Angle Hysteresis on Ice Adhesion and Growth over Superhydrophobic Surfaces under Dynamic Flow Conditions

In this paper, the iceophobic properties of superhydrophobic surfaces are investigated under dynamic flow conditions by using a closed loop low-temperature wind tunnel. Superhydrophobic surfaces were prepared by coating the substrates of aluminum and steel plates with nano-structured hydrophobic particles. The superhydrophobic plates along with uncoated control ones were exposed to an air flow of 12 m/s and 20 F accompanying micron-sized water droplets in the icing wind tunnel and the ice formation and accretion were probed by high-resolution CCD cameras. Results show that the superhydrophobic coatings significantly delay the ice formation and accretion even under the dynamic flow condition of the highly energetic impingement of accelerated super-cooled water droplets. It is found that there is a time scale for this phenomenon (delay of the ice formation) which has a clear correlation with the contact angle hysteresis and the length scale of surface roughness of the superhydrophobic surface samples, being the highest for the plate with the lowest contact angle hysteresis and finer surface roughness. The results suggest that the key parameter for designing iceophobic surfaces is to retain a low contact angle hysteresis (dynamic property) and the non-wetting superhydrophobic state under the hydrodynamic pressure of impinging droplets, rather than to only have a high contact angle (static property), in order to result in efficient anti-icing properties under dynamic conditions such as forced flows.

36 MATERIALS SCIENCE↗

3D mesh regularization within an ALE code using a weighted line sweeping method

The Lagrangian formalism is widely used to simulate hydrodynamic responses in complex engineering applications, particularly those involving strong shock waves. However, as the mesh moves with the fluid, it can become highly distorted, requiring a regularization step. This involves constructing a new grid and remapping conservative quantities onto it to restore mesh quality. This work introduces a regularization method for block-structured meshes within a 3D ALE (Arbitrary Lagrangian-Eulerian) code. The proposed approach prevents mesh tangling while preserving the anisotropic features of the initial Lagrangian mesh. This regularization technique incorporates aspect ratio-based weights to control mesh smoothing. Unlike uniform rezoning techniques, this weighted approach maintains proximity to the Lagrangian mesh while improving mesh quality. Here, the method effectively handles concave geometries by mitigating the grid attraction phenomenon, which typically leads to mesh concentration along concave edges. Numerical experiments demonstrate its efficiency in regularizing severely deformed meshes, and its integration within the ALE framework is validated on challenging hydrodynamic test cases, including the triple point problem.

42 ENGINEERING↗

Brownian dynamics of confined rigid bodies

We introduce numerical methods for simulating the diffusive motion of rigid bodies of arbitrary shape immersed in a viscous fluid. We parameterize the orientation of the bodies using normalized quaternions, which are numerically robust, space efficient, and easy to accumulate. We construct a system of overdamped Langevin equations in the quaternion representation that accounts for hydrodynamic effects, preserves the unit-norm constraint on the quaternion, and is time reversible with respect to the Gibbs-Boltzmann distribution at equilibrium. We introduce two schemes for temporal integration of the overdamped Langevin equations of motion, one based on the Fixman midpoint method and the other based on a random finite difference approach, both of which ensure that the correct stochastic drift term is captured in a computationally efficient way. We study several examples of rigid colloidal particles diffusing near a no-slip boundary and demonstrate the importance of the choice of tracking point on the measured translational mean square displacement (MSD). We examine the average short-time as well as the long-time quasi-two-dimensional diffusion coefficient of a rigid particle sedimented near a bottom wall due to gravity. For several particle shapes, we find a choice of tracking point that makes the MSD essentially linear with time, allowing us to estimate the long-time diffusion coefficient efficiently using a Monte Carlo method. However, in general, such a special choice of tracking point does not exist, and numerical techniques for simulating long trajectories, such as the ones we introduce here, are necessary to study diffusion on long time scales.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Formulation and Validation of an Efficient Computational Model for a Dilute, Settling Suspension Undergoing Rotational Mixing

Designing processing equipment for the mixing of settling suspensions is a challenging problem. Achieving low-cost mixing is especially difficult for the application of slowly reacting suspended solids because the cost of impeller power consumption becomes quite high due to the long reaction times (batch mode) or due to large-volume reactors (continuous mode). Further, the usual scale-up metrics for mixing, e.g., constant tip speed and constant power per volume, do not apply well for mixing of suspensions. As an alternative, computational fluid dynamics (CFD) can be useful for analyzing mixing at multiple scales and determining appropriate mixer designs and operating parameters. We developed a mixture model to describe the hydrodynamics of a settling cellulose suspension. The suspension motion is represented as a single velocity field in a computationally efficient Eulerian framework. The solids are represented by a scalar volume-fraction field that undergoes transport due to particle diffusion, settling, fluid advection, and shear stress. A settling model and a viscosity model, both functions of volume fraction, were selected to fit experimental settling and viscosity data, respectively. Simulations were performed with the open-source Nek5000 CFD program, which is based on the high-order spectral-finite-element method. Simulations were performed for the cellulose suspension undergoing mixing in a laboratory-scale vane mixer. The settled-bed heights predicted by the simulations were in semi-quantitative agreement with experimental observations. Further, the simulation results were in quantitative agreement with experimentally obtained torque and mixing-rate data, including a characteristic torque bifurcation. In future work, we plan to couple this CFD model with a reaction-kinetics model for the enzymatic digestion of cellulose, allowing us to predict enzymatic digestion performance for various mixing intensities and novel reactor designs.

09 BIOMASS FUELS↗

LOW-MASS STAR FORMATION TRIGGERED BY EARLY SUPERNOVA EXPLOSIONS

We study the formation of low-mass and extremely metal-poor stars in the early universe. Our study is motivated by the recent discovery of a low-mass (M {sub *} {<=} 0.8 M {sub Sun }) and extremely metal-poor (Z {<=} 4.5 Multiplication-Sign 10{sup -5} Z {sub Sun }) star in the Galactic halo by Caffau et al. We propose a model that early supernova (SN) explosions trigger the formation of low-mass stars via shell fragmentation. We first perform one-dimensional hydrodynamic simulations of the evolution of an early SN remnant. We show that the shocked shell undergoes efficient radiative cooling and then becomes gravitationally unstable to fragment and collapse in about a million years. We then follow the thermal evolution of the collapsing fragments using a one-zone code. Our one-zone calculation treats chemistry and radiative cooling self-consistently in low-metallicity gas. The collapsing gas cloud evolves roughly isothermally, until it cools rapidly by dust continuum emission at the density 10{sup 13}-10{sup 14} cm{sup -3}. The cloud core then becomes unstable and fragments again. We argue that early SNe can trigger the formation of low-mass stars in the extremely metal-poor environment as Caffau et al. discovered recently.

79 ASTRONOMY AND ASTROPHYSICS↗

RANS Simulation of the Heave Response of a Two-Body Floating Point Wave Absorber: Preprint

A preliminary study on a two-body floating wave absorbers is presented in this paper. A Reynolds-Averaged Navier-Stokes computational method is applied for analyzing the hydrodynamic heave response of the absorber in operational wave conditions. The two-body floating wave absorber contains a float section and a submerged reaction section. For validation purposes, our model is first assumed to belocked. The two sections are forced to move together with each other. The locked single body model is used in a heave decay test, where the RANS result is validated with the experimental measurement. For the two-body floating point absorber simulation, the two sections are connected through a mass-spring-damper system, which is applied to simulate the power take-off mechanism under design waveconditions. Overall, the details of the flow around the absorber and its nonlinear interaction with waves are investigated, and the power absorption efficiency of the two-body floating wave absorber in waves with a constant value spring-damper system is examined.

17 WIND ENERGY↗

High efficiency proton beam generation through target thickness control in femtosecond laser-plasma interactions

Bright proton beams with maximum energies of up to 30 MeV have been observed in an experiment investigating ion sheath acceleration driven by a short pulse (<50 fs) laser. The scaling of maximum proton energy and total beam energy content at ultra-high intensities of ∼10{sup 21} W cm{sup −2} was investigated, with the interplay between target thickness and laser pre-pulse found to be a key factor. While the maximum proton energies observed were maximised for μm-thick targets, the total proton energy content was seen to peak for thinner, 500 nm, foils. The total proton beam energy reached up to 440 mJ (a conversion efficiency of 4%), marking a significant step forward for many laser-driven ion applications. The experimental results are supported by hydrodynamic and particle-in-cell simulations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗