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 181 records · Page 10
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Sensitivity of thermodynamic profiles retrieved from ground-based microwave and infrared observations to additional input data from active remote sensing instruments and numerical weather prediction models

Accurate and continuous estimates of the thermodynamic structure of the lower atmosphere are highly beneficial to meteorological process understanding and its applications, such as weather forecasting. In this study, the Tropospheric Remotely Observed Profiling via Optimal Estimation (TROPoe) physical retrieval is used to retrieve temperature and humidity profiles from various combinations of input data collected by passive and active remote sensing instruments, in situ surface platforms, and numerical weather prediction models. Among the employed instruments are microwave radiometers (MWRs), infrared spectrometers (IRSs), radio acoustic sounding systems (RASSs), ceilometers, and surface sensors. TROPoe uses brightness temperatures and/or radiances from MWRs and IRSs, as well as other observational inputs (virtual temperature from the RASS, cloud-base height from the ceilometer, pressure, temperature, and humidity from the surface sensors) in a physical iterative retrieval approach. This starts from a climatologically reasonable profile of temperature and water vapor, with the radiative transfer model iteratively adjusting the assumed temperature and humidity profiles until the derived brightness temperatures and radiances match those observed by the MWR and/or IRS instruments within a specified uncertainty, as well as within the uncertainties of the other observations, if used as input. In this study, due to the uniqueness of the dataset that includes all the above-mentioned sensors, TROPoe is tested with different observational input combinations, some of which also include information higher than 4 km above ground level (a.g.l.) from the operational Rapid Refresh numerical weather prediction model. These temperature and humidity retrievals are assessed against independent collocated radiosonde profiles under non-cloudy conditions to assess the sensitivity of the TROPoe retrievals to different input combinations.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRSDate not supplied

Optimal Angle of Attack Control for Aerocapture Maneuvers

Recent investigations of aerocapture maneuvers have considered direct force control as a promising alternative to bank angle control. By modulating the angle of attack and sideslip angle independently, the longitudinal and lateral channels can be decoupled and thus solved separately. However, the optimal angle of attack profile for minimizing post-atmospheric propellant expenditure was previously unknown. This research applies principles of optimal control theory to numerically solve for the optimal angle of attack profile for control of the longitudinal channel. Direct methods are first employed to discretize the angle of attack profile and show that an optimal solution follows a nearly bang-bang (full lift up then full lift down) structure under a variety of conditions along the flight path angle corridor. To corroborate these results, indirect methods were also applied to scope the Optimal Control Problem (OCP) as a Two Point Boundary Value Problem (2P-BVP) which could then be solved numerically. The indirect optimization results aligned closely with the direct results, providing further evidence that a bang-bang structure provides a good approximation for the optimal angle of attack profile for aerocapture maneuvers using direct force control.

aerocapture↗
NASA NTRS1998

Numerical Simulation of the Generation of Axisymmetric Mode Jet Screech Tones

An imperfectly expanded supersonic jet, invariably, radiates both broadband noise and discrete frequency sound called screech tones. Screech tones are known to be generated by a feedback loop driven by the large scale instability waves of the jet flow. Inside the jet plume is a quasi-periodic shock cell structure. The interaction of the instability waves and the shock cell structure, as the former propagates through the latter, is responsible for the generation of the tones. Presently, there are formulas that can predict the tone frequency fairly accurately. However, there is no known way to predict the screech tone intensity. In this work, the screech phenomenon of an axisymmetric jet at low supersonic Mach number is reproduced by numerical simulation. The computed mean velocity profiles and the shock cell pressure distribution of the jet are found to be in good agreement with experimental measurements. The same is true with the simulated screech frequency. Calculated screech tone intensity and directivity at selected jet Mach number are reported in this paper. The present results demonstrate that numerical simulation using computational aeroacoustics methods offers not only a reliable way to determine the screech tone intensity and directivity but also an opportunity to study the physics and detailed mechanisms of the phenomenon by an entirely new approach.

Shen, Hao↗
NASA NTRS1980

Nitric oxide delta band absorption measurements in the lower thermosphere

