Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamics modeling”

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 631 records · Page 35

Quantitative radiography for determining density fluctuations in HED experiments

We have developed a method to extract density fluctuation measurements from x-ray radiographs of high-energy density (HED) instability growth and turbulence experiments. We use this information to calculate density fluctuation statistics for constraining the performance of turbulent mix models in HED systems. The density calculation combines image filtering, removal of systemic effects such as backlighter variation, calculation of transmission across multiple materials, and use of tracer materials to generate an approximate single-material density field. From the density map, we calculate both average density and a variance-like moment b (density-specific-volume covariance), which we compare to our models. We infer both quantities from a single image, which is significantly more information than the historic single scalar mix width measurements. We also develop a method of analyzing simulation outputs that incorporate both the density fluctuation metric from a turbulence model and the bulk material maps from the hydrodynamic code. This analysis helps address the question of how to initialize the simulations for best comparison to data from systems with large separations of scale in the mixing perturbation initial condition. We find that our data analysis method yields 1D average density and b curves with similar morphology and amplitudes as those from preliminary simulation comparisons.

47 OTHER INSTRUMENTATION↗

WEC-Sim Modeling of Laminar Scientific's Patented Seesaw Wave Energy Converter: Preprint

Laminar Scientific's patented seesaw wave energy converter was modeled in WEC-Sim to predict performance. The device operates by utilizing ocean surface waves to rotate a truss in pitch about a pivot. The pivot is located at the top of two pylons, which are embedded in the seafloor. The seesaw has a float on either end, and the buoyancy forces from each float cause the system to rise or fall with passing waves. Device performance relies upon seesaw length and ocean wavelength creating an antiphase effect. The seesaw truss has an adjustable length intended to achieve this effect. The operation method enforces a narrow band of wavelengths which induce the largest rotational motion from the device. The hydrodynamic analysis of the device was performed using Capytaine, and the results confirmed that the device operates best in a narrow frequency band. Four float-to-float spacing cases and three pylon radii were examined. The hydrodynamic results indicate a match between the model and the physical expectations for the device, and that varying the pylon radii by 0.1-m increments for three instances creates minimal changes in hydrodynamic properties. Power matrices for three float spacing cases of the device were simulated with Joint North Sea Wave Project spectra waves and optimal power take-off damping in WEC-Sim. The maximum average power production for the 15-m spacing case was 14.1 kW with a 5.0-s peak wave period and 4-m significant wave heights. Plots of capture-width ratios indicated that the device performance was linear and confirmed that the device is optimal in a narrow frequency band. The maximum percentage of the available wave power produced by the 15-m device was approximately 16%. Simulations of the device in regular waves were used to produce plots of average power compared to a ratio of float spacing to wavelength. These plots indicate that the power production is maximized at a ratio of 0.5, and further confirm that the device has a narrow frequency response. The device was simulated at an example field location, where the device produced an annual average power rating of 1.6 kW given an average omnidirectional wave climate of 10.3 kW m-1 and an optimal, linearized power take-off model. While the maximum predicted device performance is reliant upon a narrow band of wave frequencies, the conducted analysis provides an opportunity to improve device design prior to prototyping and testing. Modifying the design to respond to a broader frequency range would improve device performance.

numerical modeling↗

The Progenitors of Superluminous Type Ia Supernovae

Abstract Recent observations of type Ia supernovae (SNe Ia) have discovered a subclass of “super-Chandrasekhar” SNe Ia (SC SNe Ia) whose high luminosities and low ejecta velocities suggest that they originate from the explosions of white dwarfs (WDs) with masses that exceed the Chandrasekhar mass limit. Different models have been proposed to explain the progenitors of these explosions, including a “magnetized WD” model and a “WD merger” model. To test the robustness of these models, we conduct a 1D numerical parameter survey of WD explosions using these models as initial conditions. We follow the explosions using the hydrodynamics code Castro and then use the radiation transport code SuperNu to create light curves and spectra for the models. We find that while both classes of models fall within the range of SC SNe Ia observations on the light-curve width–luminosity relation, only the WD merger models reproduce the observed low ejecta velocities. The light curves of our merger models are more similar photometrically to observations than our magnetized models. Given this, we discuss possible explanations for the brightest SC SNe Ia observations that cannot be reproduced with our WD merger models. This study provides the basis for future SC SNe Ia observations and higher-dimensional numerical models.

