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,783 records · Page 99

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↗

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.↗

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↗

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↗

Small-scale flow topologies in decaying isotropic turbulence laden with finite-size droplets

The topology of the fine-scale motions in decaying isotropic turbulence laden with droplets of super-Kolmogorov size is investigated using results from direct numerical simulations. The invariants of the velocity-gradient, rate-of-strain, and rate-of-rotation tensors are computed in the carrier phase. The joint probability density functions of the invariants are calculated and conditioned on different distances from the droplet surface. The results show that outside the viscous region near the interface, the flow topologies favor stable focus/stretching and unstable node/saddle/saddle structures, which is in agreement with those found in canonical homogeneous isotropic turbulence. Inside the viscous layer at the droplet surface, the flow topologies shift from a preference for high-enstrophy/low-dissipation motions to favoring low-enstrophy/high-dissipation. At the droplet surface, there is a strong tendency for boundary-layer-like and vortex-sheet flow topologies in which the strain and rotation rates are positively correlated. Furthermore, an interesting observation is that the shapes of the invariant distributions at the droplet surface are remarkably similar to those reported in the viscous sublayer of turbulent wall-bounded flows. Also, the results show that the smallest hydrodynamic length scale of the carrier fluid turbulence is located at the droplet interface and that this length scale is two to three times smaller than that of the surrounding bulk flow. From a computational viewpoint, this suggests a more stringent spatial resolution requirement for the direct numerical simulation of finite-size droplets in isotropic turbulence than its single-phase counterpart.

42 ENGINEERING↗

Two-dimensional numerical simulation of a Stirling engine heat exchanger

The first phase of an effort to develop multidimensional models of Stirling engine components is described; the ultimate goal is to model an entire engine working space. More specifically, parallel plate and tubular heat exchanger models with emphasis on the central part of the channel (i.e., ignoring hydrodynamic and thermal end effects) are described. The model assumes: laminar, incompressible flow with constant thermophysical properties. In addition, a constant axial temperature gradient is imposed. The governing equations, describing the model, were solved using Crank-Nicloson finite-difference scheme. Model predictions were compared with analytical solutions for oscillating/reversing flow and heat transfer in order to check numerical accuracy. Excellent agreement was obtained for the model predictions with analytical solutions available for both flow in circular tubes and between parallel plates. Also the heat transfer computational results are in good agreement with the heat transfer analytical results for parallel plates.

Ibrahim, Mounir B.↗

Two-dimensional numerical simulation of a Stirling engine heat exchanger

The first phase of an effort to develop multidimensional models of Stirling engine components is described. The ultimate goal is to model an entire engine working space. Parallel plate and tubular heat exchanger models are described, with emphasis on the central part of the channel (i.e., ignoring hydrodynamic and thermal end effects). The model assumes laminar, incompressible flow with constant thermophysical properties. In addition, a constant axial temperature gradient is imposed. The governing equations describing the model have been solved using the Crack-Nicloson finite-difference scheme. Model predictions are compared with analytical solutions for oscillating/reversing flow and heat transfer in order to check numerical accuracy. Excellent agreement is obtained for flow both in circular tubes and between parallel plates. The computational heat transfer results are in good agreement with the analytical heat transfer results for parallel plates.

Ibrahim, Mounir↗

Coupled Aero-Hydro-Mechanical Hybrid Simulation Testing of Offshore Wind Turbines Subjected to Operational and Extreme Loading Conditions

Understanding the response of the Offshore Wind Turbine (OWT) subjected to realistic applied loads requires modeling the whole structure including its soil-foundation system. This requires unique and innovative testing facilities. OWT systems experience cyclic and dynamic loading due to wind, wave, current, rotor vibrations (i.e., 1P load) and vibrations caused by the blade shadowing effects (2P/3P loads). These loads are complicated in nature and have varying amplitudes, frequencies, and directions. Investigating the response of the entire OWT system including the soil-foundation system under these complex loading conditions, requires: (1) full understanding of the loading characteristics including: the power take-off mechanical load (1P and 3P), and areo- and hydrodynamic loads that the OWT system is subjected to; (2) testing facility with unique multidirectional loading capabilities that allows for simultaneous application of realistic wind, wave and machine loads, axial gravity loads, and induced overturning moments; and (3) unique and cost-effective testing techniques that allow for accurate analysis of the overall response of the OWT system under realistic conditions such as: Real-Time Hybrid Simulation (RTHS).

17 WIND ENERGY↗

Possible Early Linear Acceleration of Proto-neutron Stars via Asymmetric Neutrino Emission in Core-collapse Supernovae