The first successful measurements of nitric oxide band absorption in the earth's atmosphere have been made by using the high-resolution ultraviolet spectrometers on OSO 8. These measurements were conducted between January and August 1977 and between 30 deg north and south latitudes. The nitric oxide absorption data, at the (1,0) delta band wavelength of 1829.41 A, are numerically inverted to produce nitric oxide density profiles between 85 and 145 km. The major finding of the observations is that for altitudes below 95 km there is more nitric oxide at sunrise than at sunset. A numerical model, incorporating time-dependent ionospheric chemistry and vertical diffusive transport (both eddy and molecular), is used to simulate the OSO 8 nitric oxide profiles. The peak of the model sunset profile lies at 95-km altitude, in agreement with observation, if the vertical component of the eddy diffusion coefficient is 10 million sq cm/s near 100 km.

Massie, S. T.↗
NASA NTRS1997

Noise from Supersonic Coaxial Jets: Mean Flow Predictions - Part 1

Recent theories for supersonic jet noise have used an instability wave noise generation model to predict radiated noise. This model requires a known mean flow that has typically been described by simple analytic functions for single jet mean flows. The mean flow of supersonic coaxial jets is not described easily in terms of analytic functions. To provide these profiles at all axial locations, a numerical scheme is developed to calculate the mean flow properties of a coaxial jet. The Reynolds-averaged, compressible, parabolic boundary layer equations are solved using a mixing length turbulence model. Empirical correlations are developed to account for the effects of velocity and temperature ratios and Mach number on the shear layer spreading. Both normal velocity profile and inverted velocity profile coaxial jets are considered. The mixing length model is modified in each case to obtain reasonable results when the two stream jet merges into a single fully developed jet. The mean flow calculations show both good qualitative and quantitative agreement with measurements in single and coaxial jet flows.

Dahl, Milo D.↗
DOE OSTI2022

Sensitivity Analysis for Solutions to Heterogeneous Nonlocal Systems. Theoretical and Numerical Studies

The paper presents a collection of results on continuous dependence for solutions to nonlocal problems under perturbations of data and system parameters. The integral operators appearing in the systems capture interactions via heterogeneous kernels that exhibit different types of weak singularities, space dependence, even regions of zero-interaction. Here, the stability results showcase explicit bounds involving the measure of the domain and of the interaction collar size, nonlocal Poincaré constant, and other parameters. In the nonlinear setting, the bounds quantify in different L p norms the sensitivity of solutions under different nonlinearity profiles. The results are validated by numerical simulations showcasing discontinuous solutions, varying horizons of interactions, and symmetric and heterogeneous kernels.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1995

Short-range precipitation forecasts using assimilation of simulated satellite water vapor profiles and column cloud liquid water amounts

These observing system simulation experiments investigate the assimilation of satellite-observed water vapor and cloud liquid water data in the initialization of a limited-area primitive equations model with the goal of improving short-range precipitation forecasts. The assimilation procedure presented includes two aspects: specification of an initial cloud liquid water vertical distribution and diabatic initialization. The satellite data is simulated for the next generation of polar-orbiting satellite instruments, the Advanced Microwave Sounding Unit (AMSU) and the High-Resolution Infrared Sounder (HIRS), which are scheduled to be launched on the NOAA-K satellite in the mid-1990s. Based on cloud-top height and total column cloud liquid water amounts simulated for satellite data a diagnostic method is used to specify an initial cloud water vertical distribution and to modify the initial moisture distribution in cloudy areas. Using a diabatic initialization procedure, the associated latent heating profiles are directly assimilated into the numerical model. The initial heating is estimated by time averaging the latent heat release from convective and large-scale condensation during the early forecast stage after insertion of satellite-observed temperature, water vapor, and cloud water formation. The assimilation of satellite-observed moisture and cloud water, together withy three-mode diabatic initialization, significantly alleviates the model precipitation spinup problem, especially in the first 3 h of the forecast. Experimental forecasts indicate that the impact of satellite-observed temperature and water vapor profiles and cloud water alone in the initialization procedure shortens the spinup time for precipitation rates by 1-2 h and for regeneration of the areal coverage by 3 h. The diabatic initialization further reduces the precipitation spinup time (compared to adiabatic initialization) by 1 h.

Wu, Xiaohua↗
NASA NTRS2004

Source Term Model for Vortex Generator Vanes in a Navier-Stokes Computer Code

