Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “interfacial velocity”

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.

169 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

High-speed impact of micron-sized diesel drop trains—Splashing dynamics, secondary droplet formation, and effects of pre-existing film thickness

This work investigates surface impingement of monodispersed diesel drop trains using computational simulations. The three-dimensional simulations include high-speed impact of micrometer-sized fuel drops onto initially dry and wetted stainless steel substrates. The drop size and impact velocity are representative of fuel injection conditions in internal combustion engines (ICEs). The drop trains serve as a simplified representation of fuel spray. To accurately capture the film that naturally develops on initially dry substrates, a dynamic diesel-stainless steel contact angle model was developed using experiments of single drop impact under ICE representative temperature. Simulations are presented for a highly splashing case and the effects of a pre-existing film on splashing dynamics were investigated, including the temporal evolution of splashed mass and film thickness. It was concluded that for thin films, the effect of pre-existing film thickness is less significant than in single drop impingement. Secondary droplet characterization was performed on simulation results, leading to instantaneous and time-averaged distributions of secondary droplet size, velocity magnitude, and trajectory angle. It was found that for each drop impingement, approximately 58% of the splashed mass is from the impinging drop itself, while the remainder of the splashed mass is composed of film liquid. From a high resolution case, details of secondary droplet formation are observed and three distinct phases of secondary droplet formation are identified. The detailed analysis of drop train impingement under engine-relevant conditions serves as a first step toward a robust understanding of fuel impingement and the development of cleaner and more efficient ICEs.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles

Heat transfer in directly-irradiated high-temperature solid–gas flows laden with polydisperse particles is investigated using a novel transient three-dimensional computational fluid dynamics model. The model couples particle–gas hydrodynamics of solid–gas flows laden with polydisperse particles, radiative heat transfer in non-grey absorbing, emitting and anisotropically-scattering multi-component participating media, conduction heat transfer in the gas phase, and interfacial convection heat transfer. The multiphase particle-in-cell method is used to predict high-fidelity solid–gas flow characteristics, such as the local discrete particle size distribution, with increased computational efficiency by combining the advantages of both Eulerian and Lagrangian methods. The multi-component radiative transfer model is implemented using an advanced collision-based Monte Carlo ray-tracing method. The number of the prescribed discrete particle components is found to be the key parameter affecting the computational accuracy and efficiency, which primarily depends on the size distribution of the particles. For the model particle–gas flow featuring free-falling Gamma-distributed ceramic particles exposed to concentrated solar irradiation, the particle volume fraction, radiative, fluid flow and thermal characteristics appear to converge with the increasing number of the discrete particle components. Five particle components are sufficient to obtain physically meaningful results. A further increase in the number of the particle components only slightly increases the accuracy of the numerical predictions at the expense of a rapidly increasing computational time. For five particle components, the particle vertical velocity at the receiver exit for particles with the diameter of 43.4 μm is 57% of that for the particles with the diameter of 202.8 μm. The temperatures of these two particle components increase from the initial ambient values by factors of 2 and 1.2, respectively, during the simulation time. In conclusion, the model developed allows for increased fidelity of particle–gas flow simulations with significant radiative effects.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Correlation of Early-Stage Growth Process Conditions with Dislocation Evolution in MOCVD-Based GaP/Si Heteroepitaxy

To identify the complex relationships between early-stage growth processes and the resultant defect microstructure in GaP/Si heteroepitaxy, a holistic study of several key metal-organic chemical vapor deposition (MOCVD) parameters was conducted, focusing on Si surface preparation and GaP atomic layer epitaxy (ALE) based nucleation processes. Here, crystalline defects related to the lattice mismatch and/or interfacial heterovalency, namely misfit dislocations (MD), threading dislocations (TD), and stacking fault pyramids (SFP), were quantitatively characterized via electron channeling contrast imaging (ECCI) and correlated against the different process variations. Choice of Si surface preparation method between the two examined (dilute SiH4 annealing versus Si2H6 based homoepitaxy) had little impact on resultant GaP film morphology and defect content, whereas differing GaP ALE nucleation conditions produced much more substantial changes. In particular, the initial precursor species (tert-butylphosphine versus triethygallium) and ALE cycle purge times both yielded significant influence over threading dislocation densities (TDD) in thin (100 nm), post-critical thickness GaP/Si films, with TDD spanning two orders of magnitude, from 6.7×10 7 cm -2 to 7.1×10 5 cm -2 , depending on the specific process conditions employed. SFP densities were also found to follow a similar trend, ranging from 2.0×10 7 cm -2 to 1.8×10 5 cm -2 , but with no apparent causal relationship between SFP density and TDD. To help explain the dramatic differences observed, detailed, large-area MD network characterization was used to provide statistically-relevant quantitative analyses of the critical dislocation dynamics (introduction rates and glide velocities) associated with the different process variants. These extracted values are then correlated against the ALE process variants to provide insight into the potential mechanistic roles of the different growth processes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Pore-Scale Dynamics of Liquid CO2–Water Displacement in 2D Axisymmetric Porous Micromodels Under Strong Drainage and Weak Imbibition Conditions: High-Speed μPIV Measurements