79 ASTRONOMY AND ASTROPHYSICS↗

A&L Annual Report: Kapton Strength

Exploding Foil Initiator (EFI) systems function by launching a polymer flyer (often referred to as a slapper) at a high explosive (HE) pellet, which is intended to shock-initiate the HE. Flyer characteristics (e.g., shape, planarity, etc.) are extremely important in this respect, governing their ability to initiate the HE. Polymer models currently used in these systems are largely unvalidated, potentially compromising the integrity of magneto-hydrodynamic (MHD) predictions of EFI function. An improved EFI MHD modelling capability is expected to: (1) expedite optimization of EFI-based initiation system design, (2) enable consideration of EFI-related aging and corresponding lifetime predictions, and (3) potentially minimize the number of costly experiments required to certify EFI designs. To address this knowledge gap, several years ago we began a computational project to develop a higher-fidelity model for Kapton, a Polyimide copolymer most commonly employed as the flyer material in EFI devices.

36 MATERIALS SCIENCE↗

Correlations in the Initial Conditions of Heavy-Ion Collisions

Ultracentral collisions of heavy nuclei, in which the impact parameter is nearly zero, are especially sensitive to the details of the initial state model and the microscopic mechanism for collective flow. In a hydrodynamic “flow” picture, the final state momentum correlations are a direct response to the fluctuating initial geometry, although models of the initial geometry differ widely. Alternatively, dynamical mechanisms based in the color glass condensate (CGC) formalism can naturally lead to many-body correlations with very different systematics. Here we present a calculation of event-by-event elliptic flow in both the hydrodynamic and CGC paradigms and show that they can be qualitatively distinguished in ultracentral collisions of deformed nuclei. Specifically, the multiplicity dependence in such collisions is qualitatively opposite, with the CGC correlations increasing with multiplicity while the hydrodynamic correlations decrease. The consistency of the latter with experimental data on UU collisions appears to rule out a CGC-mediated explanation. We find that these qualitative features also persist in small deformed systems and can therefore be a valuable test of the microscopic physics in that regime.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Parallel Implicit Hydrodynamics for High Explosive Burn Calculations

High explosives in hostile environments will require calculational capabilities that model processes, which evolve on timescales from minutes to nanoseconds. Eventually the HE will begin to move metal. To handle this temporal evolution an implicit hydrodynamics coupled to the chemical release of the HE energy is required. In addition, the use of chemical kinetics to model the transition from, the initially, slow heating of a confined high explosive through to deflagration and on to detonation requires many computational zones to model high explosive engineering systems. This requirement means that a fully parallel implicit hydrodynamics is essential. In this paper we present the calculation of a nonlinear matrix equation for the advanced particle pressure that has been made parallel and implemented in our AMR code, BABBO. This new parallel implicit hydrodynamics has been applied to a cookoff problem, as well as, one and two dimensional shock problems. Results are presented and discussed.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Hydrodynamic theory of scrambling in chaotic long-range interacting systems

The Fisher-Kolmogorov-Petrovsky-Piskunov (FKPP) equation provides a mean-field theory of out-of-time-ordered commutators in locally interacting quantum chaotic systems at high energy density. In systems with power-law interactions, the corresponding fractional-derivative FKPP equation provides an analogous mean-field theory. However, the fractional FKPP description is potentially subject to strong quantum fluctuation effects, so it is not clear a priori if it provides a suitable effective description for generic chaotic systems with power-law interactions. Here, in this work, we study this problem using a model of coupled quantum dots with interactions decaying as 1/r α , where each dot hosts N degrees of freedom. The large-N limit corresponds to the mean-field description, while quantum fluctuations contributing to the OTOC can be modeled by 1/N corrections consisting of a cutoff function and noise. Within this framework, we show that the parameters of the effective theory can be chosen to reproduce the butterfly light cone scalings previously found for N=1 and generic finite N. In order to reproduce these scalings, the fractional index μ in the FKPP equation needs to be shifted from the naïve value of μ=2⁢α–1 to a renormalized value μ=2⁢α–2. We provide supporting analytic evidence for the cutoff model and numerical confirmation for the full fractional FKPP equation with cutoff and noise.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Development, verification, and validation of comprehensive acoustic fluid-structure interaction capabilities in an open-source computational platform

