Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numeric profiling”

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 1,729 records · Page 96
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Synthetic Absorption Lines from Simulations of Multiphase Gas in Galactic Winds

Supernova-driven galactic winds are multiphase streams of gas that are often observed flowing at a range of velocities out of star-forming regions in galaxies. In this study, we use high-resolution 3D simulations of multiphase galactic winds modeled with the hydrodynamics code Cholla to investigate the connection between numerical studies and observations. Using a simulated interaction between a hot $T$ ~ 10 4 K supernova-driven wind and a cool $T$ ~ 10 7 K cloud of interstellar material, we create mock observables, including the optical depth ($\tau$) and covering fraction ($C$ $f$ ) of six commonly observed ions (Si II , C II , Si IV , C IV , N V , and O VI ) as a function of gas velocity. We compare our mock observables to surveys of galactic winds in the literature, finding good agreement with velocities and profiles of the low ions. We then compute "empirical" values for the optical depth and covering fraction following observational techniques and compare them to the values calculated directly from the simulation data. We find that the empirically computed values tend to underestimate the "true" value of $\tau$ for ions with high optical depth and overestimate the "true" value of $\tau$ for ions with low optical depth relative to the simulated data. In conclusion, the empirically estimated covering fractions match our direct calculations very well for the low-ionization ions; for the high-ionization ions, the empirical covering fractions underestimate the directly computed values by up to ~40%.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1973

Preliminary study of a wall stabilized constructed arc

An iterative, implicit, finite-difference numerical technique is described which is suitable for obtaining solutions to the governing equations for a gas flowing in an axially symmetric constricted-arc heater. The method is shown to provide adequate solutions for three cases of simple pipe flows found in the literature, and for flow in a constricted-arc heater by direct comparison with experimental data. The comparison with arc-heater data includes static pressure, arc voltage, and wall heat flux, all as a function of axial location, and a radial temperature profile at one axial station. The arc-heater data were taken with air as the test gas at a heater inlet pressure approximately 0.40 atm and two flow rates of 2.2 and 4.8 g/sec. The arc currents investigated were between 377 and 584 amperes.

Graves, R. A.↗
NASA NTRS1991

Infrared cooling rate calculations in operational general circulation models - Comparisons with benchmark computations

The performance of several parameterized models is described with respect to numerical prediction and climate research at GFDL, NCAR, and GISS. The radiation codes of the models were compared to benchmark calculations and other codes for the intercomparison of radiation codes in climate models (ICRCCM). Cooling rates and fluxes calculated from the models are examined in terms of their application to established general circulation models (GCMs) from the three research institutions. The newest radiation parameterization techniques show the most significant agreement with the benchmark line-by-line (LBL) results. The LBL cooling rates correspond to cooling rate profiles from the models, but the parameterization of the water vapor continuum demonstrates uncertain results. These uncertainties affect the understanding of some lower tropospheric cooling, and therefore more accurate parameterization of the water vapor continuum, as well as the weaker absorption bands of CO2 and O3 is recommended.

Kiehl, J. T.↗
NASA NTRS1992

Aeroloads and secondary flows in a transonic mixed flow turbine stage

A numerical simulation of a transonic mixed flow turbine stage has been carried out using an average passage Navier-Stokes analysis. The mixed flow turbine stage considered here consists of a transonic nozzle vane and a highly loaded rotor. The simulation was run at the design pressure ratio and is assessed by comparing results with those of an established throughflow design system. The three-dimensional aerodynamic loads are studied as well as the development and migration of secondary flows and their contribution to the total pressure loss. The numerical results indicate that strong passage vortices develop in the nozzle vane, mix out quickly, and have little impact on the rotor flow. The rotor is highly loaded near the leading edge. Within the rotor passage, strong spanwise flows and other secondary flows exist along with the tip leakage vortex. The rotor exit loss distribution is similar in character to that found in radial inflow turbines. The secondary flows and non-uniform work extraction also tend to significantly redistribute a non-uniform inlet total temperature profile by the exit of the stage.

Kirtley, K. R.↗
NASA NTRS2008

CALIPSO Satellite Lidar Identification Of Elevated Dust Over Australia Compared With Air Quality Model PM60 Forecasts

