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 649 records · Page 36

On the nature of models in remote sensing

An explicit framework can provide a better understanding of remote sensing models and their interrelationships. This framework distinguishes between the scene, which is real and exists on the ground, and the image, which is a collection of spatially arranged masurements drawn from the scene. The scene model generalizes and parameterizes the essential qualities of the scene. Scene models may be discrete, in which the scene model consists of discrete elements with boundaries, or continuous, in which matter and energy flows are taken to be continuous and there are no clear or sharp boundaries in the scene. In the discrete case, there are two possibilities for models: H- and L-resolution. In the H-resolution case, the resolution cells of the image are smaller than the elements, and thus the elements may be individually resolved. In the L-resolution case, the resolution cells are larger than the elements and cannot be resolved. Most canopy models are L-resolution, deterministic, and noninvertible in nature; image processing models, however, tend to be H-resolution, empirical, and invertible. This taxonomy helps add insight to the development of remote sensing theory and point the way to new, productive areas of research.

Strahler, A. H.↗

Numerical analysis of three-dimensional viscous internal flows

A 3-D Navier-Stokes code has been developed for analysis of turbomachinery blade rows and other internal flows. The Navier-Stokes equations are written in a Cartesian coordinate system rotating about the x-axis, and then mapped to a general body-fitted coordinate system. Streamwise viscous terms are neglected using the thin-layer assumption, and turbulence effects are modeled using the Baldwin-Lomax turbulence model. The equations are discretized using finite differences on stacked C-type grids and are solved using a multistage Runge-Kutta algorithm with a spatially-varying time step and implicit residual smoothing. Calculations have been made of a horseshoe vortex formed in front of a flat plate with a round leading edge standing in a turbulent endwall boundary layer. Comparisons are made with experimental data taken by Eckerle and Langston for a circular cylinder under similar conditions. Computer and measured results are compared in terms of endwall flow visualization pictures and total pressure loss contours and vector plots on the symmetry plane. Calculated details of the primary vortex show excellent agreement with the experimental data. The calculations also show a small secondary vortex not seen experimentally.

Chima, Rodrick V.↗

Numerical simulation of three-dimensional unsteady vortex flow using a compact vorticity-velocity algorithm

A numerical algorithm is presented which is used to solve the unsteady, fully three-dimensional, incompressible Navier-Stokes equations in vorticity-velocity variables. A discussion of the discrete approximation scheme is presented as well as the solution method used to solve the resulting algebraic set of difference equations. Second order spatial and temporal accuracy is verified through solution comparisons with exact results obtained for steady three-dimensional stagnation point flow and unsteady axisymmetric vortex spin-up. In addition, results are presented for the problem of unsteady bubble-type vortex breakdown with emphasis on internal bubble dynamics and structure.

Gatski, T. B.↗

Numerical analysis of three-dimensional viscous internal flows

A 3-D Navier-Stokes code has been developed for analysis of turbomachinery blade rows and other internal flows. The Navier-Stokes equations are written in a Cartesian coordinate system rotating about the x-axis, and then mapped to a general body-fitted coordinate system. Streamwise viscous terms are neglected using the thin-layer assumption, and turbulence effects are modeled using the Baldwin-Lomax turbulence model. The equations are discretized using finite differences on stacked C-type grids and are solved using a multistage Runge-Kutta algorithm with a spatially-varying time step and implicit residual smoothing. Calculations have been made of a horseshoe vortex formed in front of a flat plate with a round leading edge standing in a turbulent endwall boundary layer. Comparisons are made with experimental data taken by Eckerle and Langston for a circular cylinder under similar conditions. Computer and measured results are compared in terms of endwall flow visualization pictures and total pressure loss contours and vector plots on the symmetry plane. Calculated details of the primary vortex show excellent agreement with the experimental data. The calculations also show a small secondary vortex not seen experimentally.

Chima, Rodrick V.↗

Structure of an energetic narrow discrete arc

