Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Artificial diffusivity”

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 37 records · Page 2

Leading-edge receptivity for blunt-nose bodies

This research program investigates boundary-layer receptivity in the leading-edge region for bodies with blunt leading edges. Receptivity theory provides the link between the unsteady distrubance environment in the free stream and the initial amplitudes of the instability waves in the boundary layer. This is a critical problem which must be addressed in order to develop more accurate prediction methods for boundary-layer transition. The first phase of this project examines the effects of leading-edge bluntness and aerodynamic loading for low Mach number flows. In the second phase of the project, the investigation is extended to supersonic Mach numbers. Singular perturbation techniques are utilized to develop an asymptotic theory for high Reynolds numbers. In the first year, the asymptotic theory was developed for leading-edge receptivity in low Mach number flows. The case of a parabolic nose is considered. Substantial progress was made on the Navier-Sotkes computations. Analytical solutions for the steady and unsteady potential flow fields were incorporated into the code, greatly expanding the types of free-stream disturbances that can be considered while also significantly reducing the the computational requirements. The time-stepping algorithm was modified so that the potential flow perturbations induced by the unsteady pressure field are directly introduced throughout the computational domain, avoiding an artificial 'numerical diffusion' of these from the outer boundary. In addition, the start-up process was modified by introducing the transient Stokes wave solution into the downstream boundary conditions.

Kerschen, Edward J.↗

Some Basic Laws of Isotropic Turbulent Flow

An Investigation is made of the diffusion of artificially produced turbulence behind screens or other turbulence producers. The method is based on the author's concept of disturbance moment as a certain theoretically well-founded measure of turbulent disturbances.

Loitsianskii, L. G.↗

Snow Depth from AMSR-2 Using Multispectral Satellite Data in an Artificial Neural Network

By using diffusion theory and Monte Carlo lidar radiative transfer simulations, Hu et al. (2022b) has derived snow depth from the first-, second- and third-order moments of the lidar backscattering pathlength distribution. Lu et al. (2022) calculated the snow depth by applying the methods to the satellite ICESat-2 lidar measurements over the Arctic sea ice, as well as land surfaces of Northern Hemisphere. In this paper, an artificial neural network (ANN) algorithm, employing several channels from Advanced Microwave Scanning Radiometer 2 (AMSR-2) and the humidity vertical profiles from Global Modeling and Assimilation Office (GMAO) Goddard Earth Observing System for Instrument Teams (GEOS-IT) product, is trained to determine snow depth identified by time and geolocation matched 2019 ICESat-2 snow-depth data during winter months over the Arctic sea ice. The trained ANN snow-depth was applied to 2018 AMSR-2 clear pixel data, although the algorithms perform reasonably well in thinner clouds. The validation data (different from the training set) of ANN snow depth from AMSR-2 showed a good agreement with time matched and co-located snow-depth values from ICESat-2. The bias was near zero, with mean absolute error (MAE) 0.05 cm and a root-mean-square-error (RMSE) 0.08 cm. Prior applying the trained ANN snow depth to AMSR-2 data, a cloud screening algorithm was developed with a similar approach. A separate ANN cloud mask was trained to determine an AMSR-2 pixel is clear or cloudy with time and geolocation matched 2015 CALIOP Vertical Feature Mask (VFM) over Arctic sea ice. The ANN cloud mask from AMSR-2 under-estimated cloud fraction by 3-6% compared to CALIOP . The additional research is needed to conclusively evaluate the ANN cloud mask accuracy. Finally, this paper will lay the foundation for a sustained long-term snowfall and snow-storm monitoring system. The future Cloud Aerosol LIdar for Global scale Observations of the ocean-Land Atmosphere system (CALIGOLA) mission will provide a means to calculate snow depth from the lidar backscattering pathlength distribution, benefiting from the UV, visible and infrared pulses. With the calculated snow depth as the truth one could develop a machine learning algorithm, as it was done in this paper, using a passive microwave instrument available at that time to generate a wide range of snow depth data, covering extensive spatial areas in the cross-orbit direction.

Neural Network↗

Propagation of Secondary Antiprotons and Cosmic Rays in the Galaxy

Recent more accurate antiproton data obtained by the BESS team during the last solar minimum pose a challenge to conventional propagation models of cosmic rays. In particular, the diffusive reacceleration model, which matches well key secondary/primary isotope ratios in cosmic rays, fails to reproduce the secondary antiproton spectrum. Tuning both secondary/primary isotope ratios and antiprotons is possible, but requires artificial breaks in the diffusion coefficient and the injection spectrum of primaries. We will discuss some possibilities to overcome these difficulties in the propagation models. We will present new results of our calculation of CR propagation in the Galaxy using the GALPROP code.

