Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow conservation”

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
All sourcesNASA NTRSDOE OSTI
DOE OSTI2026

Viscous Gubser flow with conserved charges to benchmark fluid simulations

We present semi-analytical solutions for the evolution of both the temperature and chemical potentials for viscous Gubser flow with conserved charges. Such a solution can be especially useful in testing numerical codes intended to simulate relativistic fluids with large chemical potentials. The freeze-out hypersurface profiles for constant energy density are calculated, along with the corresponding normal vectors, and presented as a new unit test for numerical codes. We also compare the influence of the equation of state on the semi-analytical solutions. We benchmark the newly developed smoothed particle hydrodynamics code ccake that includes both shear viscosity and three conserved charges. Here, the numerical solutions are in excellent agreement with the semi-analytical solution and are also able to accurately reproduce the hypersurface at freeze-out.

Hydrodynamic models↗
DOE OSTI2024

The MOOSE fluid properties module

The Fluid Properties module within the Multiphysics Object-Oriented Simulation Environment (MOOSE) is used to compute fluid properties for numerous applications, ranging from nuclear reactor thermal hydraulics to geothermal energy. Those applications drove the development of the module to enable numerous different fluid equations of states, property lookups with primitive and conserved flow variable to cater to pressure and density-driven solvers, and an object-oriented design facilitating expansion and maintenance. Each fluid property is implemented in its own class but inherits capabilities such as automatic differentiation, automated out-of-bounds handling or variable conversion capabilities. Here, this paper presents the module, its design, its user and developer interface, its content in terms of fluids and properties, and several of its applications showing its major role in the MOOSE simulation ecosystem.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2022

A consistent and conservative Phase-Field method for multiphase incompressible flows

In the present study, a consistent and conservative Phase-Field method, including both the model and scheme, is developed for multiphase flows with an arbitrary number of immiscible and incompressible fluid phases. The consistency of mass conservation and the consistency of mass and momentum transport are implemented to address the issue of physically coupling the Phase-Field equation, which locates different phases, to the hydrodynamics. These two consistency conditions, as illustrated, provide the “optimal” coupling because (i) the new momentum equation resulting from them is Galilean invariant and implies the kinetic energy conservation, regardless of the details of the Phase-Field equation, and (ii) failures of satisfying the second law of thermodynamics or the consistency of reduction of the multiphase flow model only result from the same failures of the Phase-Field equation but are not due to the new momentum equation. Physical interpretation of the consistency conditions and their formulations are first provided, and general formulations that are obtained from the consistency conditions and independent of the interpretation of the velocity are summarized. Then, the present consistent and conservative multiphase flow model is completed by selecting a reduction consistent Phase-Field equation. Several novel techniques are developed to inherit the physical properties of the multiphase flows after discretization, including the gradient-based phase selection procedure, the momentum conservative method for the surface force, and the general theorems to preserve the consistency conditions on the discrete level. Equipped with those novel techniques, a consistent and conservative scheme for the present multiphase flow model is developed and analyzed. The scheme satisfies the consistency conditions, conserves the mass and momentum, and assures the summation of the volume fractions to be unity, on the fully discrete level and for an arbitrary number of phases. All those properties are numerically validated. Finally, numerical applications demonstrate that the present model and scheme are robust and effective in studying complicated multiphase dynamics, especially for those with large-density ratios.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Path-Integrated X-Ray Images for Multi-Surface Digital Image Correlation (PI-DIC)