Particle distributions, waves, dc electric fields, and magnetic fields were measured by two sounding rockets at altitudes of 950 and 430 km through an energetic (greater than 5 keV) narrow (about 10 km) stable discrete arc. Although the payloads' magnetic footprints were separated by only 50 km, differences in the arc's structure were observed including the spatial width, peak energy, and characteristic spectra. The energetic electron precipitation included both slowly varying isotropic fluxes that formed an inverted-V energy-time signature and rapidly varying field-aligned fluxes at or below the isotropic spectral peak. The isotropic precipitation had a flux discontinuity inside the arc indicating the arc was present on a boundary between two different magnetospheric plasmas. Dispersive and nondispersive bursts of field-aligned electrons were measured throughout the arc, appearing over broad energy ranges or as monoenergetic beams. Dispersive bursts gave variable source distances less than 8000 km. Plateauing of some of the most intense bursts suggests that waves stabilized these electrons. During the lower altitude arc crossing, the field-aligned component formed a separate inverted-V energy-time signature whose peak energy was half the isotropic peak energy.

Mcfadden, J. P.↗

The Atmospheres of Saturn and Titan in the Near-Infrared: First Results of Cassini/VIMS

The wide spectral coverage and extensive spatial, temporal, and phase-angle mapping capabilities of the Visual Infrared Mapping Spectrometer (VIMS) onboard the Cassini-Huygens Orbiter are producing fundamental new insights into the nature of the atmospheres of Saturn and Titan. For both bodies, VIMS maps over time and solar phase angles provide information for a multitude of atmospheric constituents and aerosol layers, providing new insights into atmospheric structure and dynamical and chemical processes. For Saturn, salient early results include evidence for phosphine depletion in relatively dark and less cloudy belts at temperate and mid-latitudes compared to the relatively bright and cloudier Equatorial Region, consistent with traditional theories of belts being regions of relative downwelling. Additional Saturn results include (1) the mapping of enhanced trace gas absorptions at the south pole, and (2) the first high phase-angle, high-spatial-resolution imagery of CH4 fluorescence. An additional fundamental new result is the first nighttime near-infrared mapping of Saturn, clearly showing discrete meteorological features relatively deep in the atmosphere beneath the planet's sunlit haze and cloud layers, thus revealing a new dynamical regime at depth where vertical dynamics is relatively more important than zonal dynamics in determining cloud morphology. Zonal wind measurements at deeper levels than previously available are achieved by tracking these features over multiple days, thereby providing measurements of zonal wind shears within Saturn's troposphere when compared to cloudtop movements measured in reflected sunlight. For Titan, initial results include (1) the first detection and mapping of thermal emission spectra of CO, CO2, and CH3D on Titan's nightside limb, (2) the mapping of CH4 fluorescence over the dayside bright limb, extending to approximately 750 km altitude, (3) wind measurements of approximately 0.5 ms(exp -1), favoring prograde, from the movement of a persistent (multiple months) south polar cloud near 88 deg S latitude, and (4) the imaging of two transient mid-southern-latitude cloud features.

Cassini Mission↗

Entropy Stable Staggered Grid Spectral Collocation for the Burgers' and Compressible Navier-Stokes Equations

Staggered grid, entropy stable discontinuous spectral collocation operators of any order are developed for Burgers' and the compressible Navier-Stokes equations on unstructured hexahedral elements. This generalization of previous entropy stable spectral collocation work [1, 2], extends the applicable set of points from tensor product, Legendre-Gauss-Lobatto (LGL) to a combination of tensor product Legendre-Gauss (LG) and LGL points. The new semi-discrete operators discretely conserve mass, momentum, energy and satisfy a mathematical entropy inequality for both Burgers' and the compressible Navier-Stokes equations in three spatial dimensions. They are valid for smooth as well as discontinuous flows. The staggered LG and conventional LGL point formulations are compared on several challenging test problems. The staggered LG operators are significantly more accurate, although more costly to implement. The LG and LGL operators exhibit similar robustness, as is demonstrated using test problems known to be problematic for operators that lack a nonlinearly stability proof for the compressible Navier-Stokes equations (e.g., discontinuous Galerkin, spectral difference, or flux reconstruction operators).

Carpenter, Mark H.↗

Maven Observations of Electron-Induced Whistler Mode Waves in the Martian Magnetosphere

