SEARCH · Engineering Papers

Results for “Smooth approximation”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 793 records · Page 44

All sources NASA NTRS DOE OSTI

NASA NTRS2015

Spatiotemporal Prediction of Fine Particulate Matter Using High-Resolution Satellite Images in the Southeastern US 2003-2011

Numerous studies have demonstrated that fine particulate matter (PM(sub 2.5), particles smaller than 2.5 micrometers in aerodynamic diameter) is associated with adverse health outcomes. The use of ground monitoring stations of PM(sub 2.5) to assess personal exposure, however, induces measurement error. Land-use regression provides spatially resolved predictions but land-use terms do not vary temporally. Meanwhile, the advent of satellite-retrieved aerosol optical depth (AOD) products have made possible to predict the spatial and temporal patterns of PM(sub 2.5) exposures. In this paper, we used AOD data with other PM(sub 2.5) variables, such as meteorological variables, land-use regression, and spatial smoothing to predict daily concentrations of PM(sub 2.5) at a 1 sq km resolution of the Southeastern United States including the seven states of Georgia, North Carolina, South Carolina, Alabama, Tennessee, Mississippi, and Florida for the years from 2003 to 2011. We divided the study area into three regions and applied separate mixed-effect models to calibrate AOD using ground PM(sub 2.5) measurements and other spatiotemporal predictors. Using 10-fold cross-validation, we obtained out of sample R2 values of 0.77, 0.81, and 0.70 with the square root of the mean squared prediction errors of 2.89, 2.51, and 2.82 cu micrograms for regions 1, 2, and 3, respectively. The slopes of the relationships between predicted PM2.5 and held out measurements were approximately 1 indicating no bias between the observed and modeled PM(sub 2.5) concentrations. Predictions can be used in epidemiological studies investigating the effects of both acute and chronic exposures to PM(sub 2.5). Our model results will also extend the existing studies on PM(sub 2.5) which have mostly focused on urban areas because of the paucity of monitors in rural areas.

PM2.5↗

DOE OSTI2024

Ensemble variational Fokker-Planck methods for data assimilation

Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Itˆo process. This work introduces the Variational Fokker-Planck (VFP) framework for data assimilation, a general approach that includes previously known particle flow filters as special cases. The McKean-Vlasov-Itˆo process that transforms particles is defined via an optimal drift that depends on the selected diffusion term. It is established that the underlying probability density - sampled by the ensemble of particles - converges to the Bayesian posterior probability density. For a finite number of particles the optimal drift contains a regularization term that nudges particles toward becoming independent random variables. Based on this analysis, we derive computationally-feasible approximate regularization approaches that penalize the mutual information between pairs of particles, and avoid particle collapse. Moreover, the diffusion plays a role akin to a particle rejuvenation approach that aims to alleviate particle collapse. The VFP framework is very flexible. Different assumptions on prior and intermediate probability distributions can be used to implement the optimal drift, and localization and covariance shrinkage can be applied to alleviate the curse of dimensionality. A robust implicit-explicit method is discussed for the efficient integration of stiff McKean- Vlasov-Itˆo processes. Here, the effectiveness of the VFP framework is demonstrated on three progressively more challenging test problems, namely the Lorenz ’63, Lorenz ’96 and the quasi-geostrophic equations.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1994

Magnetocentrifugally driven flows from young stars and disks. 3: Numerical solution of the sub-Alfvenic region

We construct steady, axisymmetric, numerical models of the sub-Alfvenic regions of cool winds driven by the X-celerator mechanism. We find that smooth acceleration to speeds of the order of the escape speed accompanied by substantial collimation of the flow can be achieved within the Alfven surface, located characteristically at several stellar radii. We apply the nondimensional results to winds which emerge from the equator of magnetized YSOs rotating at breakup and from circumstellar disks interacting with the magnetospheres of slowly rotating T Tauri stars. Stellar magnetic fields of kilogauss strength can drive wind mass-loss rates of 10(exp -6) solar mass/year from rapidly accreting YSOs and 10(exp -8) solar mass/year from slowly accreting T Tauri stars. The X-celerator mechanism can accelerate winds from these systems to velocities of hundreds of km/sec within a few stellar radii in the case of embedded YSOs and approximately 10 stellar radii in the case of T Tauri stars.

