Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “spatial discretization”

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 613 records · Page 34
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Verification of TROPOMI NO2 Product Using OMI NO2 algorithm

We evaluated S5P TROPOspheric Monitoring Instrument (TROPOMI) operational nitrogen dioxide (NO2) product by comparing with the NO2 retrievals from Ozone Monitoring Instrument (OMI) onboard NASA’s Aura satellite. We compared spatially matched NO2 vertical column density (VCD) data from OMI and TROPOMI to identify any discrepancies between the two operational products. We also applied OMI Air-Mass Factor (AMF) algorithm to TROPOMI Slant Column Density (SCD) data (v02.03.01) aiming to create long-term merged NO2 Earth Science Data Record (https://disc.gsfc.nasa.gov/datasets/TROPOMI_MINDS_NO2_1.1/summary ). The algorithm (version 4.0) uses new Geometry-dependent Lambertian Equivalent surface Reflectivity (GLER) product available for each TROPOMI pixel. GLER is pre-calculated using the vector linearized discrete ordinate radiative transfer (VLIDORT) model, which uses as input high-resolution bidirectional reflectance distribution function (BRDF) information from NASA's Aqua Moderate Resolution Imaging Spectroradiometer (MODIS) over land and the wind-dependent Cox–Munk wave-facet slope distribution over water, the latter with a contribution from the water-leaving radiance based on MODIS gap-filled in-water chlorophyll-a data. The GLER data, combined with consistently retrieved cloud parameters, provide improved information for the calculation of the scattering weight profiles. Additional AMF improvements result from using high resolution (0.25o latitude x 0.25o longitude) a priori NO2 profile shapes and other auxiliary information from the Global Modeling Initiative (GMI) Replay simulation sampled using S5P orbital simulator. Our TROPOMI_MINDS_NO2 algorithm employs the stratosphere-troposphere separation scheme, de-striping, and surface snow/ice treatment consistently with the latest OMI NO2 re-processing (version 4). We evaluate the TROPOMI_MINDS_NO2 product by comparing with the S5P operational NO2 product as well as with independent NO2 observations from ground-based Pandora and aircraft observations. Our results indicate that the new TROPOMI_MINDS_NO2 retrievals are generally higher than the S5P operational NO2 product over polluted regions and show improved agreement with independent validation data.

TROPOMI↗
NASA NTRS2008

Direct-View Multi-Point Two-Component Interferometric Rayleigh Scattering Velocimeter

This paper describes an instantaneous velocity measurement system based on the Doppler shift of elastically scattered laser light from gas molecules (Rayleigh scattering) relative to an incident laser. The system uses a pulsed laser as the light source, direct-viewing optics to collect the scattered light, an interferometer to analyze spectrally the scattered light mixed with the incident laser light, and a CCD camera to capture the resulting interferogram. The system is capable of simultaneous, spatially (approximately 0.2 mm(exp 3)) and temporally (approximately 40 ns) resolved, multiple point measurements of two orthogonal components of flow velocity in the presence of background scattered light, acoustic noise and vibrations, and flow particulates. Measurements in a large-scale axi-symmetric Mach 1.6 H2-air combustion-heated jet running at a flow sensible enthalpy specific to Mach 5.5 hypersonic flight are performed to demonstrate the technique. The measurements are compared with CFD calculations using a finite-volume discretization of the Favre-averaged Navier-Stokes equations (VULCAN code).

Bivolaru, Daniel↗
NASA NTRS2024

Physics of the Seasonal Sea Ice Zone

The seasonal sea ice zone encompasses the region between the winter maximum and summer minimum sea ice extent. In both the Arctic and Antarctic, the majority of the ice cover can now be classified as seasonal. Here, we review the sea ice physics that governs the evolution of seasonal sea ice in the Arctic and Antarctic, spanning sea ice growth, melt, and dynamics and including interactions with ocean surface waves as well as other coupled processes. The advent of coupled wave–ice modeling and discrete-element modeling, together with improved and expanded satellite observations and field campaigns, has yielded advances in process understanding. Many topics remain in need of further investigation, including rheologies appropriate for seasonal sea ice, wave-induced sea ice fracture, welding for sea ice freeze-up, and the distribution of snow on seasonal sea ice. Future research should aim to redress biases (such as disparities in focus between the Arctic and Antarctic and between summer and winter processes) and connect observations to modeling across spatial scales.

sea ice↗
NASA NTRS2021

Prompt Response of the Dayside Magnetosphere to Discrete Structures within the Sheath Region of a Coronal Mass Ejection

A sequence of discrete solar wind structures within the sheath region of an interplanetary coronal mass ejection (ICME) on Nov 06, 2015 caused a series of compressions and releases of the dayside magnetosphere. Each compression resulted in a brief adiabatic enhancement of ions (electrons) driving bursts of electromagnetic ion cyclotron (whistler mode chorus) wave growth across the dayside magnetosphere. Fine-structured rising tones were observed in the EMIC wave bursts, resulting in nonlinear scattering of relativistic electrons in the outer radiation belt. Multipoint observations allow us to study the spatial structure and evolution of these sheath structures as they propagate Earthward from L1 as well as the spatio-temporal characteristics of the magnetospheric response. This event highlights the importance of fine-scale solar wind structure, in particular within complex sheath regions, in driving dayside phenomena within the inner magnetosphere.

L. W. Blum↗
NASA NTRS1981

A Quasi-Parabolic technique for computation of three-dimensional viscous flows

A computational technique is presented for obtaining flowfield solutions to a parabolic form of the Navier-Stokes equations. The point of departure is the General Interpolant Method (GIM) which provides a discretization for partial differential equations on arbitrary three-dimensional geometries. The new scheme, termed Quasi-Parabolic, treats the parabolized equations but with 'time-like' terms appended. Addition of these extra terms, which are relaxed by iteration, avoid many of the singularities inherent in classical parabolic Navier-Stokes methods. Streamwise derivatives are approximated by three-point backward differences and the cross plane operators use an alternating forward-backward sweep. A two-step sequence is used to implement the difference scheme in the spatial dimensions and a time-like relaxation converges the quasi-parabolic procedure at each plane. Solutions are presented for flows in two and three dimensions. Inviscid flows are solved for internal and external applications and viscous flows in boundary layers and free shear layers are also computed with the GIM/Quasi-Parabolic scheme.

Spradley, L. W.↗
NASA NTRS1994

Features of steady magnetospheric convection

The large-scale patterns of ionospheric convection and particle precipitation are described during two intervals of steady magnetospheric convection (SMC) on November 24, 1981. The unique data set used in the analysis includes recordings from the worldwide network of magnetometers and all-sky cameras, global auroral images from the Dynamics Explorer (DE) 1 spacecraft, and particle precipitation data from low-altitude National Oceanic and Atmospheric Administration (NOAA) 6 and NOAA 7 spacecraft. The data show that intense magnetospheric convection continued during more than 10 hours under the steady southward interplanetary magnetic field without any distinct substorm signatures. All data sets available confirmed the stable character of the large-scale magnetospheric configuration during this period. In particular, the magnetic flux threading the polar cap was stable (within 10%) during 3.5 hours of continued DE 1 observations. The dayside cusp was located at an unusually low latitude (70 deg CGL). The nightside auroral pattern consisted of two distinct regions. The diffuse aurora in the equatorward half of the expanded (10 deg wide) auroral oval was well-separated from the bright, active auroral forms found in the vicinity of the poleward boundary of the oval. The twin-vortex convection pattern had no signature of the Harang discontinuity; its nightside 'convection throat' was spatially coincident with the poleward active auroras. This region of the auroral oval was identified as the primary site of the short-lived transient activations during the SMC intervals. The energetic particle observations show that the auroral precipitation up to its high-latitude limit is on closed field lines and that particle acceleration up to greater than 30-keV energy starts close to this limit. The isotropic boundaries of the greater than 30-keV protons and electrons were found close to each other, separating regions of discrete and diffuse precipitation. This suggests that these precipitation types originate on the very taillike and very dipolelike field lines, respectively.

Yahnin, A.↗
NASA NTRS1995

A deterministic particle method for one-dimensional reaction-diffusion equations

We derive a deterministic particle method for the solution of nonlinear reaction-diffusion equations in one spatial dimension. This deterministic method is an analog of a Monte Carlo method for the solution of these problems that has been previously investigated by the author. The deterministic method leads to the consideration of a system of ordinary differential equations for the positions of suitably defined particles. We then consider the time explicit and implicit methods for this system of ordinary differential equations and we study a Picard and Newton iteration for the solution of the implicit system. Next we solve numerically this system and study the discretization error both analytically and numerically. Numerical computation shows that this deterministic method is automatically adaptive to large gradients in the solution.

Mascagni, Michael↗
NASA NTRS2002

High-Efficiency High-Resolution Global Model Developments at the NASA Goddard Data Assimilation Office

The Data Assimilation Office (DAO) has been developing a new generation of ultra-high resolution General Circulation Model (GCM) that is suitable for 4-D data assimilation, numerical weather predictions, and climate simulations. These three applications have conflicting requirements. For 4-D data assimilation and weather predictions, it is highly desirable to run the model at the highest possible spatial resolution (e.g., 55 km or finer) so as to be able to resolve and predict socially and economically important weather phenomena such as tropical cyclones, hurricanes, and severe winter storms. For climate change applications, the model simulations need to be carried out for decades, if not centuries. To reduce uncertainty in climate change assessments, the next generation model would also need to be run at a fine enough spatial resolution that can at least marginally simulate the effects of intense tropical cyclones. Scientific problems (e.g., parameterization of subgrid scale moist processes) aside, all three areas of application require the model's computational performance to be dramatically improved as compared to the previous generation. In this talk, I will present the current and future developments of the "finite-volume dynamical core" at the Data Assimilation Office. This dynamical core applies modem monotonicity preserving algorithms and is genuinely conservative by construction, not by an ad hoc fixer. The "discretization" of the conservation laws is purely local, which is clearly advantageous for resolving sharp gradient flow features. In addition, the local nature of the finite-volume discretization also has a significant advantage on distributed memory parallel computers. Together with a unique vertically Lagrangian control volume discretization that essentially reduces the dimension of the computational problem from three to two, the finite-volume dynamical core is very efficient, particularly at high resolutions. I will also present the computational design of the dynamical core using a hybrid distributed-shared memory programming paradigm that is portable to virtually any of today's high-end parallel super-computing clusters.

Lin, Shian-Jiann↗
NASA NTRS2002

High-Efficiency High-Resolution Global Model Developments at the NASA Goddard Data Assimilation Office

The Data Assimilation Office (DAO) has been developing a new generation of ultra-high resolution General Circulation Model (GCM) that is suitable for 4-D data assimilation, numerical weather predictions, and climate simulations. These three applications have conflicting requirements. For 4-D data assimilation and weather predictions, it is highly desirable to run the model at the highest possible spatial resolution (e.g., 55 kin or finer) so as to be able to resolve and predict socially and economically important weather phenomena such as tropical cyclones, hurricanes, and severe winter storms. For climate change applications, the model simulations need to be carried out for decades, if not centuries. To reduce uncertainty in climate change assessments, the next generation model would also need to be run at a fine enough spatial resolution that can at least marginally simulate the effects of intense tropical cyclones. Scientific problems (e.g., parameterization of subgrid scale moist processes) aside, all three areas of application require the model's computational performance to be dramatically improved as compared to the previous generation. In this talk, I will present the current and future developments of the "finite-volume dynamical core" at the Data Assimilation Office. This dynamical core applies modem monotonicity preserving algorithms and is genuinely conservative by construction, not by an ad hoc fixer. The "discretization" of the conservation laws is purely local, which is clearly advantageous for resolving sharp gradient flow features. In addition, the local nature of the finite-volume discretization also has a significant advantage on distributed memory parallel computers. Together with a unique vertically Lagrangian control volume discretization that essentially reduces the dimension of the computational problem from three to two, the finite-volume dynamical core is very efficient, particularly at high resolutions. I will also present the computational design of the dynamical core using a hybrid distributed- shared memory programming paradigm that is portable to virtually any of today's high-end parallel super-computing clusters.

Lin, Shian-Jiann↗
NASA NTRS1997

Bio-Optical Properties of the Arabian Sea as Determined by In Situ and Sea WiFS Data

The overall objective of this work was to characterize optical and fluorescence properties in the euphotic zone during two British Ocean Flux Study (BOFS) Arabian Sea cruises. This was later expanded in 1995 to include three U.S. JGOFS Arabian Sea Cruises. The region was to be divided into one or more "bio-optical provinces," within each of which a single set of regression models was to be developed to relate the vertical distribution of irradiance attenuation and normalized fluorescence (SF and NF) to remote sensing reflectance and diffuse attenuation coefficient. The working hypothesis was that over relatively large spatial and temporal scales, the vertical profiles of bio-optical properties were predictable. The specific technical objectives were: (1) To characterize the vertical distribution of the inherent and apparent optical properties by measuring downwelling and upwelling irradiances, upwelling radiances, scalar irradiance of PAR, and beam transmissions at each station - from these data, spectral diffuse attenuation coefficients, irradiance reflectances, remote sensing reflectances, surface-leaving radiances and beam attenuation coefficients were determined; (2) To characterize the spectral absorption of total particulate, detrital, and dissolved organic material at each station from discrete water samples; (3) To describe the vertical distribution of photoadaptive properties in the water column by measuring profiles of stimulated (SF) and natural (NF) fluorescence and examining relationships between SF and NF as a function of diffuse optical depth, pigment biomass and primary productivity; and (4) To establish locally derived, in-water algorithms relating remote sensing reflectance spectra to diffuse attenuation coefficients, phytoplankton pigment concentrations and primary productivity, through intercomparisons with in situ measurements, for application to SeaWiFS data.

Trees, Charles C.↗
NASA NTRS1975

Relation of lineaments to sulfide deposits: Bald Eagle Mountain, Centre County, Pennsylvania

The author has identified the following significant results. Discrete areas of finely-fractured and brecciated sandstone float are present along the crest of Bald Mountain and are commonly sites of sulfide mineralization, as evidenced by the presence of barite and limonite gossans. The frequency distributions of the brecciated float as the negative binomial distribution supports the interpretation of a separate population of intensely fractured material. Such zones of concentrated breccia float have an average width of one kilometer with a range from 0.4 to 1.6 kilometers and were observed in a quarry face to have subvertical dips. Direct spatial correlation of the Landsat-derived lineaments to the fractured areas on the ridge is low; however, the mineralized and fracture zones are commonly assymetrical to the lineament positions. Such a systematic dislocation might result from an inherent bias in the float population or could be the product of the relative erosional resistance of the silicified material in the mineralized areas in relation to the erosionally weak material at the stream gaps.

Mcmurtry, G. J.↗
NASA NTRS2001

Studying Earth's Environment From Space: Classroom and Laboratory Activities with Instructor Resources

Standard, text-book based learning for earth, ocean, and atmospheric sciences has been limited by the unavailability of quantitative teaching materials. While a descriptive presentation, in a lecture format, of discrete satellite images is often adequate for high school classrooms, this is seldom the case at the undergraduate level. In order to address these concerns, a series of numerical exercises for the Macintosh was developed for use with satellite-derived Sea Surface Temperature, pigment and sea ice concentration data. Using a modified version of NIH Image, to analyze actual satellite data, students are able to better understand ocean processes, such as circulation, upwelling, primary production, and ocean/atmosphere coupling. Graphical plots, image math, and numerical comparisons are utilized to substantiate temporal and spatial trends in sea surface temperature and ocean color. Particularly for institutions that do not offer a program in remote sensing, the subject matter is presented as modular units, each of which can be readily incorporated into existing curricula. These materials have been produced in both CD-ROM and WWW format, making them useful for classroom or lab setting. Depending upon the level of available computer support, graphics can be displayed directly from the CD-ROM, or as a series of color view graphs for standard overhead projection.

Smith, Elizabeth A.↗
NASA NTRS2003

Charge Loss and Charge Sharing Measurements for Two Different Pixelated Cadmium-Zinc-Telluride Detectors

As part of ongoing research at Marshall Space Flight Center, Cadmium-Zinc- Telluride (CdZnTe) pixilated detectors are being developed for use at the focal plane of the High Energy Replicated Optics (HERO) telescope. HERO requires a 64x64 pixel array with a spatial resolution of around 200 microns (with a 6m focal length) and high energy resolution (< 2% at 60keV). We are currently testing smaller arrays as a necessary first step towards this goal. In this presentation, we compare charge sharing and charge loss measurements between two devices that differ both electronically and geometrically. The first device consists of a 1-mm-thick piece of CdZnTe that is sputtered with a 4x4 array of pixels with pixel pitch of 750 microns (inter-pixel gap is 100 microns). The signal is read out using discrete ultra-low-noise preamplifiers, one for each of the 16 pixels. The second detector consists of a 2-mm-thick piece of CdZnTe that is sputtered with a 16x16 array of pixels with a pixel pitch of 300 microns (inter-pixel gap is 50 microns). Instead of using discrete preamplifiers, the crystal is bonded to an ASIC that provides all of the front-end electronics to each of the 256 pixels. what degree the bias voltage (i.e. the electric field) and hence the drift and diffusion coefficients affect our measurements. Further, we compare the measured results with simulated results and discuss to

Gaskin, Jessica↗
NASA NTRS1999

Data Structure and Parallel Decomposition Considerations on a Fibonacci Grid

The Fibonacci grid, proposed by Swinbank and Purser (see companion abstract), provides attractive properties for global numerical atmospheric prediction by offering an optimally homogeneous, geometrically regular, and approximately isotropic discretization, with only the polar regions requiring special numerical treatment. It is a mathematical idealization, applied to the sphere, of the multi-spiral patterns often found in botanical structures, such as in pine cones and sunflower heads. Computationally, it is natural to organize the domain, into zones, in each of which the same pair, or triple, of "Fibonacci spirals" dominate. But the further subdivision of such zones into "tiles" of a shape and size suitable for distribution to the processors of a massively parallel computer requires very careful consideration if the subsequent spatial computations along the respective spirals, especially those computations (such as compact differencing schemes) that involve recursion, can be implemented in an efficient "load-balanced "manner without requiring excessive amounts of inter-processor communications. In this paper we show how certain "number theoretic" properties of the Fibonacci sequence (whose numbers prescribe the multiplicity of successive spirals) may be exploited in the decomposition of grid zones into tidy arrangements of triangular grid tiles, each tile possessing one side approximately parallel to the constant-latitude zone boundary. We also describe how the spatially recursive processes may be decomposed across such a tiling, and the directionality of the recursions reversed on alternate grid lines, to ensure a very high degree of load balancing throughout the execution of the computations required for one time step of a global model.

Michalakes, John↗
NASA NTRSDate not supplied

Computational Prediction of Broadband Noise from a Representative Unmanned Aerial System Rotor

This work details the application of a lattice-Boltzmann method–very-large-eddy simulation (LBM-VLES) employed by the software suite, PowerFLOW. This LBM-VLES simulation predicted the aeroacoustic noise emanating from a representative, small unmanned aircraft system rotor, namely, the DJI-9450 in a hover condition. Predicted total aerodynamic loading as well as 2D aerodynamic loading along discrete spanwise sections of a rotor blade were compared to lower fidelity predictions and experimental results acquired in the Structural Acoustic Loads and Transmission anechoic chamber facility at the NASA Langley Research Center. The total acoustic spectra were decomposed into tonal and broadband components, which showed that broadband noise was a dominant contributor above 1 kHz for this rotor. These data were then compared to experimentally acquired data, showing good agreement up to approximately 11 kHz. Above 11 kHz, however, a grid sensitivity study showed dependency of the highest resolvable frequency on the spatial resolution of the computational domain, explaining the roll off in predicted data. Individual broadband noise sources were further investigated by calculating one-third octave sound pressure levels of the unsteady pressure fluctuations acting on the rotor, providing evidence that blade self-noise was the prominent noise source. Using these results, blade wake interaction noise was seen to be negligible for this particular rotor, which was further validated by calculating blade vortex miss distances and comparing to theory.

Christopher S. Thurman↗
NASA NTRS1978

The correlation function for density perturbations in an expanding universe. III The three-point and predictions of the four-point and higher order correlation functions

Higher-order correlation functions for the large-scale distribution of galaxies in space are investigated. It is demonstrated that the three-point correlation function observed by Peebles and Groth (1975) is not consistent with a distribution of perturbations that at present are randomly distributed in space. The two-point correlation function is shown to be independent of how the perturbations are distributed spatially, and a model of clustered perturbations is developed which incorporates a nonuniform perturbation distribution and which explains the three-point correlation function. A model with hierarchical perturbations incorporating the same nonuniform distribution is also constructed; it is found that this model also explains the three-point correlation function, but predicts different results for the four-point and higher-order correlation functions than does the model with clustered perturbations. It is suggested that the model of hierarchical perturbations might be explained by the single assumption of having density fluctuations or discrete objects all of the same mass randomly placed at some initial epoch.

Mcclelland, J.↗
NASA NTRS1989

Inverse dynamics of a 3 degree of freedom spatial flexible manipulator

A technique is presented for solving the inverse dynamics and kinematics of 3 degree of freedom spatial flexible manipulator. The proposed method finds the joint torques necessary to produce a specified end effector motion. Since the inverse dynamic problem in elastic manipulators is closely coupled to the inverse kinematic problem, the solution of the first also renders the displacements and rotations at any point of the manipulator, including the joints. Furthermore the formulation is complete in the sense that it includes all the nonlinear terms due to the large rotation of the links. The Timoshenko beam theory is used to model the elastic characteristics, and the resulting equations of motion are discretized using the finite element method. An iterative solution scheme is proposed that relies on local linearization of the problem. The solution of each linearization is carried out in the frequency domain. The performance and capabilities of this technique are tested through simulation analysis. Results show the potential use of this method for the smooth motion control of space telerobots.

Bayo, Eduardo↗
NASA NTRS1990

Low-frequency ionization-driven instability of a discrete auroral arc model

The low-frequency (time scales of tens of seconds) dynamics of the auroral ionosphere differs from that of the nonauroral ionosphere by the presence of strong, unstable space- and time-dependent ionospheric ionization produced by precipitating auroral electrons. If recombination is relatively unimportant (as at high ionospheric heights), it is shown that, in general, transport processes cannot remove this ionization as fast as it is created, and no equilibrium is possible. These nonequilibrium phenomena are investigated in the context of a nonlinear adiabatic auroral model, which has previously been studied in static situations. A linearized local perturbation analysis is given of what amounts to a current-driven E x B gradient-drift instability with an ionization source, as well as some exact nonlinear solutions valid in a finite but limited spatial region. These solutions show continuing motion of auroral potential and plasma density, as the aurora tries to shift its ionization problems from one place to another. The analysis gives clues to the possibility of generation of chaos and of fine-scale spatial structure.

Cornwall, John M.↗