Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “flow field”

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 415 records · Page 23

Recovering the effects of subgrid heterogeneity in simulations of radionuclide transport through fractured media

Groundwater flow and contaminant transport through fractured media can be simulated using Discrete Fracture Network (DFN) models which provide a natural description of structural heterogeneity. However, this approach is computationally expensive, with the large number of intersecting fractures necessitated by many real-world applications requiring modeling simplifications to be made for calculations to be tractable. Upscaling methods commonly used for this purpose can result in some loss of local-scale variability in the groundwater flow velocity field, resulting in underestimation of particle travel times, transport resistance and retention in transport calculations. In this paper, a transport downscaling algorithm to recover the transport effects of heterogeneity is tested on a synthetic Brittle Fault Zone model, motivated by the problem of large safety assessment calculations for geological repositories of spent nuclear fuel. We show that the variability in the local-scale velocity field which is lost by upscaling can be recovered by sampling from a library of DFN transport paths, accurately reproducing DFN transport statistic distributions and radionuclide breakthrough curves in an upscaled model.

54 ENVIRONMENTAL SCIENCES↗

An Improved Model For Determining Salinity Recharge Time for Deep Borehole Disposal - 20405

In a previous publication (K. P. Travis, D. Burley and F. G. F. Gibb, WM2017 Conference, Phoenix Arizona paper no. 17480) we introduced a numerical model to estimate the time it would take for a pressure perturbation to subside following its creation from the construction of a deep geological borehole in water-saturated rock. The model was based on the notion of a point source of momentum - a solution of the time-dependent pressure-diffusion equation. Using the model, we estimated that it would take on the order of 10 k years for physical equilibrium to become re-established following the sinking of a 5 km borehole in granite. The significance of such a model is that it places an upper bound on the lifetime required of a borehole sealing system - a necessary input to a borehole post-closure safety case assessment. No engineered sealing system has ever been devised which is capable of retaining its sealing properties for 10 half-lives of long-lived radioisotopes in spent fuel or high-level waste. One of the key advantages of disposing of nuclear waste via Deep Borehole Disposal (DBD) is the natural sealing provided by density stratified groundwater. Upon sinking of a borehole, and subsequent filling with fresh water, brine or drilling mud, this natural barrier may be temporarily damaged. Over time, fresh brine from the far-field will flow towards or away from the hole (driven by a pressure gradient) and will eventually re-establish the original salinity gradient. It follows that the engineered seals need last only as long as the time required for this salinity gradient to reset itself. One of the limitations of our previous model was the use of a static boundary condition on the borehole wall. The model used a boundary pressure which varied quadratically with depth (arising from differences between the pressure of a column of fresh water in the borehole and that of a column of brine in the host rock). However as brine replaces fresh water in the borehole (driven by a pressure gradient), the boundary function must change with time. We now introduce an improved model which takes this time dependent boundary condition into account. The model also takes into account the time taken for a mixture of brine and fresh water to re-establish chemical equilibrium through the process of diffusion. The paper contains the mathematical details of our new iterative model as well as results showing the time taken to reach steady state, and concentration profiles for the components of the brine-filled borehole as a function of time together with a discussion on the implications of the results for developing a post-closure safety assessment for DBD. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Transported PDF Modeling of Compressible Turbulent Reactive Flows by using the Eulerian Monte Carlo Fields Method

Although the transported probability density function (PDF) method has been developed for decades, its application has been mainly focused on the low-Mach number flow problems. This work extends the transported PDF method to compressible flow problems. The Eulerian Monte Carlo fields (EMCF) solution method is employed to solve the transported PDF equation for compressible flow problems. A pseudo stagnation enthalpy is introduced and its stochastic partial differential equation is derived to ensure total energy conservation numerically. A new mixing model called interaction by partial exchange with mean (IPEM) is introduced to expand the available choices of mixing models for the EMCF method. The consistency of the EMCF method is examined for solving the transported PDF equation. Numerical implementation details are discussed, such as the density coupling between the compressible flow solver and the EMCF solver, discretization schemes for the mixing terms and the stochastic terms. The implemented compressible flow solver coupled with the EMCF solver is verified and validated in a series of test cases with increasing level of complexity, ranging from a statistically one-dimensional turbulent mixing layer to a self-excited resonance model rocket combustor. It is observed that in general with the increase of compressibility, there is an increase in the sensitivity of the modeling results to the different models and algorithms. This makes it necessary to develop a thorough understanding of the model sensitivity in order to develop a robust and accurate simulation solver for highly compressible turbulent reactive flows. The thermo-acoustic instability inside the model rocket combustor case is captured reasonably, which demonstrates the overall capability of the developed compressible turbulent combustion solver based on the transported PDF method.