We report on narrowband electromagnetic waves at frequencies between the local electron cyclotron and lower hybrid frequencies observed by the Mars Atmosphere and Volatile EvolutioN (MAVEN) spacecraft in the Martian induced magnetosphere. The peaked electric field wave spectra below the electron cyclotron frequency were first observed by Phobos-2 in the Martian magnetosphere, but the lack of magnetic field wave data prevented definitive identification of the wave mode and their generation mechanisms remain unclear. Analysis of electric and magnetic field wave spectra obtained by MAVEN demonstrates that the observed narrowband waves have properties consistent with the whistler mode. Linear growth rates computed from the measured electron velocity distributions suggest that these whistler mode waves can be generated by cyclotron resonance with anisotropic electrons. Large electron anisotropy in the Martian magnetosphere is caused by absorption of parallel electrons by the collisional atmosphere. The narrowband whistler mode waves and anisotropic electrons are observed on both open and closed field lines and have similar spatial distributions in MSO and planetary coordinates. Some of the waves on closed field lines exhibit complex frequency-time structures such as discrete elements of rising tones and two bands above and below half the electron cyclotron frequency. These MAVEN observations indicate that whistler mode waves driven by anisotropic electrons, which are commonly observed in intrinsic magnetospheres and at unmagnetized airless bodies, are also present at Mars. The wave-induced electron precipitation into the Martian atmosphere should be evaluated in future studies.

martian↗

Performance Assessment of the Eastern Range Mosquito Lagoon 915-MHz Doppler Radar Wind Profiler

Please find enclosed the performance assessment for the Mosquito Lagoon 915-MHz Doppler Radar Wind Profiler (DRWP). This assessment was conducted following the installation of the new DRWP systems at the Eastern Range. DRWP profiles were compared to Automated Meteorological Profiling System balloons launched from Cape Canaveral Space Force Station. Launching balloons from the DRWP site was discussed, but was determined to be unfeasible. Therefore, these results include errors due to the spatial separation of winds. Mean wind component deltas generally ranged between -1.4-2.0 m/s between the balloons and DRWP systems. Data availability at discrete altitudes ranged from 7% to 98% with the least amount of available profiles being in the upper altitude ranges between 4000 m and 5000 m as well as below 300 m. The effective vertical resolution was not able to be estimated in this study, likely due to the relatively small dataset. The results from this analysis were requested by the United States Space Force for use of the acceptance of the DRWP system and for internal EV44 reference.

meteorology↗

OpenSn: A massively parallel, open-source simulation environment for discrete ordinates radiation transport

OpenSn is an open-source, massively parallel deterministic radiation transport code for solving the discrete-ordinates ( S N ) form of the Boltzmann transport equation on unstructured, arbitrary polyhedral meshes. It supports high-fidelity simulations involving steady-state, eigenvalue, and adjoint problems for neutral particles (e.g., neutrons, photons, multi-particles), using the multigroup approximation in energy. OpenSn combines angular discretization via discrete ordinates with a discontinuous Galerkin finite element method (DGFEM) in space, enabling accurate resolution of transport physics on arbitrary polyhedral cells, included locally refined spatial grids. It includes multiple angular quadrature types, including locally refined angular quadratures. Written in modern C++ with a Python API, OpenSn runs efficiently on platforms ranging from laptops to supercomputers. The transport sweep algorithm is implemented using a task-based, directed-acyclic-graph (DAG) approach for each angle and supports asynchronous parallelism across thousands of MPI ranks. Group-set aggregation improves compute intensity, and synthetic acceleration techniques (e.g., diffusion synthetic acceleration, second-moment method) enhance solver convergence. OpenSn has been verified on reactor physics problems and demonstrated excellent weak and strong scaling performance on more than 32,768 processes, making it a versatile and robust platform for large-scale transport simulations in complex geometries.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Impact of Numerical Hydrodynamics in Turbulent Mixing Transition Simulations