Najita, Joan R.↗

NASA NTRS1994

Nonaxisymmetric evolution in protostellar disks

We present a two-dimensional, multigridded hydrodynamical simulation of the collapse of an axisymmetric, rotating, 1 solar mass protostellar cloud, which forms a resolved, hydrotastic disk. The code includes the effects of physical viscosity, radiative transfer and radiative acceleration but not magnetic fields. We examine how the disk is affected by the inclusion of turbulent viscosity by comparing a viscous simulation with an inviscid model evolved from the same initial conditions, and we derive a disk evolutionary timescale on the order of 300,000 years if alpha = 0.01. Effects arising from non-axisymmetric gravitational instabilities in the protostellar disk are followed with a three-dimensional SPH code, starting from the two-dimensional structure. We find that the disk is prone to a series of spiral instabilities with primary azimulthal mode number m = 1 and m = 2. The torques induced by these nonaxisymmetric structures elicit material transport of angular momentum and mass through the disk, readjusting the surface density profile toward more stable configurations. We present a series of analyses which characterize both the development and the likely source of the instabilities. We speculate that an evolving disk which maintains a minimum Toomre Q-value approximately 1.4 will have a total evolutionary span of several times 10(exp 5) years, comparable to, but somewhat shorter than the evolutionary timescale resulting from viscous turbulence alone. We compare the evolution resulting from nonaxisymmetric instabilities with solutions of a one-dimensional viscous diffusion equation applied to the initial surface density and temperature profile. We find that an effective alpha-value of 0.03 is a good fit to the results of the simulation. However, the effective alpha will depend on the minimum Q in the disk at the time the instability is activated. We argue that the major fraction of the transport characterized by the value of alpha is due to the action of gravitational torques, and does not arise from inherent viscosity within the smoothed particle hydrodynamics method.

Laughlin, Gregory↗

NASA NTRS1990

An adaptive grid algorithm for one-dimensional nonlinear equations

Richards' equation, which models the flow of liquid through unsaturated porous media, is highly nonlinear and difficult to solve. Step gradients in the field variables require the use of fine grids and small time step sizes. The numerical instabilities caused by the nonlinearities often require the use of iterative methods such as Picard or Newton interation. These difficulties result in large CPU requirements in solving Richards equation. With this in mind, adaptive and multigrid methods are investigated for use with nonlinear equations such as Richards' equation. Attention is focused on one-dimensional transient problems. To investigate the use of multigrid and adaptive grid methods, a series of problems are studied. First, a multigrid program is developed and used to solve an ordinary differential equation, demonstrating the efficiency with which low and high frequency errors are smoothed out. The multigrid algorithm and an adaptive grid algorithm is used to solve one-dimensional transient partial differential equations, such as the diffusive and convective-diffusion equations. The performance of these programs are compared to that of the Gauss-Seidel and tridiagonal methods. The adaptive and multigrid schemes outperformed the Gauss-Seidel algorithm, but were not as fast as the tridiagonal method. The adaptive grid scheme solved the problems slightly faster than the multigrid method. To solve nonlinear problems, Picard iterations are introduced into the adaptive grid and tridiagonal methods. Burgers' equation is used as a test problem for the two algorithms. Both methods obtain solutions of comparable accuracy for similar time increments. For the Burgers' equation, the adaptive grid method finds the solution approximately three times faster than the tridiagonal method. Finally, both schemes are used to solve the water content formulation of the Richards' equation. For this problem, the adaptive grid method obtains a more accurate solution in fewer work units and less computation time than required by the tridiagonal method. The performance of the adaptive grid method tends to degrade as the solution process proceeds in time, but still remains faster than the tridiagonal scheme.

