Search NASA⌕ Search

SEARCH · Search NASA

Results for “HYDRODYNAMIC MODEL”

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 19 records

Smoothed particle hydrodynamics model for Landau-Lifshitz-Navier-Stokes and advection-diffusion equations

We propose a novel smoothed particle hydrodynamics (SPH) discretization of the fully coupled Landau-Lifshitz-Navier-Stokes (LLNS) and stochastic advection-diffusion equations. The accuracy of the SPH solution of the LLNS equations is demonstrated by comparing the scaling of velocity variance and the self-diffusion coefficient with kinetic temperature and particle mass obtained from the SPH simulations and analytical solutions. The spatial covariance of pressure and velocity fluctuations is found to be in a good agreement with theoretical models. To validate the accuracy of the SPH method for coupled LLNS and advection-diffusion equations, we simulate the interface between two miscible fluids. We study formation of the so-called “giant fluctuations” of the front between light and heavy fluids with and without gravity, where the light fluid lies on the top of the heavy fluid. We find that the power spectra of the simulated concentration field are in good agreement with the experiments and analytical solutions. In the absence of gravity, the power spectra decay as the power −4 of the wavenumber—except for small wavenumbers that diverge from this power law behavior due to the effect of finite domain size. Gravity suppresses the fluctuations, resulting in much weaker dependence of the power spectra on the wavenumber. Finally, the model is used to study the effect of thermal fluctuation on the Rayleigh-Taylor instability, an unstable dynamics of the front between a heavy fluid overlaying a light fluid. The front dynamics is shown to agree well with the analytical solutions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Modeling Hydrodynamic Instabilities and Mix in National Ignition Facility Hohlraums

Radiation hydrodynamics modeling of hohlraum experiments on the National Ignition Facility has shown a chronic overprediction of hohlraum performance compared to experiment. Depending on the particular experiment, up to several hundred kilojoules of laser energy must be artificially subtracted from the simulations in order to match experimentally measured quantities like the time of peak neutron production (bang time). Fine-scale hydrodynamic mixing of the high-Z hohlraum wall with the low-Z hohlraum gas, an effect omitted in current hohlraum simulations, has long been hypothesized as possible source of some of this degradation. A preliminary computational assessment using a spherical hohlraum model has been made. It shows that the wall-gas interface is both unstable during certain phases of the experiment and has the potential for turbulent mixing of the gas and wall material, depending on the hohlraum gas fill and pulse shape. The same simulations also indicate that mixing at this interface can delay the simulated bang time by up to 140 ps with high gas fills but appears negligible for low gas fills.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Remote Sensing-Driven Hydrodynamic Modeling in Data-Scarce Regions: Integrating ICESat-2, Sentinel-2, SWOT and Re-analysis Models for Coastal Monitoring

Hydrodynamic models in coastal and estuarine systems are typically constrained by sparse bathymetry, boundary, and validation data, especially in regions where field campaigns are costly or impractical. Here we develop and test a fully satellite-driven framework for hydrodynamic modeling in South Africa’s Langebaan Lagoon without using any local in situ measurements. Bathymetry is derived by training multispectral Sentinel-2 reflectance against ICESat-2 ATL24 photon-derived depths using an XGBoost model optimized with Bayesian search. The final satellite-derived bathymetry reproduces independent ATL24 points with RMSE = 0.45 m and R 2 = 0.97. This bathymetry was used in a depth-averaged Delft3D Flexible Mesh model driven at the open boundary by TPXO tidal harmonics and by ERA5 winds. We validate modeled water surface elevation against 16 SWOT low-rate (250 m, unsmoothed) passes in 2023. SWOT–model comparisons yield an overall RMSE of 0.11 m and R 2 = 0.61, with typical point differences <0.10 m (∼7% of the 1.5 m tidal range), and showed consistent spatial gradients in water level from the offshore boundary, through Saldanha Bay, and into the lagoon. At the offshore boundary, TPXO and SWOT sea surface heights agree closely (R 2 = 0.86). A simple phase adjustment of ∼26,min between TPXO and SWOT lowers the RMSE from 0.18,m to 0.11,m, showing that phase offset accounts for some of the discrepancy, with additional errors likely linked to non-tidal signals. Our results demonstrate that combining passive optical, photon-counting LiDAR, radar interferometry, and global tidal/atmospheric models enables robust, transferrable hydrodynamic modeling in data-scarce coastal systems, offering a cost-effective pathway for monitoring.