In this Letter we present the result of an axisymmetric core-collapse supernovae simulation conducted with appropriate treatments of neutrino transport and proper motions of proto-neutron stars (PNSs), in which a remarkable PNS acceleration is observed in association with asymmetric neutrino emissions 300 ms after bounce. We find that these asymmetric neutrino emissions play important roles in the acceleration of PNSs in this phase. The correlation between the PNS proper motion and the asymmetric ejecta is similar to that in a neutron star (NS) kick of hydrodynamic origin. Both electron-type neutrinos (ν e ) and their anti-particles (${\bar{\nu }}_{{\rm{e}}}$) have a ~10% level of asymmetry between the northern and southern hemispheres, while other heavy-leptonic neutrinos (ν x ) have much a smaller asymmetry of ~1%. The emissions of ${\bar{\nu }}_{{\rm{e}}}$ and ν x are higher in the hemisphere of stronger shock expansion, whereas the ν e emission is enhanced in the opposite hemisphere: in total, the neutrinos carry some linear momentum to the hemisphere of the stronger shock expansion. This asymmetry is attributed to the non-spherical distribution of electron-fraction (Y e ) in the envelope of PNS. Although it is similar to lepton-emission self-sustained asymmetry, the Y e asymmetry seems to be associated with the PNS motion: the latter triggers lateral circular motions in the envelope of PNS by breaking the symmetry of the matter distribution there, which is then sustained by a combination of convection, lateral neutrino diffusion, and matter-pressure gradient. Our findings may have an influence on the current theories on the NS kick mechanism, although long-term simulations are required to assess their impact on later evolution.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling hydrodynamic and biomass pyrolysis effects of recycled product gases in a bubbling fluidized bed reactor

Fast pyrolysis of biomass in a fluidized bed reactor is typically conducted in a nitrogen gas environment. Recycling product gas can improve the economics of operating such a system by reducing reliance on pure process streams, but much less is known about how recycling pyrolysis product gas may affect fluidization behavior and pyrolysis kinetics. Therefore, gas effects in a fluidized bed biomass pyrolysis reactor were investigated using engineering correlations, low-order models, and CFD simulations for N 2 , H 2 , CO, CO 2 , and CH 4 carrier gas mixtures. Here, our findings reveal viscosity of a gas mixture can be significantly underestimated depending on the model and correlation. Furthermore, fluidization characteristics such as U mf and gas-solid convective heat transfer can be greatly affected by the gas properties. By utilizing H 2 as the fluidizing gas (instead of N 2 ), while maintaining a constant fluidization ratio (U s /U mf ), the bio-oil yields can be increased ~5%. This is due to the lower density H 2 producing similar hydrodynamics as N 2 at higher gas flow rates. These higher flow rates result in shorter gas residence times, and as a result, less secondary reactions that convert bio-oil to light gases and char. Model results also suggest that bio-oil yield is not significantly affected by the type of carrier gas used, with bio-oil yield varying by ~2% across different gas mixtures while maintaining constant flow rate. In conclusion, this indicates that recycled pyrolytic gas can be used as the carrier gas for biomass pyrolysis.

09 BIOMASS FUELS↗

TEAMER: Crossflow Turbine Fairing Geometry Optimization - Report and CFD Modeling Files

The dataset includes computational fluid dynamics (CFD) models and simulation files for crossflow turbines as well as a detailed project report. The report documents the project undertaken by the Ocean Renewable Power Company (ORPC) to design and optimize a modular fairing for the Modular RivGen Marine Hydrokinetic (MHK) turbine, which enhances the efficient deployment and operation of turbine arrays. The project focused on optimizing the hydrodynamic performance of the fairing using CFD, with an emphasis on two key geometric parameters: the fairing's cross-sectional shape and the spacing between the rotor and the fairing. The analysis aimed to maximize net power output while also assessing discretized loading to evaluate ultimate and fatigue loads on the turbine components. The numerical modeling was conducted using both the commercial CFD software Star-CCM+ and the open-source code openFOAM, with the latter utilizing the actuator line library, turbinesFOAM. This dual-code approach was intended to increase confidence in the results and demonstrate the viability of using open-source tools for high-fidelity marine energy modeling. This dataset includes all necessary files for actuator line simulations in openFOAM, as well as 2D blade-resolved CFD results, along with Python and Java scripts for setting up and post-processing simulations.

16 TIDAL AND WAVE POWER↗

Machine learning assisted bayesian inference of mix and hot-spot conditions in NIF implosions