42 ENGINEERING↗

Mechanism for sequestering magnetic energy at large scales in shear-flow turbulence

Straining of magnetic fields by large-scale shear flow, which is generally assumed to lead to intensification and generation of small scales, is reexamined in light of the persistent observation of large-scale magnetic fields in astrophysics. It is shown that, in magnetohydrodynamic turbulence, unstable shear flows have the unexpected effect of sequestering magnetic energy at large scales due to counteracting straining motion of nonlinearly excited large-scale stable eigenmodes. This effect is quantified via dissipation rates, energy transfer rates, and visualizations of magnetic field evolution by artificially removing the stable modes. In conclusion, these analyses show that predictions based upon physics of the linear instability alone miss substantial dynamics, including those of magnetic fluctuations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Implementing contact angle boundary conditions for second-order Phase-Field models of wall-bounded multiphase flows

In the present work, a general formulation is proposed to implement the contact angle boundary conditions for the second-order Phase-Field models, which is applicable to N-phase (N ≥ 2) moving contact line problems. To remedy the issue of mass change due to the contact angle boundary condition, a source term or Lagrange multiplier is added to the original second-order Phase-Field models, which is determined by the consistent and conservative volume distribution algorithm so that the summation of the order parameters and the consistency of reduction are not influenced. To physically couple the proposed formulation to the hydrodynamics, especially for large-density-ratio problems, the consistent formulation is employed. The reduction-consistent conservative Allen-Cahn models are chosen as examples to illustrate the application of the proposed formulation. The numerical scheme that preserves the consistency and conservation of the proposed formulation is employed to demonstrate its effectiveness. Results produced by the proposed formulation are in good agreement with the exact and/or asymptotic solutions. The proposed method captures complex dynamics of moving contact line problems having large density ratios.

97 MATHEMATICS AND COMPUTING↗

Flow Sensor Test Article (F-STAr) (Design and Fabrication Status Report)

The Flow Sensor Test Article (F-STAr) is a test article currently under construction for the Mechanisms Engineering Test Loop (METL). F-STAr was designed to provide high sodium flowrate capabilities for sensor calibration, component testing, and fluid studies. Figure 1 shows two solid model views of the new test article. F-STAr includes a high-capacity pump that can provide a nominal flowrate of 120 GPM; a test section support structure that can accommodate a wide array of sub-test articles and their instrumentation; and finally, a heating and cooling system to aid in controlling the testing environment. Initially, F-STAr will be configured to test liquid metal flow sensors, specifically field shift sensors like the Eddy Current Flow Sensors (ECFS) based on the RDT C4-7T standard. To test these sensors, two test sections were designed that attempt to model Sodium Fast Reactor (SFR) outlet conditions. Figure 2 shows models of the “Full-Scaled Test Section” (FSTS) that represents a generic SFR fuel handling socket and a “Pseudo-Scaled Test Section” (PSTS) that represents a generic array of scaled fuel handling sockets. While F-STAr will be configured to test ECFS’s, it can also be outfitted to meet other experimental needs. For example, F-STAr could be setup with a test section that includes a fluidic diode or a component test investigating the performance of hydrodynamic bearings. Other test sections include studies of sodium thermal hydraulics like thermal striping. Overall, F-STAr is a flexible test article designed to accommodate many needs. This report will provide a status update on the design and construction of F-STAr. Manufacturing of all components has commenced, and several components have already been completed. These components include the pump and test section assembly stand. Other components, such as the heater, have been completed but were rejected due to the vendor not meeting the requirements of the purchase order. Lastly, the submersible flowmeter and main flange components are under construction, and their status will be reviewed in this report. Overall, most of the F-STAr components and parts will be completed by the end of September 2022.

42 ENGINEERING↗

Observations of Stably Stratified Flow through a Microscale Gap