ICESat-2↗

Radiation hydrodynamics modeling of the highest compression inertial confinement fusion ignition experiment from the National Ignition Campaign

The recently completed National Ignition Campaign (NIC) on the National Ignition Facility (NIF) showed significant discrepancies between post-shot simulations of implosion performance and experimentally measured performance, particularly in thermonuclear yield. This discrepancy between simulation and observation persisted despite concerted efforts to include all of the known sources of performance degradation within a reasonable two-dimensional (2-D), and even three-dimensional (3-D), simulation model, e.g., using measured surface imperfections and radiation drives adjusted to reproduce observed implosion trajectories [Clark et al., Phys. Plasmas 20, 056318 (2013)]. Since the completion of the NIC, several effects have been identified that could explain these discrepancies and that were omitted in previous simulations. In particular, there is now clear evidence for larger than anticipated long-wavelength radiation drive asymmetries and a larger than expected perturbation seeded by the capsule support tent. This paper describes an updated suite of one-dimensional (1-D), 2-D, and 3-D simulations that include the current best understanding of these effects identified since the NIC, as applied to a specific NIC shot. The relative importance of each effect on the experimental observables is compared. In combination, these effects reduce the simulated-to-measured yield ratio from 125:1 in 1-D to 1.5:1 in 3-D, as compared to 15:1 in the best 2-D simulations published previously. While the agreement with the experimental data remains imperfect, the comparison to the data is significantly improved and suggests that the largest sources for the previous discrepancies between simulation and experiment are now being included.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Development of Kinetic and Radiation Hydrodynamics Modeling of Thermonuclear Burn Propagation in Isochoric p - 11 B Through the Support of the INFUSE Program

The report summarizes DOE INFUSE-supported work between HB11 Energy and the University of Rochester’s TriForce Institute to improve computational modeling of advanced fusion fuels, especially proton–boron-11 (p- 11 B). The project extended the TriForce particle-in-cell/Monte Carlo collision code to include physics needed for dense, high-temperature p- 11 B burn studies, including p- 11 B fusion reactions, three-alpha-particle reaction products, relativistic Coulomb collisions, large-angle nuclear scattering, bremsstrahlung radiation, inverse bremsstrahlung absorption, and photon transport. The upgraded models were verified against focused physics tests and against known deuterium–tritium burn behavior. The study then used one-dimensional spherical simulations to estimate the conditions required for thermonuclear burn propagation in isochoric p- 11 B fuel. The calculations found that burn propagation is possible in the model, but only under very extreme hot-spot conditions, such as about 7000 g/cm 3 at 500 keV or 9000 g/cm 3 at 300 keV for a 20-micron hot spot. These conditions are much more demanding than current demonstrated inertial confinement fusion hot spots. The report concludes that the INFUSE collaboration successfully advanced kinetic and radiation modeling capabilities for p- 11 B fusion and provided useful estimates of ignition requirements. However, the simulated fuel gains remain below what would be needed for practical inertial fusion energy, and further work is needed to reconcile differences among kinetic, radiation-hydrodynamic, and analytic models and to identify more achievable target designs.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Validation of the High-Resolution Salish Sea Tidal Hydrodynamic Model

In this study, a tidal hydrodynamic model was developed and validated to simulate tidal currents in Puget Sound, Washington, to support tidal energy resource characterization using the unstructured-grid, Finite Volume Community Ocean Model (FVCOM). The Salish Sea tidal hydrodynamic model was driven by tides along two open boundaries at the entrance of the Strait of Juan de Fuca and north end of Georgia Strait, and river flows from 19 major rivers in the Salish Sea. To simulate the tidal current in Puget Sound, a high-resolution model grid is required to accurately represent the complex coastlines and bathymetry. The spatial resolution of the model grid varies from ~10 m near river boundaries and ~30 m in small tidal channels and estuaries to near 1000 m inside Georgia Strait and at the open boundaries. Model validation was carried out by comparing simulated and observed water levels at 12 tidal stations and currents at 135 Acoustic Doppler Current Profiler stations in the model domain. A set of model performance metrics, including root mean square error, scatter index, bias, and linear correlation coefficient, were used to quantify the model skills in simulating the tidal hydrodynamics in Puget Sound. Error statistics showed an overall good agreement between simulated and observed tidal elevations and currents, which demonstrated that the Puget Sound tidal model can be used to accurately characterize the tidal stream energy resource in Puget Sound.