Gutierrez, William E.↗

NASA NTRS1998

Horizontal Radiative Fluxes in Clouds at Absorbing Wavelengths

We discuss the effect of horizontal fluxes on the accuracy of a conventional plane-parallel radiative transfer calculation for a single pixel, known as the Independent Pixel Approximation (IPA) at absorbing wavelengths. Vertically integrated horizontal fluxes can be represented as a sum of three components; each component is the IPA accuracy on a pixel-by-pixel basis for reflectance, transmittance and absorptance, respectively. We show that IPA accuracy for reflectance always improves with more absorption, while the IPA accuracy for transmittance is less sensitive to the changes in absorption: with respect to the non-absorbing case, it may first deteriorate for weak absorption and then improve again for strongly absorbing wavelengths. EPA accuracy for absorptance always deteriorates with more absorption. As a result, vertically integrated horizontal fluxes, as a sum of IPA accuracies for reflectance, transmittance and absorptance, increase with more absorption. Finally, the question of correlations between horizontal fluxes, IPA uncertainties and radiative smoothing is addressed using wavenumber spectra of radiation fields reflected from or transmitted through fractal clouds.

Marshak, A.↗

NASA NTRS2021

ARTEMIS Observations of Plasma Waves in Laminar and Perturbed Interplanetary Shocks

The ‘Acceleration, Reconnection, Turbulence and Electrodynamics of the Moon's Interaction with the Sun’ (ARTEMIS) mission provides a unique opportunity to study the structure of interplanetary shocks and the associated generation of plasma waves with frequencies between ~50-8000 Hz due to its long duration electric and magnetic field burst waveform captures. We compare wave properties and occurrence rates at 11 quasi-perpendicular interplanetary shocks with burst data within 10 minutes (~3200 proton gyroradii upstream, ~1900 downstream) of the shock ramp. A perturbed shock is defined as possessing a large amplitude whistler precursor in the quasi-static magnetic field with an amplitude greater than ⅓ the difference between the upstream and downstream average magnetic field magnitudes; laminar shocks lack these large precursors and have a smooth, step function-like transition. In addition to wave modes previously observed, including ion acoustic, whistler, and electrostatic solitary waves, waves in the ion acoustic frequency range that show rapid temporal frequency change are common. Three shocks had burst captures in the ramp; of these, the two laminar shocks with burst data in the ramp contained a wide range of large amplitude wave modes in the ramp whereas the one perturbed shock contained no such waves. Thus, energy dissipation through wave-particle interactions is more prominent in these two laminar shocks than in the perturbed shock. Based on observations from all 11 shocks, the The wave occurrence rates for laminar shocks are higher in the transition region, especially the ramp, than downstream. In contrast, perturbed shocks have approximately 2-3 times the wave occurrence rate downstream than laminar shocks.

L. A. Davis↗

DOE OSTI2024

Sprain energy consequences for damage localization and fracture mechanics

The 2023 smooth Lagrangian Crack-Band Model (slCBM), inspired by the 2020 invention of the gap test, prevented spurious damage localization during fracture growth by introducing the second gradient of the displacement field vector, named the “sprain,” as the localization limiter. The key idea was that, in the finite element implementation, the displacement vector and its gradient should be treated as independent fields with the lowest ( C 0 ) continuity, constrained by a second-order Lagrange multiplier tensor. Coupled with a realistic constitutive law for triaxial softening damage, such as microplane model M7, the known limitations of the classical Crack Band Model were eliminated. Here, we show that the slCBM closely reproduces the size effect revealed by the gap test at various crack-parallel stresses. To describe it, we present an approximate corrective formula, although a strong loading-path dependence limits its applicability. Except for the rare case of zero crack-parallel stresses, the fracture predictions of the line crack models (linear elastic fracture mechanics, phase-field, extended finite element method (XFEM), cohesive crack models) can be as much as 100% in error. We argue that the localization limiter concept must be extended by including the resistance to material rotation gradients. We also show that, without this resistance, the existing strain-gradient damage theories may predict a wrong fracture pattern and have, for Mode II and III fractures, a load capacity error as much as 55%. Finally, we argue that the crack-parallel stress effect must occur in all materials, ranging from concrete to atomistically sharp cracks in crystals.

