Search NASA⌕ Search

SEARCH · Search NASA

Results for “Stokes parameters”

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.

At least 343 records · Page 19

Optimization of Plug Nozzles Using CFD Informed Design Space Querying

This paper presents an axisymmetric CFD approach based on Reynolds-averaged Navier-Stokes (RANS) calculations for the shape optimization of plug nozzle contours for both supersonic cruise and landing/takeoff (LTO) conditions. Three key external plug design parameters were defined as independent variables in the optimization scheme, while nozzle performance was evaluated using thrust coefficient (Cfg) and discharge coefficient (Cd) as the primary dependent variables. By varying these parameters, a wide range of axisymmetric conic section plug profiles can be achieved, simplifying the problem. The paper also discusses the supersonic plug nozzle optimization approach, design space seeding, meshing, and CFD simulations for the optimization process. The results demonstrate that the approach is effective in exploring the parameter space and reducing computational expenses, making it a viable option for optimizing nozzle plug contours. Using this optimization toolchain, the study yielded preliminary designs for LTO, cruise, and LTO-cruise balanced plug nozzles.

CFD↗

Large Eddy Simulation of Lean Mixed-Mode Combustion Assisted by Partial Fuel Stratification in a Spark-Ignition Engine

Partial fuel stratification (PFS) is a promising fuel injection strategy to improve the stability of lean combustion by applying a small amount of pilot injection right before spark timing. Mixed-mode combustion, which makes use of end-gas autoignition following conventional deflagration-based combustion, can be further utilized to speed up the overall combustion. In this study, PFS-assisted mixed-mode combustion in a lean-burn direct injection spark-ignition (DISI) engine is numerically investigated using multi-cycle large eddy simulation (LES). A previously developed hybrid G-equation/well-stirred reactor combustion model for the well-mixed operation is extended to the PFS-assisted operation. The experimental spray morphology is employed to derive spray model parameters for the pilot injection. The LES-based model is validated against experimental data and is further compared with the Reynolds-averaged Navier–Stokes (RANS)-based model. Overall, both RANS and LES predict the mean pressure and heat release rate traces well, while LES outperforms RANS in capturing the cycle-to-cycle variation (CCV) and the combustion phasing in the mass burned space. Liquid and vapor penetrations obtained from the simulations agree reasonably well with the experiment. Detailed flame structures predicted from the simulations reveal the transition from a sooting diffusion flame to a lean premixed flame, which is consistent with experimental findings. LES captures more wrinkled and stretched flames than RANS. Finally, the LES model is employed to investigate the impacts of fuel properties, including heat of vaporization (HoV) and laminar burning speed (SL). Combustion phasing is found more sensitive to SL than to HoV, with a larger fuel property sensitivity of the heat release rate from autoignition than that from deflagration. Moreover, the combustion phasing in the PFS-assisted operation is shown to be less sensitive to SL compared with the well-mixed operation.

fuel property sensitivity↗

On the penetration of a hot diapir through a strongly temperature-dependent viscosity medium

The ascent of a hot spherical body through a fluid with a strongly temperature-dependent viscosity has been studied using an axisymmetric finite element method. Numerical solutions range over Peclet numbers of 0.1 - 1000 from constant viscosity up to viscosity variations of 100,000. Both rigid and stress-free boundary conditions were applied at the surface of the sphere. The dependence of drag on viscosity variation was shown to have no dependence on the stress boundary condition except for a Stokes flow scaling factor. A Nusselt number parameterization based on the stress-free constant viscosity functional dependence on the Peclet number scaled by a parameter depending on the viscosity structure fits both stress-free and rigid boundary condition data above viscosity variations of 100. The temperature scale height was determined as a function of sphere radius. For the simple physical model studied in this paper pre-heating is required to reduce the ambient viscosity of the country rock to less than 10 to the 22nd sq cm/s in order for a 10 km diapir to penetrate a distance of several radii.

Daly, S. F.↗

Computational analysis of the flowfield of a two-dimensional ejector nozzle

