Engineering Papers⌕ Search

SEARCH · Engineering Papers

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 451 records · Page 25

Geometric Scale-up Experiments on Fluidization of Geldart B Glass Beads

The objective of this work is to provide a valuable database from controlled experiments for validating computational models. Recently, coarse-grained techniques such as particle-in-cell (PIC) or coarse-grained discrete element modeling (DEM) have gained popularity due to their computational efficiency while modeling large-scale systems; however, the influence of model parameters and their sensitivities at different geometric scales and flow conditions remain to be analyzed. These datasets are critical for the multiphase flow research community to assess predictive capability of modeling techniques as well as elucidate the hydrodynamic behavior in these systems. This study performed fluidization experiments using three different test sections with internal diameters of 2.5, 4, and 6 in. The operating conditions, bed material, and range of flow velocities at the inlet were constant in all the units, which were not hydrodynamically scaled. Glass beads having a Sauter Mean Diameter of 332 μm were used. Superficial velocity was varied from 2.97 to 5.35 times the minimum fluidization velocity and the initial static bed height was 0.1524 m. The order in which the experiments were performed was randomized and replicates were included to provide uncertainty in measurements. Statistics of differential pressure and bed height from these tests were reported. Future plans include validating PIC methodology in the open-source software, MFiX (Multiphase Flow with Interphase Exchanges) using results from this study. This could further be extended to determine optimal model parameters using inverse techniques such as deterministic calibration or Bayesian inference.

42 ENGINEERING↗

TURBULENCE, TRANSPORT, AND WAVES IN OHMIC DEAD ZONES

We use local numerical simulations to study a vertically stratified accretion disk with a resistive mid-plane that damps magnetohydrodynamic (MHD) turbulence. This is an idealized model for the dead zones that may be present at some radii in protoplanetary and dwarf novae disks. We vary the relative thickness of the dead and active zones to quantify how forced fluid motions in the dead zone change. We find that the residual Reynolds stress near the mid-plane decreases with increasing dead zone thickness, becoming negligible in cases where the active to dead mass ratio is less than a few percent. This implies that purely Ohmic dead zones would be vulnerable to episodic accretion outbursts via the mechanism of Martin and Lubow. We show that even thick dead zones support a large amount of kinetic energy, but this energy is largely in fluid motions that are inefficient at angular momentum transport. Confirming results from Oishi and Mac Low, the perturbed velocity field in the dead zone is dominated by an oscillatory, vertically extended circulation pattern with a low frequency compared to the orbital frequency. This disturbance has the properties predicted for the lowest order r mode in a hydrodynamic disk. We suggest that in a global disk similar excitations would lead to propagating waves, whose properties would vary with the thickness of the dead zone and the nature of the perturbations (isothermal or adiabatic). Flows with similar amplitudes would buckle settled particle layers and could reduce the efficiency of pebble accretion.

79 ASTRONOMY AND ASTROPHYSICS↗

Data Assimilation and Transport Modeling in Terrestrial and Planetary Atmospheres

Data assimilation is a blanket term used to describe a number of techniques for retrieving important physical parameters from observational data, subject to constraints imposed by prior knowledge (such as, in the case of meteorology, the primitive equations that govern atmospheric motion). Since these newly developed methods make efficient use of computational resources, they are of great importance in the interpretation of the voluminous datasets that are now produced by satellite missions. As proposed, these techniques have been applied to the study of the Martian and terrestrial atmospheres based on available satellite observations. In addition, a sophisticated hydrodynamic model (non-hydrostatic, and therefore applicable to the study of the interiors of the giant planets) has also been developed and successfully applied to the study of tidally induced motions in Jupiter.

Houben, Howard C.↗

Design and evaluation of a dilute flow particle-to-air heat exchanger for energy storage applications