Global measurements of the vertical distribution of clouds and aerosols have been recorded by the lidar on board the CALIPSO (Cloud Aerosol Lidar Infrared Pathfinder Satellite Observations) satellite since June 2006. Such extensive, height-resolved measurements provide a rare and valuable opportunity for developing, testing and validating various atmospheric models, including global climate, numerical weather prediction, chemical transport and air quality models. Here we report on the initial results of an investigation into the performance of the Australian Air Quality Forecast System (AAQFS) model in forecasting the distribution of elevated dust over the Australian region. The model forecasts of PM60 dust distribution are compared with the CALIPSO lidar Vertical Feature Mask (VFM) data product. The VFM classifies contiguous atmospheric regions of enhanced backscatter as either cloud or aerosols. Aerosols are further classified into six subtypes. By comparing forecast PM60 concentration profiles to the spatial distribution of dust reported in the CALIPSO VFM, we can assess the model s ability to predict the occurrence and the vertical and horizontal extents of dust events within the study area.

Young, Stuart A.↗
NASA NTRSDate not supplied

Extinction of the Stagnation Point Diffusion Flame: Effect of Conductive Heat Loss into Solid Interior

The planned microgravity experiment Growth and Extinction Limit (GEL) is designed to study the material flammability and extinction of thermally thick solid with an emphasis to quantify the effect of conductive heat loss into the solid interior by measuring the in-depth temperature profile. To obtain a more detailed understanding of the extinction phenomena, a previously developed one-dimensional axisymmetric stagnation-point diffusion flame model is expanded to account for the excessive heat loss into the solid beyond the steady burning state. This is important because for thick solids, this is important because the burning period can have a long and slow solid transient with a quasi-steady gas -phase. A non-dimensional excessive conductive heat loss parameter Ψ is defined in the paper and extensive numerical computations have been performed to determine the extinction boundary and the characteristics of near-limit flames. The classical U-shaped flammability boundary in the oxygen vs. stretch rate plot now becomes a series of U-shaped boundaries as a function of Ψ. For a given oxygen ambient, extinction boundary is a function of Ψ. The boundary consists of a blowoff branch and a radiative quenching branch. This is similar to the steady burning flammability boundary (Ψ =0) in the oxygen vs stretch rate plot. Extensive computed data including flame and pyrolysis temperatures, burning rate, flame standoff distance, and reaction rate are presented along the extinction boundaries. In particular, we note that the non-dimensional gas reaction rate is the ratio of diffusion reaction time to the reaction diffusion time, so it has the classical Damhkoler number character. By examining the maximum non-dimensional reaction rates, we observe that along the extinction boundaries, they are approximate constant across the entire range of stretch rates (from quenching to blow-off). These critical values are of the order unity but may vary in different oxygen ambient. The implication will be discussed in the presentation together with the near-limit flame structures.

Chengyao Li↗
DOE OSTI2021

Anisotropic fluid dynamical simulations of heavy-ion collisions

We present VAH, a (3+1)–dimensional simulation that evolves the far-from-equilibrium quark-gluon plasma produced in ultrarelativistic heavy-ion collisions with anisotropic fluid dynamics. We solve the hydrodynamic equations on an Eulerian grid using the Kurganov–Tadmor algorithm in combination with a new adaptive Runge–Kutta method. Our numerical scheme allows us to start the simulation soon after the nuclear collision, largely avoiding the need to integrate it with a separate pre-equilibrium dynamics module. We test the code's performance by simulating on the Eulerian grid conformal and non-conformal Bjorken flow as well as conformal Gubser flow, whose (0+1)–dimensional solutions are precisely known. Finally, we compare non-conformal anisotropic hydrodynamics to second-order viscous hydrodynamics in central Pb+Pb collisions and find that the former's longitudinal flow profile responds more consistently to the fluid's gradients along the spacetime rapidity direction.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Deprojection of X-ray data in galaxy clusters: confronting simulations with observations

ABSTRACT Numerical simulations with varying realism indicate an emergent principle−multiphase condensation and large cavity power occur when the ratio of the cooling time to the free-fall time (tcool/tff) falls below a threshold value close to 10. Observations indeed show cool-core signatures when this ratio falls below 20–30, but the prevalence of cores with tcool/tff ratio below 10 is rare as compared to simulations. In X-ray observations, we obtain projected spectra from which we have to infer radial gas density and temperature profiles. Using idealized models of X-ray cavities and multiphase gas in the core and 3D hydro jet-ICM simulations, we quantify the biases introduced by deprojection based on the assumption of spherical symmetry in determining tcool/tff. We show that while the used methods are able to recover the tcool/tff ratio for relaxed clusters, they have an uncertainty of a factor of 2−3 in systems containing large cavities (≳ 20 kpc). We also show that the mass estimates from these methods, in the absence of X-ray spectra close to the virial radius, suffer from a degeneracy between the virial mass (M200) and the concentration parameter (c) in the form of M200c2 ≈ constant. Additionally, the lack of soft-X-ray (≲ 0.5 keV) coverage and poor spatial resolution makes us overestimate min(tcool/tff) by a factor of few in clusters with min(tcool/tff) ≲ 5. This bias can largely explain the lack of cool-core clusters with min(tcool/tff) ≲ 5.