A time-iterative full Navier-Stokes code, PARC, is used to analyze the flowfield of a two-dimensional ejector nozzle system. A parametric study was performed for two controlling parameters, duct to nozzle area ratio and nozzle pressure ratio. Results show that there is an optimum area ratio for the efficient pumping of secondary flow. At high area ratios, a freestream flow passes directly through the mixing duct without giving adequate pumping. At low area ratios, the jet boundary blocks the incoming flow. The nozzle pressure ratio variation shows that the pumping rate increases as the pressure ratio increases, provided there is no interaction between the shroud wall and the shock cell structure.

Choi, Y. H.↗

Embedded adaptation for time dependent real fluid/gas flows

An embedded adaptation algorithm is used to illustrate applications and requirements for Navier-Stokes and nonequilibrium unsteady two-dimensional flow fields. Both grid refinement and equation modification are included for unstructured distributions of quadrilateral cells. Multiple refinement parameters for distinct physical phenomena serve to monitor the evolving nonuniform field and to recognize flow features in isolated and overlapping subdomains. Examples from subsonic, hypersonic, internal, and external fields show substantial advantages for the adaptive approach.

Baron, Judson R.↗

Computational analysis of the flowfield of a two-dimensional ejector nozzle

A time-iterative full Navier-Stokes code, PARC, is used to analyze the flowfield of a two-dimensional ejector nozzle system. A parametric study was performed for two controlling parameters, duct to nozzle area ratio and nozzle pressure ratio. Results show that there is an optimum area ratio for the efficient pumping of secondary flow. At high area ratios, a freestream flow passes directly through the mixing duct without giving adequate pumping. At low area ratios, the jet boundary blocks the incoming flow. The nozzle pressure ratio variation shows that the pumping rate increases as the pressure ratio increases, provided there is no interaction between the shroud wall and the shock cell structure.

Choi, Y. H.↗

Design and Modeling of an 80 bar Oxy-Combustor for Direct Fired Supercritical CO2 Applications

As supercritical CO2 power cycles for fossil energy power generation continue to generate interest, there are significant issues and unanswered questions regarding injector design, flame stabilization, wall heat transfer, CO emissions, combustion dynamics and other combustion phenomenon. For natural gas, direct-fired cycles with carbon capture it is believed that Computational Fluid Dynamics (CFD) modeling will play an essential role in the combustor design process. To accurately model turbulent reacting flows at these unique conditions, experimental data is needed to validate CFD codes and sub-models and is currently lacking at conditions relevant for these cycles. This paper presents the conceptual design and CFD simulations of an experimental 80 bar oxy-combustion facility and test article currently under construction at NETL. The facility is targeted towards the testing of a single injector, direct-fired sCO2 combustor at the 100 kW thermal output level. While these conditions do not reflect the actual Allam cycle operating conditions (300 bar) they are viewed as a stepping stone in the model validation process at supercritical conditions. Reynolds Averaged Navier-Stokes as well as Large Eddy Simulation are used to model the turbulent combustion process and aid in the design of the combustor and injector. Process parameters including oxidizer preheat temperature and combustor flowrate are investigated.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Assessing the Heat Transfer Modeling Capabilities of CFD Software for Involute-Shaped Plate Research Reactors

The ongoing efforts to convert High-Performance Research Reactors (HPRRs) using Highly Enriched Uranium (HEU) to Low-Enriched Uranium (LEU) fuel require reliable thermal–hydraulic assessments of modified core designs. The involute-shaped fuel plates used in several major HPRRs present unique modeling challenges due to their compact core geometries and high heat flux conditions. This study evaluates the capability of three commercial CFD tools, STAR-CCM+, COMSOL, and ANSYS CFX, to predict cladding-to-coolant heat transfer using Reynolds-Averaged Navier–Stokes (RANS) methods within the thermal–hydraulic regimes of involute-shaped plate reactors. Broad sensitivity analysis was conducted across a range of reactor-relevant parameters using two turbulence models (k−ϵ and k−ω SST) and different near-wall treatment strategies. The results were benchmarked against the Sieder–Tate correlation and experimental data from historic studies. The codes produced consistent results, showing good agreement with the empirical correlation of Sieder–Tate and the experimental measurements. The findings support the use of these commercial CFD codes as effective tools for assessing the thermal–hydraulic performance of involute-shaped plate HPRRs and guide future LEU core development.

CFD↗

Similitude requirements for hypersonic, rarefied, nonequilibrium flow