The use of inert and redox-active particles for high-temperature energy storage requires the development of components that can efficiently transfer energy to high-pressure working fluids like supercritical carbon dioxide (sCO 2 ). Dilute flow reactors can enable high working fluid outlet temperatures and minimal parasitic losses compared to moving packed bed and fluidized bed reactors. This research uses both computational and experimental methods to explore the design trade-offs and practical challenges of a novel component for transferring energy from dilute flows of hot, reduced metal oxide (MO x ) particles to sCO 2 in tubes. A discretized thermal resistance network model, which accounts for particle hydrodynamics, multi-mode heat transfer, and reaction equilibrium, guides the design of a prototype device. This device is experimentally tested with a surrogate heat transfer fluids and inert particle temperatures up to 400°C and a heat duty exceeding 1 kW. The data are used to validate the thermal hydraulic sub-models, allowing for the simulation of reacting particle scenarios. Under nominal design conditions, the flow rate of reactive particles is predicted to be 30% lower than that of inert particles for the same energy recovered, with over 70% of the stored particle energy transferred to the sCO 2 . Furthermore, these findings can inform the design of more efficient energy recovery reactors for particle-based systems and can be integrated into system-level concentrated solar power models with thermal storage to optimize operating conditions.

14 SOLAR ENERGY↗

Low Mach number fluctuating hydrodynamics for electrolytes

Here, we formulate and study computationally the low Mach number fluctuating hydrodynamic equations for electrolyte solutions. We are also interested in studying transport in mixtures of charged species at the mesoscale, down to scales below the Debye length, where thermal fluctuations have a significant impact on the dynamics. Continuing our previous work on fluctuating hydrodynamics of multicomponent mixtures of incompressible isothermal miscible liquids (A. Donev, et al., Physics of Fluids, 27, 3, 2015), we now include the effect of charged species using a quasielectrostatic approximation. Localized charges create an electric field, which in turn provides additional forcing in the mass and momentum equations. Our low Mach number formulation eliminates sound waves from the fully compressible formulation and leads to a more computationally efficient quasi-incompressible formulation. Furthermore, we demonstrate our ability to model saltwater (NaCl) solutions in both equilibrium and nonequilibrium settings. We show that our algorithm is second-order in the deterministic setting, and for length scales much greater than the Debye length gives results consistent with an electroneutral/ambipolar approximation. In the stochastic setting, our model captures the predicted dynamics of equilibrium and nonequilibrium fluctuations. We also identify and model an instability that appears when diffusive mixing occurs in the presence of an applied electric field.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Toward connecting core-collapse supernova theory with observations. I. Shock revival in a 15 M {sub ☉} blue supergiant progenitor with SN 1987A energetics

We study the evolution of the collapsing core of a 15 M {sub ☉} blue supergiant supernova progenitor from the core bounce until 1.5 s later. We present a sample of hydrodynamic models parameterized to match the explosion energetics of SN 1987A. We find the spatial model dimensionality to be an important contributing factor in the explosion process. Compared to two-dimensional (2D) simulations, our three-dimensional (3D) models require lower neutrino luminosities to produce equally energetic explosions. We estimate that the convective engine in our models is 4% more efficient in 3D than in 2D. We propose that the greater efficiency of the convective engine found in 3D simulations might be due to the larger surface-to-volume ratio of convective plumes, which aids in distributing energy deposited by neutrinos. We do not find evidence of the standing accretion shock instability or turbulence being a key factor in powering the explosion in our models. Instead, the analysis of the energy transport in the post-shock region reveals characteristics of penetrative convection. The explosion energy decreases dramatically once the resolution is inadequate to capture the morphology of convection on large scales. This shows that the role of dimensionality is secondary to correctly accounting for the basic physics of the explosion. We also analyze information provided by particle tracers embedded in the flow and find that the unbound material has relatively long residency times in 2D models, while in 3D a significant fraction of the explosion energy is carried by particles with relatively short residency times.

79 ASTRONOMY AND ASTROPHYSICS↗

Drivers and impacts of sediment deposition in Amazonian floodplains