A source term model for an array of vortex generators was implemented into a non-proprietary Navier-Stokes computer code, OVERFLOW. The source term models the side force created by a vortex generator vane. The model is obtained by introducing a side force to the momentum and energy equations that can adjust its strength automatically based on the local flow. The model was tested and calibrated by comparing data from numerical simulations and experiments of a single low profile vortex generator vane on a flat plate. In addition, the model was compared to experimental data of an S-duct with 22 co-rotating, low profile vortex generators. The source term model allowed a grid reduction of about seventy percent when compared with the numerical simulations performed on a fully gridded vortex generator on a flat plate without adversely affecting the development and capture of the vortex created. The source term model was able to predict the shape and size of the stream-wise vorticity and velocity contours very well when compared with both numerical simulations and experimental data. The peak vorticity and its location were also predicted very well when compared to numerical simulations and experimental data. The circulation predicted by the source term model matches the prediction of the numerical simulation. The source term model predicted the engine fan face distortion and total pressure recovery of the S-duct with 22 co-rotating vortex generators very well. The source term model allows a researcher to quickly investigate different locations of individual or a row of vortex generators. The researcher is able to conduct a preliminary investigation with minimal grid generation and computational time.

Waithe, Kenrick A.↗
NASA NTRS2002

Numerical Simulations of Vortex Generator Vanes and Jets on a Flat Plate

Numerical simulations of a single low-profile vortex generator vane, which is only a small fraction of the boundary-layer thickness, and a vortex generating jet have been performed for flows over a flat plate. The numerical simulations were computed by solving the steady-state solution to the Reynolds-averaged Navier-Stokes equations. The vortex generating vane results were evaluated by comparing the strength and trajectory of the streamwise vortex to experimental particle image velocimetry measurements. From the numerical simulations of the vane case, it was observed that the Shear-Stress Transport (SST) turbulence model resulted in a better prediction of the streamwise peak vorticity and trajectory when compared to the Spalart-Allmaras (SA) turbulence model. It is shown in this investigation that the estimation of the turbulent eddy viscosity near the vortex core, for both the vane and jet simulations, was higher for the SA model when compared to the SST model. Even though the numerical simulations of the vortex generating vane were able to predict the trajectory of the stream-wise vortex, the initial magnitude and decay of the peak streamwise vorticity were significantly under predicted. A comparison of the positive circulation associated with the streamwise vortex showed that while the numerical simulations produced a more diffused vortex, the vortex strength compared very well to the experimental observations. A grid resolution study for the vortex generating vane was also performed showing that the diffusion of the vortex was not a result of insufficient grid resolution. Comparisons were also made between a fully modeled trapezoidal vane with finite thickness to a simply modeled rectangular thin vane. The comparisons showed that the simply modeled rectangular vane produced a streamwise vortex which had a strength and trajectory very similar to the fully modeled trapezoidal vane.

Allan, Brian G.↗
DOE OSTI2021

Logarithmic profile of temperature in sheared and unstably stratified atmospheric boundary layers

The impact of buoyancy on the mean velocity, temperature, and scalar concentration profiles in the lower atmosphere is typically investigated within the framework of Monin-Obukhov similarity theory (MOST). MOST is the theoretical foundation for parametrizing surface-atmosphere exchanges in nearly all weather, climate, and hydrological models.According to MOST, the classic logarithmic profiles of mean velocity and temperature break down as the buoyancy effects become important. However, recent studies on turbulent Rayleigh-Bénard convection and natural convection along vertical walls suggest that the mean temperature in the near-surface region still follows a logarithmic profile.Motivated by these new results, we study the mean potential temperature profile in sheared and unstably stratified atmospheric boundary layers using direct numerical simulations and field observations. We find that the mean potential temperature profile remains logarithmic across a wide range of stability parameters, which characterizes the relative importance of buoyancy versus shear effects. Compared to MOST, our results suggest that the buoyancy force does not modify the logarithmic nature of the mean potential temperature profile, but instead modulates its slope, which is no longer universal and differs from 1/κ,whereκisthe von Kármán constant. Furthermore, this study provides another perspective on scalar turbulence in the atmospheric boundary layer

58 GEOSCIENCES↗
NASA NTRS1972

Numerical solution of laminar jet mixing with and without free stream.