Science & Technology - Other Topics↗

DOE OSTI2026

Process-oriented evaluation of quasi-stationary Rossby waves and their impact on surface air temperature extremes in dynamical downscaling over North America

Quasi-stationary Rossby waves are a crucial component of the general circulation and play a significant role in regional water and energy cycles, as well as in extreme events. However, process-oriented evaluation for Rossby waves is rarely performed for dynamical downscaling simulations. To close this gap, we evaluate three classes of dynamical downscaling approaches, with a focus on quasi-stationary Rossby waves and their impact on surface air temperature over North America during Northern Hemisphere summer. The three classes of models differ in the way large-scale forcing is provided: a limited-area model (LAM) constrained only by lateral boundary conditions, represented by RegCM4 from the North American branch of the Coordinated Regional Downscaling Experiment (NA-CORDEX), a LAM with spectral nudging to maintain consistency in large-scale dynamics with the forcing data, represented by the Weather Research and Forecasting (WRF) model simulation in NA-CORDEX, and a global variable-resolution model with smoothly varying grid spacings, represented by the Community Atmosphere Model version 5.4, with the Model for Prediction Across Scales (MPAS) as its dynamical core (CAM-MPAS). With no constraints on the atmospheric dynamics, CAM-MPAS exhibits several mean biases in the upper-level circulations over the Pacific Coast region: a weaker subtropical jet, a northward-shifted mid-latitude jet, and an overestimated southerly flow. With the lateral boundary constraint alone, RegCM4 also exhibits weaker jets and overestimated southerly winds off the West Coast. Rossby ray theory reveals that those wind biases direct incoming Rossby waves northward. The erroneously routed Rossby waves distort the relationship between the accumulation of wave activity over the US West Coast and surface temperature anomalies over the Southern Great Plains, which emerges approximately 4 d after the convergence of wave-activity flux in the ERA-Interim reanalysis. Furthermore, the response of heatwaves to the extreme wave activity flux is not reproduced by the two models, a serious drawback as a dynamical downscaling framework is expected to connect large-scale forcing to local-scale phenomena. The WRF model employing spectral nudging is largely free from the aforementioned problems. A pair of sensitivity simulations suggests that spectral nudging is the key to improving the dynamics of quasi-stationary Rossby waves and their impact on surface air temperature. Our results also demonstrate the effectiveness of Rossby wave diagnostics that allow for realistic background flows for assessing the credibility of dynamical downscaling over North America, where incoming Rossby waves propagate through complex circulation patterns before traveling across the continent.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1994

The expanding photosphere method applied to SN 1992am AT cz = 14 600 km/s

We present photometry and spectroscopy of Supernova (SN) 1992am for five months following its discovery by the Calan Cerro-Tololo Inter-American Observatory (CTIO) SN search. These data show SN 1992am to be type II-P, displaying hydrogen in its spectrum and the typical shoulder in its light curve. The photometric data and the distance from our own analysis are used to construct the supernova's bolometric light curve. Using the bolometric light curve, we estimate SN 1992am ejected approximately 0.30 solar mass of Ni-56, an amount four times larger than that of other well studied SNe II. SN 1992am's; host galaxy lies at a redshift of cz = 14 600 km s(exp -1), making it one of the most distant SNe II discovered, and an important application of the Expanding Photsphere Method. Since z = 0.05 is large enough for redshift-dependent effects to matter, we develop the technique to derive luminosity distances with the Expanding Photosphere Method at any redshift, and apply this method to SN 1992am. The derived distance, D = 180(sub -25) (sup +30) Mpc, is independent of all other rungs in the extragalactic distance ladder. The redshift of SN 1992am's host galaxy is sufficiently large that uncertainties due to perturbations in the smooth Hubble flow should be smaller than 10%. The Hubble ratio derived from the distance and redshift of this single object is H(sub 0) = 81(sub -15) (sup +17) km s(exp -1) Mpc(exp -1). In the future, with more of these distant objects, we hope to establish an independent and statistically robust estimate of H(sub 0) based solely on type II supernovae.

