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.
Short duration technique providing simulation of thermodynamic properties and composition of exhaust products of liquid and solid propellant rocket engines
A large-eddy simulation (LES) of transition in plane channel flow was carried out. The LES results were compared with those of a fine direct numerical simulation (DNS), and with those of a coarse DNS that uses the same mesh as the LES, but does not use a residual stress model. While at the early stages of transition, LES and coarse DNS give the same results: the presence of the residual stress model was found to be necessary to predict accurately mean velocity and Reynolds stress profiles during the late stages of transition (after the second spike stage). The evolution of single Fourier modes is also predicted more accurately by the LES than by the DNS. As small scales are generated, the dissipative character of the residual stress starts to reproduce correctly the energy cascade. As transition progresses, the flow approaches its fully developed turbulent state, the subgrid scales tend towards equilibrium, and the model becomes more accurate.
Here, we study the application of a finite element numerical upscaling technique to the incompressible two-phase porous media total velocity formulation. Specifically, an element-agglomeration-based algebraic multigrid (AMGe) technique with improved approximation properties is used, for the first time, to generate upscaled and accurate coarse systems for the reservoir simulation equations. The upscaling technique is applied to both the mixed system for velocity and pressure and to the hyperbolic transport equations, providing fully upscaled systems. By introducing additional degrees of freedom associated with nonplanar interfaces between agglomerates, the coarse velocity space has guaranteed approximation properties. The employed AMGe technique provides coarse spaces with desirable local mass conservation and stability properties analogous to the original pair of Raviart--Thomas and piecewise discontinuous polynomial spaces, resulting in strong mass conservation for the upscaled systems. Due to the guaranteed approximation properties and the generic nature of the AMGe method, recursive multilevel upscaling is automatically obtained. Furthermore, this technique works for both structured and unstructured meshes. Multiscale mixed finite elements exhibit accuracy for general unstructured meshes but do not in general lead to nested hierarchy of spaces. Multiscale multilevel mimetic finite differences generate nested spaces but lack the adaptivity of the flux representation on coarser levels that the proposed AMGe approach offers. Thus, the proposed approach can be seen as a rigorous bridge that merges the best properties of these two existing methods. The accuracy and stability of the studied multilevel AMGe upscaling technique is demonstrated on two challenging test cases.
Efforts to improve the characteristics of fuel-oxidizer mixing in liquid rocket combustors have lead to a swirl element design for a liquid oxygen injector which is being considered for use on the STME. For the design which is the subject of this investigation, the oxygen enters the injector element perpendicular to the injector axis and nearly tangent to the circular injector wall. This swirl element is at one end of a tube and the injector exit is at the other. This geometric configuration creates a plume in the shape of a conical sheet. This sheet is either primarily contiguous liquid or droplets depending on the pressure drop in the injector and the distance from the injector exit. Probe-based devices such as two-dimensional grid patternators have been used to investigate simulated LOX injector flow fields (Hulka). The primary work described herein is an effort to use optical techniques to investigate the plume of a swirl injector element. For this investigation, a high pressure (500 psig) cold flow test facility was constructed. Water was used as the LOX simulate and air pressure was used to drive the injector flow field. Laser-induced fluorescence (LIF) from dye seeded into the water was used to obtain quantitative measurements of the time-averaged water concentration distribution in the plume. Scattered laser light and LIF were used for time averaged plume visualization and scattered light from a strobe with a 1 microsecond pulse was used for time-resolved plume visualization. During the Summer Faculty Fellowship for which this report was developed, an additional effort, unrelated to the swirl injector investigation, was made to resolve fluctuations in the combustion product composition in the exhaust of a hybrid rocket motor. A brief description of this effort is included herein.
The ability to model subsurface processes at the field-scale is essential to site closure and remedial action decisions. Field-scale models often require high-performance computing to meet the computational demands associated with complex systems, model calibration, and uncertainty quantification. eSTOMP, a parallel version of STOMP, is an analytical tool that meets the need for investigating coupled processes in the subsurface involving regional multi-fluid flow and biogeochemistry. The scalability of eSTOMP makes it ideal for simulating contaminant and geochemical reactive transport when long run times result from execution of a serial code. To demonstrate the impact of high-performance computation on subsurface analyses, the application of eSTOMP is presented for two sites. The first site involves uranium surface complexation within a highly dynamic, variably saturated flow field. The use of this model in conjunction with continued field monitoring is expected to provide a rigorous basis for developing operational strategies for field remediation and for defining defensible remediation endpoints. The second site involves the long-term dissolution of low-level radioactive glass emplaced within a near-surface repository. The waste form release rate is evaluated by modeling the basic physical and chemical processes that are known to control the waste form dissolution behavior, including pH and composition of the fluid contacting the glass. Simulations demonstrate that glass dissolution rates vary both in time and as a function of position within the disposal system, and that a 'leach rate' or radionuclide release rate parameter can only be estimated after identifying these impacts in a dynamic system. (authors)
Numerical simulation of flow over complex geometry needs accurate and efficient computational methods. Different techniques are available to handle complex geometry. The unstructured grid and multi-block body-fitted grid techniques have been widely adopted for complex geometry in engineering applications. In atmospheric applications, terrain fitted single grid techniques have found common use. Although these are very effective techniques, their implementation, coupling with the flow algorithm, and efficient parallelization of the complete method are more involved than a Cartesian grid method. The grid generation can be tedious and one needs to pay special attention in numerics to handle skewed cells for conservation purposes. Researchers have long sought for alternative methods to ease the effort involved in simulating flow over complex geometry.
Simulation of blood flows in the pulmonary artery provides some insight into certain diseases by examining the relationship between some continuum metrics, eg, the wall shear stress acting on the vascular endothelium, which responds to flow-induced mechanical forces by releasing vasodilators/constrictors. V. Kheyfets, in his previous work, studies numerically a patient-specific pulmonary circulation to show that decreasing wall shear stress is correlated with increasing pulmonary vascular impedance. We develop a scalable parallel algorithm based on domain decomposition methods to investigate an unsteady model with patient-specific pulsatile waveforms as the inlet boundary condition. The unsteady model offers tremendously more information about the dynamic behavior of the flow field, but computationally speaking, the simulation is a lot more expensive since a problem which is similar to the steady-state problem has to be solved many times, and therefore, the traditional sequential approach is not suitable anymore. We show computationally that simulations using the proposed parallel approach with up to 10 000 processor cores can be obtained with much reduced compute time. Finally, this makes the technology potentially usable for the routine study of the dynamic behavior of blood flows in the pulmonary artery, in particular, the changes of the blood flows and the wall shear stress in the spatial and temporal dimensions.
Flow-through time resolved analysis (FT-TRA) involves subjecting small mineral samples (< 10 mg) inserted in a miniature flow-through cell (50 μL) to controlled flows of eluent analyzed on-line by ICP-MS. In this study, FT-TRA is used to empirically determine the dissolution regimes for the two well-studied minerals forsterite and calcite, representing minerals with relatively slow and fast dissolution kinetics. A proportional increase in steady-state effluent [Mg, Si] concentrations with increasing flow-through cell eluent residence times confirms a dominantly surface-controlled dissolution regime for a powdered forsterite sample at pH 2.3, implying that transport limitations are negligible. In contrast, the relationship between flow rates and dissolution rates for single grain calcite samples at pH 2.3-4 reveals that transport limitations affect the rate of calcite dissolution. In order to provide a quantitative and process-based assessment of the effect of diffusive transport limitations, simulations of the calcite experiments were performed with a high resolution, pore-scale model that considers the geometry of the calcite grain and the FT-TRA flow-through reactor. The pore-scale model reproduces the observed effluent [Ca] concentrations for all experimental conditions using a single set of surface kinetic parameters, by accounting for the formation of a diffusive boundary layer (DBL) that varies in thickness as a function of flow rates. These results demonstrate that combining FT-TRA with pore-scale modeling makes it possible to obtain unprecedented insights not achievable by either method separately, including quantification of DBL thicknesses and the determination of transport controls as a function of pH, flow velocity and residence times.
Multidimensional simulations have been conducted to simulate atmospheric pressure, flat-flame/McKenna-burner-driven-flow tube experiments targeted to obtain NO x speciation data for predicting/analyzing syngas combustion emissions. In a prior work, we demonstrated the impacts of multidimensional transport on post flame region prediction departures from those assuming unidimensional flow/transport conditions. In this work, we develop and utilize a multidimensional laminar reacting flow solver to simulate the fully coupled flame and post flame regions to further elucidate the impacts of the earlier unidimensional modeling assumptions on interpreting post flame NO x experimental data. The model is used to simulate a lean, premixed syngas/air flame and its associated post flame regions within a cylindrical flow-tube-like arrangement. The combustion process takes place under atmospheric condition with trace amount of NO x seeding fed into the inlet gas stream. The spatial evolution of NO x species (NO and NO 2 ) in the flame and in the post-combustion zone suggests two distinct regions: 1) a region encompassing the flame structure itself; and 2) a post flame region in which the temperature decays due to both axial and radial transport processes. The predictions show that for the conditions studied, a pulsatile flow field exists due to the formation of an expanding and contracting recirculation zone in the outer periphery of the flow tube. By resolving the nature of the flow, the resulting time-averaged temperature and species concentrations show improved agreement with existing experimental measurements. The flow-field interaction results in radial inhomogeneities in the NO 2 profiles with the maximum concentration offset from the flow centerline. The location of the peak in NO 2 is coupled with radial temperature gradients from wall cooling effects and their significant influence on NO/NO 2 interconversion kinetics, producing notable NO 2 accumulation in regions near the wall. Geometrical configurations capable of suppressing/minimizing the pulsatile nature are also investigated and the results are compared. Other experimental configurations could be considered in parametric simulations to determine the optimal configuration that would minimize non-idealities in the observations. The work shows the value in performing such computations in advance of settling on a particular design for flow tube/flow reactor experiments.
In this paper, we developed a machine learning-based Bayesian approach to inversely quantify and reduce the uncertainties of multiphase computational fluid dynamics (MCFD) simulations for bubbly flows. The proposed approach is supported by high-resolution two-phase flow measurements, including those by double-sensor conductivity probes, high-speed imaging, and particle image velocimetry. Local distributions of key physical quantities of interest (QoIs), including the void fraction and phasic velocities, are obtained to support the Bayesian inference. In the process, the epistemic uncertainties of the closure relations are inversely quantified while the aleatory uncertainties from stochastic fluctuations of the system are evaluated based on experimental uncertainty analysis. The combined uncertainties are then propagated through the MCFD solver to obtain uncertainties of the QoIs, based on which probability-boxes are constructed for validation. The proposed approach relies on three machine learning methods: feedforward neural networks and principal component analysis for surrogate modeling, and Gaussian processes for model form uncertainty modeling. The whole process is implemented within the framework of an open-source deep learning library PyTorch with graphics processing unit (GPU) acceleration, thus ensuring the efficiency of the computation. The results demonstrate that with the support of high-resolution data, the uncertainties of MCFD simulations can be significantly reduced. The proposed approach has the potential for other applications that involve numerical models with empirical parameters.
An efficient particle simulation technique for hypersonic rarefied flows is presented at an algorithmic and implementation level. The implementation is for a vector computer architecture, specifically the Cray-2. The method models an ideal diatomic Maxwell molecule with three translational and two rotational degrees of freedom. Algorithms are designed specifically for compatibility with fine grain parallelism by reducing the number of data dependencies in the computation. By insisting on this compatibility, the method is capable of performing simulation on a much larger scale than previously possible. A two-dimensional simulation of supersonic flow over a wedge is carried out for the near-continuum limit where the gas is in equilibrium and the ideal solution can be used as a check on the accuracy of the gas model employed in the method. Also, a three-dimensional, Mach 8, rarefied flow about a finite-span flat plate at a 45 degree angle of attack was simulated. It utilized over 10 to the 7th particles carried through 400 discrete time steps in less than one hour of Cray-2 CPU time. This problem was chosen to exhibit the capability of the method in handling a large number of particles and a true three-dimensional geometry.
In nuclear thermal-hydraulic studies, counter-current flow limitation (CCFL) typically refers to steam rising at a fast rate such that it prevents coolant from draining down within a confined channel. CCFL is a crucial issue in nuclear reactor safety analysis. This study investigates CCFL in debris bed channels using high-resolution interface-capturing simulations. A novel proportional-integral-derivative flow rate controller is developed to efficiently achieve the CCFL conditions. Verification studies confirm that CCFL occurs under the same conditions with or without the controller, demonstrating that PID control ensures accurate prediction. Three debris bed channel geometries were examined: a cylindrical channel, a channel with small obstacles, and a channel with large obstacles. Results show that obstacles significantly impact flow behavior, interfacial shear, wall shear, and pressure gradients required for CCFL. Furthermore, the comparison with experimental data confirmed that simulations incorporating geometric complexities align more closely with experimental CCFL conditions. A pressure gradient correlation was also developed for CCFL prediction.
The present algorithm for the numerical simulation of flows about complex configurations (whose multiple and nonsimilar components have arbitraty geometries) employs a hybridization of the domain decomposition techniques for grid generation as well as to reduce computer-memory requirements. A fully vectorized, finite-volume, upwind-biased, approximately factored multigrid method is used to solve 3D Reynolds-averaged unsteady and compressible Navier-Stokes equations simulating supersonic flows past an ogive-nose-cylinder near or within a cavity. The time-averaged surface pressures obtained compare favorably with the wind tunnel data.
This paper presents computational results obtained with the direct simulation Monte Carlo (DSMC) method for several biconic test cases in which shock interactions and flow separation-reattachment are key features of the flow. Recent ground-based experiments have been performed for several biconic configurations, and surface heating rate and pressure measurements have been proposed for code validation studies. The present focus is to expand on the current validating activities for a relatively new DSMC code called DS2V that Bird (second author) has developed. Comparisons with experiments and other computations help clarify the agreement currently being achieved between computations and experiments and to identify the range of measurement variability of the proposed validation data when benchmarked with respect to the current computations. For the test cases with significant vibrational nonequilibrium, the effect of the vibrational energy surface accommodation on heating and other quantities is demonstrated.
This paper describes the computational work performed on the simulation of a 16-in shock-tunnel facility. The numerical problems encountered during the computation of these flows are discussed along with the validity of some approximations used, notably concerning the reduction of the problem into problems of smaller dimensionality. Quasi-1D simulations can be used to help design experiments, or to better understanding the characteristics of the facility. An application to the design of a nonintrusive diagnostic is shown. The multidimensional flow transients computed include the shock reflection at the end of the driven tube, the shock propagation down the nozzle, and the breaking of the main diaphragm.
A 'smart' wind tunnel concept for V/STOL transition regime studies is presented and performance is projected with numerical simulations. The wind tunnel would provide a flow that faithfully represents free air flow and thereby permits accurate lift calculations for V/STOL aircraft. The design concept is based on free-air tunnels used for planar wings. Modifications are necessary to account for the highly deflected, vortical wake produced in low-speed take-off and landing situations. A three-step numerical procedure covering the potential flow of a high lift system in free-air, controlled flow data with a model in a wind tunnel and addition of the wall conditions to the simulation as a feedback loop to correct the closed-tunnel is described. The initial simulation is for three-dimensional full span jet flapped wings, and is derived from extended two-dimensional calculations. The results indicate that only partial control of the tunnel flow is sufficient to serve as a valid environment for V/STOl studies.
Detached Eddy Simulation (DES) of flap side-edge flow was performed with a wing and half-span flap configuration used in previous experimental and numerical studies. The focus of the study is the unsteady flow features responsible for the production of far-field noise. The simulation was performed at a Reynolds number (based on the main wing chord) of 3.7 million. Reynolds Averaged Navier-Stokes (RANS) simulations were performed as a precursor to the DES. The results of these precursor simulations match previous experimental and RANS results closely. Although the present DES simulations have not reached statistical stationary yet, some unsteady features of the developing flap side-edge flowfield are presented. In the final paper it is expected that statistically stationary results will be presented including comparisons of surface pressure spectra with experimental data.
This paper extends the particle Green’s function approach of Gallis et al. [Phys. Fluids 13, 3482–3492 (2001)] to compute the force and heat transfer on nonspherical, sublimating particles in rarefied gas flows. Green’s functions are derived for rods and plates, which are representative of ice particle shapes, and a numerical procedure is developed for computing the Green’s functions for a general discretized convex shape. Additionally, a sublimation model is developed to simulate the evolution of ice particles in rarefied gas flows. These models are implemented in a direct simulation Monte Carlo code and demonstrated for the case of ice-particle-laden hypersonic flow over a ramp.