16 TIDAL AND WAVE POWER↗

Puget Sound Dissolved Oxygen Modeling Study: Development of an Intermediate-Scale Hydrodynamic Model

The Washington State Department of Ecology contracted with Pacific Northwest National Laboratory to develop an intermediate-scale hydrodynamic and water quality model to study dissolved oxygen and nutrient dynamics in Puget Sound and to help define potential Puget Sound-wide nutrient management strategies and decisions. Specifically, the project is expected to help determine 1) if current and potential future nitrogen loadings from point and non-point sources are significantly impairing water quality at a large scale and 2) what level of nutrient reductions are necessary to reduce or dominate human impacts to dissolved oxygen levels in the sensitive areas. In this study, an intermediate-scale hydrodynamic model of Puget Sound was developed to simulate the hydrodynamics of Puget Sound and the Northwest Straits for the year 2006. The model was constructed using the unstructured Finite Volume Coastal Ocean Model. The overall model grid resolution within Puget Sound in its present configuration is about 880 m. The model was driven by tides, river inflows, and meteorological forcing (wind and net heat flux) and simulated tidal circulations, temperature, and salinity distributions in Puget Sound. The model was validated against observed data of water surface elevation, velocity, temperature, and salinity at various stations within the study domain. Model validation indicated that the model simulates tidal elevations and currents in Puget Sound well and reproduces the general patterns of the temperature and salinity distributions.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Impact of the numerical solution approach of a plant hydrodynamic model (v0.1) on vegetation dynamics

Abstract. Numerous plant hydrodynamic models have started to be implemented in vegetation dynamics models, reflecting the central role of plant hydraulic traits in driving water, energy, and carbon cycles, as well as plant adaptation to climate change. Different numerical approximations of the governing equations of the hydrodynamic models have been documented, but the numerical accuracy of these models and its subsequent effects on the simulated vegetation function and dynamics have rarely been evaluated. Using different numerical solution methods (including implicit and explicit approaches) and vertical discrete grid resolutions, we evaluated the numerical performance of a plant hydrodynamic module in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES-HYDRO version 0.1) based on single-point and global simulations. Our simulation results showed that when near-surface vertical grid spacing is coarsened (grid size >10 cm), the model significantly overestimates aboveground biomass (AGB) in most of the temperate forest locations and underestimates AGB in the boreal forest locations, as compared to a simulation with finer vertical grid spacing. Grid coarsening has a small effect on AGB in the tropical zones of Asia and South America. In particular, coarse surface grid resolution should not be used when there are large and prolonged water content differences among soil layers at depths due to long dry-season duration and/or well-drained soil or when soil evaporation is a dominant fraction of evapotranspiration. Similarly, coarse surface grid resolution should not be used when there is lithologic discontinuity along the soil depth. This information is useful for uncertainty quantification, sensitivity analysis, or the training of surrogate models to design the simulations when computational cost limits the use of ensemble simulations.

54 ENVIRONMENTAL SCIENCES↗

Quantum corrected hydrodynamic models for charge transport in graphene

A quantum hydrodynamic model for charge transport in graphene is derived from a moment expansion of the Wigner–Boltzmann equation. The needed closure relations are obtained by adding quantum corrections based on the equilibrium Wigner function to the semiclassical model formulated in Camiola and Romano (2014) and Luca and Romano (2018) [2,3] by exploiting the Maximum Entropy Principle. The expression of the equilibrium Wigner function which takes into account the energy band of graphene has been obtained by solving the corresponding Bloch equation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Integrating Deep Learning and Hydrodynamic Modeling to Improve the Great Lakes Forecast