The acoustic fluid-structure interaction (FSI) formulation is a practical numerical approach for the seismic analysis of fluid-filled tanks. However, there are no verification and validation studies reported in the literature that demonstrate the ability of an acoustic FSI numerical model to predict responses important to structural and mechanical design for intense translational and rotational earthquake inputs. Herein, an acoustic FSI formulation is implemented in the open-source Multiphysics Object-Oriented Simulation Environment (MOOSE), and is formally verified and validated using analytical solutions and code-to-code verification, and experimental data, respectively. The analytical solutions are for small amplitude, unidirectional seismic inputs. The code-to-code verification utilizes a previously verified and validated Arbitrary Lagrangian-Eulerian (ALE) numerical model in the commercial finite element code LS-DYNA. The validation studies utilize a comprehensive data set assembled from results of 3D earthquake-simulator tests of a fluid-filled vessel. The acoustic numerical model in MOOSE is verified and validated for hydrodynamic pressures and support reactions except for cases that involve significant convective response. For small amplitude inputs, numerically predicted wave heights match those of the analytical solutions. The numerical model is not verified and validated for wave height calculations under intense 3D seismic inputs. The run times for the acoustic FSI simulations in MOOSE are an order of magnitude, or more, shorter than for the corresponding ALE simulations in LS-DYNA. The utility of the MOOSE acoustic FSI implementation is demonstrated by seismic analysis of a building equipped with a fluid-filled, advanced nuclear reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Analysis into Asymptotic Convergence to Full Nonlinear Solutions and Exploration of the Implication of Numerical Operator Mutation of Differential Systems

A robust, sufficiently accurate and practical hydrodynamic simulation toolset is required as a key component of the modeling and simulation of air-gap electrostatic discharge events. This work was performed to complement these ongoing efforts. In particular, hydrodynamic simulations must be vetted to ensure they are robust and sufficiently accurate over relevant characteristic scales. Verification models were generated in order to cultivate the technical knowledge and expertise needed to properly create, implement and execute numerical simulations. Furthermore, this effort was utilized extensively to educate students on the mathematical and numerical principles underlying hydrodynamic simulations. This education opportunity, provided in a holistic and rigorous manner, has greatly benefited developing scientists and engineers with the necessary understandings and toolsets required to excel at accomplishing the task at hand, and, more generally, it has enabled them to generate key programmatic deliverables. This report articulates several subtilties; specifically, how perturbations, nonlinear behavior, and dissipative mechanisms influence numerical stability, how to properly structure mathematical and numerical solutions, and how to properly generate error estimation/assignment. A more rigorous discussion of the consequences of such topics can be found in the body of this report in Chapters 2 and 3 with qualitative findings discussed in Chapter 4.

97 MATHEMATICS AND COMPUTING↗

Validation of Combined Analytical Methods to Predict Slip in Cylindrical Roller Bearings

This paper presents a combination of models that together calculate the cage and roller speeds of a cylindrical roller bearing. The models consider elastohydrodynamic lubrication and contact elasticity between the roller and raceway, roller centrifugal forces, hydrodynamic lubrication at the cage pocket, and frictional forces. Using these models, the predicted cage and roller speeds and the extent of slip are compared to measurements acquired on cylindrical roller bearings in a commercial gearbox in steady-state and transient operating conditions of a wind turbine. In steady-state conditions at low wind speeds and low lubricant temperatures, cage and roller slip up to 60% occur in the loaded zone of the bearing. In the unloaded zone, up to 80% roller slip occurs. Cage and roller slip then decrease as the lubricant temperature and wind speed increases. In general, the analytical model results match experimental measurements within 10% for lubricant temperatures above 40 degrees C and wind speeds over 10 meters per second. The analytical model is further evaluated during a transient start-up event and highly dynamic emergency stop event and also used to examine changes in the bearing design or lubricant properties. Furthermore, roller and cage slip in cylindrical bearings is a combined effect of the bearing design, applied load, shaft speed, and lubricant properties and temperature, and can be quickly evaluated with the combined analytical models.