Schmidt, Brian P.↗

NASA NTRS1953

Rocket-Model Investigation of the Longitudinal Stability, Drag, and Duct Performance Characteristics of the North American MX-770 (X-10) Missile at Mach Numbers from 0.80 to 1.70

A free-flight 0.12-scale rocket-boosted model of the North American MX-770 (X-10) missile has been tested in flight by the Pilotless Aircraft Research Division of the Langley Aeronautical Laboratory. Drag, longitudinal stability, and duct performance data were obtained at Mach numbers from 0.8 to 1.7 covering a Reynolds number range of about 9 x 10(exp 6) to 24 x 10(exp 6) based on wing mean aerodynamic chord. The lift-curve slope, static stability, and damping-in-pitch derivatives showed similar variations with Mach number, the parameters increasing from subsonic values in the transonic region and decreasing in the supersonic region. The variations were for the most part fairly smooth. The aerodynamic center of the configuration shifted rearward in the transonic region and moved forward gradually in the supersonic region. The pitching effectiveness of the canard control surfaces was maintained throughout the flight speed range, the supersonic values being somewhat greater than the subsonic. Trim values of angle of attack and lift coefficient changed abruptly in the transonic region, the change being associated with variations in the out-of-trim pitching moment, control effectiveness, and aerodynamic-center travel in this speed range. Duct total-pressure recovery decreased with increase in free-stream Mach number and the values were somewhat less than normal-shock recovery. Minimum drag data indicated a supersonic drag coefficient about twice the subsonic drag coefficient and a drag-rise Mach number of approximately 0.90. Base drag was small subsonically but was about 25 percent of the minimum drag of the configuration supersonically.

Bond, Aleck C.↗

NASA NTRS1999

Space Weathering in the Fine Size Fractions of Lunar Soils: Soil Maturity Effects

