Study of radiative heat transfer and flow physics from medium-scale methanol pool fire simulations
Not Available
SEARCH · Engineering Papers
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.
Not Available
Explore the source record for details and available documents.
The issues and components supporting the engine power stream are reviewed. It is essential that companies pay close attention to engine sealing issues, particularly on the high-pressure spool or high-pressure pumps. Small changes in these systems are reflected throughout the entire engine. Although cavity, platform, and tip sealing are complex and have a significant effect on component and engine performance, computational tools (e.g., NASA-developed INDSEAL, SCISEAL, and ADPAC) are available to help guide the designer and the experimenter. Gas turbine engine and rocket engine externals must all function efficiently with a high degree of reliability in order for the engine to run but often receive little attention until they malfunction. Within the open literature statistically significant data for critical engine components are virtually nonexistent; the classic approach is deterministic. Studies show that variations with loading can have a significant effect on component performance and life. Without validation data they are just studies. These variations and deficits in statistical databases require immediate attention.
Explore the source record for details and available documents.
The combination of interdigitated flow fields (IDFF) with porous electrodes offers lower pressure drop and better performance than conventional flow‐through porous electrodes in redox flow batteries. Comprehensive three‐dimensional and simplified one‐dimensional + two‐dimensional models describing flow uniformity and pressure losses within flow through, parallel, and interdigitated flow fields were developed and used to demonstrate the benefits of IDFF. Analytical solutions for IDFF that compare favorably with computational fluid dynamics quantify the trade between pressure loss and velocity maldistribution both along the channels and within the electrode. © 2018 American Institute of Chemical Engineers AIChE J , 64: 3746–3755, 2018
Dissolution trapping is one of the most important mechanisms for geological carbon storage (GCS). Recent laboratory and field experiments have shown non-equilibrium dissolution of supercritical CO2 (scCO2) and coupled scCO2 dissolution and water flow, i.e., scCO2 dissolution at local pores/pore throats creating new water-flow paths, which in turn enhance dissolution by increased advection and interfacial area. However, the impacts of pore-scale characteristics on these coupled processes have not been investigated. In this study, imbibition and dissolution experiments were conducted under 40°C and 9?MPa using a homogeneous/isotropic hexagonal micromodel, two homogeneous elliptical micromodels with low or high anisotropy, and a heterogeneous sandstone-analog micromodel. The four micromodels, initially saturated with deionized (DI)-water, were drained by injecting scCO2 to establish a stable scCO2 saturation. DI water was then injected at different rates with logCa (the capillary number) ranging from -6.56 to -4.34. Results show that bypass of scCO2 by displacing water is the dominant mechanism contributing to the residual CO2 trapping, triggered by heterogeneity in pore characteristics or pore-scale scCO2-water distribution. Bypass can be enhanced by pore heterogeneity or reduced by increasing transverse permeability, resulting in relatively low (<2% of CO2 solubility) or high (9–13% of CO2 solubility) dissolved CO2 concentration in displacing water. The overall dissolution of residual scCO2 increases with decreasing Ca, and approaches to their solubility at low Ca value with sufficient residence time. This main trend is similar to a capillary desaturation curve that represents the relationship between the residual saturation and Ca. Spatially, dissolution initiates along the boundary of bypassed scCO2 cluster(s) in a non-equilibrium manner, and the coupling of water flow and dissolution occurs which fragments the bypassed scCO2 clusters and enhance scCO2 dissolution.
A new direct numerical simulation (DNS) code for multi-component gaseous reacting flows has been developed at KAUST, with the state-of-the-art programming model for next generation high performance computing platforms. The code, named KAUST Adaptive Reacting Flows Solver (KARFS), employs the MPI+X programming, and relies on Kokkos for “X” for performance portability to multi-core, many-core and GPUs, providing innovative software development while maintaining backward compatibility with established parallel models and legacy code. The capability and potential of KARFS to perform DNS of reacting flows with large, detailed reaction mechanisms is demonstrated with various model problems involving ignition and turbulent flame propagations with varying degrees of chemical complexities.
Water exchange between the surface and subsurface is important for both water resource management and environmental protection. In this paper, we develop coupled surface and subsurface flow simulation capability in a parallel subsurface flow and reactive transport code PFLOTRAN. We sequentially couple the diffusion wave-based surface flow with the subsurface flow governedby the Richards equation in PFLOTRAN. These two flow domains are linked with a boundary condition switching method that ensures continuity of pressure and flux at the surface-subsurface interface. We verify the coupled code against other existing hydrologic models and observation data using a number of numerical experiments. The coupled hydrological model exhibits good performance in strong parallel scaling tests. The new coupled surface and subsurface simulator significantly advance community simulation capability towards improving integrated hydrologic and biogeochemical understanding of complex systems such as watersheds and river corridors. Keywords: Surface flow, Integrated hydrological modeling, Boundary condition switching, Parallel computing
We investigate the development of the directed, v 1 , and elliptic flow, v 2 , in heavy ion collisions in mid-central Au+Au reactions at E lab = 1.23A GeV. We demonstrate that the elliptic flow of hot and dense matter is initially positive (v 2 >0) due to the early pressure gradient. This positive v 2 transfers its momentum to the spectators, which leads to the creation of the directed flow v 1 . In turn, the spectator shadowing of the in-plane expansion leads to a preferred decoupling of hadrons in the out-of-plane direction and results in a negative v 2 for the observable final state hadrons. We propose a measurement of v 1 −v 2 flow correlations and of the elliptic flow of dileptons as methods to pin down this evolution pattern. The elliptic flow of the dileptons allows then to determine the early-state EoS more precisely, because it avoids the strong modifications of the momentum distribution due to shadowing seen in the protons. This opens the unique opportunity for the HADES and CBM collaborations to measure the Equation-of-State directly at 2-3 times nuclear saturation density.
We propose a fully spatio-temporal approach for identifying spatially varying modes of oscillation in fluid dynamics simulation output by means of multitaper frequency–wavenumber spectral analysis. One-dimensional spectrum estimation has proven to be a valuable tool in the analysis of turbulence data applied spatially to determine the rate of energy transport between spatial scales, or temporally to determine frequencies of oscillatory flows. It also allows for the quantitative comparison of flow characteristics between two scenarios using a standard basis. It has the limitation, however, that it neglects coupling between spatial and temporal structures. Two-dimensional frequency–wavenumber spectral analysis allows one to decompose waveforms into standing or travelling variety. The extended higher-dimensional multitaper method proposed here is shown to have improved statistical properties over conventional non-parametric spectral estimators, and is accompanied by confidence intervals which estimate their uncertainty. Multitaper frequency–wavenumber analysis is applied to a canonical benchmark problem, namely, a direct numerical simulation of von Kármán vortex shedding off a square wall-mounted cylinder with two inflow scenarios with matching momentum-thickness Reynolds numbers $Re_{\unicode[STIX]{x1D703}}\approx 1000$ at the obstacle. Frequency–wavenumber analysis of a two-dimensional section of these data reveals that although both the laminar and turbulent inflow scenarios show a turbulent $-5/3$ cascade in wavenumber ( $\unicode[STIX]{x1D708}$ ) and frequency ( $f$ ), the flow characteristics differ in that there is a significantly more prominent discrete harmonic oscillation near $(f,\unicode[STIX]{x1D708})=(0.2,0.21)$ in wavenumber and frequency in the laminar inflow scenario than the turbulent scenario. This frequency–wavenumber pair corresponds to a travelling wave with velocity near one near the centre path of the vortex street.
Abstract Multiscale heterogeneity and insufficient characterization data for a specific subsurface formation of interest render predictions of multi‐phase fluid flow in geologic formations highly uncertain. Quantification of the uncertainty propagation from the geomodel to the fluid‐flow response is typically done within a probabilistic framework. This task is computationally demanding due to, for example, the slow convergence of Monte Carlo simulations (MCS), especially when computing the tails of a distribution that are necessary for risk assessment and decision‐making under uncertainty. The frozen streamlines method (FROST) accelerates probabilistic predictions of immiscible two‐phase fluid flow problems; however, FROST relies on MCS to compute the travel‐time distribution, which is then used to perform the transport (phase saturation) computations. To alleviate this computational bottleneck, we replace MCS with a deterministic equation for the cumulative distribution function (CDF) of travel time. The resulting CDF‐FROST approach yields the CDF of the saturation field without resorting to sampling‐based strategies. Our numerical experiments demonstrate the high accuracy of CDF‐FROST in computing the CDFs of both saturation and travel time. For the same accuracy, it is about 5 and 10 times faster than FROST and MCS, respectively.
Abstract Environmental water allocation in California is a complex legal process involving various government agencies and stakeholders. E-flow requirements can be based on annual runoff typologies called water year types (WYTs), which dictate water volume, timing, and duration. In this study, we examined hydropower licensing documents of the major water and power projects in the Central Sierra Nevada to catalog e-flow requirements by WYT. In this study case, we identify how WYT classification systems and categories vary across and within different basins. Additionally, we assessed the impacts of climate change on hydrology, the frequency of WYTs identified, and the reliability and resilience of e-flows using future projections (2031–2060) of 10 Global Circulation Models (GCMs). We then propose a potential adaptation strategy using a 30 year moving percentiles approach to recalculate WYTs. We identified eight WYT classifications systems were identified, and their WYT distributions statistically significantly changes across all GCMs, even though most GCMs indicate no statistically significant change in hydrology. Disparities in future impacts are observed among and within hydropower projects, with some river reaches showing negative impacts on reliability and resilience. The adaptation strategy can generally boost resilience and improve reliability, but simply updating existing WYT thresholds without flexible regulatory frameworks reconsidering WYTs and e-flows thresholds, may not yield substantial improvements. Challenges in managing e-flows in California within regulatory and hydroclimatic contexts are intricate due to the lack of standardized approaches, leading to inconsistencies and potential conflicts among stakeholders, that will likely be exacerbated by climate change. Thus, we emphasize that targeted, site-specific, and adaptive management strategies are crucial, besides the need for a harmonized and consistent approach to defining and applying WYT categories and methods and/or e-flow assessments.
Particle-to-supercritical carbon dioxide (sCO 2 ) heat exchanger is a critical component in next-generation concentrating solar power (CSP) plants. The inherently low heat transfer between falling particles and sCO 2 imposes a challenge toward economic justification of levelized cost of electricity produced through solar energy. Introduction of integrated porous media with the walls bounding particle flow has the potential to enhance the overall particle-to-sCO 2 heat exchanger performance. Here, this paper presents an experimental study on heat transfer characterization of additively manufactured lattice frame material based on Octet-shaped unit cell with particles and air as working fluids. The lattice structures were additively manufactured in stainless steel (SS) 316L and SS420 (with 40% bronze infiltration) via Binder jetting process, where the lattice porosities were varied between 0.75 and 0.9. The mean particle diameters were varied from 266 μm to 966 μm. The effective thermal conductivity and averaged heat transfer coefficient were determined through steady-state experiments. It was found that the presence of lattice enhances the effective thermal conductivity by 2–4 times when compared to packed bed of particles alone. Furthermore, for gravity-assisted particle flow through lattice panel, significantly high convective heat transfer coefficients ranging from 200 W/m2K to 400 W/m2K were obtained for the range of particle diameters tested. The superior thermal transport properties of Octet-shape-based lattice frame for particle flow makes it a very promising candidate for particle-to-sCO 2 heat exchanger for CSP application.
We are investigating a microwave resonant cavity transducer for flow sensing in the vessel of a high temperature fluid advanced reactor (AR), such as a molten salt cooled reactor (MSCR) or a sodium fast reactor (SFR). This transducer is a hollow metallic cylindrical cavity, with one of the flat walls of the cylinder flexible enough to undergo microscopic deflection due to dynamic fluid pressure. Membrane deflection leads to a shift in the microwave resonant frequency, which can be detected with a spectrum analyzer. We have developed a continuum electromechanics model of the microwave resonant cavity transducer performance. The deflection of the membrane is calculated using equations from the theory of plates and shells for deflection of radially constrained thin circular disk under uniform load. The microwave frequency shift is calculated using the equations from the microwave cavity perturbation theory. The transducer model was validated with proof-of-concept experimental data of water flow sensing with a Brass 360 right circular cylinder with 22.2mm diameter and 203µm thick wall. The cavity was excited through a subwavelength hole in TE 011 mode with resonant frequency f ≈ 17.8GHz. Subsequently, we have performed a preliminary proof-of-principle test of flow sensing in high temperature liquid sodium in environment. For this test, we have developed a cylindrical resonator with the same dimensions as for the Brass 360 cavity. The resonator for high temperature liquid sodium test was machined from stainless steel 316 and electroplated with silver on the interior surfaces. We have also developed and insertion probe consisting of a 50cm WR-42 brass waveguide enclosed in a protective SS316 tube. The cavity was excited through a subwavelength hole on the side of the wall of the cylinder in the TE 011 mode with resonant frequency f ≈ 17.8GHz. The liquid sodium setup consists of a cylindrical vessel with a center feed line, where a transducer inserted through the top cap of the vessel measures velocity of the impinging liquid jet. The sodium flow sensing study was performed in a liquid sodium vessel at 340°C temperature and ambient pressure. The flow rate was changed by varying sodium pump power, and a frequency shift of several hundred kHz was observed. The transducer was removed after 70 days of testing. No visual damage was observed, and the electromagnetic response remained the same. As a calibration experiment, we have assembled a water vessel with a center feed, and with dimensions similar to those of the liquid sodium setup. Water flow sensing was performed before and after liquid sodium test. The responses in both cases are similar and agree with COMSOL computer simulations. Because long-term (multi-year) experimental tests of transducer resilience to harsh environment are not practical, we have developed a probabilistic model of creep to estimate transducer resilience to the harsh environment. The probabilistic model considers diffusion creep under the condition of high temperature and low stress, where the stress and temperature are allowed to be random variables with Gaussian distributions. Using the probabilistic model, we estimate inelastic membrane deflections due to creep for several temperature ranges. We conclude that for temperatures less than 650°C, creep has negligible long-term effect on the transducer performance.
The $O(3)$ non-linear sigma model (NLSM) is a prototypical field theory for QCD and ferromagnetism, featuring topological qualities. Though the topological susceptibility should vanish in physical theories, lattice simulations of the NLSM find that the topological susceptibility diverges in the continuum limit. We introduce a $\theta$-term and analyze the topological charge as a function of $\theta$ under the gradient flow. Our results show that divergence persists in the presence of the flow, even at non-zero $\theta$.
A visualization study of the flow through a three ft dia model of a four bladed inducer, which is operated in air at a flow coefficient of 0.065, is reported in this paper. The flow near the blade surfaces, inside the rotating passages, downstream and upstream of the inducer is visualized by means of smoke, tufts, ammonia filament, and lampblack techniques. Flow is found to be highly three dimensional, with appreciable radial velocity throughout the entire passage. The secondary flows observed near the hub and annulus walls agree with qualitative predictions obtained from the inviscid secondary flow theory.
A computer program written in conjunction with the numerical solution of the flow of chemically reacting gas-particle mixtures was documented. The solution to the set of governing equations was obtained by utilizing the method of characteristics. The equations cast in characteristic form were shown to be formally the same for ideal, frozen, chemical equilibrium and chemical non-equilibrium reacting gas mixtures. The characteristic directions for the gas-particle system are found to be the conventional gas Mach lines, the gas streamlines and the particle streamlines. The basic mesh construction for the flow solution is along streamlines and normals to the streamlines for axisymmetric or two-dimensional flow. The analysis gives detailed information of the supersonic flow and provides for a continuous solution of the nozzle and exhaust plume flow fields. Boundary conditions for the flow solution are either the nozzle wall or the exhaust plume boundary.
A flow model has been developed to study the flow development after reattachment with supersonic external streams. Special attention is given to the pressure difference across the viscous layer, and it is suggested that such a flow redevelopment can be treated as a relaxation of this pressure difference. Upon correlating the pressure difference with a slope parameter of the velocity profile, the system of equations governing the flow would produce a saddle point singularity corresponding to the fully rehabilitated asymptotic flow condition. A method of calculation for this flowfield, in conjunction with the matching of the upstream flow, has been derived and is discussed. Samples of calculations are also presented. Reasonably good agreement with experimental data has also been observed.