Resolving pore-scale transient flow dynamics is crucial to understanding the physics underlying multiphase flow in porous media and informing large-scale predictive models. Surface properties of the porous matrix play an important role in controlling such physics, yet interfacial mechanisms remain poorly understood, in part due to a lack of direct observations. This study reports on an experimental investigation of the pore-scale flow dynamics of liquid CO 2 and water in two-dimensional (2D) circular porous micromodels with different surface characteristics employing high-speed microscopic particle image velocimetry (μPIV). The design of the micromodel minimized side boundary effects due to the limited size of the domain. The high-speed μPIV technique resolved the spatial and temporal dynamics of multiphase flow of CO 2 and water under reservoir-relevant conditions, for both drainage and imbibition scenarios. When CO 2 displaced water in a hydrophilic micromodel (i.e., drainage), unstable capillary fingering occurred and the pore flow was dominated by successive pore-scale burst events (i.e., Haines jumps). When the same experiment was repeated in a nearly neutral wetting micromodel (i.e., weak imbibition), flow instability and fluctuations were virtually eliminated, leading to a more compact displacement pattern. Energy balance analysis indicates that the conversion efficiency between surface energy and external work is less than 30%, and that kinetic energy is a disproportionately smaller contributor to the energy budget. This is true even during a Haines jump event, which induces velocities typically two orders of magnitude higher than the bulk velocity. These novel measurements further enabled direct observations of the meniscus displacement, revealing a significant alteration of the pore filling mechanisms during drainage and imbibition. While the former typically featured burst events, which often occur only at one of the several throats connecting a pore, the latter is typically dominated by a cooperative filling mechanism involving simultaneous invasion of a pore from multiple throats. This cooperative filling mechanism leads to merging of two interfaces and releases surface energy, causing instantaneous high-speed events that are similar, yet fundamentally different from, burst events. Finally, pore-scale velocity fields were statistically analyzed to provide a quantitative measure of the role of capillary effects in these pore flows.

Li, Yaofa↗
NASA NTRS1994

Convective flow effects on protein crystal growth

A high-resolution microscopic interferometric setup for the monitoring of protein morphologies has been developed. Growth or dissolution of a crystal can be resolved with a long-term depth resolution of 200 A and a lateral resolution of 2 microns. This capability of simultaneously monitoring the interfacial displacement with high local depth resolution has yielded several novel results. We have found with lysozyme that (1) the normal growth rate is oscillatory, and (2) depending on the impurity content of the solution, the growth step density is either greater or lower at the periphery of a facet than in its center. The repartitioning of Na plus and Cl minus ions between lysozyme solutions and crystals was studied for a wide range of crystallization conditions. A nucleation-growth-repartitioning model was developed, to interpret the large body of data in unified way. The results strongly suggest that (1) the ion to lysozyne ratio in the crystal depends mostly on kinetic rather than crystallographic parameters, and (2) lysozyme crystals possess a salt-rich core with a diameter electron microscopy results appear to confirm this finding, which could have far-reaching consequences for x-ray diffraction studies. A computational model for diffusive-convective transport in protein crystallization has been applied to a realistic growth cell geometry, taking into account the findings of the above repartitioning studies and our kinetics data for the growth of lysozyme. The results show that even in the small cell employed, protein concentration nonuniformities and gravity-driven solutal convection can be significant. The calculated convection velocities are of the same order to magnitude as those found in earlier experiments. As expected, convective transport, i.e., at Og, lysozyme crystal growth remains kinetically limited. The salt distribution in the crystal is predicted to be non-uniform at both 1g and 0g, as a consequence of protein depletion in the solution. Static and dynamic light scattering studies in undersaturated and supersaturated solutions have been performed. Diffusivities in undersaturated solutions, were found to vary with lysozyme concentrations. Depending on the salt concentration, the diffusivities either increase or decrease. Interestingly, the corresponding static scattering intensities behave oppositely, Our current analysis indicates that these changes are inconsistent with aggregation in undersaturated solutions. However, the data are compatible with concentration-dependent changes of the interactions between protein and salt.

Rosenberger, Franz↗
DOE OSTI2025

Modeling and Simulation of Fuel Dispersal During the Loss-of-Coolant Accident

