Search NASA⌕ Search

SEARCH · Search NASA

Results for “hydrodynamics methods”

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 955 records · Page 53

The Role of Atmospheric Rivers on Groundwater: Lessons Learned From an Extreme Wet Year

Abstract In the coastal regions of the western United States, atmospheric rivers (ARs) are associated with the largest precipitation generating storms and contribute up to half of annual precipitation, but the impact of ARs on the integrated hydrologic cycle, specifically on groundwater storage and hydrodynamics, is largely unknown. To better explore the hydrologic behavior of AR versus non‐AR event precipitation, we present a novel combination of two water tracking methods (one in the atmosphere and one in the subsurface) to explicitly track the full lifecycle of water parcels generated by ARs. Simulations of northern California's Cosumnes River watershed during the record wet 2017 water year are performed via the coupling of a high‐resolution regional climate model and a land surface‐groundwater model accounting for lateral groundwater flow. Despite ARs contributing more precipitation than non‐AR storms, we find less AR water is preferentially stored in aquifers by year end. Fractionally, ARs result in 300% less snow derived groundwater‐recharged compared to non‐AR precipitation. Rain‐on‐snow (RoS) plays an important role in AR‐driven discharge, where over 50% of total discharge from ARs snow is from RoS events. Finally, despite record‐breaking annual precipitation, simulated groundwater depletion occurs by year end due to estimates of groundwater pumping activities. The results from these simulations serve as a partial analogue of future hydrologic conditions where ARs are expected to intensify and provide a greater fraction of annual precipitation due to climate change.

54 ENVIRONMENTAL SCIENCES↗

Traveling waves in Hall-magnetohydrodynamics and the ion-acoustic shock structure

Hall-magnetohydrodynamics (HMHD) is a mixed hyperbolic-parabolic partial differential equation that describes the dynamics of an ideal two fluid plasma with massless electrons. We study the only shock wave family that exists in this system (the other discontinuities being contact discontinuities and not shocks). We study planar traveling wave solutions and we find solutions with discontinuities in the hydrodynamic variables, which arise due to the presence of real characteristics in Hall-MHD. We introduce a small viscosity into the equations and use the method of matched asymptotic expansions to show that solutions with a discontinuity satisfying the Rankine-Hugoniot conditions and also an entropy condition have continuous shock structures. The lowest order inner equations reduce to the compressible Navier-Stokes equations, plus an equation which implies the constancy of the magnetic field inside the shock structure. We are able to show that the current is discontinuous across the shock, even as the magnetic field is continuous, and that the lowest order outer equations, which are the equations for traveling waves in inviscid Hall-MHD, are exactly integrable. We show that the inner and outer solutions match, which allows us to construct a family of uniformly valid continuous composite solutions that become discontinuous when the diffusivity vanishes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Time-dependent density-functional-theory calculations of the nonlocal electron stopping range for inertial confinement fusion applications

Nonlocal electron transport is important for understanding laser-target coupling for laser-direct-drive (LDD) inertial confinement fusion (ICF) simulations. Current models for the nonlocal electron mean free path in radiation-hydrodynamic codes are based on plasma-physics models developed decades ago; improvements are needed to accurately predict the electron conduction in LDD simulations of ICF target implosions. Here we utilized time-dependent density functional theory (TD-DFT) to calculate the electron stopping power (SP) in the so-called conduction-zone plasmas of polystyrene in a wide range of densities and temperatures relevant to LDD. Compared with the modified Lee-More model, the TD-DFT calculations indicated a lower SP and a higher stopping range for nonlocal electrons. We fit these electron SP calculations to obtain a global analytical model for the electron stopping range as a function of plasma conditions and the nonlocal electron kinetic energy. This model was implemented in the one-dimensional radiation-hydrodynamic code LILAC to perform simulations of LDD ICF implosions, which are further compared with simulations by the standard modified Lee-More model. In conclusion, results from these integrated simulations are discussed in terms of the implications of this TD-DFT-based mean-free-path model to ICF simulations.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Three-dimensional rotational plasma flows near solid surfaces in an axial magnetic field