Similitude requirements for hypersonic, rarefied flow with nonequilibrium chemistry and vibration are presented. The full Navier-Stokes equations with catalytic or noncatalytic walls and with or without slip conditions are nondimensionalized. The heat transfer coefficient is written in terms of fourteen dimensionless parameters and reduced to four by making the binary scaling assumption. Duplication of blunt and sharp nose heat transfer requires the use of air over a geometrically similar model with the same free stream velocity, wall temperature and product of free stream density and characteristic length. Estimates of this heat transfer coefficient are also presented.

Hendricks, W. L.↗

Supersonic Wind Tunnel Tests of a Half-axisymmetric 12 Deg-spike Inlet to a Rocket-based Combined-cycle Propulsion System

Results of an isolated inlet test for NASA's GTX air-breathing launch vehicle concept are presented. The GTX is a Vertical Take-off/ Horizontal Landing reusable single-stage-to-orbit system powered by a rocket-based combined-cycle propulsion system. Tests were conducted in the NASA Glenn 1- by 1-Foot Supersonic Wind Tunnel during two entries in October 1998 and February 1999. Tests were run from Mach 2.8 to 6. Integrated performance parameters and static pressure distributions are reported. The maximum contraction ratios achieved in the tests were lower than predicted by axisymmetric Reynolds-averaged Navier-Stokes computational fluid dynamics (CFD). At Mach 6, the maximum contraction ratio was roughly one-half of the CFD value of 16. The addition of either boundary-layer trip strips or vortex generators had a negligible effect on the maximum contraction ratio. A shock boundary-layer interaction was also evident on the end-walls that terminate the annular flowpath cross section. Cut-back end-walls, designed to reduce the boundary-layer growth upstream of the shock and minimize the interaction, also had negligible effect on the maximum contraction ratio. Both the excessive turning of low-momentum comer flows and local over-contraction due to asymmetric end-walls were identified as possible reasons for the discrepancy between the CFD predictions and the experiment. It is recommended that the centerbody spike and throat angles be reduced in order to lessen the induced pressure rise. The addition of a step on the cowl surface, and planar end-walls more closely approximating a plane of symmetry are also recommended. Provisions for end-wall boundary-layer bleed should be incorporated.

DeBonis, J. R.↗

Learning Optimal Aerodynamic Designs

This project created a framework for efficient, accurate, and scalable deep neural network representations of design optimization problem solutions. The inputs to these DNN representations are the vector of design requirement parameters, the outputs are the optimal design variables, and the goal is to learn the map from inputs to outputs (i.e., inverse design). The team addressed the problem of the optimal shape design of aerodynamic lifting surfaces—in particular aircraft wings—using a Reynolds-Average Navier Stokes model to govern the CFD-based aerodynamic shape optimization. The inverse design map for such problems is very complex and high-dimensional, involving inputs and outputs on the order of 1000s. To approximate this inverse design map, the team developed algorithms to construct parsimonious DNN architectures, which automatically identify low-dimensional manifolds in which design requirements affect optimal shape parameters, and trained these architectures with multifidelity optimization methods. The resulting methodology accurately and automatically designs optimal aerodynamic lifting surfaces with very high accuracy (99%) at interactive speeds, of the order of milliseconds, resulting in factors of one million or more speedup relative to CFD-based design optimization.

97 MATHEMATICS AND COMPUTING↗

An AOTV aeroheating and thermal protection study

The aerothermodynamics and thermal protection of aerobraking orbital transfer vehicles are studied with a view toward design methodology and technology requirements. A particular class of geometries is investigated: the ellipsoidally blunted raked-off elliptic cone. Discussed are the influence of various geometrical parameters on the trajectory and the heat flux as well as heating prediction methods ranging from an engineering correlation formula to computational solutions of the Navier-Stokes equations. Finite-rate catalytic thermal protection materials and gas radiation are considered. Estimates of the specific mass of possible thermal protection systems indicate that support structure is a large contributor to the aerobrake mass as compared with the thermal protection insulation itself. Sizing curves are included which enable the designer to estimate the size aerobrake required from heating considerations for different mass vehicles having aerobrakes of several specific masses.

Scott, C. D.↗

Integrated Oxygen Flow Meter/Heat Exchanger for Portable Life Support Systems