The effects of space weathering on the optical properties of lunar materials have been well documented. These effects include a reddened continuum slope, lowered albedo, and attenuated absorption features in reflectance spectra of lunar soils as compared to finely comminuted rocks from the same Apollo sites. However, the regolith processes that cause these effects are not well known, nor is the petrographic setting of the products of these processes fully understood. A Lunar Soil Characterization Consortium has been formed with the purpose of systematically integrating chemical and mineralogical data with the optical properties of lunar soils. Understanding space-weathering effects is critical in order to fully integrate the lunar sample collection with remotely-sensed data from recent robotic missions (e.g., Lunar Prospector, Clementine, and Galileo) We have shown that depositional processes (condensation of impact-derived vapors, sputter deposits, accreted impact material, e.g., splash glass, spherules, etc.) are a major factor in the modification of the optical surfaces of lunar regolith materials. In mature soils, it is the size and distribution of the nanophase metal in the soil grains that has the major effect on optical properties. In this report, we compare and contrast the space-weathering effects in an immature and a mature soil with similar elemental compositions. For this study, we analyzed <10 micron sieve fractions of two Apollo 17 soils, 79221 (mature, Is/FeO = 81) and 71061 (immature, Is/FeO = 14). Details of the sieving procedures and allocation scheme are given else where. The results of other detailed chemical, mineralogical, and spectroscopic analyses of these soil samples are reported elsewhere. A representative sample of each soil was embedded in low-viscosity epoxy, and thin sections (about 70nm thick) were obtained through ultra microtomy. The thin sections used for these analyses typically contained cross sections of up to 500 individual grains. The thin sections were studied using a JEOL 2010 transmission electron microscope (TEM) equipped with a thin window energy-dispersive X-ray (EDX) spectrometer. An individual thin section was selected from each soil, and for each grain in the section we determined (1) the elemental composition by EDX; (2) whether the grain was crystalline or glassy using electron diffraction and darkfield imaging; (3) the presence or absence of rims and accreted material; and (4) the distribution of nanophase Fe where present. Most of the categories are self-evident; however, we divide the agglutinate derived material into agglutinitic glass (glass with approximately the same composition as the bulk soil that contains nanophase Fe with or without vesicles) and agglutinate fragments, which are composed of crystalline grains and agglutinitic glass. Lithic fragments are defined as polymineralic grains with no glass. Pyroxene grains have been divided into high- and low-Ca groups. As expected, there are a number of differences in the petrography of the <10-microns fractions of 79221 and 71061 given the great difference in their respective maturities, but we focus here on two major distinctions: agglutinate content and the number of grains with micropatina. Slightly over 50% of the particles in 79221 consist of agglutinitic glass and agglutinate fragments, while the remainder are predominantly crystalline mineral grains. The agglutinic glass particles contain abundant nanophase Fe and vesicles. Angular particles are rare, with most showing smooth, rounded exteriors, Of the mineral grains analyzed thus far, over 90% of the grains have amorphous rims that contain nanophase Fe (these rims are believed to have formed by vapor deposition and irradiation effects). The nanophase Fe in these rims probably accounts for a significant fraction of the increase in Is/FeO measured in these size fractions. In addition to the rims, the majority of particles also show abundant accreted material in the form of glass splashes and spherules that also contain nanophase Fe. In stark contrast, the surfaces of the mineral grains in the 71061 sample are relatively prisitine, as only about 14% of the mineral grains in the sample exhibited amorphous rims. Furthermore, the mineral particles are more angular and show greater surface roughness than in the mature sample. Accreted material on particle surfaces is rare. Agglutinitic material is a major component of the 71061 sample; however, nanophase Fe and vesicles are not as well developed as in the 79221 sample. It is now recognized that nanophase Fe is probably the main agent in modifying the optical properties of lunar soil grains. The most important result of this study is the observation that in the fine size fractions of mature soils, nearly every grain has nanophase Fe within 100 run of the particle surface. (Additional Information contained in original)

Keller, L. P.↗

NASA NTRS1994

Investigation of a Laminar Flow Leading Edge

The recent resurgence of interest in utilizing laminar flow on aircraft surfaces for reduction in skin friction drag has generated a considerable amount of research in natural laminar flow (NLF) and hybrid laminar flow control (HLFC) on transonic aircraft wings. This research has focused primarily on airfoil design and understanding transition behavior with little concern for the surface imperfections and manufacturing variations inherent to most production aircraft. In order for laminar flow to find wide-spread use on production aircraft, techniques for constructing the wings must be found such that the large surface imperfections present in the leading edge region of current aircraft do not occur. Toward this end, a modification to existing leading edge construction techniques was devised such that the resulting surface did not contain large gaps and steps as are common on current production aircraft of this class. A lowspeed experiment was first conducted on a simulation of the surface that would result from this construction technique. Preston tube measurements of the boundary layer downstream of the simulated joint and flow visualization using sublimation chemicals validated the literature on the effects of steps on a laminar boundary layer. These results also indicated that the construction technique was indeed compatible with laminar flow. In order to fully validate the compatibility of this construction technique with laminar flow, thus proving that it is possible to build wings that are smooth enough to be used on business jets and light transports in a manner compatible with laminar flow, a flight experiment is being conducted. In this experiment Mach number and Reynolds number will be matched in a real flight environment. The experiment is being conducted using the NASA Dryden F-104 Flight Test Fixture (FTF). The FTF is a low aspect ratio ventral fin mounted beneath an F-104G research aircraft. A new nose shape was designed and constructed for this experiment. This nose shape provides an accelerating pressure gradient in the leading edge region. By flying the aircraft at appropriate Mach numbers and altitudes, this nose shape simulates the leading edge region of a laminar flow wing for a business jet or light transport. Manufactured into the nose shape is a spanwise slot located approximately four inches downstream of the leading edge. The slot, which is an inch wide and one-eighth of an inch deep allows the simulation of surface imperfections, such as gaps and steps at skin joints, which will occur on aircraft using this new construction technique. By placing strips of aluminum of various sizes and shapes in the slot, the effect on the boundary layer of different sizes and shapes of steps and gaps will be examined. It is planned to use five different configurations, differing primarily in the size and number of gaps. Downstream of the slot, the state of the boundary layer is determined using hot film gages and Stanton gages. Agreement between these two very different techniques of measuring boundary layer properties is considered important to being able to state with confidence the effects on the boundary layer of the simulated manufacturing imperfections. To date, the aircraft has not flown. First flights of the aircraft are on schedule to begin October 4, 1993. Low-speed, preliminary experiments at matching Reynolds numbers have been completed.

