Linear and Nonlinear Solvers for Simulating High-Temperature Multiphase Flow Within Large-Scale Engineered Subsurface Systems.
Abstract not provided.
SEARCH · Search NASA
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.
Abstract not provided.
A systematic theoretical investigation has been carried out on the role of dust charging dynamics on the nature and stability of DIA (Dust Ion Acoustic) mode in complex plasma. The study has been made for both linear and non-linear scale regime of DIA mode. The observed results have been characterized in terms of background plasma responses towards dust surface responsible for dust charge fluctuation, invoking important dusty plasma parameters, especially the ion flow speed and dust size. The linear analyses confirm the nature of instability in DIA mode in presence of dust charge fluctuation. The instability shows a damping of DIA mode in subsonic flow regime followed by a gradual growth in instability in supersonic limit of ion flow. The strength of non-linearity and their existence domain is found to be driven by different dusty plasma parameters. As dust is ubiquitous in interstellar medium with plasma background, the study also addresses the possible effect of dust charging dynamics in gravito-electrostatic characterization and the stability of dust molecular clouds especially in proto-planetary disc. The observations are influential and interesting towards the understanding of dust settling mechanism and formation of dust environments in different regions in space.
Optimal design of hybrid renewable mini-grids requires both economic and power quality considerations. Existing modeling approaches address these considerations via separate or loosely coupled models. Here, we extend REopt—a techno-economic optimization model developed at the National Renewable Energy Laboratory—to consider both within a single model. REopt formulates the design problem as a mixed-integer linear program that solves for a site's optimal technology mix, sizing, and operation to minimize life cycle cost. REopt has traditionally assumed a single node system. In the work presented here, we expand the REopt platform to consider multiple connected nodes with associated voltage constraints. In order to do this, we model power flow using a fixed-point linear approximation method. Additionally, we then use the model to explore design considerations of mini-grids in Sub-Saharan Africa. Specifically, we evaluate under what combinations of transmission line distance and capacity it is technically viable and economically preferable to build multiple isolated mini-grids versus an interconnected, centralized system.
A weakly nonlinear wave-interaction theory is used to explain the decay of linearly unstable oblique waves in rotating channel flow at low Reynolds number. At the low Reynolds number studied, the waves tend to stabilize each other and therefore when no waves are preferentially excited, the linearly most unstable mode reaches a nonlinear amplitude first and causes the stabilization of all competing modes.
Two approaches for the identification of internal gravity waves in sheared and unsheared homogeneous stratified turbulence are investigated. First, the phase angle between the vertical velocity and density fluctuations is considered. It was found, however, that a continuous distribution of the phase angle is present in weakly and strongly stratified flow. Second, a projection onto the solution of the linearized inviscid equations of motion of unsheared stratified flow is investigated. It was found that a solution of the fully nonlinear viscous Navier-Stokes equations can be represented by the linearized inviscid solution. The projection yields a decomposition into vertical wave modes and horizontal vortical modes.
Two approaches for the identification of internal gravity waves in sheared and unsheared homogeneous stratified turbulence are investigated. First, the phase angle between the vertical velocity and density fluctuations is considered. It is found, however, that a continuous distribution of the phase angle is present in weakly and strongly stratified flow. Second, a projection onto the solution of the linearized inviscid equations of motion of unsheared stratified flow is investigated. It is found that a solution of the fully nonlinear viscous Navier-Stokes equations can be represented by the linearized inviscid solution. The projection yields a decomposition into vertical wave modes and horizontal vortical modes.
Isothermal, pressure driven rarefied gas flows through long channels with octagonal cross section shapes are analyzed computationally. The capillary is between inlet and outlet reservoirs. The cross section is constant along the axial direction. The boundary condition at the solid-gas interface is assumed to be diffuse reflection. Since the channel is long, the gaseous velocity is small compared to the average molecular speed. Consequently, a linearized description can be used. The flow is described by the linearized Bhatnagar-Gross-Krook kinetic model. The solution of the problem is divided into two stages. First, the local flow field is determined by assuming the local pressure gradient. Secondly, the global flow behavior is deduced by the consideration of the conservation of the mass along the axis of the capillary. The kinetic equation is solved by the discrete velocity method on the cross section. Both spatial and velocity spaces are discretized. A body fitted rectangular grid is used for the spatial space. Near the boundary, first-order, while in the interior part of the flow domain, second-order finite-differences are applied to approximate the spatial derivatives. This combination results into an efficient and straightforward numerical treatment. The velocity space is represented by a Gauss-Legendre quadrature. The kinetic equation is solved in an iterative manner. The local dimensionless flow rate is calculated and tabulated for a wide range of the gaseous rarefaction for octagonal cross sections with various geometrical parameters. It exhibits the Knudsen minimum phenomenon. The flow rates in the octagonal channel are compared to those through capillaries with circular and square cross sections. Typical velocity profiles are also shown. The mass flow rate and the distribution of the pressure are determined and presented for global pressure driven flows.
An investigation is conducted into the determination of the credibility of interacting boundary layers in predicting compressible subsonic flows over smooth surface imperfections. The case of smooth backward-facing steps is considered. The predicted mean flows are compared with those obtained using a Navier-Stokes solver. Moreover, the linear 2-D compressible stability characteristics of both mean flows are compared. The results show that the interacting boundary-layer formulation produces accurate mean flows that yield accurate linear stability characteristics, such as growth rates and amplification factors.
An investigation is conducted into the determination of the credibility of interacting boundary layers in predicting compressible subsonic flows over smooth surface imperfections. The case of smooth backward-facing steps is considered. The predicted mean flows are compared with those obtained using a Navier-Stokes solver. Moreover, the linear 2-D compressible stability characteristics of both mean flows are compared. The results show that the interacting boundary-layer formulation produces accurate mean flows that yield accurate linear stability characteristics, such as growth rates and amplification factors.
The linearized problem of free-molecular gas flow in a long elliptic channel with a circular cylindrical element inside has been solved in the kinetic approximation. The flow in the channel is caused by temperature and pressure drops between its ends. The Boltzmann kinetic equation for collisionless gas has been used as a basic equation, and the boundary condition has been set in terms of the diffuse reflection model. The distribution of the mass velocity of the gas over the cross section of the channel has been obtained. The mass flow rate of the gas in the channel versus the temperature and pressure drops between its ends has been calculated. It has been found that the mass flow of the gas substantially depends on the radius of the inner cylinder.
A novel reduced model for ion temperature gradient (ITG) turbulent transport in helical plasmas is presented. The model enables one to predict nonlinear gyrokinetic simulation results from linear gyrokinetic analyses. It is shown from nonlinear gyrokinetic simulations of the ITG turbulence in helical plasmas that the transport coefficient can be expressed as a function of the turbulent fluctuation level and the averaged zonal flow amplitude. Then, the reduced model for the turbulent ion heat diffusivity is derived by representing the nonlinear turbulent fluctuations and zonal flow amplitude in terms of the linear growth rate of the ITG instability and the linear response of the zonal flow potentials. It is confirmed that the reduced transport model is in a good agreement with nonlinear gyrokinetic simulation results for high ion temperature plasmas in the large helical device.
A Galerkin-weighted residuals formulation is employed to establish an implicit finite element solution algorithm for generally nonlinear initial-boundary value problems. Solution accuracy, and convergence rate with discretization refinement, are quantized in several error norms, by a systematic study of numerical solutions to several nonlinear parabolic and a hyperbolic partial differential equation characteristic of the equations governing fluid flows. Solutions are generated using selective linear, quadratic and cubic basis functions. Richardson extrapolation is employed to generate a higher-order accurate solution to facilitate isolation of truncation error in all norms. Extension of the mathematical theory underlying accuracy and convergence concepts for linear elliptic equations is predicted for equations characteristic of laminar and turbulent fluid flows at nonmodest Reynolds number. The nondiagonal initial-value matrix structure introduced by the finite element theory is determined intrinsic to improved solution accuracy and convergence. A factored Jacobian iteration algorithm is derived and evaluated to yield a consequential reduction in both computer storage and execution CPU requirements while retaining solution accuracy.
An analysis has been developed to predict the unsteady force and moment generated by the passage of a timewise periodic total pressure distortion through an arbitrary cascade of airfoils. The mathematical formulation of this analysis is based on the assumption that the magnitudes of the timewise fluctuations of the variables which describe the flow field are small compared to their time average values. This assumption permits the development of a linear unsteady perturbation analysis about a steady flow field. In addition to this linearization assumption the fluid medium is assumed to be incompressible and inviscid. The mathematical development begins by decomposing the velocity field surrounding an infinite cascade of airfoils into its irrotational and rotational components. The rotational component is associated with an upstream unsteady total pressure distortion and is defined in terms of the vorticity field associated with the distortion pattern. The irrotational component is further decomposed into a steady and unsteady part. A combined analytical and numerical procedure has been developed to solve the field equations which govern the rotational and irrotational velocity fields. Results of this analysis show a strong influence of mean loading on the unsteady force generated by the passage of a one dimensional gust through a cascade of compressor blades.
Background: Assessment of infarct size after myocardial infarction is predictive of subsequent morphological changes and clinical outcome. This study aimed to assess subacute post-intervention Rubidium-82 ({sup 82}Rb)-PET imaging in predicting left ventricle ejection fraction, regional wall motion, and final infarct size by CMR at 3-months after STEMI. Methods: STEMI patients undergoing percutaneous coronary intervention were included prospectively. Rest-only {sup 82}Rb-PET perfusion imaging was performed at median 36 hours [IQR: 22 to 50] after the treatment. The extent of hypoperfusion and absolute blood flow (mL·min·g) were estimated on a global and a 17-segment model with dedicated software. At 3-months follow-up patients completed the CMR functional and late gadolinium enhancement imaging. Results: 42 patients were included, but only 35 had follow-up CMR and constituted the study population. Absolute blood flow was significantly lower in the infarct-related territory compared to remote myocardium, P < .005. Extent of hypoperfusion correlated with final infarct size, r = 0.58, P < .001, while blood flow correlated with ejection fraction, r = 0.41, P < .05. In linear mixed models, higher subacute absolute blood flow (β = 4.6, confidence interval [3.5; 5.2], P < .001, R{sup 2} = 0.67) was associated with greater wall motion. Segmental extent of subacute hypoperfusion (β = 0.43 [0.38; 0.49], P < .001, R{sup 2} = 0.58) was associated with the degree of late gadolinium enhancement at 3-months. Conclusions: Subacute rest-only {sup 82}Rb-PET is feasible following STEMI and seems predictive of myocardial function and infarct size at 3-months.
A ground test was performed to determine the effects of compressor bleed flow extraction on the performance of F404-GE-400 afterburning turbofan engines. The two engines were installed in the F/A-18 High Alpha Research Vehicle at the NASA Dryden Flight Research Facility. A specialized bleed ducting system was installed onto the aircraft to control and measure engine bleed airflow while the aircraft was tied down to a thrust measuring stand. The test was conducted on each engine and at various power settings. The bleed air extraction levels analyzed included flow rates above the manufacturer's maximum specification limit. The measured relationship between thrust and bleed flow extraction was shown to be essentially linear at all power settings with an increase in bleed flow causing a corresponding decrease in thrust. A comparison with the F404-GE-400 steady-state engine simulation showed the estimation to be within +/- 1 percent of measured thrust losses for large increases in bleed flow rate.
A ground test was performed to determine the effects of compressor bleed flow extraction on the performance of F404-GE-400 afterburning turbofan engines. The two engines were installed in the F/A-18 High Alpha Research Vehicle at the NASA Dryden Flight Research Facility. A specialized bleed ducting system was installed onto the aircraft to control and measure engine bleed airflow while the aircraft was tied down to a thrust measuring stand. The test was conducted on each engine and at various power settings. The bleed air extraction levels analyzed included flow rates above the manufacturer's maximum specification limit. The measured relationship between thrust and bleed flow extraction was shown to be essentially linear at all power settings with an increase in bleed flow causing a corresponding decrease in thrust. A comparison with the F404-GE-400 steady-state engine simulation showed the estimation to be within +/- 1 percent of measured thrust losses for large increases is bleed flow rate.
The inviscid, linear, nonaxisymmetric, temporal stability of the boundary layer associated with the supersonic flow past axisymmetric bodies (with particular emphasis on long thin, straight circular cylinders), subject to heated or cooled wall conditions is investigated. The eigenvalue problem is computed in some detail for a particular Mach number or 3.8, revealing that the effect of curvature and the choice of wall conditions both have a significant effect on the stability of the flow. Both the asymptotic, large azimuthal wavenumber solution and the asymptotic, far downstream solution are obtained for the stability analysis and compared with numerical results. Additionally, asymptotic analyses valid for large radii of curvature with cooled/heated wall conditions, are presented. In general, important differences were found to exist between the wall temperature conditions imposed and the adiabatic wall conditions considered previously.
The inviscid, linear, nonaxisymmetric, temporal stability of the boundary layer associated with the supersonic flow past axisymmetric bodies (with particular emphasis on long thin, straight circular cylinders), subject to heated or cooled wall conditions is investigated. The eigenvalue problem is computed in some detail for a particular Mach number or 3.8, revealing that the effect of curvature and the choice of wall conditions both have a significant effect on the stability of the flow. Both the asymptotic, large azimuthal wavenumber solution and the asymptotic, far downstream solution are obtained for the stability analysis and compared with numerical results. Additionally, asymptotic analyses valid for large radii of curvature with cooled/heated wall conditions are presented. In general, important differences were found to exist between the wall temperature conditions imposed and the adiabatic wall conditions considered previously.