Underresolved simulations are unavoidable in high Reynolds (Re) and Mach (Ma) number turbulent flow applications at scale. Implicit large-Eddy simulation (ILES) often becomes the effective strategy to capture the dominating effects of convectively driven flow instabilities. We evaluate the impact of three distinct numerical strategies in simulations of transition and turbulence decay with ILES: the Harten–Lax–van Leer (HLL) Riemann solver applying Strang splitting and a Lagrange-plus-Remap formalism to solve the directional sweep—denoted split; the Harten–Lax–Van Leer-Contact (HLLC) Riemann solver using a directionally unsplit strategy and parabolic reconstruction—denoted unsplit; and the HLLC Riemann solver using unsplit and a low-Ma correction (LMC)—denoted unsplit*. Three case studies are considered: (1) a shock tube problem prototyping shock-driven turbulent mixing, (2) the Taylor–Green Vortex (TGV) prototyping transition to turbulence, and, (3) an homogeneous isotropic turbulence (HIT) case, focusing on the impact of discretization on transition and decay from fixed well-characterized initial conditions. Significantly more accurate predictions are provided by the unsplit schemes, in particular, when augmented with the LMC. For given resolution, only the unsplit schemes predict the turbulent mixing transition after reshock observed in the shock tube experiments. Relevant comparisons of ILES based on Euler and Navier–Stokes equations addressing potential occurrence of low-Re regimes in the applications are presented. Unsplit schemes are instrumental in allowing to capture the spatial development of the TGV flow and its validation at prescribed Re with significantly less resolution. HIT analysis confirms higher simulated turbulence Re and increased small-scale content associated with the unsplit discretizations.

74 ATOMIC AND MOLECULAR PHYSICS↗

Update and review of accuracy assessment techniques for remotely sensed data

Research performed in the accuracy assessment of remotely sensed data is updated and reviewed. The use of discrete multivariate analysis techniques for the assessment of error matrices, the use of computer simulation for assessing various sampling strategies, and an investigation of spatial autocorrelation techniques are examined.

Congalton, R. G.↗

A minimum entropy principle in the gas dynamics equations

Let u(x bar,t) be a weak solution of the Euler equations, governing the inviscid polytropic gas dynamics; in addition, u(x bar, t) is assumed to respect the usual entropy conditions connected with the conservative Euler equations. We show that such entropy solutions of the gas dynamics equations satisfy a minimum entropy principle, namely, that the spatial minimum of their specific entropy, (Ess inf s(u(x,t)))/x, is an increasing function of time. This principle equally applies to discrete approximations of the Euler equations such as the Godunov-type and Lax-Friedrichs schemes. Our derivation of this minimum principle makes use of the fact that there is a family of generalized entrophy functions connected with the conservative Euler equations.

Tadmor, E.↗

Sensor/actuator placement for flexible space structures

A new approach for the placement of sensors and actuators in the active control of flexible space structures is developed. The approach converts the discrete nature of the sensor and actuator positioning problem to a nonlinear programming optimization through approximation of the control forces and output measurements by spatially continuous functions. The locations of the sensors and actuators are optimized in order to move the transmission zeros of the system farther to the left of the imaginary axis. The criterion for sensor/actuator placement can be quite useful for optimal regulation and tracking problems, as well as for low-authority controller designs. Two performance metrics are considered for the optimization and are applied to the sensor/actuator positioning of a large-order flexible space structure.

Maghami, P. G.↗

Illinois Disasters: Utilizing NASA Earth Observations to Enhance Drought Monitoring in Illinois

Drought and flooding in Illinois have severe impacts on the communities and ecosystems of the state. Soil moisture is a valuable indicator of drought and flood vulnerability but can be difficult to measure since in situ monitoring is limited to discrete stations throughout the state. The team created a framework to compare in situ, modeled, and NASA satellite soil moisture measurements to increase the spatial coverage of soil moisture monitoring. The team partnered with the Illinois State Weather Survey, USDA Midwest Climate Hub, NOAA Regional Climate Services of the Central Region, NOAA National Integrated Drought Information System’s Midwest Drought Early Warning System, and the NOAA North Central River ForecastCenter. The team standardized and compared soil moisture data from NASA’s Soil Moisture Active Passive (SMAP) mission, modeled soil moisture outputs from NASA’s SPoRT Land Information System (SPoRT-LIS), and in situ measurements from the Illinois State Weather Survey’s Water and Atmospheric Resources Monitoring (WARM) program. Compared to the WARM data, the satellite and modeled data showed seasonally variable differences and bias. The difference was highest in the winter months and lowest in the late summer and early fall months for the SMAP and SPoRT-LIS data products. SPoRT-LIS produced lower seasonal variability and SMAP demonstrated higher correlation values and lower differences. These analyses suggest that both SMAP and SPoRT-LIS products offer unique strengths and limitations when used for soil moisture monitoring