17 WIND ENERGY↗

Impact of moment-based, energy integrated neutrino transport on microphysics and ejecta in binary neutron star mergers

We present an extensive study of the effects of neutrino transport in three-dimensional general relativistic radiation hydrodynamics (GRHD) simulations of binary neutron star (BNS) mergers using our moment-based, energy-integrated neutrino radiation transport (M1) scheme. Here, we consider a total of eight BNS configurations, while varying equation of state models, mass ratios, and grid resolutions, for a total of 16 simulations. We find that M1 neutrino transport is crucial in modeling the local absorption of neutrinos and the deposition of lepton number throughout the medium. We provide an in-depth look at the effects of neutrinos on the fluid dynamics and luminosity during the late inspiral and postmerger phases, the properties of ejecta and outflow, and the postmerger nucleosynthesis. The simulations presented in this work comprise an extensive study of the combined effect of the equation of state and M1 neutrino transport in GRHD simulations of BNS mergers, and establish that the solution provided by our M1 scheme is robust across system properties.

150 ≤ A ≤ 18959 ≤ A ≤ 8990 ≤ A ≤ 149↗

Dwarf Galaxies at Cosmic Noon: New JWST Constraints on Satellite Models and Subhalo Tidal Evolution

The advent of JWST has revolutionized the study of faint satellite galaxies at z ≳ 1, enabling statistical constraints on galaxy evolution and the galaxy–halo connection in a previously unexplored mass and redshift regime. We compare satellite abundances at 1 < z < 3.5 from recent JWST observations with predictions from cosmological dark-matter-only zoom-in simulations. We identify and quantify several sources of biases that can impact theoretical satellite counts, finding that assumptions about subhalo tidal evolution introduce the largest uncertainty in predictions for the satellite mass function. Using a flexible galaxy disruption model, we explore a range of disruption scenarios, spanning hydrodynamically motivated and idealized prescriptions, to bracket plausible physical outcomes. We show that varying galaxy durability can change the predicted satellite mass functions by a factor of ∼3.5. The JWST data and our fiducial model are consistent within 1σ–2σ across the full redshift (1 < z < 3.5) and stellar mass (M⋆ > 10 7 M ⊙ ) range probed. We find evidence that subhalos are at least as long-lived as predicted by hydrodynamic simulations. Our framework will enable robust constraints on the tidal evolution of subhalos with future observations. This work presents the first direct comparison between cosmological models and observations of the high-redshift satellite population in this low-mass regime. These results showcase JWST’s emerging power to test structure formation in the first half of the Universe in a new domain and to constrain the physical processes driving the evolution of low-mass galaxies across cosmic time.

79 ASTRONOMY AND ASTROPHYSICS↗

Development of a multi-layer canopy model for E3SM Land Model with support for heterogeneous computing

The vertical structure of vegetation canopies creates micro-climates. However, the land components of most Earth System Models, including the Energy Exascale Earth System Model (E3SM), typically neglect vertical canopy structure by using a single layer big-leaf representation to simulate water, CO 2 , and energy exchanges between the land and the atmosphere. In this study, we developed a Multi-Layer Canopy Model for the E3SM Land Model to resolve the micro-climate created by vegetation canopies. The model developed in this study re-implements the CLM-ml_v1 to support heterogeneous computing architectures consisting of CPUs and GPUs and includes three additional optimization-based stomatal conductance models. The use of Portable, Extensible Toolkit for Scientific Computation provides a speedup of 25–50 times on a GPU relative to a CPU. The numerical implementation of the model was verified against CLM-ml_v1 for a month-long simulation using data from the Ameriflux US-University of Michigan Biological Station site. Model structural uncertainty was explored by performing control simulations for five stomatal conductance models that exclude and include the control of plant hydrodynamics (PHD) on photosynthesis. The bias in simulated sensible and latent heat fluxes was lower when PHD was accounted for in the model. Additionally, six idealized simulations were performed to study the impact of three environmental variables (i.e. air temperature, atmospheric CO 2 , and soil moisture) on canopy processes (i.e. net CO 2 assimilation, leaf temperature, and leaf water potential). Increasing air temperature reduced net CO 2 assimilation and increased air temperature. Net CO 2 assimilation increased at higher atmospheric CO 2 , while decreasing soil moisture resulted in lower leaf water potential.