X-ray imaging offers unique possibilities for Digital Image Correlation (DIC), opening the door for full-field deformation measurements of a test article in complex environments where optical DIC suffers severe biases or is impossible. While X-ray DIC has been performed in the past with standard DIC codes designed for optical images, the path-integrated nature of X-ray images places constraints on the experimental setup, predominantly that only a single surface of interest moves/deforms. These requirements are difficult to realize for many practical situations and limit the amount of information that can be garnered in a single test. Other X-ray based diagnostics such as Digital Volume Correlation (DVC) and Projection DVC (P-DVC) overcome these obstacles, but DVC is limited to quasi-static tests, and both DVC and P-DVC necessitate high-resolution computed tomography (CT) scan(s) and often require a potentially invasive pattern throughout the volume of the specimen. Here this work presents a novel approach to measure time-resolved displacements and strains on multiple surfaces from a single series of 2D, path-integrated (PI) X-ray images, called PI-DIC. The principle of optical flow or conservation of intensity—the foundation of DIC—was reframed for path-integrated images, for an exemplar setup comprised of two plates moving and deforming independently. Synthetic images were generated for rigid translations, rigid rotations, and uniform stretches, where each plate underwent a unique motion/deformation. Experimental specimens were fabricated (either an aluminum plate with tantalum features or a plastic plate with steel features) and the two specimens were independently translated. PI-DIC was successfully demonstrated with the synthetic images and validated with the experimental images. Prescribed displacements were recovered for each plate from the single set of path-integrated, deformed images. Errors were approximately 0.02 px for the synthetic images with 1.5% image noise, and 0.05 px for the experimental images. These results provide the foundation for PI-DIC to measure motion and deformation of multiple, independent surfaces with subpixel accuracy from a single series of path-integrated X-ray images.

47 OTHER INSTRUMENTATION↗
DOE OSTI2021

A consistent and conservative model and its scheme for N -phase- M -component incompressible flows

Here, we propose a consistent and conservative model for multiphase and multicomponent incompressible flows, where there can be arbitrary numbers of phases and components. Each phase has a background fluid called the pure phase, each pair of phases is immiscible, and components are dissolvable in some specific phases. The model is developed based on the multiphase Phase-Field model including the contact angle boundary condition, the diffuse domain approach, and the analyses on the proposed consistency conditions for multiphase and multicomponent flows. The model conserves the mass of individual pure phases, the amount of each component in its dissolvable region, and thus the mass of the fluid mixture, and the momentum of the flow. It ensures that no fictitious phases or components can be generated and that the summation of the volume fractions from the Phase-Field model is unity everywhere so that there is no local void or overfilling. It satisfies a physical energy law and it is Galilean invariant. A corresponding numerical scheme is developed for the proposed model, whose formal accuracy is 2nd-order in both time and space. It is shown to be consistent and conservative and its solution is demonstrated to preserve the Galilean invariance and energy law. Numerical tests indicate that the proposed model and scheme are effective and robust to study various challenging multiphase and multicomponent flows.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Degradation of performance in ICF implosions due to Rayleigh–Taylor instabilities: A Hamiltonian perspective

The Rayleigh–Taylor instability (RTI) is an ubiquitous phenomenon that occurs in inertial-confinement-fusion (ICF) implosions and is recognized as an important limiting factor of ICF performance. To analytically understand the RTI dynamics and its impact on ICF capsule implosions, we develop a first-principle variational theory that describes an imploding spherical shell undergoing RTI. The model is based on a thin-shell approximation and includes the dynamical coupling between the imploding spherical shell and an adiabatically compressed fluid within its interior. Using a quasilinear analysis, we study the degradation trends of key ICF performance metrics (e.g., stagnation pressure, residual kinetic energy, and areal density) as functions of initial RTI parameters (e.g., the initial amplitude and Legendre mode), as well as the 1D implosion characteristics (e.g., the convergence ratio). We compare analytical results from the theory against nonlinear results obtained by numerically integrating the governing equations of this reduced model. Our findings emphasize the need to incorporate polar flows in the calculation of residual kinetic energy and demonstrate that higher convergence ratios in ICF implosions lead to significantly greater degradation of key performance metrics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Relativistic Dynamics of Fluctuations and QCD Critical Point

Herein, we review the recent development on a systematic deterministic formalism describing dynamics of both bulk and fluctuations in an arbitrary relativistic hydrodynamic flow carrying conserved charges. In particular, we discuss the implementation of such general formalism near the QCD critical point.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Modeling and simulations of hydrodynamic shocks in a plasma flowing across randomized ICF scale laser beams