A rotational flow of a conducting viscous medium near an extended dielectric disk in a uniform axial magnetic field is analyzed in the magnetohydrodynamic (MHD) approach. An analytical solution to the system of nonlinear differential MHD equations of motion in the boundary layer for the general case of different rotation velocities of the disk and medium is obtained using a modified Slezkin–Targ method. A particular case of a medium rotating near a stationary disk imitating the end surface of a laboratory device is considered. The characteristics of a hydrodynamic flow near the disk surface are calculated within the model of a finite-thickness boundary layer. The influence of the magnetic field on the intensity of the secondary flow is studied. Calculations are performed for a weakly ionized dense plasma flow without allowance for the Hall effect and plasma compressibility. An MHD flow in a rotating cylinder bounded from above by a retarding cap is considered. The results obtained can be used to estimate the influence of the end surfaces on the main azimuthal flow, as well as the intensities of circulating flows in various devices with rotating plasmas, in particular, in plasma centrifuges and laboratory devices designed to study instabilities of rotating plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Evolution of Mass and Velocity Field in the Cosmic Web: Comparison Between Baryonic and Dark Matter

We investigate the evolution of the cosmic web since z = 5 in grid-based cosmological hydrodynamical simulations, focusing on the mass and velocity fields of both baryonic and cold dark matter. The tidal tensor of density is used as the main method for web identification, with λ {sub th} = 0.2–1.2. The evolution trends in baryonic and dark matter are similar, although moderate differences are observed. Sheets appear early, and their large-scale pattern may have been set up by z = 3. In terms of mass, filaments supersede sheets as the primary collapsing structures from z ∼ 2–3. Tenuous filaments assembled with each other to form prominent ones at z < 2. In accordance with the construction of the frame of the sheets, the cosmic divergence velocity, v {sub div}, was already well-developed above 2–3 Mpc by z = 3. Afterwards, the curl velocity, v {sub curl}, grew dramatically along with the rising of filaments, becoming comparable to v {sub div}, for <2–3 Mpc at z = 0. The scaling of v {sub curl} can be described by the hierarchical turbulence model. The alignment between the vorticity and the eigenvectors of the shear tensor in the baryonic matter field resembles that in the dark matter field, and is even moderately stronger between ω and e{sub 1}, and ω and e{sub 3}. Compared with dark matter, there is slightly less baryonic matter found residing in filaments and clusters, and its vorticity developed more significantly below 2–3 Mpc. These differences may be underestimated because of the limited resolution and lack of star formation in our simulation. The impact of the change of dominant structures in overdense regions at z ∼ 2–3 on galaxy formation and evolution is shortly discussed.

79 ASTRONOMY AND ASTROPHYSICS↗

Computation of turbulent flows

Substantial advances made over the past decade in the prediction of turbulent flows are discussed. There was extensive work in the development of turbulence models, particularly for use in boundary layer calculations. Basic aspects of several important methods based on partial differential equations for the mean velocity field and turbulence quantities, including the relationship between the methods and suggestions for future development were reviewed. Work on three-dimensional time-dependent large eddy simulations is discussed. The emphasis is on the hydrodynamics of incompressible flows, but sources for consideration of heat transfer and compressibility are mentioned.

Reynolds, W. C.↗

A Fully Nonlinear, Dynamically Consistent Numerical Model for Solid-Body Ship Motion. I. Ship Motion with Fixed Heading