Moskalenko, I. V.↗

An experimental investigation of two large annular diffusers with swirling and distorted inflow

Two annular diffusers downstream of a nacelle-mounted fan were tested for aerodynamic performance, measured in terms of two static pressure recovery parameters (one near the diffuser exit plane and one about three diameters downstream in the settling duct) in the presence of several inflow conditions. The two diffusers each had an inlet diameter of 1.84 m, an area ratio of 2.3, and an equivalent cone angle of 11.5, but were distinguished by centerbodies of different lengths. The dependence of diffuser performance on various combinations of swirling, radially distorted, and/or azimuthally distorted inflow was examined. Swirling flow and distortions in the axial velocity profile in the annulus upstream of the diffuser inlet were caused by the intrinsic flow patterns downstream of a fan in a duct and by artificial intensification of the distortions. Azimuthal distortions or defects were generated by the addition of four artificial devices (screens and fences). Pressure recovery data indicated beneficial effects of both radial distortion (for a limited range of distortion levels) and inflow swirl. Small amounts of azimuthal distortion created by the artificial devices produced only small effects on diffuser performance. A large artificial distortion device was required to produce enough azimuthal flow distortion to significantly degrade the diffuser static pressure recovery.

Eckert, W. T.↗

ULTRA-SHARP nonoscillatory convection schemes for high-speed steady multidimensional flow

For convection-dominated flows, classical second-order methods are notoriously oscillatory and often unstable. For this reason, many computational fluid dynamicists have adopted various forms of (inherently stable) first-order upwinding over the past few decades. Although it is now well known that first-order convection schemes suffer from serious inaccuracies attributable to artificial viscosity or numerical diffusion under high convection conditions, these methods continue to enjoy widespread popularity for numerical heat transfer calculations, apparently due to a perceived lack of viable high accuracy alternatives. But alternatives are available. For example, nonoscillatory methods used in gasdynamics, including currently popular TVD schemes, can be easily adapted to multidimensional incompressible flow and convective transport. This, in itself, would be a major advance for numerical convective heat transfer, for example. But, as is shown, second-order TVD schemes form only a small, overly restrictive, subclass of a much more universal, and extremely simple, nonoscillatory flux-limiting strategy which can be applied to convection schemes of arbitrarily high order accuracy, while requiring only a simple tridiagonal ADI line-solver, as used in the majority of general purpose iterative codes for incompressible flow and numerical heat transfer. The new universal limiter and associated solution procedures form the so-called ULTRA-SHARP alternative for high resolution nonoscillatory multidimensional steady state high speed convective modelling.

Leonard, B. P.↗

Chemical and physical effects in the bulk of cometary analogs

KOSI comet simulation experiments were designed as a macroscopic tester for the studies of physicochemical problems inherent to comet bodies. The analog samples consist of H2O and CO2 ice, organic admixtures, mineral dust, and carbon. Two of the fundamental changes the analogs undergo when submitted to 'insolation' by artificial sunlight, i.e., the diffusion of frozen gases and subsequent crust formation and the natural isotopic fractionation, are reported.

Roessler, K.↗

Convergence properties of finite-difference hydrodynamics schemes in the presence of shocks

We investigate the asymptotic convergence of finite-difference schemes for the Euler equations when the limiting solution contains shocks. The Lax-Wendroff theorem guarantees that certain conservative schemes converge to correct, physically valid solutions. We focus on two one-dimensional operator-split schemes with explicit artificial-viscosity terms. One, an internal-energy scheme, does not satisfy the assumptions of Lax-Wendroff; the other, a conservative total-energy scheme, does. With viscous lengths chosen proportional to the grid size, we find that both schemes converge to their zero-grid-size limits at the theoretically expected rate, but only the conversative scheme converges toward correct solutions of the inviscid fluid equations. We show that the difference in their behaviors results directly from the presence of shocks in the limiting solution. Empirically, we find that when the viscous lenghts tend toward zero more slowly than the grid size, however the nonconservative scheme also converges toward correct solutions. We characterize the asymptotic behavior of the total-energy scheme in a particular problem in which a shock forms. As the grid is refined, a Cauchy error approaches the expected rate of change slowly. We show that the changes in the artificial viscosity alter the diffusion of small-amplitude waves. The differences associated with such waves make the dominant contribution to the Cauchy error. We formulate an analytic model to relate the rate of approach to the effect of varying diffusion in waves and find quantitative agreement with our numerical results.