Abstract The Amazon River carries enormous amounts of sediment from the Andes mountains, much of which is deposited in its floodplains. However, accurate quantification of the sediment sink at fine spatiotemporal scales is still challenging. Here, we present a high-resolution hydrodynamic-sediment model to simulate sediment deposition in a representative Amazon/Solimões floodplain. The process is found to be jointly driven by inundation, suspended sediment concentration in the Amazon River, and floodplain hydrodynamics and only weakly correlated with inundation level. By upscaling the sediment deposition rate (1.33 ± 0.24 kg m −2 yr −1 ), we estimate the trapping of 77.3 ± 13.9 Mt (or 6.1 ± 1%) of the Amazon River sediment by the Amazon/Solimões floodplains every year. Widespread deforestation would reduce the trapping efficiency of the floodplains over time, exacerbating downstream river aggradation. Additionally, we show that the deposition of sediment-associated organic carbon plays a minor role in fueling carbon dioxide and methane emissions in the Amazon.

54 ENVIRONMENTAL SCIENCES↗

Advancing material modeling in hydrocodes using a concurrent finite-element and molecular dynamics multiscale framework

We present a multiscale simulation framework that couples the finite-element method with molecular dynamics. Bypassing traditional equations of state (EOS) by using in-line atomistic simulations, the method offers the advantage of incorporating detailed microscale physics not easily represented with coarse-grained models. Coupling consistency with the continuum code is ensured through the use of lifting and restriction operators, in line with heterogeneous multiscale methods. The concurrent continuum-atomistic framework is validated through comparison with experimental results and conventional EOS models, and demonstrated in a shock-driven hydrodynamic flow simulation under extreme conditions. We further evaluate the framework's usability by comparing it to state-of-the-art EOS models of deuterium. A computational performance study reveals that the atomistic EOS evaluation is a feasible alternative to conventional approaches, and demonstrates a weak scaling of 99% efficiency. These results highlight the framework's potential for large-scale multiscale modeling across a broad range of materials and conditions.

Computer science↗

CAN STAR–DISK COLLISIONS EXPLAIN THE MISSING RED GIANTS PROBLEM IN THE GALACTIC CENTER?

Observations have revealed a relative paucity of red giant (RG) stars within the central 0.5 pc in the Galactic Center (GC). Motivated by this finding we investigate the hypothesis that collisions of stars with a fragmenting accretion disk are responsible for the observed dearth of evolved stars. We use three-dimensional hydrodynamic simulations to model a star with radius 10 R {sub ⊙} and mass 1 M {sub ⊙}, representative of the missing population of RGs, colliding with high density clumps. We find that multiple collisions with clumps of column density ≳10{sup 8} g cm{sup −2} can strip a substantial fraction of the star’s envelope and in principle render it invisible to observations. Simulations confirm that repeated impacts are particularly efficient in driving mass loss as partially stripped RGs expand and have increased cross sections for subsequent collisions. Because the envelope is unbound on account of the kinetic energy of the star, any significant amount of stripping of the RG population in the GC should be mirrored by a systematic decay of their orbits and possibly by their enhanced rotational velocity. To be viable, this scenario requires that the total mass of the fragmenting disk has been several orders of magnitude higher than that of the early-type stars which now form the stellar disk in the GC.

79 ASTRONOMY AND ASTROPHYSICS↗

A comparison of two- and three-dimensional neutrino-hydrodynamics simulations of core-collapse supernovae

We present numerical results on two- (2D) and three-dimensional (3D) hydrodynamic core-collapse simulations of an 11.2 M {sub ☉} star. By changing numerical resolutions and seed perturbations systematically, we study how the postbounce dynamics are different in 2D and 3D. The calculations were performed with an energy-dependent treatment of the neutrino transport based on the isotropic diffusion source approximation scheme, which we have updated to achieve a very high computational efficiency. All of the computed models in this work, including nine 3D models and fifteen 2D models, exhibit the revival of the stalled bounce shock, leading to the possibility of explosion. All of them are driven by the neutrino-heating mechanism, which is fostered by neutrino-driven convection and the standing-accretion-shock instability. Reflecting the stochastic nature of multi-dimensional (multi-D) neutrino-driven explosions, the blast morphology changes from model to model. However, we find that the final fate of the multi-D models, whether an explosion is obtained or not, is little affected by the explosion stochasticity. In agreement with some previous studies, higher numerical resolutions lead to slower onset of the shock revival in both 2D and 3D. Based on the self-consistent supernova models leading to the possibility of explosions, our results systematically show that the revived shock expands more energetically in 2D than in 3D.