In this paper, we describe the details of our numerical model for simulating ship solidbody motion in a given environment. In this model, the fully nonlinear dynamical equations governing the time-varying solid-body ship motion under the forces arising from ship wave interactions are solved with given initial conditions. The net force and moment (torque) on the ship body are directly calculated via integration of the hydrodynamic pressure over the wetted surface and the buoyancy effect from the underwater volume of the actual ship hull with a hybrid finite-difference/finite-element method. Neither empirical nor free parametrization is introduced in this model, i.e. no a priori experimental data are needed for modelling. This model is benchmarked with many experiments of various ship hulls for heave, roll and pitch motion. In addition to the benchmark cases, numerical experiments are also carried out for strongly nonlinear ship motion with a fixed heading. These new cases demonstrate clearly the importance of nonlinearities in ship motion modelling.

Lin, Ray-Quing↗

Trigger of Standard and Blowout Solar Jets: I. Parametric Studies of the Magnetic Field Geometry

Context. Jets are dynamic, impulsive, well collimated plasma events developing at many different scales and in differentlayers of the solar atmosphere.Aims. Jets are believed to be induced by magnetic reconnection a process central to many astrophysical phenomena.Studying their dynamics can help us to better understand the processes acting in larger eruptive events (e.g. flares andcoronal mass ejection) as well as mass, magnetic helicity and energy transfer at all scale in the solar atmosphere. Therelative simplicity of their magnetic geometry and topology, compared with larger solar active events, makes jets idealcandidate for studying the fundamental role of reconnection in energetic events. Methods. In this study, using our state-of-the-art numerical solver ARMS, we present several parametric studies of anumerical three-dimensional magneto-hydrodynamic model of solar jet-like events. We study the impact of the magneticfield inclination and photospheric field distribution on the generation and properties of two types of observed solar jets:the so-called standard and blowout jets.Results. The present parametric studies validate our model of jets for different geometric properties of the magneticconfiguration. On one hand, a blowout jets is eventually always trigger for the range of parameters that we tested. The3D magnetic null-point configuration is a very robust structure for the energy storage and impulsive release characteristicof blowout jets. On the other hand, the existence of the standard jet depends on the magnetic geometry. We show thatthe reconnection occurring during the standard jet phase influences the properties of the trigger of the blowout jet.Conclusions. Our results allow us to better understand the energization, triggering and driving processes of standardand blowout jets. Our model allows us to predict the impulsiveness and energetics of jets in terms of the surroundingmagnetic field configuration. Finally we discuss the interpretation of standard jets and blowout jets and the physicalfactors that determine which type of jets will occur.

Pariat, E.↗

Surface spin-electron acoustic waves in magnetically ordered metals

Degenerate plasmas with motionless ions show existence of three surface waves: the Langmuir wave, the electromagnetic wave, and the zeroth sound. Applying the separated spin evolution quantum hydrodynamics to half-space plasma, we demonstrate the existence of the surface spin-electron acoustic wave (SSEAW). We study dispersion of the SSEAW. We show that there is hybridization between the surface Langmuir wave and the SSEAW at rather small spin polarization. In the hybridization area, the dispersion branches are located close to each other. In this area, there is a strong interaction between these waves leading to the energy exchange. Consequently, generating the Langmuir waves with the frequencies close to hybridization area we can generate the SSEAWs. Thus, we report a method of creation of the spin-electron acoustic waves.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Efficiently resolving the terrestrial-aquatic interface in E3SM with sub-grid methods to improve coastal simulations (Final Technical Report)

The broad goal of this project was to extend the capabilities of MPAS‐Ocean, the ocean component of DOE’s Energy Exascale Earth System Model (E3SM), to better resolve the hydrodynamics of the terrestrial‐aquatic interface without significantly increasing computing resource requirements. Specifically, our goal was to better resolve the small-scale features across the terrestrial‐aquatic interface which includes the low lying coastal floodplain and a dendritic network of hydraulically efficient channels connected to estuarine systems and back bays, barrier islands and high energy inlets. The processes driving flow across the terrestrial‐aquatic interface include upland and floodplain hydrology, tides, wind and atmospheric pressure (including coastally generated storm surge across continental shelves and within estuaries and bays), wind waves and water temperature and salinity gradients.