Systematic numerical solutions of two-dimensional and axisymmetrical laminar jet of an incompressible fluid with and without free stream have been obtained. For the case without free stream, the exact numerical solutions have been obtained for different initial velocity profiles and compared with experimental results and similarity solutions. For the case with free stream, the numerical solutions have been compared with the linearized analytical solutions. An approximate numerical solution for the axial velocity distribution of the three-dimensional laminar jet is proposed. The accuracy of the method has been determined from the exact solutions of the two limiting cases of three-dimensional jets, i.e., the two-dimensional and axisymmetric cases.

Pai, S. I.↗
NASA NTRS1993

Correcting for diffraction in the far-infrared reflectance measurement of rough surfaces

The manner in which diffraction can increase the reflectance measured from very rough surfaces is shown by observations of diffracted light within the instrument profile of a far-IR reflectometer system. A correction to the calibration signal based on numerical integration of the diffracted part of the instrument profile is described. Diffraction correction factors as large as 2.94 have been found with small optics at long wavelength (630 microns). The effect of diffraction on diffuse reflectance measurements of a very rough perfect reflector is shown at wavelengths from 56 to 200 microns.

Smith, Sheldon M.↗
NASA NTRS2017

Simple Analytical Expressions for Electron Pitch Angle Diffusion Coefficients

In the present paper, we have calculated electron pitch angle diffusion coefficients due to resonant interactions with whistler mode lower band chorus (LBC), upper band chorus (UBC), and electrostatic electron cyclotron harmonic (ECH) waves. Calculations have been performed at two values of the ratio of electron plasma frequency to gyro-frequency and thirteen representative values of electron energies for the plasma sheet electrons. The numerical data of diffusion coefficients have been fitted to simple analytical expressions for each wave mode. These analytical expressions allow for simple evaluation of pitch angle diffusion coefficients for the arbitrary pitch angle, energy, the ambient magnetic field, the wave amplitude, and the ratio of plasma frequency to gyro-frequency. In the case of LBC waves, the analytical coefficients are generally within a factor of two to the numerical coefficients, except at higher pitch angles where the numerical coefficients drop to show negligible values. Likewise, also for UBC waves, the analytical coefficients are generally within a factor of two to the numerical coefficients. The analytical coefficients for ECH waves are generally in agreement with numerical coefficients. However, the analytical expressions developed for ECH waves do not reproduce sharp fluctuations, dips, and gaps, the so-called banded structure observed in the data of numerical coefficients. The effect of the shape variation on the profile of the wave spectral intensity and the effect of cold temperature on the numerical coefficients are also investigated. Applications of the analytical expressions of diffusion coefficients are discussed.

coefficients↗
NASA NTRS1975

Calculation of three-dimensional compressible laminar and turbulent boundary layers. An implicit finite-difference procedure for solving the three-dimensional compressible laminar, transitional, and turbulent boundary-layer equations

An implicit finite-difference procedure is presented for solving the compressible three-dimensional boundary-layer equations. The method is second-order accurate, unconditionally stable (conditional stability for reverse cross flow), and efficient from the viewpoint of computer storage and processing time. The Reynolds stress terms are modeled by (1) a single-layer mixing length model and (2) a two-layer eddy viscosity model. These models, although simple in concept, accurately predicted the equilibrium turbulent flow for the conditions considered. Numerical results are compared with experimental wall and profile data for a cone at an angle of attack larger than the cone semiapex angle. These comparisons clearly indicate that the numerical procedure and turbulence models accurately predict the experimental data with as few as 21 nodal points in the plane normal to the wall boundary.

Harris, J. E.↗
NASA NTRS1991

Shot model parameters for Cygnus X-1 through phase portrait fitting

Shot models for systems having about 1/f power density spectrum are developed by utilizing a distribution of shot durations. Parameters of the distribution are determined by fitting the power spectrum either with analytic forms for the spectrum of a shot model with a given shot profile, or with the spectrum derived from numerical realizations of trial shot models. The shot fraction is specified by fitting the phase portrait, which is a plot of intensity at a given time versus intensity at a delayed time and in principle is sensitive to different shot profiles. These techniques have been extensively applied to the X-ray variability of Cygnus X-1, using HEAO 1 A-2 and an Exosat ME observation. The power spectra suggest models having characteristic shot durations lasting from milliseconds to a few seconds, while the phase portrait fits give shot fractions of about 50 percent. Best fits to the portraits are obtained if the amplitude of the shot is a power-law function of the duration of the shot. These fits prefer shots having a symmetric exponential rise and decay. Results are interpreted in terms of a distribution of magnetic flares in the accretion disk.