Drake, Aaron↗

NASA NTRS2012

Modifications Of Discrete Ordinate Method For Computations With High Scattering Anisotropy: Comparative Analysis

A numerical accuracy analysis of the radiative transfer equation (RTE) solution based on separation of the diffuse light field into anisotropic and smooth parts is presented. The analysis uses three different algorithms based on the discrete ordinate method (DOM). Two methods, DOMAS and DOM2+, that do not use the truncation of the phase function, are compared against the TMS-method. DOMAS and DOM2+ use the Small-Angle Modification of RTE and the single scattering term, respectively, as an anisotropic part. The TMS method uses Delta-M method for truncation of the phase function along with the single scattering correction. For reference, a standard discrete ordinate method, DOM, is also included in analysis. The obtained results for cases with high scattering anisotropy show that at low number of streams (16, 32) only DOMAS provides an accurate solution in the aureole area. Outside of the aureole, the convergence and accuracy of DOMAS, and TMS is found to be approximately similar: DOMAS was found more accurate in cases with coarse aerosol and liquid water cloud models, except low optical depth, while the TMS showed better results in case of ice cloud.

Korkin, Sergey V.↗

NASA NTRS2011

A Sensitivity Study on the Effects of Particle Chemistry, Asphericity and Size on the Mass Extinction Efficiency of Mineral Dust in the Earth's Atmosphere: From the Near to Thermal IR

To determine a plausible range of mass extinction efficiencies (MEE) of terrestrial atmospheric dust from the near to thermal IR, sensitivity analyses are performed over an extended range of dust microphysical and chemistry perturbations. The IR values are subsequently compared to those in the near-IR, to evaluate spectral relationships in their optical properties. Synthesized size distributions consistent with measurements, model particle size, while composition is defined by the refractive indices of minerals routinely observed in dust, including the widely used OPAC/Hess parameterization. Single-scattering properties of representative dust particle shapes are calculated using the T-matrix, Discrete Dipole Approximation and Lorenz-Mie light-scattering codes. For the parameterizations examined, MEE ranges from nearly zero to 1.2 square meters per gram, with the higher values associated with non-spheres composed of quartz and gypsum. At near-IR wavelengths, MEE for non-spheres generally exceeds those for spheres, while in the thermal IR, shape-induced changes in MEE strongly depend on volume median diameter (VMD) and wavelength, particularly for MEE evaluated at the mineral resonant frequencies. MEE spectral distributions appear to follow particle geometry and are evidence for shape dependency in the optical properties. It is also shown that non-spheres best reproduce the positions of prominent absorption peaks found in silicates. Generally, angular particles exhibit wider and more symmetric MEE spectral distribution patterns from 8-10 micrometers than those with smooth surfaces, likely due to their edge-effects. Lastly, MEE ratios allow for inferring dust optical properties across the visible-IR spectrum. We conclude the MEE of dust aerosol are significant for the parameter space investigated, and are a key component for remote sensing applications and the study of direct aerosol radiative effects.