Joshua Green↗

X-ray Polarization Detection of Cassiopeia A with IXPE

We report on a ∼5σ detection of polarized 3–6 keV X-ray emission from the supernova remnant Cassiopeia A (Cas A) with the Imaging X-ray Polarimetry Explorer (IXPE). The overall polarization degree of 1.8% ± 0.3% is detected by summing over a large region, assuming circular symmetry for the polarization vectors. The measurements imply an average polarization degree for the synchrotron component of ∼2.5%, and close to 5% for the X-ray synchrotron-dominated forward shock region. These numbers are based on an assessment of the thermal and nonthermal radiation contributions, for which we used a detailed spatial-spectral model based on Chandra X-ray data. A pixel-by-pixel search for polarization provides a few tentative detections from discrete regions at the ∼ 3σ confidence level. Given the number of pixels, the significance is insufficient to claim a detection for individual pixels, but implies considerable turbulence on scales smaller than the angular resolution. Cas A's X-ray continuum emission is dominated by synchrotron radiation from regions within ≲1017 cm of the forward and reverse shocks. We find that (i) the measured polarization angle corresponds to a radially oriented magnetic field, similar to what has been inferred from radio observations; (ii) the X-ray polarization degree is lower than in the radio band (∼5%). Since shock compression should impose a tangential magnetic-field structure, the IXPE results imply that magnetic fields are reoriented within ∼1017 cm of the shock. If the magnetic-field alignment is due to locally enhanced acceleration near quasi-parallel shocks, the preferred X-ray polarization angle suggests a size of 3 × 1016 cm for cells with radial magnetic fields.

Polarimetry↗

Final Report for Grant

This is the final report for this project. The major objectives were to develop full models and standard cells of CMOS silicon photomultipliers (SiPM) for neutron detection applications. Specifically, the research seeks to replace photomultiplier tubes, used in neutron detection applications, with SiPMs containing integrated readout electronics fabricated in standard CMOS processes. The detectors improve on the detection efficiency, readout speed, resolution (sensitivity), and spatial resolution of current systems. In particular, this work involved collaboration with the instrument and detector group at ORNL who build detectors based on discrete commercially available PMTs and SiPMs. The SiPM detectors developed in this work are based on Perimeter Gated Single Photon Avalanche Diodes (PGSPAD). These devices have tunable breakdown voltage which affects other performance characteristics of the detectors. The ability of the detectors to have integrated readout and tune the noise floor and dynamic range in addition to offering a method for non-uniformity compensation in a compact form significantly impacts the imaging research field and includes the nuclear sciences field, circuits and systems field in addition to varied fields such as geography and biology/medicine.

42 ENGINEERING↗

Nonlinear soliton confinement in weakly coupled antiferromagnetic spin chains

In this work, we analyze the low-energy dynamics of quasi-one-dimensional, large-$\textit{S}$ quantum antiferromagnets with easy-axis anisotropy, using a semiclassical nonlinear sigma model. The saddle point approximation leads to a sine-Gordon equation which supports soliton solutions. These correspond to the movement of spatially extended domain walls. Long-range magnetic order is a consequence of a weak interchain coupling. Below the ordering temperature, the coupling to nearby chains leads to an energy cost associated with the separation of two domain walls. From the kink-antikink two-soliton solution, we compute the effective confinement potential. At distances large compared to the size of the solitons the potential is linear, as expected for pointlike domain walls. At small distances the gradual annihilation of the solitons weakens the effective attraction and renders the potential quadratic. From numerically solving the effective one-dimensional Schrödinger equation with this nonlinear confinement potential we compute the soliton bound state spectrum. We apply the theory to CaFe 2 O 4 , an anisotropic $\textit{S}$ = 5/2 magnet based upon antiferromagnetic zigzag chains. Using inelastic neutron scattering, we are able to resolve seven discrete energy levels for spectra recorded slightly below the Néel temperature $T_N$ ≈ 200 K. These modes are well described by our nonlinear confinement model in the regime of large spatially extended solitons.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