Lochner, James C.↗
DOE OSTI2026

Impact of in situ nuclear networks and atomic opacities on neutron star merger ejecta dynamics, nucleosynthesis, and kilonovae

Context. Binary neutron star merger (BNSM) ejecta are key sites of rapid neutron capture (r-process) nucleosynthesis and they produce kilonovae powered by the radioactive decay of freshly synthesized nuclei. Modeling their evolution requires multi-physics simulations involving hydrodynamics, nuclear reactions, and radiative processes. The impact of nuclear burning and atomic opacity is poorly understood and often treated with simplified prescriptions. Aims. We systematically investigate different treatments of nuclear heating, particle thermalization, and atomic opacities in radiation-hydrodynamics simulations of BNSM ejecta and kilonova light curves. Methods. Ejecta profiles from long-term numerical-relativity simulations of asymmetric neutron star binaries with a massive neutron star remnant were evolved to ∼30 days using a 2D ray-by-ray approach. We compared simplified heating-rate and thermalization prescriptions with in situ Nuclear reaction Network (NN) calculations that track nuclear energy deposition and include a composition-dependent thermalization scheme. We also contrasted various gray opacity models with a frequency-dependent treatment based on atomic calculations. Results. Coupling NN and hydrodynamics significantly affects nucleosynthesis and kilonova emission. Assuming homologous expansion alters abundance evolution and produces a narrower, less populated second r-process peak and a third peak shifted to higher mass numbers. The back-reaction of nuclear heating affects the temperature evolution enough to delay and redden the early (t∼ hours) kilonova peaks. A constant thermalization efficiency underestimates and reddens the early emission while overestimating the late-time luminosity compared to the composition-dependent treatment. Analytical opacity prescriptions yield a more extended, colder photosphere, resulting in dimmer, redder kilonovae at early times (t≲ hour), while the delayed recession of the photosphere prolongs the red emission at t ≳ 5 days. Conclusions. Coupling hydrodynamics to an in situ NN is crucial for reliable nucleosynthesis and kilonova predictions. Resolving the first several hundred milliseconds of the hydrodynamics is essential for robust nucleosynthesis calculations. Composition-dependent thermalization and frequency-dependent, atomic-physics-based opacities are needed to accurately capture the temperature evolution of the ejecta and the brightness and color evolution of the kilonova. Calibrated analytic nuclear-power fits with simplified thermalization and opacity prescriptions can still reproduce the density and temperature evolution of the ejecta.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2024

On the accuracy of compressibility transformations

This study highlights the importance of satisfying the eddy viscosity equivalence below the logarithmic layer, to deriving accurate compressibility transformations. First, we analyze the ability of known transformations to satisfy the eddy viscosity equivalence and show that the accuracy of these transformations is strongly dependent on this ability. Second, in a step-by-step manner, we devise new transformations that satisfy this hypothesis. An approach based on curve fitting of the incompressible Direct Numerical Simulation data for eddy viscosity profiles below the logarithmic layer provides an extremely accurate transformation, which motivates self-contained methods, making use of mixing length formulas in the inner region. It is shown that the accuracy of existing transformations can be significantly improved by applying these ideas, below the logarithmic layer. Motivated by the effectiveness of the formulations derived from eddy viscosity equivalence, we introduce a new integral transformation based on Reynolds number equivalence between compressible and incompressible flows. This approach is based on defining a new compressible velocity scale, which affects the accuracy of transformations. Several choices for the velocity scale are tested, and in each attempt, it is shown that the eddy viscosity equivalence plays a very important role for the accuracy of compressibility transformations.

42 ENGINEERING↗
NASA NTRS1972

6300 A quantum efficiency of the recombination mechanism in the night-time F layer.

Simultaneous airglow and electron content measurements made at Hawaii are used to infer the number of 6300 and 6364 A quanta produced per electron lost in the nighttime F layer of the ionosphere. The equation of continuity of electrons is then solved numerically to obtain the electron density profile, and the amount of quenching is estimated. This leads to the number of excitations of O(super-1 D) per O2(+) recombination (epsilon). We find, for an exospheric temperature of 1100 K, epsilon is equal to 1.1 plus or minus 0.6, in good agreement with Zipf's laboratory measurement at 300 K.

Brown, W. E.↗