The Laurentian Great Lakes, one of the world’s largest surface freshwater systems, pose a modeling challenge in seasonal forecast and climate projection. While physics-based hydrodynamic modeling is a fundamental approach, improving the forecast accuracy remains critical. In recent years, machine learning (ML) has quickly emerged in geoscience applications, but its application to the Great Lakes hydrodynamic prediction is still in its early stages. This work is the first one to explore a deep learning approach to predicting spatiotemporal distributions of the lake surface temperature (LST) in the Great Lakes. Our study shows that the Long Short-Term Memory (LSTM) neural network, trained with the limited data from hypothetical monitoring networks, can provide consistent and robust performance. The LSTM prediction captured the LST spatiotemporal variabilities across the five Great Lakes well, suggesting an effective and efficient way for monitoring network design in assisting the ML-based forecast. Furthermore, we employed an explainable artificial intelligence (XAI) technique named SHapley Additive exPlanations (SHAP) to uncover how the features impact the LSTM prediction. Our XAI analysis shows air temperature is the most influential feature for predicting LST in the trained LSTM. The relatively large bias in the LSTM prediction during the spring and fall was associated with substantial heterogeneity of air temperature during the two seasons. In contrast, the physics-based hydrodynamic model performed better in spring and fall yet exhibited relatively large biases during the summer stratification period. Finally, we developed a statistical integration of the hydrodynamic modeling and deep learning results based on the Best Linear Unbiased Estimator (BLUE). The integration further enhanced prediction accuracy, suggesting its potential for next-generation Great Lakes forecast systems.

Xue, Pengfei (ORCID:000000025702421X)↗

Hydrodynamic models of a Cepheid atmosphere. I - Deep envelope models

The implicit hydrodynamic code of Kutter and Sparks has been modified to include radiative transfer effects. This modified code has been used to compute deep envelope models of a classical Cepheid with a period of 12 days. It is shown that in this particular model the hydrogen ionization region plays only a small role in producing the observed phase lag between the light and velocity curves. The cause of the bumps on the model's light curve is examined, and a mechanism is presented to explain those Cepheids with two secondary features on their light curves. This mechanism is shown to be consistent with the Hertzsprung sequence only if the evolutionary mass-luminosity law is used.

Karp, A. H.↗

Recent Hydrodynamic Modeling Enhancements in OpenFAST

To improve the hydrodynamic modeling capabilities of the offshore wind design and modeling tool OpenFAST, a restructuring of the HydroDyn hydrodynamics module was undertaken with several new features implemented. The generation of the wave field is now separated from HydroDyn into a new module called SeaState. Unlike previous versions of HydroDyn, which precomputed the wave kinematics at the undisplaced positions of the hydrodynamic nodes, the new SeaState module computes the complete time history of the wave field at the vertices of a user-defined wave grid. During the simulation, the wave kinematics at any point within the grid can be efficiently interpolated, allowing the strip-theory wave loads to be evaluated based on the wave kinematics at the displaced structure position. A phase correction was also implemented for the potential-flow wave-exciting loads to account for large structure displacements. Several wave-stretching methods were implemented for the strip-theory solution along with a wave-load redistribution method that ensures smooth variation of the nodal loads as the nodes of a discrete hydrodynamic mesh enter and exit the water. The load smoothing is included to avoid the excitation of unphysical high-frequency structural vibration during hydroelastic simulations. The MacCamy-Fuchs diffraction correction for the strip-theory solution and the capability to insert a constrained NewWave into extreme stochastic sea states with directional spreading were also added. The present article documents the formulations of the new improvements to HydroDyn with example applications and numerical results.

constrained wave↗

Plasmon excitation in metal slab by fast point charge: The role of additional boundary conditions in quantum hydrodynamic model

We study the wake effect in the induced potential and the stopping power due to plasmon excitation in a metal slab by a point charge moving inside the slab. Nonlocal effects in the response of the electron gas in the metal are described by a quantum hydrodynamic model, where the equation of electronic motion contains both a quantum pressure term and a gradient correction from the Bohm quantum potential, resulting in a fourth-order differential equation for the perturbed electron density. Thus, besides using the condition that the normal component of the electron velocity should vanish at the impenetrable boundary of the metal, a consistent inclusion of the gradient correction is shown to introduce two possibilities for an additional boundary condition for the perturbed electron density. We show that using two different sets of boundary conditions only gives rise to differences in the wake potential at large distances behind the charged particle. On the other hand, the gradient correction in the quantum hydrodynamic model is seen to cause a reduction in the depth of the potential well closest to the particle, and a reduction of its stopping power. Even for a particle moving in the center of the slab, we observe nonlocal effects in the induced potential and the stopping power due to reduction of the slab thickness, which arise from the gradient correction in the quantum hydrodynamic model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Assessment of Energy Removal Impacts on Physical Systems: Hydrodynamic Model Domain Expansion and Refinement, and Online Dissemination of Model Results