Hansell, R. A., Jr.↗

NASA NTRS2016

VIIRS Reflective Solar Band Radiometric and Stability Evaluation Using Deep Convective Clouds

This work takes advantage of the stable distribution of deep convective cloud (DCC) reflectance measurements to assess the calibration stability and detector difference in Visible Infrared Imaging Radiometer Suite (VIIRS) reflective bands. VIIRS Sensor Data Records (SDRs) from February 2012 to June 2015 are utilized to analyze the long-term trending, detector difference, and half angle mirror (HAM) side difference. VIIRS has two thermal emissive bands with coverage crossing 11 microns for DCC pixel identification. The comparison of the results of these two processing bands is one of the indicators of analysis reliability. The long-term stability analysis shows downward trends (up to approximately 0.4 per year) for the visible and near-infrared bands and upward trends (up to 0.5per year) for the short- and mid-wave infrared bands. The detector difference for each band is calculated as the difference relative to the average reflectance overall detectors. Except for the slightly greater than 1 difference in the two bands at 1610 nm, the detector difference is less than1 for other solar reflective bands. The detector differences show increasing trends for some short-wave bands with center wavelengths from 400 to 600 nm and remain unchanged for the bands with longer center wavelengths. The HAM side difference is insignificant and stable. Those short-wave bands from 400 to 600 nm also have relatively larger HAM side difference, up to 0.25.Comparing the striped images from SDR and the smooth images after the correction validates the analyses of detector difference and HAM side difference. These analyses are very helpful for VIIRS calibration improvement and thus enhance product quality

VIIRS↗

NASA NTRS1993

Numerical simulations of galactic wakes

We have performed a series of numerical simulations which examine the simple but realistic case of a galaxy with a smooth potential corresponding to a King profile, and with no interstellar medium. This corresponds to the case of an early-type cluster galaxy that has previously been stripped of its gas. The simulations provide a numerical basis for future work involving galaxies with interstellar media. We use a time-explicit finite difference code (Lufkin & Hawley 1992, ApJ submitted) to obtain solutions to the nonlinear equations for adiabatic flow in axisymmetry. The computational grid covers a region extending 500 kpc from a stationary galaxy radially and along the symmetry axis. We use a graded mesh, allowing for full resolution near the galaxy. Because the assumed potential has a finite depth, no artificial inner boundary condition is necessary; reflection symmetry at the axis is assumed. With a total grid size of 64 radial x 128 vertical, each run requires approximately 5 minutes of cpu time on the NCSA Cray-2.

Lufkin, Eric A.↗

DOE OSTI2025

High-order limiting methods using maximum principle bounds derived from the Boltzmann equation I: Euler equations

The use of limiting methods for high-order numerical approximations of hyperbolic conservation laws generally requires defining an admissible region/bounds for the solution. In this work, we present a novel approach for computing solution bounds and limiting for the Euler equations through the kinetic representation provided by the Boltzmann equation, which allows for extending limiters designed for linear advection directly to the Euler equations. Given an arbitrary set of solution values to compute bounds over (e.g., numerical stencil) and a desired linear advection limiter, the proposed approach yields an analytic expression for the admissible region of particle distribution function values, which may be numerically integrated to yield a set of bounds for the density, momentum, and total energy. Further, these solution bounds are shown to preserve positivity of density/pressure/internal energy and, when paired with a limiting technique, can robustly resolve strong discontinuities while recovering high-order accuracy in smooth regions without any ad hoc corrections (e.g., relaxing the bounds). This approach is demonstrated in the context of an explicit unstructured high-order discontinuous Galerkin/flux reconstruction scheme for a variety of difficult problems in gas dynamics, including cases with extreme shocks and shock-vortex interactions. Furthermore, this work presents a foundation for limiting techniques for more complex macroscopic governing equations that can be derived from an underlying kinetic representation for which admissible solution bounds are not well-understood.

42 ENGINEERING↗