54 ENVIRONMENTAL SCIENCES↗

Mechanistic mass transfer in hollow fiber membrane solvent extraction for bio-based isobutanol

Membrane solvent extraction (MSE) has emerged as a promising method for selectively recovering bioproducts from complex aqueous streams. Bio-isobutanol, a next-generation feedstock for biofuel, remains challenging to recover because of its low concentration and the presence of inhibitory substances. This study explores the potential of hollow fiber (HF) MSE for bio-isobutanol recovery and systematically examines the coupled effects of fiber packing, shell-side flow dynamics, and aqueous chemistry on performance. A resistance-in-series model is applied to understand mass transfer in the HF MSE modules, quantify local resistances, and validate overall performance. The results show that increasing the fiber packing provides a larger interfacial area but induces poor flow distribution and channeling, hindering effective isobutanol transport. Meanwhile, increasing the shell-side velocity improves isobutanol recovery due to reductions in the boundary layer thickness. The presence of salts, added to mimic fermentation broth, increases the partition coefficient through salting-out effects, further improving isobutanol flux. A modified correlation for the shell-side mass transfer coefficient (k s,ϕ+v ), integrating geometric and hydrodynamic effects, was developed and validated. The proposed model achieves highly predictive accuracy (r 2 = 0.9808) across a wide range of conditions, outperforming previous models. The findings provide mechanistic insight into the interaction of geometric packing, hydrodynamics, and chemistry in governing mass transfer in HF MSE. Overall, this work demonstrates the potential of HF MSE for efficient bio-isobutanol recovery and also provides practical guidelines on critical factors (packing fraction, partition coefficient, and shell-side velocity), aiding in the design and scaling of MSE systems for resource recovery.

Aqueous chemistry↗

Atomic kinetics of laboratory photoionized plasmas relevant to astrophysics. Final report

The major goal of this project was to study the atomic and radiation physics of photoionized plasmas driven by a broadband x-ray flux and its impact on plasma x-ray heating and ionization, emissivity and opacity, and radiation-hydrodynamics. This problem is relevant to the physics of accretion disks surrounding black holes, x-ray binaries, and active galactic nuclei in astrophysics. In addition, it is also important for understanding non-equilibrium laboratory plasmas driven by a distribution of photons. To this end, and motivated by previous and ongoing gas-cell experiments performed by the PI at the Z facility of Sandia National Laboratories, we developed a new experimental platform at the 1MA Zebra pulsed power accelerator of the University of Nevada, Reno based on using a gas jet produced by a supersonic nozzle. The gas jet is turned into a photoionized plasma driven by the x- ray flux from a wire-array z-pinch. Establishing this new experimental platform at Zebra was a major goal of this project. In addition, we have also dedicated effort to continue complementary experiments at the Z facility using the gas-cell photoionized plasma set up with support from the Fundamental Science Program at Z as well as to perform data processing and analysis, and theory and modeling simulations, including electron kinetics and radiation-hydrodynamics. Since we seek to test and establish what physics models are needed to describe photoionized plasma guided by experimental observation, these complementary experimental efforts were central to achieving the goals of this project.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Minimal one-dimensional model of bad metal behavior from fast particle-hole scattering