58 GEOSCIENCES↗

Hydrodynamic Analysis of a Suspended Cylinder Under Regular Wave Loading Based on Computational Fluid Dynamics

An investigation into the computation of hydrodynamic loads on a suspended cylinder in regular waves is presented. The primary goal was to perform a three-way validation of the loads between experimental measurements and simulations from two computational methods. Experimental measurements of the longitudinal in-line force on a cylinder suspended at a fixed position were available from the Offshore Code Comparison Collaboration, Continued, with Correlation (OC5) project, Phase Ia. These measurements were compared to computational fluid dynamics (CFD) simulations based on the solution of Reynolds-averaged Navier-Stokes (RANS) equations, as implemented in STAR-CCM+. The study encompassed a sensitivity analysis of the loads computed in STAR-CCM+ based on wave modeling, boundary conditions, turbulence modeling, and spatial and temporal discretization. The analysis was supplemented by results generated with the offshore wind turbine engineering software OpenFAST, based on a hybrid combination of second-order potential flow and viscous drag from Morison's equation. The focus of the investigation was on the assessment of the accuracy of the computation of first- and higher-order hydrodynamic loads. Substantial differences were observed in the numerical prediction of the second and third harmonic force contribution. Local flow field analysis with CFD was applied to study the physics of wave run-up and diffraction dynamics to identify the causes.

16 TIDAL AND WAVE POWER↗

Hydrodynamic Analysis of a Suspended Cylinder Under Regular Wave Loading Based on Computational Fluid Dynamics: Preprint

An investigation into the computation of hydrodynamic loads on a suspended cylinder in regular waves is presented. The primary goal was to perform a three-way-validation of the loads between experimental measurements and simulations from two computational methods. Experimental measurements of the longitudinal inline force on a cylinder suspended at a fixed po- sition were available from the Offshore Code Comparison Col- laboration, Continued, with Correlation (OC5) project, Phase Ia. These measurements were compared to Computational Fluid Dynamics (CFD) simulations based on the solution of Reynolds- averaged Navier-Stokes (RANS) equations, as implemented in STAR-CCM plus. The study encompassed a sensitivity analysis of the loads computed in STAR-CCM plus based on wave modeling, boundary conditions, turbulence modeling, and spatial and tem- poral discretization. The analysis was supplemented by results generated with the offshore wind turbine engineering software OpenFAST, based on a hybrid combination of second-order po- tential flow and viscous drag fromMorison's equation. The focus of the investigation was on the assessment of the accuracy of the computation of first- and higher-order hydrodynamic loads. Sub- stantial differences were observed in the numerical prediction of the second and third harmonic force contribution. Local flow field analysis with CFD was applied to study the physics of wave run-up and diffraction dynamics to identify the causes.

16 TIDAL AND WAVE POWER↗

Bottomonium suppression in an open quantum system using the quantum trajectories method

We solve the Lindblad equation describing the Brownian motion of a Coulombic heavy quark-antiquark pair in a strongly coupled quark-gluon plasma using the highly efficient Monte Carlo wave-function method. The Lindblad equation has been derived in the framework of pNRQCD and fully accounts for the quantum and non-Abelian nature of the system. The hydrodynamics of the plasma is realistically implemented through a 3+1D dissipative hydrodynamics code. We compute the bottomonium nuclear modification factor and compare with the most recent LHC data. The computation does not rely on any free parameter, as it depends on two transport coefficients that have been evaluated independently in lattice QCD. Our final results, which include late-time feed down of excited states, agree well with the available data from LHC 5.02 TeV PbPb collisions.

79 ASTRONOMY AND ASTROPHYSICS↗

Oxygen reduction reaction catalyzed by platinum nanonetwork prepared by template free one step synthesis for polymer electrolyte membrane fuel cells