79 ASTRONOMY AND ASTROPHYSICS↗

Growth of First Galaxies: Impacts of Star Formation and Stellar Feedback

Recent observations have detected galaxies at high-redshift z∼6--11, and revealed the diversity of their physical properties, from normal star-forming galaxies to starburst galaxies. To understand the properties of these observed galaxies, it is crucial to understand the star formation (SF) history of high-redshift galaxies under the influence of stellar feedback. In this work, we present the results of cosmological hydrodynamic simulations with zoom-in initial conditions, and investigate the formation of the first galaxies and their evolution toward observable galaxies at z∼6. We focus on three different galaxies that end up in halos with masses M{sub h}=2.4×10{sup 10} h{sup −1} M{sub ⊙} (Halo-10), 1.6×10{sup 11} h{sup −1} M{sub ⊙} (Halo-11), and 0.7×10{sup 12} h{sup −1} M{sub ⊙} (Halo-12) at z = 6. Our simulations also probe the impacts of different subgrid assumptions, i.e., SF efficiency and cosmic reionization, on SF histories in the first galaxies. We find that SF occurs intermittently due to supernova (SN) feedback at z≳10, and then it proceeds more smoothly as the halo mass grows at lower redshifts. Galactic disks are destroyed due to SN feedback, while galaxies in simulations with no feedback or lower SF efficiency models can sustain a galactic disk for long periods ≳10 Myr. The expulsion of gas at the galactic center also affects the inner dark matter density profile for a short period. Our simulated galaxies in Halo-11 and Halo-12 reproduce the SF rates and stellar masses of observed Lyα emitters at z∼7--8 fairly well given the observational uncertainties.

79 ASTRONOMY AND ASTROPHYSICS↗

Carbon Storage in Fold‐and‐Thrust Belts: An Overlooked Gigatonne Storage Opportunity

This study presents numerical investigations of the trapping characteristics of fold-and-thrust belt structures, defining three carbon capture and storage (CCS) play types that could be used to store commercial volumes (millions of tonnes) of CO 2 . Specifically, we present simulations of CO 2 storage in three fold-and-thrust belt models comprising a thrust-ramp, duplex, and thrust-fold geometry. To constrain these play types in realistic geology, each model is based on a study site, including a novel investigation of a greenfield saline reservoir in Virginia, USA, being considered for commercial carbon storage and two well-characterized petroleum fields: the Wilburton field in Oklahoma, USA, and the Incahuasi field in Bolivia. Our results provide insight into several key parameters, such as the long-term security of injected CO 2 in these geologies and injection strategies for maximizing storage efficiency while reducing pressure-related risk. These results improve the understanding of CCS in fold-and-thrust belt storage sites globally by describing general storage parameters that may be applied to site-specific projects. We find that thrust-ramp geometries may securely trap CO 2 through solubility and hydrodynamic trapping under suitable reservoir conditions, duplex structures may store some quantities of CO 2 but are pressure-constrained, and that thrust-ramp structures may store large quantities of CO 2 by maximizing fetch volume, which simultaneously lowers geomechanical risk by reducing pressure buildup along zones of weakness.

42 ENGINEERING↗

Heating the intergalactic medium by X-rays from population III binaries in high-redshift galaxies