High-energy laser beams interacting with flowing plasmas can produce a plasma response that leads to deflection of the beam, beam bending. Such beams have usually a speckle structure generated by optical smoothing techniques that reduce the spatial and temporal coherence in the laser field pattern. The cumulative plasma response from laser speckles slows down the velocity of the incoming flow by momentum conservation. For slightly super-sonic flow the cumulative plasma response to the ponderomotive force exerted by the beam speckle ensemble is the strongest, such that slowing down the flow to subsonic velocities leads eventually to the generation of a shock around the cross section of the beam. This scenario has been predicted theoretically and is confirmed here by our hydrodynamic simulations in two dimensions with speckled beams and in one dimension with a reduced model. The conditions of shock generation are given in terms of the ponderomotive pressure, speckle size and the flow velocity. The nonlinear properties of the shocks are analyzed using Rankine–Hugoniot relations. According to linear theory, temporally smoothed laser beams exhibit a higher threshold for shock generation. Numerical simulations with beams that are smoothed by spectral dispersion compare well with the linear theory results, diverging from those produced by beams with only a random phase plates in the nonlinear regime. The conditions necessary for shock generation and their effects on the laser plasma coupling in inertial confinement fusion (ICF) experiments are also discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Adaptive clipping‐and‐redistribution algorithms for bounded and conservative high‐order interpolations applied to discontinuous and reactive flows

Abstract A new adaptive clipping‐and‐redistribution method is presented which provides bounds‐preservation for multidimensional interpolation in the context of high‐order finite‐volume discretizations with adaptive mesh refinement (AMR). The underlying finite‐volume method (FVM) for the computational fluid dynamics applications is fourth‐order accurate for smooth solutions and utilizes AMR for computational efficiency in solving multiscale problems involving turbulence and combustion. High‐order interpolation between different AMR levels is required. However, this operation often leads to numerical issues because combustion species must have physical bounds preserved. The present study overcomes two major challenges in the development of the high‐order interpolation method. First, the method needs to be bound‐preserving near extrema or discontinuities to prevent the emergence of unphysical oscillations while maintaining fourth‐order accuracy in smooth flows. Second, the method needs to satisfy the conservation requirement in multiple dimensions, particularly in the context of curvilinear coordinate transformations. Additionally, the method is designed to be localized and computationally inexpensive. The new interpolation scheme is demonstrated by solving reacting flows, which are extremely sensitive to unphysical overshoots in conserved quantities. The test problems are shock‐induced ‐ combustion and a ‐air flame in a practical bluff‐body combustor. Results show the method prevents new extrema near discontinuities while maintaining high‐order accuracy in smooth regions. In particular, the method is extremely beneficial for combustion with stiff chemistry. With the proposed new method, even if flame fronts cross AMR interfaces or new grids are created in the vicinity of the flame, solution stability is retained.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Quantifying the Impact of Lagged Hydrological Responses on the Effectiveness of Groundwater Conservation

Many irrigated agricultural areas seek to prolong the lifetime of their groundwater resources by reducing pumping. However, it is unclear how lagged responses, such as reduced groundwater recharge caused by more efficient irrigation, may impact the long-term effectiveness of conservation initiatives. Here, we use a variably saturated, simplified surrogate groundwater model to: 1) analyze aquifer responses to pumping reductions, 2) quantify time lags between reductions and groundwater level responses, and 3) identify the physical controls on lagged responses. We explore a range of plausible model parameters for an area of the High Plains Aquifer (USA) where stakeholder-driven conservation has slowed groundwater depletion. We identify two types of lagged responses that reduce the long-term effectiveness of groundwater conservation, recharge-dominated and lateral-flow-dominated, with vertical hydraulic conductivity (K Z ) the major controlling variable. When high K z allows percolation to reach the aquifer, more efficient irrigation reduces groundwater recharge. By contrast, when low K z impedes vertical flow, short term changes in recharge are negligible, but pumping reductions alter the lateral flow between the groundwater conservation area and the surrounding regions (lateral-flow-dominated response). For the modeled area, we found that a pumping reduction of 30% resulted in median usable lifetime extensions of 20 or 25 years, depending on the dominant lagged response mechanism (recharge- vs. lateral-flow-dominated). These estimates are far shorter than estimates that do not account for lagged responses. To conclude, results indicate that conservation-based pumping reductions can extend aquifer lifetimes, but lagged responses can create a sizable difference between the initially perceived and actual long-term effectiveness.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Hanford Waste Treatment Plant Low Activity Waste Facility Stack Effluent Monitoring: Sampling Probe Location Qualification Evaluation (Rev.1)