Space suits for future exploration missions will have multi-mission goals with new and challenging requirements for the portable life support system (PLSS). In particular, the space suit ventilation loop requires cooling and flow measurement components that must meet specifications that go well beyond the capabilities of the components used for the existing Extravehicular Mobility Unit. The flow meter must have high measurement accuracy over a wide flow range, compatibility with pure oxygen, low pressure losses, and very compact size. The heat exchanger that cools the ventilation loop must be built from materials that are compatible with the liquid cooling loop, and it must provide efficient gas cooling in a small package across conditions ranging from normal suit pressure to elevated pressure. This paper describes the development of a novel device that combines the flow measurement and cooling functions in a single, compact flow meter/heat exchanger (FMHX). We have developed design methods that enable us to assess trade-offs, optimize performance, and specify the design of an FMHX that meets the requirements and constraints for operation in future PLSSs. We used empirical design correlations combined with computational fluid dynamics analysis to design the FMHX design. Data from tests of a proof-of-concept FMHX validate the design methods and show that the device meets all design requirements. We used the results from these tests to refine the design parameters and predict performance of an optimized, prototype FMHX.

Izenson, Michael G.↗

Direct Measurement of the Selective Microwave-Induced Heating of Agglomerates of Dipolar Molecules: The Origin of and Parameters Controlling a Microwave Specific Superheating Effect

Agglomerates of polar molecules in nonpolar solvents are selectively heated by microwave radiation. The magnitude of the selective heating was directly measured by using the temperature dependence of the intensities of the Stokes and anti-Stokes bands in the Raman spectra of p-nitroanisole (pNA) and mesitylene. Under dynamic heating conditions, a large apparent temperature difference (ΔT) of over 100 °C was observed between the polar pNA solute and the nonpolar mesitylene solvent. This represents the first direct measurement of the selective microwave heating process. The magnitude of the selective microwave heating was affected by the properties of the agglomerated pNA. As the concentration of the pNA increases, the magnitude of the selective heating of the pNA was observed to decrease. This is explained by the tendency of the pNA dipoles to orient in an antiparallel fashion in the aggregates as measured by the Kirkwood g value, which decreased with increasing concentration. This effect reduces the net dipole moment of the agglomerates, which decreases the microwave absorption. After the radiation was terminated, the effective temperature of the dipolar molecules returned slowly to that of the medium. The slow heat transfer was modeled successfully by treating the solutions as a biphasic solvent/solute system. Based on modeling and the fact that the agglomerate can be heated above the boiling temperature of the solvent, an insulating layer of solvent vapor is suggested to form around the heated agglomerate, slowing convective heat transfer out of the agglomerate. Furthermore, this is an effect unique to microwave heating.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Global Sensitivity Analysis of Simulated Remote Sensing Polarimetric Observations Over Snow

This study presents a detailed theoretical assessment of the information content of passive polarimetric observations over snow scenes, using a global sensitivity analysis (GSA) method. Conventional sensitivity studies focus on varying a single parameter while keeping all other parameters fixed. In contrast, the GSA correctly addresses the covariance of state parameters across their entire parameter space, hence favoring a more correct interpretation of inversion algorithms and the optimal design of their state vectors. The forward simulations exploit a vector radiative transfer model to obtain the Stokes vector emerging at the top of the atmosphere for different solar zenith angles, when the bottom boundary consists of a vertically resolved snowpack of non-spherical grains. The presence of light-absorbing particulates (LAPs), either embedded in the snow or aloft in the atmosphere above in the form of aerosols, is also considered. The results are presented for a set of wavelengths spanning the visible (VIS), near-infrared (NIR), and shortwave infrared (SWIR) region of the spectrum. The GSA correctly captures the expected, high sensitivity of the reflectance to LAPs in the VIS–NIR and to grain size at different depths in the snowpack in the NIR–SWIR. With adequate viewing geometries, mono-angle measurements of total reflectance in the VIS–SWIR (akin to those of the Moderate Resolution Imaging Spectroradiometer, MODIS) resolve grain size in the top layer of the snowpack sufficiently well. The addition of multi-angle polarimetric observations in the VIS–NIR provides information on grain shape and microscale roughness. The simultaneous sensitivity in the VIS–NIR to both aerosols and snow-embedded impurities can be disentangled by extending the spectral range to the SWIR, which contains information on aerosol optical depth while remaining essentially unaffected when the same particulates are mixed with the snow. Multi-angle polarimetric observations can therefore (i) effectively partition LAPs between the atmosphere and the surface, which represents a notorious challenge for snow remote sensing based on measurements of total reflectance only and (ii) lead to better estimates of grain shape and roughness and, in turn, the asymmetry parameter, which is critical for the determination of albedo. The retrieval uncertainties are minimized when the degree of linear polarization is used in place of the polarized reflectance. The Sobol indices, which are the main metric for the GSA, were used to select the state parameters in retrievals performed on data simulated for multiple instrument configurations. Improvements in retrieval quality with the addition of measurements of polarization, multi-angle views, and different spectral channels reflect the information content, identified by the Sobol indices, relative to each configuration. The results encourage the development of new remote sensing algorithms that fully leverage multi-angle and polarimetric capabilities of modern remote sensors. They can also aid flight planning activities, since the optimal exploitation of the information content of multi-angle measurements depends on the viewing geometry. The better characterization of surface and atmospheric parameters in snow-covered regions advances research opportunities for scientists of the cryosphere and ultimately benefits albedo estimates in climate models.