Due to their long mean free path, X-rays are expected to have an important impact on cosmic reionization by heating and ionizing the intergalactic medium (IGM) on large scales, especially after simulations have suggested that Population III (Pop III) stars may form in pairs at redshifts as high as 20-30. We use the Pop III distribution and evolution from a self-consistent cosmological radiation hydrodynamic simulation of the formation of the first galaxies and a simple Pop III X-ray binary model to estimate their X-ray output in a high-density region larger than 100 comoving (Mpc){sup 3}. We then combine three different methods—ray tracing, a one-zone model, and X-ray background modeling—to investigate the X-ray propagation, intensity distribution, and long-term effects on the IGM thermal and ionization state. The efficiency and morphology of photoheating and photoionization are dependent on the photon energies. The sub-kiloelectronvolt X-rays only impact the IGM near the sources, while the kiloelectronvolt photons contribute significantly to the X-ray background and heat and ionize the IGM smoothly. The X-rays just below 1 keV are most efficient in heating and ionizing the IGM. We find that the IGM might be heated to over 100 K by z = 10 and the high-density source region might reach 10{sup 4} K, limited by atomic hydrogen cooling. This may be important for predicting the 21 cm neutral hydrogen signals. On the other hand, the free electrons from X-ray ionizations are not enough to contribute significantly to the optical depth of the cosmic microwave background to the Thomson scattering.

79 ASTRONOMY AND ASTROPHYSICS↗

HIGH-TEMPERATURE PROCESSING OF SOLIDS THROUGH SOLAR NEBULAR BOW SHOCKS: 3D RADIATION HYDRODYNAMICS SIMULATIONS WITH PARTICLES

A fundamental, unsolved problem in solar system formation is explaining the melting and crystallization of chondrules found in chondritic meteorites. Theoretical models of chondrule melting in nebular shocks have been shown to be consistent with many aspects of thermal histories inferred for chondrules from laboratory experiments; but, the mechanism driving these shocks is unknown. Planetesimals and planetary embryos on eccentric orbits can produce bow shocks as they move supersonically through the disk gas, and are one possible source of chondrule-melting shocks. We investigate chondrule formation in bow shocks around planetoids through three-dimensional radiation hydrodynamics simulations. A new radiation transport algorithm that combines elements of flux-limited diffusion and Monte Carlo methods is used to capture the complexity of radiative transport around bow shocks. An equation of state that includes the rotational, vibrational, and dissociation modes of H{sub 2} is also used. Solids are followed directly in the simulations and their thermal histories are recorded. Adiabatic expansion creates rapid cooling of the gas, and tail shocks behind the embryo can cause secondary heating events. Radiative transport is efficient, and bow shocks around planetoids can have luminosities ∼few× 10{sup –8} L{sub ☉}. While barred and radial chondrule textures could be produced in the radiative shocks explored here, porphyritic chondrules may only be possible in the adiabatic limit. We present a series of predicted cooling curves that merit investigation in laboratory experiments to determine whether the solids produced by bow shocks are represented in the meteoritic record by chondrules or other solids.

79 ASTRONOMY AND ASTROPHYSICS↗

Bulge-driven fueling of seed black holes

We examine radiation-regulated accretion onto intermediate-mass and massive black holes (BHs) embedded in a bulge component. Using spherically symmetric one-dimensional radiation-hydrodynamics simulations, we track the growth of BHs accreting from a cold, neutral gas reservoir with temperature T{sub ∞}=10{sup 4} K. We find that the accretion rate of BHs embedded in bulges is proportional to r{sub B,eff}/r{sub B}, where r{sub B,eff} is the increased effective Bondi radius that includes the gravitational potential of the bulge, and r{sub B} is the Bondi radius of the BH. The radiative feedback from the BH suppresses the cold accretion rate to ∼1% of the Bondi rate when a bulge is not considered. However, we find that the BH fueling rate increases rapidly when the bulge mass M{sub bulge} is greater than the critical value of ∼10{sup 6} M{sub ☉} and is proportional to r{sub B,eff}/r{sub B}≃M{sub bulge}/M{sub BH}, where M{sub BH} is the BH mass. Since the critical bulge mass is independent of the central BH mass, the growth rate of BHs with masses M{sub BH}=10{sup 2}, 10{sup 4}, and 10{sup 6} M{sub ☉} exhibits distinct dependencies on the bulge-to-BH mass ratio. Our results imply that light seed BHs (≲10{sup 2} M{sub ☉}), which might be the remnants of the Pop III stars, cannot grow through accretion coevally with the early assembly of the bulge of the host galaxies until the bulge reaches the critical mass. However, massive BH seeds (≳10{sup 5} M{sub ☉}), which may form via direct collapse, are more likely to be embedded in a supercritical bulge, and thus can grow efficiently coupling to the host galaxies and driving the early evolution of the M{sub BH}–σ relationship.