This paper reports the findings of a comprehensive field investigation on flow through a mountain gap subject to a range of stably stratified environmental conditions. This study was embedded within the Perdigão field campaign, which was conducted in a region of parallel double-ridge topography with ridge-normal wind climatology. One of the ridges has a well-defined gap (col) at the top, and an array of in situ and remote sensors, including a novel triple Doppler lidar system, was deployed around it. The experimental design was mostly guided by previous numerical and theoretical studies conducted with an idealized configuration where a flow (with characteristic velocity U 0 and buoyancy frequency N ) approaches normal to a mountain of height h with a gap at its crest, for which the governing parameters are the dimensionless mountain height G= (Nh)/U 0 and various gap aspect ratios. Modified forms of G were proposed to account for real-world atmospheric variability, and the results are presented in terms of a gap-averaged value G c . The nature of gap flow was highly dependent on G c , wherein a nearly neutral flow regime (G c <1), a transitional mountain wave regime (G c ~O(1)), and a gap-jetting regime ( G c >O(1)) were identified. The measurements were in broad agreement with previous numerical and theoretical studies on a single ridge with a gap or double ridge topography, although details vary. This is the first ever detailed field study reported on microscale (O(100 m)) gap flows, and provides useful data and insights for future theoretical and numerical studies.

Microscale, Gap, environmental fluid dynamics, win↗

Effects of an External Magnetic Field on Keyhole Mode Laser Melting of 316 Stainless Steel

Keyhole-mode laser melting is an efficient method for joining or cutting large, thick components, but controlling keyhole depth and fluctuations has remained challenging. Applying an external magnetic field can control melt pool flows and indirectly influence keyhole morphology and dynamics. The induced Lorentz force, comprising Seebeck and damping components, plays a crucial role in the melt pool dynamics, depending on temperature gradient, flow rates, and magnetic field orientation and magnitude. This research investigates the effects of an external magnetic field on keyhole behavior during laser spot melting of 316 stainless steel using synchronized high-speed synchrotron X-ray and thermal imaging. Findings revealed that a longitudinal magnetic field (120 mT) increased keyhole depth but exacerbated lateral fluctuations, resulted in a 20% increase in the melt pool temperature gradient and a 27% decrease in cooling rate. Conversely, a transverse magnetic field (760 mT) reduced keyhole depth and improved porosity formation. The findings suggest that a decrease in keyhole depth correlates with a decrease in fluctuations, and vice versa. These insights enhance understanding of external magnetic fields’ impact on laser melting, with implications for improving part quality.