Experiments on the National Ignition Facility (NIF) have provided clear evidence of ablator material mixing into the Hot-Spot, leading to degraded performance. However, inferring the amount of mix and Hot-Spot conditions from typical experimental observations (e.g. x-ray spectra and images) is highly challenging. Here, we have developed an analysis method that utilizes machine learning assisted Bayesian inference to find the probability distributions of the Hot-Spot and mix conditions. This approach uses a neural network, trained on an idealized 2-dimensional representation of the Hot-Spot and mix distribution, and Bayesian inference to find the statistical distributions of Hot-Spot conditions that provide a match with observations. We have tested this method with synthetic data from simulations.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Diffusivity and Structure of Room Temperature Ionic Liquid in Various Organic Solvents

Room-temperature ionic liquids (RTILs) hold promise for applications in electric double layer capacitors (EDLCs), owing to a much wider potential window, lower vapor pressure, and better thermal and chemical stabilities compared to conventional aqueous and organic electrolytes. However, because the low diffusivity of ions in neat RTILs negates the EDLCs’ advantage of high power density, the ionic liquids are often used in mixture with organic solvents. In this study, we measured the diffusivity of cations and anions in RTIL, 1-butyl-3-methylimidazolium bis(trifluoromethylsulfonyl) ([BMIM + ][TFSI – ]), mixed with 10 organic solvents, by using the pulsed-field gradient NMR method. The ion diffusivity was found to follow that of neat solvents and in most studied solvents showed an excellent agreement with the predicted values reported in the recent molecular dynamics (MD) study [Thompson, M. W.; J. Phys. Chem. B 2019, 123, 1340-1347]. In two solvents consisting of long-chain molecules, however, the MD simulations predictions slightly underestimated the ionic diffusivities. The degree of ion dissociation was also estimated for each solvent by comparing the ionic conductivity with the molar conductivity derived from the diffusion measurements. The degree of ion dissociation and the hydrodynamic radius of ions suggest that the ions are coordinated by ~1 solvent molecule. Finally, the scarcity of solvent–ion interactions explains the fact that the diffusivity of ions in the mixture significantly depends on the viscosity of the solvent.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Microscopic Theory of Fluctuating Hydrodynamics in Nonlinear Lattices

The theory of fluctuating hydrodynamics has been an important tool for analyzing macroscopic behavior in nonlinear lattices. However, despite its practical success, its microscopic derivation is still incomplete. In this work, we provide the microscopic derivation of fluctuating hydrodynamics, using the coarse-graining and projection technique; the equivalence of ensembles turns out to be critical. The Green-Kubo (GK)-like formula for the bare transport coefficients are presented in a numerically computable form. Our numerical simulations show that the bare transport coefficients exist for a sufficiently large but finite coarse-graining length in the infinite lattice within the framework of the GK-like formula. Finally, this demonstrates that the bare transport coefficients uniquely exist for each physical system.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Benchmarking numerical methods for impact and cratering applications

Large scale computational models are important for studying impact cratering events that are prevalent both on Earth and, more broadly, in this solar system. To address these problems, models must reliably account for both large length scales (e.g., kilometers) and relatively long time scales (hundreds of seconds). This work benchmarks two such approaches, a more traditional hydrodynamics approach and a finite-discrete element method (FDEM), for impact cratering applications. Both 2D and 3D results are discussed for two different impact velocities, 5 km/s and 20 km/s, striking normal to the target and, for 3D simulations, 45° from vertical. In addition, comparisons to previously published data are presented. Finally, differences in how these methods model damage are discussed. Ultimately, both approaches show successful modeling of several different impact scenarios.

58 GEOSCIENCES↗

Core Physics and Kinetics Calculations for the Fissioning Plasma Core Reactor

Highly efficient, compact nuclear reactors would provide high specific impulse spacecraft propulsion. This analysis and numerical simulation effort has focused on the technical feasibility issues related to the nuclear design characteristics of a novel reactor design. The Fissioning Plasma Core Reactor (FPCR) is a shockwave-driven gaseous-core nuclear reactor, which uses Magneto Hydrodynamic effects to generate electric power to be used for propulsion. The nuclear design of the system depends on two major calculations: core physics calculations and kinetics calculations. Presently, core physics calculations have concentrated on the use of the MCNP4C code. However, initial results from other codes such as COMBINE/VENTURE and SCALE4a. are also shown. Several significant modifications were made to the ISR-developed QCALC1 kinetics analysis code. These modifications include testing the state of the core materials, an improvement to the calculation of the material properties of the core, the addition of an adiabatic core temperature model and improvement of the first order reactivity correction model. The accuracy of these modifications has been verified, and the accuracy of the point-core kinetics model used by the QCALC1 code has also been validated. Previously calculated kinetics results for the FPCR were described in the ISR report, "QCALC1: A code for FPCR Kinetics Model Feasibility Analysis" dated June 1, 2002.