Sarkar, Kartick C. (ORCID:0000000277678472)↗
DOE OSTI2021

Computational Fluid Dynamics Simulations to Support Efficiency Improvements in Aluminum Smelting Process

Smelting is broadly described as the extraction of a metal from its ore. In the United States, aluminum is commonly produced by smelting alumina in bauxite using the Hall-Héroult process. Optimization of equipment and processes in conventional smelting is crucial to enhancing process efficiency and productivity, is necessary for improving the techno-economic feasibility, which directly manifests as the growth of the American economy. To achieve optima, insightful data on the multiphysics phenomena that are inherent to the process must be obtained through physical investigation or high-fidelity numerical simulations. The resolution of relevant scales in time and space for smelting operations requires intensive, high-performance computing (HPC) simulations. Hostile operating conditions limit physical data acquisition to specific techniques; therefore, these data do not describe the multiscale interaction of simultaneous effects. Fortunately, in recent decades, significant advancements in computing hardware and computational methods have made the numerical resolution of such a complex process possible. In this study, a high-fidelity simulation of aluminum smelting was performed using an open-source tool, OpenFOAM, which analyzed many parameters characteristic to underlying phenomena. A multiphysics model based on the Eulerian-Eulerian multifluid approach was adopted. This model can resolve critical issues in the electrolytic smelting of aluminum, such as bubbling of carbon dioxide from the anode(s), magnetohydrodynamics from electromagnetic effects, ionic dissolution of the alumina in the electrolyte, and the evolution of thermal profiles. This study provides valuable connectivity for characteristic data that can direct the future designs of efficient smelters. A basic framework to model and simulate the smelting process using OpenFOAM is presented for user modification in keeping with process development. Of relevance to the flow field, a detailed investigation of vortices produced by bubble motion and electromagnetics is discussed, along with their impact on the evolution of thermal profiles. The predictions show small-scale vortices in the clearance between the anode and cathode caused by magnetic forces. Predictions also indicate relatively large-scale vortices in the inter-anode space resulting from carbon dioxide rising through the electrolytic flow field. The formation of vortices at the edges of anodes was shown to direct alumina charged by the feeder to the bottom of the anodes, thus preventing the entrapment of gas bubbles in the periphery of the bottom of the anode. Symmetry was observed in the location of cold spots in the electrolytic mixture in the vicinity of the feeder. Cold spots were also observed in the clearance between the anode and cathode due to the flow’s transmission of unconverted alumina to this region.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Can reanalysis products outperform mesoscale numerical weather prediction models in modeling the wind resource in simple terrain?

Mesoscale numerical weather prediction (NWP) models are generally considered more accurate than reanalysis products in characterizing the wind resource at heights of interest for wind energy, given their finer spatial resolution and more comprehensive physics. However, advancements in the latest ERA-5 reanalysis product motivate an assessment on whether ERA-5 can model wind speeds as well as a state-of-the-art NWP model – the Weather Research and Forecasting (WRF) Model. We consider this research question for both simple terrain and offshore applications. Specifically, we compare wind profiles from ERA-5 and the preliminary WRF runs of the Wind Integration National Dataset (WIND) Toolkit Long-term Ensemble Dataset (WTK-LED) to those observed by lidars at a site in Oklahoma, United States, and in a United States Atlantic offshore wind energy area. We find that ERA-5 shows a significant negative bias (~-1ms-1) at both locations, with a larger bias at the land-based site. WTK-LED-predicted wind speed profiles show a limited negative bias (~-0.5ms-1) offshore and a slight positive bias (~+0.5ms-1) at the land-based site. On the other hand, we find that ERA-5 outperforms WTK-LED in terms of the centered root-mean-square error (cRMSE) and correlation coefficient, for both the land-based and offshore cases, in all atmospheric stability conditions. We find that WTK-LED's higher cRMSE is caused by its tendency to overpredict the amplitude of the wind speed diurnal cycle. At the land-based site, this is partially caused by wind plant wake effects not being accurately captured by WTK-LED.