Kimoto, Paul A.↗

Challenging Cosmic Ray Propagation with Antiprotons: Evidence for a "Fresh" Nuclei Component?

Recent measurements of the cosmic ray (CR) antiproton flux have been shown to challenge existing CR propagation models. It was shown that the reacceleration models designed to match secondary to primary nuclei ratio (e.g., Boron/Carbon) produce too few antiprotons, while the traditional non-reacceleration models can reproduce the antiproton flux but fall short of explaining the low-energy decrease in the secondary to primary nuclei ratio. Matching both the secondary to primary nuclei ratio and antiproton flux requires artificial breaks in the diffusion coefficient and the primary injection spectrum suggesting the need for other approaches. In the present paper we discuss one possibility to overcome these difficulties. Using the measured antiproton flux to fix the diffusion coefficient, we show that the spectra of primary nuclei as measured in the heliosphere may contain a fresh local unprocessed component at low energies, thus decreasing the measured secondary to primary nuclei ratio. A model reproducing antiprotons, B/C ratio, and abundances up to Ni is presented.

Moskalenko, Igor V.↗

Some Observations on Transitory Stall in Conical Diffusers

Results from an experimental investigation on the flow through conical diffusers are presented. The mean and fluctuating velocity fields are compared for three diffusers with total diffusion angles of 16, 20 and 24 degrees, in the throat Mach number (M sub t) range of 0.05 to 0.95. Each of the diffusers were 14 cm long and had a 5.08 cm inlet diameter, and the flow exited into the ambient. The boundary layer at the throat was thin with the throat diameter (D sub t) to momentum thickness (O) ratio being as high as 800 at M(sub t) = 0.4. While the 16 deg diffuser flow exited with a top-hat mean velocity profile, increasing losses due to increasing separation resulted in fuller profiles for the 20 and 24 degree cases. A detailed flow field study was conducted for the 16 deg. diffuser. The u'-spectrum, measured at the exit plane, exhibited a peak apparently due to the ensuing jet column instability throughout the M(sub t) range covered. In addition, a much lower frequency spectral peak also occurred in the M(sub t) range of 0.3 to 0.7. Both of the spectral peaks were due to axisymmetric flow fluctuations. A self-sustaining flow oscillation occurred in the M(sub t) range of 0.6 to 0.85, emitting a loud tone, when the jet column instability frequency matched the resonance frequency of the diffuser. Limited data showed that artificial acoustic excitation was effective in reducing the flow fluctuations, with a resultant increase in the pressure recovery, at low M(sub t).

Zaman, K. B. M. Q.↗

Some observations on transitory stall in conical diffusers

Results from an experimental investigation on the flow through conical diffusers are presented. The mean and fluctuating velocity fields are compared for three diffusers with total diffusion angles of 16, 20 and 24 degrees, in the throat Mach number (M sub t) range of 0.05 to 0.95. Each of the diffusers were 14 cm long and had a 5.08 cm inlet diameter, and the flow exited into the ambient. The boundary layer at the throat was thin with the throat diameter (D sub t) to momentum thickness (O) ratio being as high as 800 at M(sub t) = 0.4. While the 16 deg diffuser flow exited with a top-hat mean velocity profile, increasing losses due to increasing separation resulted in fuller profiles for the 20 and 24 degrees cases. A detailed flow field study was conducted for the 16 deg diffuser. The u'-spectrum, measured at the exit plane, exhibited a peak apparently due to the ensuing jet column instability throughout the M(sub t) range covered. In addition, a much lower frequency spectral peak also occurred in the M(sub t) range of 0.3 to 0.7. Both of the spectral peaks were due to axisymmetric flow fluctuations. A self-sustaining flow oscillation occurred in the M(sub t) range of 0.6 to 0.85, emitting a loud tone, when the jet column instability frequency matched the resonance frequency of the diffuser. Limited data showed that artificial acoustic excitation was effective in reducing the flow fluctuations, with a resultant increase in the pressure recovery, at low M(sub t).

Zaman, K. B. M. Q.↗

Euler equation computations for the flow over a hovering helicopter rotor