Butler, C.↗

Investigation of the Role of the Role of Nuclear Physics in Heavy Element Nucleosynthesis, through the Study of Key Reactions, and the Improvement of Theoretical Reaction Rates (Final Report)

Neutrino-driven winds in core-collapse supernovae have been identified as a possible site for the production of elements heavier than iron. Traditionally, these neutrino-driven winds have been proposed as the site of the main r-process. Recent simulations fail to reproduce the conditions required for the main r-process. while they remain a promising site for producing the lightest elements beyond iron, e.g., Sr, Y, and Zr through the νp process. The efficiency of the νp process depends on the hydrodynamical conditions, the electron fraction (which is related to the neutrino properties), and the nuclear reactions on many short-lived nuclei with limited (if at all) experimental information. The reaction rates on these nuclei are based on theoretical predictions using the Hauser-Feshbach model. Recent sensitivity studies have highlighted the importance of neutron-induced reactions on these nuclei along the νp process path. This work aimed to experimentally constrain reaction rates that are known to play a key role in the neutrino-p process nucleosynthesis. A secondary subsequently-added objective was to start the implementation of techniques that improve the description of nuclear properties in the Hauser-Feshbach model by extending the microscopic nuclear level density description offered via the shell model to high excitation energies without using experiment-based renormalizations. The main objective of this work was the experimental constraint of the 56 Ni(n,p) 56 Co reaction rate via a measurement of the inverse reaction 56 Co(p,n) 56 Ni at the National Superconducting Cyclotron Laboratory (NSCL) and later the Facility for Rare Isotope Beams (FRIB). This reaction is considered the key one for determining the yields possible by the neutrino-p process. A technique for this type of measurement in inverse kinematics at low energies did not exist before this work. The work also had two secondary objectives. First, to contribute to efforts to measure the same reaction in direct kinematics using a radioactive target at Los Alamos National Laboratory (LANL), and second, to advance work to implement shell-model-deduced microscopic level densities in Hauser-Feshbach calculations. The project has resulted in the development of the first technique to perform (p,n) cross-section measurements in relevant-for-astrophysics low energies in inverse kinematics using a magnetic spectrometer or separator, and neutron detectors for neutron tagging. It has also resulted in the precise measurement of the cross-section of the 40 Ar(p,n) 40 K reaction in a proof-of-principle experiment realized by using a beam-line quadrupole of the ReA3 accelerator of NSCL/FRIB. As part of this project the technique was successfully adapted to make use of the superior acceptance of the Separator for Capture Reactions (SECAR) at FRIB. In this project, the required experimental setup simulations and beam optics were developed and tested with the measurement of the 58 Fe(p,n) 58 Cu reaction cross-section. Additionally, this project contributed with simulation work to the development of a technique to measure (n,p) reactions with radioactive targets at LANL, and which resulted in the measurement of the key 56 Ni(n,p) 56 Co reaction cross-section at neutron energies above ≈1 MeV. Last, the project initiated work in the development of shell model based level densities using the moments method.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

OC6 Phase Ia: CFD Simulations of the Free-Decay Motion of the DeepCwind Semisubmersible

Currently, the design of floating offshore wind systems is primarily based on mid-fidelity models with empirical drag forces. The tuning of the model coefficients requires data from either experiments or high-fidelity simulations. As part of the OC6 (Offshore Code Comparison Collaboration, Continued, with Correlation, and unCertainty (OC6) is a project under the International Energy Agency Wind Task 30 framework) project, the present investigation explores the latter option. A verification and validation study of computational fluid dynamics (CFD) models of the DeepCwind semisubmersible undergoing free-decay motion is performed. Several institutions provided CFD results for validation against the OC6 experimental campaign. The objective is to evaluate whether the CFD setups of the participants can provide valid estimates of the hydrodynamic damping coefficients needed by mid-fidelity models. The linear and quadratic damping coefficients and the equivalent damping ratio are chosen as metrics for validation. Large numerical uncertainties are estimated for the linear and quadratic damping coefficients; however, the equivalent damping ratios are more consistently predicted with lower uncertainty. Some difference is observed between the experimental and CFD surge-decay motion, which is caused by mechanical damping not considered in the simulations that likely originated from the mooring setup, including a Coulomb-friction-type force. Overall, the simulations and the experiment show reasonable agreement, thus demonstrating the feasibility of using CFD simulations to tune mid-fidelity models.

17 WIND ENERGY↗