17 WIND ENERGY↗
NASA NTRS1974

A global model of the earth's ionosphere for use in space applications

A general expression is derived for the F layer electron density profile as a function of latitude and longitude for that part of the earth which is in direct sunlight including dawn and dusk. Furthermore, the derived model is extended to encompass the night-time ionosphere. The expressions allow determination by standard means of the range correction for arbitrary ray path directions. It is also shown that the naive application of the Chapman ionospheric model entails range correction errors which for low elevation angles (less than 20 deg) and large solar zenith angles (40 deg) cannot be tolerated. Numerical calculations are displayed showing the dependence of the range correction on the pertinent parameters.

Von Roos, O. H.↗
NASA NTRS1988

Experimental and numerical investigation of the effect of distributed suction on oblique shock wave/turbulent boundary layer interaction

A combined experimental and numerical study of the interaction of an incident oblique shock wave with a turbulent boundary layer on a rough plate and on a porous plate with suction is presented. The experimental phase involved the acquisition of mean data upstream of, within, and downstream of the interaction region at Mach numbers 2.5 and 3.0. Data were taken at unit Reynolds numbers of 1.66 E7 and 1.85 E7 m respectively, and for flow deflection angles of 0, 4, 6 and 8 degs. Measured data include wall static pressure, pitot pressure profiles, and local bleed distributions on the porous plate. On the rough plate, with no suction, the boundary layer profiles were modified near the wall, but not separated for the 4 deg flow deflection angle. For the higher deflection angles of 6 and 8 degs, the boundary layer was separated. Suction increases the strength of the incident shock required to separate the turbulent boundary layer; for all shock strengths tested, separation is completely eliminated. The pitot pressure profiles are affected throughout the whole boundary layer; they are fuller than the ones obtained on the rough plate. It is also found that the combination of suction and roughness introduces spatial perturbations.

Benhachmi, Driss↗
NASA NTRS1990

Mode-medium instability and its correction with a Gaussian reflectivity mirror

A high power CO2 laser beam is known to deteriorate after a few microseconds due to a mode-medium instability (MMI) which results from an intensity dependent heating rate related to the vibrational-to-translational decay of the upper and lower CO2 lasing levels. An iterative numerical technique is developed to model the time evolution of the beam as it is affected by the MMI. The technique is used to study the MMI in an unstable CO2 resonator with a hard-edge output mirror for different parameters like the Fresnel number and the gas density. The results show that the mode of the hard edge unstable resonator deteriorates because of the diffraction ripples in the mode. A Gaussian-reflectivity mirror was used to correct the MMI. This mirror produces a smoother intensity profile which significantly reduces the effects of the MMI. Quantitative results on peak density variation and beam quality are presented.

Webster, K. L.↗
NASA NTRS1992

Mode-medium instability and its correction with a Gaussian-reflectivity mirror

A high-power CO2 laser beam is known to deteriorate after a few microseconds due to a mode-medium instability (MMI) which results from an intensity-dependent heating rate related to the vibrational-to-translational decay of the upper and lower CO2 lasing levels. An iterative numerical technique is developed to model the time evolution of the beam as it is affected by the MMI. The technique is used to study the MMI in an unstable CO2 resonator with a hard-edge output mirror for different parameters like the Fresnel number and the gas density. The results show that the mode of the hard edge unstable resonator deteriorates because of the diffraction ripples in the mode. A Gaussian-reflectivity mirror was used to correct the MMI. This mirror produces a smoother intensity profile which significantly reduces the effects of the MMI. Quantitative results on peak density variation and beam quality are presented.

Webster, K. L.↗
DOE OSTI2025

Experimental Investigation of Vapor Formation in Liquid CO2 Flow Through a Converging-Diverging Nozzle