This document is the compilation of the milestone portion to a larger end of project NEUP report. The executive summary of the modeling portion is provided below: In the event of cladding rupture during a postulated LOCA in a pressurized water reactor, fuel particles, along with fission gases, can be expelled into the reactor core from the fractured fuel rod, a phenomenon referred to as fuel dispersal. The initial stage of fuel dispersal is strongly influenced by the high-pressure ejection of fuel fragments, the size and geometry of the ruptured cladding, and the depressurization history of the fuel rod during the postulated LOCA transient. Depending on the location of the burst orifice relative to the quench front, the dispersal event represents an intricate three-phase flow and heat transfer phenomenon, where high-temperature fuel particles carried by the fission gases interact with the coolant within the narrow subchannels of the fuel assemblies, inducing localized phase change. Given the unique multiphysics nature of this phenomena, the current study develops a dedicated computational framework to predict the mass distribution and cooling of dispersing fuel particles, facilitating post-accident assessment and management of the fuel assemblies. Considering the scale of nuclear reactor applications, a continuum three-fluid model is proposed for simulating the transport of solids within the reactor core. With high-temperature fuel fragments within the liquid media, nucleation sites inducing phase changes are dispersed within the flow domain. Coupled with the fact that the transient dispersal event occurs on different time scales than other three-phase flow applications, this study derives a time-averaged three-fluid flow model without losing generality. The assumptions regarding the continuum treatment of the solid phase and the modeling of fuel dispersal behavior are incorporated to simplify the governing equations and derive applicable closure relations. The computational validation of the model was conducted using adiabatic experimental results obtained from ongoing research at Oregon State University, focusing on characterizing fuel dispersal behavior during simulated LOCA conditions. Settlement characteristics of the solids, quantified by the probability distribution of equivalent particles, closely matched the probability density functions reported in experimental studies. The transport of fuel particles within a scaled 5 × 5 lattice of a pressurized-water reactor rod bundle geometry was modeled through a two-fluid Eulerian framework. The required boundary conditions were evaluated from the fuel performance code BISON in a postulated large-break LOCA scenario. The modeling framework considered solid fuel particles as granular matter, interacting with the gaseous dry steam phase and fission gases through the governing interfacial momentum exchange between the participating fluids. The simulation results provided the volume fraction of the solids obtained at the bottom surface of the enclosing tank geometry. Postulated LOCA leading to fuel dispersal phenomena involves the strong coupling between fuel thermomechanics, cladding deformation, thermal-hydraulics, and fuel particle transport. Incorporation of such a strong coupling in numerical simulation is performed by coupling the multiphysics solvers. In the case of fuel dispersal, a strong coupled simulation can be performed by coupling the BISON code for fuel performance, the TRACE code for system-level thermal hydraulics, and fuel particle transport in Multiphysics Object-Oriented Simulation Environment (MOOSE). For such intricate infrastructure, the MOOSE Framework eases the data transfer between codes. The recent version of MOOSE has incorporated the Navier-Stokes module for the fluid flow. An exploratory exercise was done to gain familiarity with finite volume capabilities in the MOOSE framework to incorporate the Spalart-Allmaras (SA) turbulence model. New finite-volume and auxiliary kernels were introduced to assemble the SA transport equation, compute turbulent viscosity, and evaluate wall distance and diagnostic turbulence terms, fully integrated with existing Navier-Stokes modules. A turbulent lid-driven cavity at a Reynolds number of approximately 10,000 is used for verification. MOOSE shows the robust solver convergence and produces the turbulent features. But it underpredicts the velocity profile and turbulent quantities, emphasizing the need to develop improved SA near-wall treatments (e.g., low-Re corrections or wall functions) as a key direction for future work.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
NASA NTRS2003

NASA Glenn/AADC-Rolls Royce Collaborated to Measure Erosion Resistance on Coated Polymer Matrix Composites

Polymer matrix composites (PMCs) are increasingly used in aerospace and automotive applications because of their light weight and high strength-to-weight ratio relative to metals. However, a major drawback of PMCs is poor abrasion resistance, which restricts their use, especially at high temperatures. Simply applying a hard coating on PMCs to improve abrasion and erosion resistance is not effective since coating durability is short lived (ref. 1). Generally, PMCs have higher coefficients of thermal expansion than metallic or ceramic coatings have, and coating adhesion suffers because of poor interfacial adhesion strength. One technique commonly used to improve coating adhesion or durability is the use of bond coats that are interleaved between a coating and a substrate with vastly different coefficients of thermal expansion. An example of this remedy is the use of bondcoats for ceramic thermal barrier coatings on metallic turbine components (ref. 2). Prior collaborative research between the NASA Glenn Research Center and the Allison Advanced Development Company (AADC) demonstrated that bond coats sandwiched between PMCs and high-quality plasma-sprayed, erosion-resistant coatings substantially improved the erosion resistance of PMCs (ref. 3). One unresolved problem in this earlier collaboration was that there was no easy, accurate way to measure the coating erosion wear scar. Coating wear was determined by both profilometry and optical microscopy. Both techniques are time consuming. Wear measurement by optical microscopy requires sample destruction and does not provide a comprehensive measure of the entire wear volume. An even more subtle, yet critical, problem is that these erosion coatings contain two or more materials with different densities. Therefore, simply measuring specimen mass loss before and after erosion will not provide an accurate gauge for coating and/or substrate volume loss. By using a noncontact technique called scanning optical interferometry, which was recently developed at Glenn, researchers can accurately determine the wear performance of erosion-coated PMCs while preserving the sample. An example of this interferometry technique is shown in the preceding figure for an erosion-coated inlet guide vane from a Rolls Royce AE3007 regional gas turbine jet engine. Erosion was conducted with coated and uncoated PMC vanes, with the abrasive material moving at a velocity of 229 m/s at impingement angles of 20 and 90 degrees. The coatings for PMCs remarkably reduced the erosion volume loss by a factor of approximately 10. Currently, several erosion coatings for PMCs are being compared and downselected for engine testing at Rolls Royce.

Miyoshi, Kazuhisa↗