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.
Surface property effects on effectiveness of mass transfer cooling in laminar boundary layers with hydrogen injected into free stream of nitrogen or carbon dioxide
Computational Fluid Dynamics (CFD) validation experiments at NASA Lewis Research Center are described. The material presented summarizes the research in three areas: Inlets, Ducts and Nozzles; Turbomachinery; and Chemically Reacting Flows. The specific validation activities are concerned with shock-boundary layer interactions, vortex generator effects, large low speed centrifugal compressor measurements, transonic fan shock structure, rotor/stator kinetic energy distributions, stator wake shedding characteristics, boundary layer transition, multiphase flow and reacting shear layers. These experiments are intended to provide CFD validation data for the internal flow fields within aerospace propulsion system components.
Computational Fluid Dynamics (CFD) validation experiments at NASA Lewis are described. The material presented summarized the research in 3 areas: Inlets, ducts and nozzles; Turbomachinery; and Chemically reacting flows. The specific validation activities are concerned with shock boundary layer interactions, vortex generator effects, large low speed centrifugal compressor measurements, transonic fan shock structure, rotor/stator kinetic energy distributions, stator wake shedding characteristics, boundary layer transition, multiphase flow and reacting shear layers. These experiments are intended to provide CFD validation data for the internal flow fields within aerospace propulsion system components.
The flutter characteristics of the first AGARD standard aeroelastic configuration for dynamic response, Wing 445.6, are studied using an unsteady Navier-Stokes algorithm in order to investigate a previously noted discrepancy between Euler flutter characteristics and the experimental data. The algorithm, which is a three-dimensional, implicit, upwind Euler/Navier-Stokes code (CFL3D Version 2.1), was previously modified for the time-marching, aeroelastic analysis of wings using the unsteady Euler equations. These modifications include the incorporation of a deforming mesh algorithm and the addition of the structural equations of motion for their simultaneous time integration with the governing flow equations. In this paper, the aeroelastic method is extended and evaluated for applications that use the Navier- Stokes aerodynamics. The paper presents a brief description of the aeroelastic method and presents unsteady calculations which verify this method for Navier-Stokes calculations. A linear stability analysis and a time-marching aeroelastic analysis are used to determine the flutter characteristics of the isolated 45 deg. swept-back wing. Effects of fluid viscosity, structural damping, and number of modes in the structural model are investigated. For the linear stability analysis, the unsteady generalized aerodynamic forces of the wing are computed for a range of reduced frequencies using the pulse transfer-function approach. The flutter characteristics of the wing are determined using these unsteady generalized aerodynamic forces in a traditional V-g analysis. This stability analysis is used to determine the flutter characteristics of the wing at free-stream Mach numbers of 0.96 and 1.141 using the generalized aerodynamic forces generated by solving the Euler equations and the Navier-Stokes equations. Time-marching aeroelastic calculations are performed at a free-stream Mach number of 1.141 using the Euler and Navier-Stokes equations to compare with the linear V-g flutter analysis method. The V-g analysis, which is used in conjunction with the time-marching analysis, indicates that the fluid viscosity has a significant effect on the supersonic flutter boundary for this wing while the structural damping and number of modes in the structural model have a lesser effect.
Deployment Objectives Flexibly co-locate with new industry users of nuclear energy Steam and electric cogeneration applications Direct process heat with temperature ranges from 700°C to 950°C Enabling Requirements Meet regulatory dose limits at the Exclusion Area Boundary (EAB) 25 rem Total Effective Dose Equivalent (TEDE) for duration of the release from 10 CFR 50.34 (10 CFR 52.79) at Exclusion Area Boundary (EAB) for design basis accidents EAB is estimated approximately 400 meters from the modular HTGR plant (to support co-location with industrial facilities) Meet safety goals for cumulative individual risk for normal and off-normal operation Design goal: meet EPA Protective Action Guides (PAGs) at EAB 1 rem TEDE for sheltering Design basis and beyond design basis events are considered Realistically evaluated at the EAB Emergency planning and protection
An incompressible three-dimensional laminar boundary-layer flow over a swept wing is used as a model to study both the wall-curvature and streamline-curvature effects on the stationary crossflow instability. The basic state is obtained by solving the full Navier-Stokes (N-S) equations numerically. The linear disturbance equations are cast on a fixed, body-intrinsic, curvilinear coordinate system. Those nonparallel terms which contribute mainly to the streamline-curvature effect are retained in the formulation of the disturbance equations and approximated by their local finite difference values. The resulting eigenvalue problem is solved by a Chebyshev collocation method. The present results indicate that the convex wall curvature has a stabilizing effect, whereas the streamline curvature has a destabilizing effect. A validation of these effects with an N-S solution for the linear disturbance flow is provided.
This study focuses on the effects of drizzle in a one-dimensional third-order turbulence closure model of the nocturnal stratus-topped marine boundary layer. When the simulated drizzle rate is relatively small (maximum approximately equal to 0.6 mm/day), steady-state solutions are obtained. The boundary layer stabilizes essentially because drizzle causes evaporative cooling of the subcloud layer. This stabilization considerably reduces the buoyancy flux and turbulence kinetic energy below the stratus cloud. Thus, drizzle tends to decouple the cloud from the subcloud layer in the model, as suggested by many observational studies. In addition, the evaporation of drizzle in the subcloud layer creates small scattered clouds, which are likely to represent cumulus clouds, below the solid stratus cloud in the model. The sensitivity experiments show that these scattered clouds help maintain a coupled boundary layer. When the drizzle rate is relatively large (maximum approximately equal to 0.9 mm/day), the response of the model becomes transient with bursts in turbulent fluxes. This phenomenon is related to the formation of the scattered cloud layer below the solid stratus cloud. It appears that the model is inadequate to represent the heat and moisture transport by strong updrafts covering a small fractional area in cumulus convection.
The Langley 0.3 Meter Transonic Cryogenic Tunnel has provision for boundary removal from the sidewalls to reduce sidewall interference effects on the test data. The tests carried out to determine the change in the empty test section sidewall boundary layer thickness at the model station with upstream boundary layer mass removal are described. The boundary layer measurements showed that the upstream removal region is effective in reducing the boundary layer thickness at the model station. The boundary layer displacement thickness reduced from about 1.2 percent to about .4 percent of the test section width. The boundary layer velocity profiles followed a power law variation in the outer region and showed good correlation when plotted in terms of boundary layer momentum thickness.
Avehicle-level noise assessment has been performed for theNASAD8 concept aircraft (ND8) in the NASA Advanced Air Transport Technology Project portfolio. The NASA research-level Aircraft NOise Prediction Program (ANOPP-Research) was used to predict the noise from each source component on the ND8 to build up a noise estimate for the full aircraft. The propulsion airframe aeroacoustic (PAA) effects of the ND8, namely boundary layer ingestion (BLI) with its influence on fan noise, and the noise shielding, reflection, and diffraction mechanisms of the unconventional airframe, were empirically modeled using experimental data. Noise reduction technologies appropriate to the 2025-2035 time frame were included in this study. Including all technologies and PAA effects, the ND8 is predicted to have a cumulative margin to the Stage 4 certification metric of only 7.4 EPNdB. Boundary layer ingestion is predicted to have a detrimental impact on cumulative noise levels on the order of 15 EPNdB. Fan noise is seen to be the primary noise source at all three certification points, even if the BLI noise impact could be entirely suppressed. The impact of engine noise shielding by the airframe is limited by a lack of aft shielding and the presence of horizontal tail reflections in the aft direction. The physical constraint on engine size by the pi-tail is seen as a potential barrier to engine noise reduction through the corresponding limitation on fan bypass ratio. Mildly reduced climb performance (compared to similar reference aircraft) does not provide any benefit through increased noise propagation distance. If the boundary layer ingestion noise penalty could be suppressed such that BLI would have no effect on noise, the cumulative margin to Stage 4 would increase to 22.4 EPNdB, still below the NASA Mid Term goal of 32-42 EPNdB.
Ideally, the reflection of a shock from the closed end of a shock tube provides, for laboratory study, a quantity of stationary gas at extremely high temperature. Because of the action of viscosity, however, the flow in the real case is not one-dimensional, and a boundary layer grows in the fluid following the initial shock wave. In this paper simplifying assumptions are made to allow an analysis of the interaction of the shock reflected from the closed end with the boundary layer of the initial shock afterflow. The analysis predicts that interactions of several different types will exist in different ranges of initial shock Mach number. It is shown that the cooling effect of the wall on the afterflow boundary layer accounts for the change in interaction type. An experiment is carried out which verifies the existence of the several interaction regions and shows that they are satisfactorily predicted by the theory. Along with these results, sufficient information is obtained from the experiments to make possible a model for the interaction in the most complicated case. This model is further verified by measurements made during the experiment. The case of interaction with a turbulent boundary layer is also considered. Identifying the type of interaction with the state of turbulence of the interacting boundary layer allows for an estimate of the state of turbulence of the boundary layer based on an experimental investigation of the type of interaction. A method is proposed whereby the effect of the boundary-layer interaction on the strength of the reflected shock may be calculated. The calculation indicates that the reflected shock is rapidly attenuated for a short distance after reflection, and this result compares favorably with available experimental results.
Abstract Engineering semiconductor devices requires an understanding of the effective mass of electrons and holes. Effective masses have historically been determined in metals at cryogenic temperatures estimated using measurements of the electronic specific heat. Instead, by combining measurements of the Seebeck and Hall effects, a density of states effective mass can be determined in doped semiconductors at room temperature and above. Here, a simple method to calculate the electron effective mass using the Seebeck coefficient and an estimate of the free electron or hole concentration, such as that determined from the Hall effect, is introduced here is the Seebeck effective mass, n H is the charge carrier concentration measured by the Hall effect ( n H = 1/ eR H , R H is Hall resistance) in 10 20 cm −3 , T is the absolute temperature in K, S is the Seebeck coefficient, and k B / e = 86.3 μV K −1 . This estimate of the effective mass can aid the understanding and engineering of the electronic structure as it is largely independent of scattering and the effects of microstructure (grain boundary resistance). It is particularly helpful in characterizing thermoelectric materials.
Temperature- and irradiation-assisted failure mechanisms in miscible phase boundaries are clarified via atomistic calculations. We first establish the temperature-dependent brittle-to-ductile transition in U–Zr miscible phase boundaries. Our results confirm that these boundaries are mostly brittle at low temperatures and ductile at elevated temperatures. We then investigate the changes induced by irradiation on the fracture mechanisms in such phase boundaries. The irradiation-induced defect accumulation follows a multi-stage process that starts with the accumulation of isolated small dislocation loops before transitioning to the saturation and growth of larger dislocation loops and end up with a reorganization into forest dislocations. The accumulation of loops is the primary feature to participate in the delineation between brittle and ductile interfacial fracture in irradiated phase boundaries. At low damage levels, radiation defect interactions with the crack tip are limited and U–Zr miscible boundaries fail through the emission of dislocations ahead of the crack tip followed by brittle cleavage in agreement with the classical Griffith’s criterion for crack stability. At higher damage levels, the failure mode transitions from brittle crack growth to ductile void growth. In this case, the microcrack tip is blunted by the high density of pre-existing, radiation-induced defects in the vicinity of the crack. This interaction leads to the development and growth of a cavity at the interface as opposed to interfacial cleavage.
Building-integrated photovoltaic (BIPV) systems are intrinsically designed to generate electricity and to provide at least one building-related function. When BIPV modules act as glazing products in windows, skylights or curtain walls, their ability to control the transmission of solar energy into the building must be characterised by a Solar Heat Gain Coefficient (SHGC) or g value (also known as Total Solar Energy Transmittance – TSET – or “solar factor”). For the comparison of BIPV glazing products consisting of one PV laminate and possibly further, conventional glazing layers separated by gas-filled cavities, the procedures documented in international standards for architectural glazing (e.g. ISO 9050 and EN 410) form a suitable starting point. Easily implemented modifications to these procedures are proposed to take both optical inhomogeneity (if relevant) and extraction of electricity from BIPV glazing units into account. Geometrically complex glazing and shading devices, and light-scattering glazing layers, are outside the scope of the proposed methodology; SHGC determination for obliquely incident solar radiation is also excluded. For these cases, the experimental calorimetric approach documented in [ISO 19467:2017; ISO 19467-2:2021] is recommended. The paper also presents results and conclusions from an implementation exercise and sensitivity study carried out by participants of the IEA-PVPS Task 15 on BIPV. The cell coverage ratio in the PV laminate, the thermal resistance offered by the glazing configuration, the choice of boundary conditions and the effect of extracting electricity were all identified as parameters which significantly affect the SHGC value determined for a given type of BIPV glazing. A practicable approach to accommodate the great variety of dimensions typical for BIPV glazing is also proposed. These findings should pave the way for modifying the existing component-based standards for architectural glazing to take the specific features of BIPV glazing into account.
Tritium processing facilities may release tritium oxide (HTO) to the atmosphere which poses potential health risks to exposed co-located workers and to offsite individuals. Most radiological consequence analyses determine HTO dose by applying Gaussian plume models to simulate the transport of HTO. Within these models, deposition velocity is used to assess the sum of all deposition processes acting on the plume. While this may account for vegetative and soil uptake or respiration processes, it may currently lack inclusion of the complex interactions within heterogeneous forested environments. In this complex morphology, dispersion patterns are significantly altered by changing flow regimes above and below the forest canopy and by the transfer of plume material across the canopy boundary. To determine the effects of a heterogeneous forest canopy on an airborne HTO plume, a Gaussian plume model coupled with an advection-diffusion plume model was applied to estimate transport in the free atmosphere above the forest and within the forest canopy and understory. During 2012, wind speed and wind direction measurements taken at 5 heights, ranging from 2-m to 28-m, on an instrumented meteorological tower located in a loblolly pine forest at the Department of Energy (DOE) Savannah River Site (SRS), near Aiken, SC. From these measurements, model predictions were made over a full spectrum of meteorological conditions. Deposition and resuspension velocities were calculated based on the model-predicted flux of plume material across the top of the forest canopy. Additionally, net deposition velocity of the plume material was calculated as the difference between the deposition and resuspension velocities. The 1st and 5th percentile net deposition velocities were estimated to be 0.7 cm s-1 and 1.2 cm s-1, respectively.
This work presents a method to incorporate the micro Hall-Petch equation into the crystal plasticity finite element (CPFE) framework accounting for the microstructural features to understand the coupling between grain size, texture, and loading direction in magnesium alloys. The effect of grain size and texture is accounted for by modifying the slip resistances of individual basal and prismatic systems based on the micro Hall-Petch equation. The modification based on the micro Hall-Petch equation endows every slip system at each microstructural point with a slip system-level grain size and maximum compatibility factor, which are in turn used to modify the slip resistance. While the slip-system level grain size is a measure of the grain size, the maximum compatibility factor encodes the effect of the grain boundary on the slip system resistance modification and is computed based on the Luster-Morris factor. The model is calibrated using experimental stress-strain curves of Mg-4Al samples with three different grain sizes from which the Hall-Petch coefficients are extracted and compared with Hall-Petch coefficients predicted using original parameters from previous work. The predictability of the model is then evaluated for a Mg-4Al sample with different texture and three grain sizes subjected to loading in different directions. The calibrated parameters are then used for some parametric studies to investigate the variation of Hall-Petch slope for different degrees of simulated spread in basal texture, variation of Hall-Petch slope with loading direction relative to basal poles for a microstructure with strong basal texture, and variation of yield strength with change in grain morphology. The proposed approach to incorporate the micro Hall-Petch equation into the CPFE framework provides a foundation to quantitatively model more complicated scenarios of coupling between grain size, texture and loading direction in the plasticity of Mg alloys.
Here, in this study, microstructural studies and micro-mechanical testing of an ex-plant material harvested from the decommissioned pressurized water reactor (PWR) is carried out. Irradiated 304 stainless steel (SS) components were harvested by the Electric Power Research Institute, U.S. Nuclear Regulatory Commission, and members of the José Cabrera Nuclear Power Station (Zorita) in Spain. The bi-crystalline micro-tensile specimens, atom probe tomography (APT) tips, and transmission electron microscopy (TEM) lamellae were fabricated from the baffle plate materials irradiated to 0.05, 15 and 50 dpa to elucidate microstructural, microchemical, and local mechanical properties changes as a function of dose. TEM and APT studies reveal radiation induced Ni, Cr, and Si rich precipitates in the 15 dpa sample that become smaller and partially redissolve into the matrix of the 50 dpa sample (decreasing size and number density). Dislocation loop and cavity number density and size were quantified as a function of dose as well as swelling. Intergranular Ni and Si enrichments with concomitant Fe and Cr depletion were observed in the 15 dpa sample and were more pronounced in the 50 dpa sample. Micro-tensile testing performed in the scanning electron microscope (SEM) at room temperature and 300 °C shows that the material's local yield strength and ultimate tensile strength decreases with increased dose and elevated test temperature and provides further insights via localized strain mapping. Current and previous mechanical data was compared with calculated values using the dispersed barrier hardening model with inputs of dislocation loops, precipitates, and cavities from microstructural characterization.
Microstructures of cast aluminum alloys used in automotive engine applications often consist of intermetallic particles that can impact the tensile elongation of these alloys. Here in this paper, we investigate the effect of intermetallic grain boundary particles on tensile elongation by fabricating a series of alloys with Cu content varying between 6.0 - 9.0 wt% in cast Al–Cu–Mn–Zr (ACMZ) type compositions. The tensile elongation of as-aged ACMZ alloys decreases monotonically with increase in Cu content. While the microstructure within the grains and yield stress of the alloy remains invariant with Cu content, the decrease in tensile elongation correlates well with increase in the size and volume fraction of grain boundary particles. Microstructural observations are combined with finite element simulations to explain the trend in tensile elongation with changing Cu content. Crack initiation is found to occur by brittle fracture of the grain boundary particles. Increase in particle size promotes crack initiation by reduction in size dependent particle fracture strength. Lower inter-particle spacing at higher particle volume fraction further facilitates crack initiation by increasing stress within the particles caused by the interaction between stress fields of neighboring particles. Increase in particle volume fraction also accelerates crack propagation through the formation of macro shear zones in the microstructure. The increase in Cu content of cast ACMZ alloys, therefore, decreases tensile elongation by promoting both crack initiation and crack propagation.
Nanocrystalline alloys can have exceptional strengths, yet due to limited microstructural stability it is difficult to fabricate bulk pieces through traditional processing routes that retain nanosized grains. In this study, centimeter-sized Cu-Zr alloy pellets were fabricated via a simple and improved powder metallurgy processing route. Different consolidation temperatures and times were employed to investigate the effect of amorphous grain boundary complexions on densification and the resulting mechanical properties. Bulk compression tests were carried out, with the samples that were hot pressed at 900 °C for 10 h exhibiting an excellent combination of average yield strength of 722 ± 45 MPa and average failure strain of 25.3 ± 2.4 %. Furthermore, we find that a powder processing route which enables amorphous complexion-assisted sintering leads to samples that (1) reach full density without requiring quenching treatments or other complex processing, (2) demonstrate appreciable plasticity, and (3) have strength that competes with commercially available high-strength Cu alloys.