Alamdari, Aslan Bafahm [The Ohio State Univ., Colu↗

Including the parallel mass flow in calculating the steady-state solutions and stability of the momentum balance equations for a quasisymmetric stellarator

The Helically Symmetric Experiment (HSX) is a quasisymmetric stellarator with minimal parallel viscous damping in a helical direction. The parallel flow (Vǁ) along the magnetic field is similarly weakly damped by viscosity. In this paper, the self-consistent steady-state parallel and poloidal momentum balance equations are used to show that a large Vǁ on the order of the ion thermal velocity can increase the ion resonant radial electric field (Er) beyond the value calculated using the typical approximation that Vǁ is zero. By altering the damping of Vǁ, either by degrading the quasisymmetry or varying the neutral density, the ion resonant Er can shift in a controllable fashion. It is shown explicitly that there exist stable and unstable steady-state solutions in the two-dimensional space of Vǁ and Er. A stability analysis of each solution is performed by calculating the eigenvalues and eigenvectors of the Jacobian. The unstable solution corresponds to a saddle point in which the eigenvalues have opposite signs. The analysis leads to the conclusion that unstable solutions occur when the derivative of the total poloidal damping with respect to Er is positive. A hysteresis in Er and Vǁ is observed when the radial current density is linearly increased to a maximum and then decreased back to zero. Jumps in the radial electric field and the parallel flow are observed as the radial current density drives the evolution from one stable point to the next. This result is similar to experimental data observed on several devices.

Physics↗

QES-Plume v1.0: a Lagrangian dispersion model

Low-cost simulations providing accurate predictions of transport of airborne material in urban areas, vegetative canopies, and complex terrain are demanding because of the small-scale heterogeneity of the features influencing the mean flow and turbulence fields. Common models used to predict turbulent transport of passive scalars are based on the Lagrangian stochastic dispersion model. The Quick Environmental Simulation (QES) tool is a low-computational-cost framework developed to provide high-resolution wind and concentration fields in a variety of complex atmospheric-boundary-layer environments. Part of the framework, QES-Plume, is a Lagrangian dispersion code that uses a time-implicit integration scheme to solve the generalized Langevin equations which require mean flow and turbulence fields. Here, QES-Plume is driven by QES-Winds, a 3D fast-response model that computes mass-consistent wind fields around buildings, vegetation, and hills using empirical parameterizations, and QES-Turb, a local-mixing-length turbulence model. In this paper, the particle dispersion model is presented and validated against analytical solutions to examine QES-Plume’s performance under idealized conditions. In particular, QES-Plume is evaluated against a classical Gaussian plume model for an elevated continuous point-source release in uniform flow, the Lagrangian scaling of dispersion in isotropic turbulence, and a non-Gaussian plume model for an elevated continuous point-source release in a power-law boundary-layer flow. In these cases, QES-Plume yields a maximum relative error below 6 % when compared with analytical solutions. In addition, the model is tested against wind-tunnel data for a uniform array of cubical buildings. QES-Plume exhibits good agreement with the experiment with 99 % of matched zeros and 59 % of the predicted concentrations falling within a factor of 2 of the experimental concentrations. Furthermore, results also emphasize the importance of using high-quality turbulence models for particle dispersion in complex environments. Finally, QES-Plume demonstrates excellent computational performance.

58 GEOSCIENCES↗

Observations of pressure anisotropy effects within semi-collisional magnetized plasma bubbles

Magnetized plasma interactions are ubiquitous in astrophysical and laboratory plasmas. Various physical effects have been shown to be important within colliding plasma flows influenced by opposing magnetic fields, however, experimental verification of the mechanisms within the interaction region has remained elusive. Here we discuss a laser-plasma experiment whereby experimental results verify that Biermann battery generated magnetic fields are advected by Nernst flows and anisotropic pressure effects dominate these flows in a reconnection region. These fields are mapped using time-resolved proton probing in multiple directions. Various experimental, modelling and analytical techniques demonstrate the importance of anisotropic pressure in semi-collisional, high- β plasmas, causing a reduction in the magnitude of the reconnecting fields when compared to resistive processes. Anisotropic pressure dynamics are crucial in collisionless plasmas, but are often neglected in collisional plasmas. We show pressure anisotropy to be essential in maintaining the interaction layer, redistributing magnetic fields even for semi-collisional, high energy density physics (HEDP) regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Probing interfacial momentum closures in two-phase bubbly flow with machine learning-aided methods

Computational fluid dynamics (CFD) approach has already reached a high level of maturity for single-phase flows, however the development of closure models for two-phase flow requires additional attention. Multiphase CFD (M-CFD) methods resolve the conservation equations for mass, momentum and energy while differing in the approaches and strategies adopted in the physical closure models. The most widely adopted framework for M-CFD is the Eulerian-Eulerian two-fluid approach which assumes that all phases are co-existing inside each computational cell. For each fluid, the full set of conservation equations is solved; therefore, each fluid has a different velocity field. For adiabatic two-phase flow, the mechanisms of the interfacial momentum transfer are modeled by the interfacial forces representing different physical mechanisms. One of the crucial issues in the development and application of two-fluid model is the understanding of the interfacial momentum closures which determines the bubble distribution and migration behaviors. Dedicated experiments are performed to support the physical understanding and drive the closures’ development. However, limitations exist due to the uncertainties in the experimental measurement and the simplified analytical assumptions which have difficulties on representing the complex non-linear flow fields. In this paper, a data-driven approach, Feature Similarity Measurement (FSM), is developed and proposed to resolve the challenges of modeling the interfacial forces closures. Case study is performed with two-phase flow scenarios where the high-fidelity experimental data is available. Within the Eulerian-Eulerian two-fluid framework, only momentum equations for gas and liquid phases are solved and reduced-order interfacial momentum closures are aided with FSM. Predictions of void fraction and velocity fields are analyzed and demonstrate the potential of machine learning-driven interfacial forces closures.

97 MATHEMATICS AND COMPUTING↗

Artificial intelligence velocimetry and microaneurysm-on-a-chip for three-dimensional analysis of blood flow in physiology and disease

Significance Microfluidics is an important in vitro platform to gain insights into mechanics of blood flow and mechanisms of pathophysiology of human diseases. Extraction of 3D fields in microfluidics with dense cell suspensions remains a formidable challenge. We present artificial-intelligence velocimetry (AIV) as a general platform to determine 3D flow fields and a microaneurysm-on-a-chip to simulate blood flow in microaneurysms in patients with diabetic retinopathy. AIV is built on physics-informed neural networks that integrate seamlessly 2D images from microfluidic experiments or in vivo observations with physical laws to estimate full 3D velocity and stress fields. AIV could be integrated into imaging technologies to automatically infer key hemodynamic metrics from in vivo and in vitro biomedical images.

42 ENGINEERING↗

Unraveling the Determinant Mechanisms in Flow-Mediated Crystal Growth and Phase Behaviors

To uncover the critical mechanisms responsible for mesoscopic level development during flow-mediated crystal growth, we develop a semi-two-way hydrodynamic coupled structural phase-field crystal formalism (HXPFC-s2). The new formalism, inspired by previous attempts at coupling hydrodynamic and phase-field crystal (PFC) equations, allows for studying mesoscopic flow-mediated crystallization at diffusive timescales pertinent to industrial applications. Unlike previous efforts, the devised coupling to the structural PFC (XPFC) equations allows generalization to more complex crystal structures through explicit parameterization of the direct correlation function (DCF). Utilizing the HXPFC-s2 formalism, we seek to uncover the determinant physical mechanisms in crystallization under simple shear flows by comparing temperature-driven crystallization to flow-mediated crystallization under varying flow-strengths. Parallels and deviations of under-cooling and flow-strength effects on crystal growth are drawn using the crystal cluster-size and system ordering time evolutions. In doing so, we identify scaling behaviors with a Peclet-like number, Pe∼, a critical Peclet-like number, Pe∼*, and flow-field-crystal plane-dependent interactions. Our findings may be relevant for controlling crystal growth and phase behaviors in flow applications.

Willis, L. Connor (ORCID:0009000961321848)↗

Inferring Heat Flow in Laser Absorption Regions using Diagnostic Magnetic Fields

This research was motivated by the desire to obtain a measurement of heat-flow in overdense regions of a laser ablated plasma for the first time using a high energy proton beam. For sufficiently high energies, the proton beam can pass through overdense regions and only be deflected by electric and magnetic fields. As the transport of magnetic fields in a plasma are intrinsically linked to the transport of thermal energy, it was posited that such a measurement could be used to infer a heat-flow. While the project did not prove the method to be infeasible, difficulties with understanding the proton radiographs without an imposed magnetic field made the final aim of this project unreached. Instead, however, a closer understanding of proton radiographs of laser-solid interactions without an imposed field has been obtained. 5 shot days on The OMEGA Laser Facility have been awarded (3 through ICF and 2 as ride-along) to further investigate some of these features. Extrapolating current results to hohlraum conditions has the magnetic energy being a significant (10%) fraction of the laser drive energy, which has the potential to explain the NIF hohlraum drive deficit.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The impact of magnetic fields on momentum transport and saturation of shear-flow instability by stable modes

The Kelvin–Helmholtz (KH) instability of a shear layer with an initially uniform magnetic field in the direction of flow is studied in the framework of 2D incompressible magnetohydrodynamics with finite resistivity and viscosity using direct numerical simulations. The shear layer evolves freely, with no external forcing, and thus broadens in time as turbulent stresses transport momentum across it. As with hydrodynamic KH, the instability here features a conjugate stable mode for every unstable mode in the absence of dissipation. Stable modes are shown to transport momentum up its gradient, shrinking the layer width whenever they exceed unstable modes in amplitude. In simulations with weak magnetic fields, the linear instability is minimally affected by the field, but enhanced small-scale fluctuations relative to the hydrodynamic case are observed. These enhanced fluctuations coincide with increased energy dissipation and faster layer broadening, with these features more pronounced in simulations with stronger fields. These trends result from the magnetic field reducing the effects of stable modes relative to the transfer of energy to small scales. As field strength increases, stable modes become less excited, thus transporting less momentum against its gradient. Furthermore, the energy that would otherwise transfer back to the driving shear because of the stable modes is instead allowed to cascade to small scales, where it is lost to dissipation. Approximations of the turbulent state in terms of a reduced set of modes are explored. While the Reynolds stress is well-described using just two modes per wavenumber at large scales, the Maxwell stress is not.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