In this report we describe the 1) the expansion of the PNNL hydrodynamic model domain to include the continental shelf along the coasts of Washington, Oregon, and Vancouver Island; and 2) the approach and progress in developing the online/Internet disseminations of model results and outreach efforts in support of the Puget Sound Operational Forecast System (PS-OPF). Submittal of this report completes the work on Task 2.1.2, Effects of Physical Systems, Subtask 2.1.2.1, Hydrodynamics, for fiscal year 2010 of the Environmental Effects of Marine and Hydrokinetic Energy project.

58 GEOSCIENCES↗

Hydrodynamic model of a dielectric-barrier discharge in pure chlorine

A one-dimensional hydrodynamic model of a dielectric-barrier discharge (DBD) in pure chlorine is developed, and the properties of the discharge are modeled. The discharge is excited in an 8-mm-long discharge gap between 2-mm-thick dielectric quartz layers covering metal electrodes. The DBD spatiotemporal characteristics at gas pressures of 15–100 Torr are modeled for the case in which a 100-kHz harmonic voltage with an amplitude of 8 kV is applied to the electrodes. The average power density deposited in the discharge over one voltage period is 2.5–5.8 W/cm{sup 3}. It is shown that ions and electrons absorb about 95 and 5% of the discharge power, respectively. In this case, from 67 to 97% of the power absorbed by electrons is spent on the dissociation and ionization of Cl{sub 2} molecules. Two phases can be distinguished in the discharge dynamics: the active (multispike) phase, which follows the breakdown of the discharge gap, and the passive phase. The active phase is characterized by the presence of multiple current spikes, a relatively high current, small surface charge density on the dielectrics, and large voltage drop across the discharge gap. The passive phase (with no current spikes) is characterized by a low current, large surface charge density on the dielectrics, and small voltage drop across the discharge gap. The peak current density in the spikes at all pressures is about 4 mA/cm{sup 2}. In the multispike phase, there are distinct space charge sheaths with thicknesses of 1.5–1.8 mm and a mean electron energy of 4.3–7 eV and the central region of quasineutral plasma with a weak electric field and a mean electron energy of 0.8–3 eV. The degree of ionization of chlorine molecules in the discharge is ~0.02% at a pressure of 15 Torr and ~0.01% at 100 Torr. The DBD plasma is electronegative due to the fast attachment of electrons to chlorine atoms: e + Cl{sub 2} → Cl + Cl{sup –}. The most abundant charged particles are Cl{sub 2}{sup +} and Cl{sup −} ions, and the degree of ionization during current spikes in the active phase is (4.1–5.5) × 10{sup –7}. The mechanism of discharge sustainment is analyzed. The appearance of a series of current spikes in the active phase of the discharge is explained.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Recent Hydrodynamic Modeling Enhancements in OpenFAST: Preprint

To improve the hydrodynamic modeling capabilities of the offshore wind design and modeling tool OpenFAST, a restructuring of the HydroDyn hydrodynamics module was undertaken with several new features implemented. The generation of the wave field is now separated from HydroDyn into a new module called SeaState. Unlike previous versions of HydroDyn, which precomputed the wave kinematics at the undisplaced positions of the hydrodynamic nodes, the new SeaState module computes the complete time history of the wave field at the vertices of a user-defined wave grid. During the simulation, the wave kinematics at any point within the grid can be efficiently interpolated, allowing the strip-theory wave loads to be evaluated based on the wave kinematics at the displaced structure position. A phase correction was also implemented for the potential-flow wave-exciting loads to account for large structure displacements. Several wave-stretching methods were implemented for the strip-theory solution along with a wave-load redistribution method that ensures smooth variation of the nodal loads as the nodes of a discrete hydrodynamic mesh enter and exit the water. The load smoothing is included to avoid the excitation of unphysical high-frequency structural vibration during hydroelastic simulations. The MacCamy-Fuchs diffraction correction for the strip-theory solution and the capability to insert a constrained NewWave into extreme stochastic sea states with directional spreading were also added. The present article documents the formulations of the new improvements to HydroDyn with example applications and numerical results.

constrained wave↗