Carbon dioxide is an attractive working fluid for many cycles, including for pumped thermal energy storage (PTES). A challenge with some proposed sCO2 PTES cycles is the operation of sCO2 machinery outside the typical bounds of experience, with local phase change from the liquid state being particularly unknown. Presently, there is insufficient data in the literature regarding multiphase CO2 to adequately design a multiphase-tolerant turbine, so generation of foundational data is required. This experimental study investigates the flow characteristics of sub-sonic liquid CO2 undergoing expansion and phase change in a converging-diverging nozzle. The nozzle is instrumented to measure static pressure, unsteady pressure, temperature, and density. The static pressure transducers are located at 27 axial locations to accurately characterize the pressure profile in the nozzle. High-accuracy RTDs are located at the entrance and exit of the nozzle, and three dynamic pressure transducers are strategically located to capture any unsteady phenomena. During testing, values of mass flow and nozzle inlet pressure are swept to vary the pressure drop and fluid properties. The measured total pressure drop in the nozzle is compared to a homogenous model and the Lockhart-Martinelli correlation method, with the latter predicting loss quite closely. The resulting data set is valuable for validating multiphase numerical models in a simple geometry before implementation of these models in turbomachinery design.

25 ENERGY STORAGE↗
NASA NTRS1979

Equatorial wave-mean flow interaction - A numerical study of the role of latitudinal shear

A time-dependent primitive equation model for an equatorial channel is used to assess the interaction of equatorial Kelvin and mixed Rossby-gravity waves with the mean flow. The proposed model involves a semiimplicit time-differencing scheme and a finite-difference grid in the meridional plane. It is shown that forced equatorial waves interact with mean flow to produce equatorial jets characterized by downward-moving westerly (Kelvin wave forcing) and easterly (mixed Rossby-gravity wave forcing) shear zones, respectively. For parameters characteristic of the observed waves in the equatorial stratosphere, the wave-mean flow interaction process always reduces the amplitude of any initial cross-equatorial mean wind shear. The mean flow profile tends to become symmetric about the equator as the interaction process continues.

Holton, J. R.↗
NASA NTRS1999

Analysis of X-Ray Microradiographs of Al-Au Interface Quench Profile using Modeling of Solidification Including Double-Diffusion and Convection in the Melt

Experimental data on Al-0.8Au horizontal solidification of a 1 mm thick specimen in a BN crucible shows the effect of growth rate on the solidification interface shape. For translation rates below 0.5 micron/s the interface maintains a plain and flat shape. When the translation rate is 3 to 5 micron/s or more, the interface appearance changes to two planar zones, with the zone closer to the bottom having higher inclination. The interface shapes were measured by first quenching in place during growth. X-ray microscopy shows the interface shape within the quenched sample by viewing through the side of the specimen. In order to provide theoretical explanation of the phenomena, numerical modeling was undertaken using finite element code FIDAP. Double diffusion convection in Al-0.8Au melt and crystal-melt interface curvature during directional solidification was analyzed numerically. Actual thermophysical properties of Al-0.8Au including the binary Al-Au phase diagram were used. Although convection in the sample is weak, for the slower translation rate convection and diffusion is sufficient for the redistribution of initial compositional stratification caused by gravity. When translation rate is raised, neither convection nor diffusion can provide proper mixing so that solidification temperatures differ significantly near the bottom within the bulk of the sample. As a result, the solid-liquid interface appears to have two planar zones with different inclination.

Bune, Andris V.↗
DOE OSTI2022

Quantitative interpretation of time-lapse seismic data at Farnsworth field unit: Rock physics modeling, and calibration of simulated time-lapse velocity responses

Here, this study investigates the contribution of fluid saturation variation to the time-lapse velocity response by performing fluid substitution modeling. The methodology is exemplified by the time-lapse seismic monitoring of carbon dioxide at Farnsworth field unit (FWU). In order to evaluate the fluid distribution in a matured oil reservoir, the Southwest Regional Partnership (SWP) acquired multiple vertical seismic profile (VSP) surveys at different times during the CO 2 –water alternatinggas (WAG) injection period. In this work, we present a thorough methodology for computing the elastic response of the saturated rock for different fluid saturations using a site-specific petro-elastic model (PEM). The output from the PEM was combined with results from a fluid compositional model to compute the seismic velocities at times corresponding to each VSP survey. To produce a calibrated simulated response, the measured time-lapse seismic velocities were integrated into the numerical simulation model. The mismatches between the predicted and measured time-lapse velocities were minimized through an iterative calibration process using a trained artificial neural network proxy (ANN) coupled with a particle swarm optimizer (PSO). Our study indicates that the hybrid optimization workflow can effectively perform the history matching. With an accurate prediction of the hydrodynamic properties, the migration of CO 2 within the subsurface was modeled by predicting the spatial velocity distribution for a radius of 305 m around the injection well. The technology demonstrated and the expertise gained from this study can guide similar CO 2 -WAG projects.

58 GEOSCIENCES↗