The Hanford Tank Waste Treatment and Immobilization Plant Low Activity Waste (LAW) facility stack monitor locations were qualified using scale model stacks to mitigate the risk of discovering that sampling locations do not meet the qualification criteria on the full-scale stacks. As required by the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1-1999 standard, the scale model and its sampling location were geometrically similar to the actual stack and the Reynolds numbers for both the actual and model stacks were >10,000. An additional criterion is that the product of the hydraulic diameter and mean velocity (DV) of the full-scale stack must be between 1/6 DV and 6 DV of the scale model stack tests. Verification tests of the LAW stacks were performed at normal operating conditions. The minimum 1/6 DV value, along with the maximum 6 DV value from the scale model testing, determines the range of stack flow rates for which the full-scale stack may be operated and remain in compliance with the stack verification criterion. For this analysis, the range of qualified flow rates listed is conservatively based on the average DV through 6 DV for LV-S1, LV-S2, and LV-C2, and 1/3 DV to 3 DV for LV-S3. Table S1 lists the operating flow rates along with the conservative lower and upper qualified stack flow rates for each of the LAW facility stacks. For each stack, the operating flow is below the upper qualified stack.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2024

Limit cycle oscillations in the zonal-flow-catalyzed interactions of ion-temperature-gradient turbulence

Limit-cycle oscillations are studied for ion temperature gradient turbulence, which, in the absence of large diamagnetic (mean) shear flows, saturates through energy transfer from unstable modes to large-scale stable modes via zonal-flow intermediary modes. Oscillations of zonal flow and turbulence levels are strongly constrained by the reactive, largely non-dissipative character of the zonal flows. Since existing predator–prey models for observed oscillations in experiments do not include energy transfer through zonal flows to stable modes, low-order fluid models with this physics are constructed and investigated. A simple three-wave truncation produces low-amplitude zonal flows that slowly oscillate around a zero mean, with turbulence oscillations between coupled wavenumbers that exceed linear frequencies by orders of magnitude. This inconsistency with experimental observations is caused by the weak non-linear drive of zonal flows in three-wave systems and the lack of multiple-wavenumber turbulent interactions. A more comprehensive model that preserves multiple wavenumber interactions within the context of conservative zonal-flow-mediated energy transfer to stable modes accurately reflects observed dynamics when the phase between stable and unstable modes is occasionally randomized.

Li, P. -Y. (ORCID:0000000295254171)↗
DOE OSTI2020

Flexible and Modular Simultaneous Modeling of Flow and Reactive Transport in Rivers and Hyporheic Zones

Investigations of coupled multiphysics processes in rivers and hyporheic zones have extensively used numerical models. Most existing models use a sequential, one-way coupling between the surface and subsurface domains. Such one-way coupling potentially introduces error. To overcome this, a fully coupled model, hyporheicFoam, was developed using the open-source computational platform OpenFOAM. It captures the coupled flow and multicomponent reactive transport processes within both surface and subsurface domains and across their interface. The coupling between two domains is implemented by mapping conservative flux boundary conditions at the interface through an iterative algorithm. Reactive transport is enabled by specifying a reaction network. To start, we have implemented reaction kinetics following the double Monod-type model with inhibition. The model capability is illustrated through modeling of both conservative and reactive hyporheic flow and transport through dune bedforms. With the novel coupled model, it is now possible to quantify reactions wherein the reactants and products are constantly exchanging between domains and have feedbacks. hyporheicFoam can simulate large, three-dimensional cases owing to the computational flexibility and power offered by the code structure and parallel design of OpenFOAM.

58 GEOSCIENCES↗