remote sensing↗

Comparing Global Spray Combustion Characteristics and Local Shot-to-Shot Variations in a Reacting n -Heptane Spray

In support of efforts to develop improved models of turbulent spray behavior and combustion in diesel engines, experimental data and analysis must be obtained for guidance and validation. For Reynolds-averaged Navier–Stokes (RANS)-based Computational fluid dynamics (CFD) modeling approaches, representative ensemble average experimental results are important. For high-fidelity models such as large eddy simulations (LES)-based CFD, precise individual experimental results are desirable. However, making comparisons between a given experiment and LES is a challenge since local parameters cannot be directly compared. In this work, an optically accessible constant pressure flow rig (CPFR) is utilized to acquire diesel-like fuel injection and reaction behavior simultaneously with three optical diagnostic techniques: rainbow Schlieren deflectometry (RSD), OH* chemiluminescence (OH*), and two-color pyrometry (2CP). The CPFR allows a large number of repeated injection experiments to be performed for statistical analysis and convergence using ensemble-averaging techniques, while maintaining highly repeatable test conditions. Even for stable test conditions, variations in local turbulent fuel–air mixing introduce variability, which manifests as significant differences in OH* and 2CP results. Experimental measurements of characteristic parameters including liquid and vapor jet penetration, liftoff length, soot temperature and concentration, and turbulent flame speed, along with the shot-to-shot variability of each dataset, are presented and discussed. Finally, a statistical method is utilized to analyze the extent of this variability, and to identify superlative injections within the dataset for discussion and analysis of shot-to-shot variations.

42 ENGINEERING↗

Computational and Experimental Study of Film-Cooling Effectiveness With and Without Downstream Vortex Generators

Downstream vortex generators that involve a pair of rectangular plates arranged in an open V-shape placed just downstream of each film-cooling hole were shown to create flow and vortical structures that entrain lifted film-cooling flow back to the surface and increase its lateral spreading on the surface (GT2020–14317). In this study, computations and measurements were performed to examine the flow mechanisms this vortex generator induces to improve film-cooling effectiveness of a flat plate with the cooling flow emanating from one row of inclined holes. Parameters studied include blowing ratio (BR = 0.75 and 1.0) and temperature ratio (TR = 1.07 and 1.9). Here, the computational study is based on steady Reynolds-averaged Navier–Stokes (RANS) closed by the shear-stress transport (SST) turbulence model with and without conjugate analysis. The experimental study was conducted by using a conjugate heat transfer test rig with a plenum, where cooling flow is introduced. Measurements made include velocity and temperature profiles upstream and downstream of the film-cooling holes as well as the temperature at several locations on the hot and cold sides of the film-cooled flat plate. The computational study was validated by comparing computed results with those from measurements at BR = 0.75 and 1.0 and TR = 1.9. Computational and experimental results are presented to show the effects of BR and TR on the flow structures and how those structures improve the effectiveness of film cooling with and without the downstream vortex generators and with and without conjugate heat transfer.

42 ENGINEERING↗