A strongly interacting plasma of linearly dispersing electron and hole excitations in two spatial dimensions (2D), also known as a Dirac fluid, can be captured by relativistic hydrodynamics and shares many universal features with other quantum critical systems. Here, we propose a one-dimensional (1D) model to capture key aspects of the 2D Dirac fluid while including lattice effects and being amenable to nonperturbative computation. When interactions are added to the Dirac-like 1D dispersion without opening a gap, we show that this kind of irrelevant interaction is able to preserve Fermi-liquid-like quasiparticle features while relaxing a zero-momentum charge current via collisions between particle-hole excitations, leading to resistivity that is linear in temperature via a mechanism previously discussed for large-diameter metallic carbon nanotubes. We further provide a microscopic lattice model and obtain numerical results via density-matrix renormalization group simulations, which support the above physical picture. The limits on such fast relaxation at strong coupling are of considerable interest because of the ubiquity of bad metals in experiments.

1-dimensional systems↗

Dark energy survey year 1 results: Constraining baryonic physics in the Universe

ABSTRACT Measurements of large-scale structure are interpreted using theoretical predictions for the matter distribution, including potential impacts of baryonic physics. We constrain the feedback strength of baryons jointly with cosmology using weak lensing and galaxy clustering observables (3 × 2pt) of Dark Energy Survey (DES) Year 1 data in combination with external information from baryon acoustic oscillations (BAO) and Planck cosmic microwave background polarization. Our baryon modelling is informed by a set of hydrodynamical simulations that span a variety of baryon scenarios; we span this space via a Principal Component (PC) analysis of the summary statistics extracted from these simulations. We show that at the level of DES Y1 constraining power, one PC is sufficient to describe the variation of baryonic effects in the observables, and the first PC amplitude (Q1) generally reflects the strength of baryon feedback. With the upper limit of Q1 prior being bound by the Illustris feedback scenarios, we reach $\sim 20{{\ \rm per\ cent}}$ improvement in the constraint of $S_8=\sigma _8(\Omega _{\rm m}/0.3)^{0.5}=0.788^{+0.018}_{-0.021}$ compared to the original DES 3 × 2pt analysis. This gain is driven by the inclusion of small-scale cosmic shear information down to 2.5 arcmin, which was excluded in previous DES analyses that did not model baryonic physics. We obtain $S_8=0.781^{+0.014}_{-0.015}$ for the combined DES Y1+Planck EE+BAO analysis with a non-informative Q1 prior. In terms of the baryon constraints, we measure $Q_1=1.14^{+2.20}_{-2.80}$ for DES Y1 only and $Q_1=1.42^{+1.63}_{-1.48}$ for DESY1+Planck EE+BAO, allowing us to exclude one of the most extreme AGN feedback hydrodynamical scenario at more than 2σ.

79 ASTRONOMY AND ASTROPHYSICS↗

Trends of Sediment Resuspension and Budget in Southern Lake Michigan Under Changing Wave Climate and Hydrodynamic Environment

Sediment suspension and transport driven by waves and currents play a significant role in both the ecological and physical environments of large lakes. Lake Michigan has faced a rapidly increasing water level associated with intensified wind waves in the past decade. To investigate the spatiotemporal characteristics of suspended sediment concentration (SSC) and associated coastal sediment budgets in southern Lake Michigan, a 30-year (1991–2020) hindcast was performed using a coupled wave-current-sediment model (SWAN-FVCOM-CSTMS). We found that in southern Lake Michigan, the basin-wide mean SSC increased, and the coastal sediment loss accelerated dramatically, corresponding with intensified waves, currents and lake water level rises over the past decade. The basin-wide mean SSC, coastal sediment loss, wave height, wind speed, current speed, and water level in southern Lake Michigan are highly correlated. Spatially, the results reveal decreases in coastal SSC and sediment loss in the western portion of the southern basin, while the eastern sectors show an increase in both metrics. This reflects a clear shift in the wave climate and hydrodynamic environment. The alterations in long-term coastal sediment budgets imply that considerable shoreline transformations are being influenced by modifications in the wave climate. Understanding the spatiotemporal characteristics of SSC and coastal sediment budgets is crucial for strategic water resource management and coastal infrastructure planning.

54 ENVIRONMENTAL SCIENCES↗