Highlights: • Supportless Pt nanonetwork (Pt NN) synthesized by novel template free one step method as per our earlier reported procedure. • Electrocatalytic activity of Pt NN studied taking oxygen reduction reaction in acid medium. • Kinetic and thermodynamic parameters were deduced under hydrodynamic conditions. • ORR mechanistic pathway was proposed based on kinetic rate constants. • ADT analysis found enhanced stability (5000 cycles) for Pt NN than Pt NN/VC and reported Pt/C. - Abstract: The reduction reaction of molecular oxygen (ORR) was investigated using supportless Pt nanonetwork (Pt NN) electrocatalyst in sulfuric acid medium. Pt NN was prepared by template free borohydride reduction. The transmission electron microscope images revealed a network like nano-architecture having an average cluster size of 30 nm. The electrochemical characterization of supportless and Vulcan carbon supported Pt NN (Pt NN/VC) was carried out using rotating disc and ring disc electrodes at various temperatures. Kinetic and thermodynamic parameters were estimated under hydrodynamic conditions and compared with Pt NN/VC and reported Pt/C catalysts. The accelerated durability test revealed that supportless Pt NN is quite stable for 5000 potential cycles with 22% reduction in electrochemical surface area (ECSA). While the initial limiting current density has in fact increased by 11.6%, whereas Pt NN/VC suffered nearly 55% loss in ECSA and 13% loss in limiting current density confirming an enhanced stability of supportless Pt NN morphology for ORR compared to conventional Pt/C ORR catalysts in acid medium.

36 MATERIALS SCIENCE↗

Sagdeev potential approach for quantum ion-acoustic solitary waves in an electron-positron-ion plasma

Sagdeev pseudopotential method is employed to study the arbitrary amplitude quantum ion-acoustic solitary waves in an unmagnetized electron-positron-ion plasma by using one dimensional quantum hydrodynamic model together with the Poisson equation. Sagdeev potential function is obtained in terms of electrostatic potential and analyzed with and without the effect of quantum diffraction parameter H. Effects of the parameter H on both the amplitude and width of the solitary waves have been observed. It is also observed that the positron density can affect the wave propagation.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

X-ray heating and electron temperature of laboratory photoionized plasmas

We discuss experimental and modeling results for the x-ray heating and temperature of laboratory photoionized plasmas. A method is used to extract electron temperature based on the analysis of transmission spectroscopy data that is independent of atomic kinetics modeling. Boltzmann kinetic and radiation-hydrodynamics simulations emphasized the critical interplay between atomic physics and heating in photoionized plasmas. They also demonstrated the dramatic impact of photoexcitation on excited state populations, line emissivity, and radiation cooling. Temperature calculations performed with astrophysical codes significantly overestimated the measurement.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Finite element solution algorithm for incompressible fluid dynamics

A finite element solution algorithm is established for the two-dimensional Navier-Stokes equations governing the transient motion of a viscous incompressible fluid, i.e., hydrodynamics. Dependent variable transformation renders the differential equation description uniformly elliptic. The finite element algorithm is established using the Galerkin criterion on a local basis within the Method of Weighted Residuals. It is unconstrained with respect to system linearity, computational mesh uniformity or solution domain closure regularity. The finite element matrices are established using a linear 'natural coordinate function' description. Computational solutions using the COMOC computer program illustrate the various features of the algorithm including recirculating flows.

Baker, A. J.↗

Theoretical interpretation of Traveling Interplanetary Phenomena and their solar origins

A review is made of recent theoretical studies on Traveling Interplanetary Phenomena (TIP) and their possible solar origins in the context of mathematical methods and physical processes. Flare-generated interplanetary shocks are used as examples to illustrate the following set of contrasting approaches to the phenomena: analytical vs numerical, hydrodynamic vs magnetohydrodynamic, dissipative vs non-dissipative, continuum/macroscopic vs kinetic/microscopic. It is suggested in conclusion that future research should emphasize continuum/macroscopic physical processes together with the kinetic approach to deepen understanding of microscopic interactions.

Wu, S. T.↗