A numerical solution technique is developed for computing the flow field around an isolated helicopter rotor in hover. The flow is governed by the compressible Euler equations which are integrated using a finite volume approach. The Euler equations are coupled to a free wake model of the rotary wing vortical wake. This wake model is incorporated into the finite volume solver using a prescribed flow, or perturbation, technique which eliminates the numerical diffusion of vorticity due to the artificial viscosity of the scheme. The work is divided into three major parts: (1) comparisons of Euler solutions to experimental data for the flow around isolated wings show good agreement with the surface pressures, but poor agreement with the vortical wake structure; (2) the perturbation method is developed and used to compute the interaction of a streamwise vortex with a semispan wing. The rapid diffusion of the vortex when only the basic Euler solver is used is illustrated, and excellent agreement with experimental section lift coefficients is demonstrated when using the perturbation approach; and (3) the free wake solution technique is described and the coupling of the wake to the Euler solver for an isolated rotor is presented. Comparisons with experimental blade load data for several cases show good agreement, with discrepancies largely attributable to the neglect of viscous effects. The computed wake geometries agree less well with experiment, the primary difference being that too rapid a wake contraction is predicted for all the cases.

Roberts, Thomas Wesley↗

Free-wake analysis of a rotor in hover

A numerical method based on the axisymmetric, incompressible Navier-Stokes equations is combined with a lifting surface code to predict the vortex wake of hovering rotors. The lifting surface code, AMI Hover, is used to obtain the circulation distribution on the blade. This circulation distribution is fed into the Navier-Stokes code to compute the vortex wake under this specified circulation distribution. An iteration approach is used between these two codes to converge the circulation distribution and the shape of the vortex wake. A relaxation scheme is developed to resolve the instability encountered among the tip vortices. A reconcentration scheme is used to solve the diffusion problem due to the strong artificial viscosity. The results from the present method are compared with experimental data obtained by smoke-flow visualization and hot-wire measurements for several rotor blade configurations. The comparisons show that the present method is able to predict the complex wake system shed by a hovering rotor.

Chen, C. S.↗

The Light of the Night Sky in EUSO: Duty Cycle and Background

The EUSO detector, which will be installed on the International Space Station, will detect the light produced by the EAS generated by the UHECR when entering the Earth's atmosphere. The EUSO telescope will operate only during dark night. The light background in the EUSO detector is of outmost importance. To determine the detector duty cycle, it is necessary to know the periods in which it will be in the obscurity, taking into account both sun and moon precise ephemerides and the ISS trajectory. For signal to background determination and trigger efficiency study, the photon background has to be evaluated, taking into account the various UV light sources at night: natural sources as airglow, stars and diffuse light of the night sky, and artificial sources.

Adams, J.↗

An Evaluation of Size-Resolved Cloud Microphysics Scheme Numerics for use with Radar Observations Part II: Condensation and Evaporation

Accurate numerical modeling of clouds and precipitation is essential for weather forecasting and climate change research. While size-resolved (bin) cloud microphysics models predict particle size distributions without imposing shapes, results are subject to artificial size distribution broadening owing to numerical diffusion associated with various processes. Whereas Part 1 addressed collision-coalescence, here we investigate numerical diffusion that occurs in solving condensation and evaporation. In a parcel model framework, all of the numerical schemes examined converge to one solution of condensation and evaporation as the mass grid is refined, and the advection-based schemes are recommended over the reassigning schemes. Including Eulerian vertical advection in a column limits the convergence to some extent, but that limitation occurs at a sufficiently fine mass grid, and the number of iterations in solving vertical advection should be minimized to reduce numerical diffusion. Insubstantial numerical diffusion in solving condensation can be amplified if collision-coalescence is also active, which in turn can be substantially diminished if turbulence effects on collision are incorporated. Large-eddy simulations of a drizzling stratocumulus field reveal that changes in moments of Doppler spectra obtained using different mass grids are consistent with those obtained from the simplified framework, and that spectral moments obtained using a mass grid designed to effectively reduce numerical diffusion are generally closer to observations. Notable differences between the simulations and observations still exist, and our results suggest a need to consider whether factors other than numerical diffusion in the fundamental process schemes employed can cause such differences.

Cloud microphysics↗

The Spartan-281 Far Ultraviolet Imaging Spectrograph

The U.S. Naval Research Laboratory's Far Ultraviolet Imaging Spectrograph (FUVIS), currently under development for flight as a Spartan shuttle payload, is designed to perform spectroscopy of diffuse sources in the FUV with very high sensitivity and moderate spatial and spectral resolution. Diffuse nebulae, the general galactic background radiation, and artificially induced radiation associated with the Space Shuttle vehicle are sources of particular interest. The FUVIS instrument will cover the wavelength range of 970-2000 A with selectable resolutions of 5 and 30 A. It is a slit imaging spectrograph having 3 arcmin spatial resolution along its 2.7 deg long slit.

Carruthers, George R.↗