79 ASTRONOMY AND ASTROPHYSICS↗

Sensitivity of Simulations of Double-detonation Type Ia Supernovae to Integration Methodology

Abstract We study the coupling of hydrodynamics and reactions in simulations of the double-detonation model for Type Ia supernovae. When assessing the convergence of simulations, the focus is usually on spatial resolution; however, the method of coupling the physics together as well as the tolerances used in integrating a reaction network also play an important role. In this paper, we explore how the choices made in both coupling and integrating the reaction portion of a simulation (operator/Strang splitting versus the simplified spectral deferred corrections method we introduced previously) influences the accuracy, efficiency, and nucleosynthesis of simulations of double detonations. We find no need to limit reaction rates or reduce the simulation time step to the reaction timescale. The entire simulation methodology used here is GPU-accelerated and made freely available as part of the Castro simulation code.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Cold accretion in early galaxy formation and its Lyα signatures

Lyα emission has played an important role in detecting high-redshift galaxies, including recently distant ones at redshifts z>7. It may also contain important information concerning the origin of these galaxies. Here, we investigate the formation of a typical L{sup *} galaxy and its observational signatures at the earliest stage by combining a cosmological hydrodynamic simulation with three-dimensional radiative transfer (RT) calculations using the newly improved ART{sup 2} code. Our cosmological simulation uses the Aquila initial condition, which zooms in on a Milky-Way-like halo with high resolutions, and our RT couples multi-wavelength continuum, Lyα line, and ionization of hydrogen. We find that the modeled galaxy starts to form at redshift z ∼ 24 through the efficient accretion of cold gas, which produces a strong Lyα line with a luminosity of L{sub Lyα}∼10{sup 42} erg s{sup −1} as early as z ∼ 14. The Lyα emission appears to trace the cold, dense gas. The lines exhibit asymmetric, single-peak profiles, and are shifted to the blue wing, a characteristic feature of gas inflow. Moreover, the contribution to the total Lyα luminosity from excitation cooling increases with redshift and becomes dominant at z ≳ 6. We predict that L{sup *} galaxies such as the modeled one may be detected at z ≲ 8 by the James Webb Space Telescope and Atacama Large Millimeter Array with a reasonable integration time. Beyond redshift 12, however, the Lyα line may only be observable by spectroscopic surveys. Our results suggest that the Lyα line is one of the most powerful tools to detect the first generation of galaxies and decipher their formation mechanism.

79 ASTRONOMY AND ASTROPHYSICS↗

Meshfree Process Modeling and Experimental Validation of Friction Riveting of Aluminum 5052 to Aluminum 6061

Friction riveting (Fric-riveting) is an innovative, fast, and energy-efficient process for spot-joining metal-metal structures. Although fric-riveting has been studied experimentally in recent years, its process modeling is rarely found in the literature primarily because of the associated large material deformation, extreme thermomechanical conditions, and complex contact conditions. In this work, a mesh-free smoothed particle hydrodynamics (SPH) framework that can well handle the abovementioned numerical challenges is used to simulate the fric-riveting of AA5052 to AA6061. Predicted material morphology, multi-point temperatures, and plunge force are thoroughly validated by experimental observations. The material severe plastic deformation zone in the vicinity of the riveting zone is further predicted by the SPH model, which indicates the material mixing and potential grain refinement zone. Based on the validated model, process parameters can be optimized which yields better performance over the baseline case.

Friction riveting (fric-riveting), Smoothed partic↗