Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Spatial dependence”

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 1,765 records · Page 98

Magnetic Reconnection in a Sheared Magnetic Flux Tube: Slippage Versus Tearing

The process of magnetic reconnection in a flux tube can occur in a time-stationary fashion as slippage reconnection or in a time-dependent manner based on, for example, the tearing instability. However, it is not well known under which conditions a system can sustain slippage reconnection. Likewise, it is unclear whether systems can exhibit slippage reconnection and time-dependent reconnection simultaneously. In order to investigate these questions, we employ a set of 3D magnetohydrodynamic simulation. Using these simulations, we model a twisted flux tube with a spatially localized resistive region in the center of the simulation box, and an applied driver on one end plane of the simulation box. As a result, the Poynting flux injected at the boundary propagates to the resistive region and dissipates there. If the driver is nearly circular, we find a predominance of slippage reconnection, which decouples the applied driver entirely from one half of the flux tube. Second, a tearing-like instability occurs with slippage reconnection in the resistive region when the resistivity is strong enough and the cross section of the flux tube, as defined by the shape of the applied velocity perturbation, is elliptical enough. In this case we find the surprising feature that magnetic field perturbation generated by the tearing-like instability radiates away from the resistive region to influence the entire flux tube shape.

Hidetaka Kuniyoshi↗

NASA TROPOMI Aerosol Products: Algorithmic Upgrades and Preliminary Evaluation

This poster presentation describes an expanded NASA TROPOMI (Tropospheric Monitoring Instrument)aerosol algorithm (N-TROPOMAER) that takes advantage of TROPOMI observations in the ultraviolet and visible spectral regions. The availability of the Oxygen B-band observations, and the unprecedentedly high spatial resolution (3.5 km X 5.5 km) for a hyper-spectral sensor are significant improvements for aerosol properties retrieval. The heritage N-TROPOMAER aerosol algorithm uses near-ultraviolet radiances at 354 nm and 388 nm from Sentinel 5 Precursor-TROPOMI for simultaneously retrieving aerosol optical depth (AOD), single-scattering albedo (SSA), aerosol absorption optical depth (AAOD), and above-cloud aerosol optical depth(ACAOD) at 388 nm, along with the qualitative UV aerosol index (UVAI). We have expanded the inversion capability beyond the UV, to retrieve AOD at 466 nm and 680 nm. Surface reflectance effects at466 nm are accounted for using a recently developed geometry-dependent surface Lambertian-equivalent reflectivity (GLER) product, which is derived from the top-of-atmosphere radiance computed with Rayleigh scattering and surface bidirectional reflectance distribution function (BRDF) for the exact viewing geometry at the sensor’s spatial resolution. Aerosol layer height (ALH) and 680 nm AOD are simultaneously derived from observations at 680 nm and at the Oxygen-B band (688 nm). Another important upgrade is the use of time averaged total column carbon monoxide from the NASA GEOS-CF(Global Earth Observing System Composition Forecast) as a tracer of carbonaceous aerosols.

TROPOMI↗

MPACT 4.4 Theory Manual

MPACT is a three-dimensional (3D) full-core neutron transport code capable of calculating subpin power distributions. Calculations are based on the Boltzmann transport equation for neutron fluxes for problems in which the detailed geometrical configuration of fuel components such as the pellet and cladding are explicitly retained. The cross-section data needed for the neutron transport calculation are obtained directly from a multigroup cross section library, which has traditionally been used by lattice physics codes to generate few-group homogenized cross sections for nodal core simulators. Hence, MPACT assumes neither a priori homogenization nor group condensation for the full core spatial solution. The 3D MPACT transport solution can be obtained using the method of characteristics (MOC), which employs discrete ray tracing within each fuel pin. However, for practical reactor applications, the direct application of MOC to 3D core configurations requires an excessive amount of memory and computing time due to the very large number of rays. For practical 3D full-core calculations, MPACT commonly uses an approximate “2D/1D” method that treats the radial (x and y) variables differently from the axial (z) variable. In particular, the radial dependence of the solution is calculated using transport theory, and the axial dependence is calculated using diffusion or P 3 theory. The 2D/1D method requires the core to be divided into a vertical stack of axial slices with a thickness of Δ z ≈ 5–10 cm. Each axial slice is divided radially into coarse spatial cells with boundaries that usually constitute the pin cell boundaries, for which Δ x = Δ y ≈ 1.5 cm. Then, each coarse radial cell (pin cell) is divided into 50–100 fine radial cells, which resolve the angular flux in the fuel, cladding, and moderator regions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The simulated features of heliospheric cosmic-ray modulation with a time-dependent drift model. IV - The role of heliospheric neutral sheet deformation

Previous calculations with a time-dependent drift model revealed the model to be less successful in describing time-dependent modulation during periods of moderate to large solar activity. In this paper, it is argued that a major reason for this is that the previously used wavy heliospheric neutral sheet (HNS) description was based on an idealized HNS not subject to any spatial evolution while propagating radially outward. It is suggested that the deformation and compression of HNS wave peaks will lead to significant increases in the crossfield diffusion across these peaks (short-circuiting). The cosmic rays will effectively experience reduced tilt angles and therefore a reduction in the integrated HNS modulation effect between an observer and the heliospheric boundary. During periods of moderate to large solar activity these HNS deformation processes are progressively more frequent and should lead to a significant reduction in time-dependent modulation as predicted by drift models. Calculations done with radially propagating tilt angles that effectively decrease with radial distance give the expected reduction which improves the general description of modulation from 1987-1988.

Le Roux, J. A.↗

Development and Evaluation of Ensemble Consensus Precipitation Estimates over High Mountain Asia

Precipitation estimates are highly uncertain in complex regions such as High-Mountain Asia (HMA), where ground measurements are very difficult to obtain, and atmospheric dynamics poorly understood. Though gridded products derived from satellite-based observations and/or reanalysis can provide temporally and spatially distributed estimates of precipitation, there are significant inconsistencies in these products. As such, to date, there is little agreement in the community on the best and most accurate gridded precipitation product in HMA, which is likely area dependent because of HMA’s strong heterogeneities and complex orography. Targeting these gaps, this article presents the development of a consensus ensemble precipitation product using three gridded precipitation datasets (the Integrated Multi-satellitE Retrieals for Global Precipitation Measurement IMERG, the Climate Hazards group Infrared Precipitation with Stations CHIRPS, and the ECMWF Reanalysis ERA5) with a localized probability matched mean (LPM) approach. We evaluate the performance of the LPM estimate along with a simple ensemble mean (EM) estimate to overcome the differences and disparities of the three selected constituent products on long-term averages and trends in HMA. Our analysis demonstrates that LPM reduces the high biases embedded in the ensemble members and provides more realistic spatial patterns compared to EM. LPM is also a good alternative for merging data products with different spatio-temporal resolutions. By filtering disparities among the individual ensemble members, LPM overcomes the problem of a certain product performing well only in a particular area and provides a consensus estimate with plausible temporal trends.

Fadji Z Maina↗

Pickup Ion Distributions from Three Dimensional Neutral Exospheres

Pickup ions formed from ionized neutral exospheres in flowing plasmas have phase space distributions that reflect their source's spatial distributions. Phase space distributions of the ions are derived from the Vlasov equation with a delta function source using three.dimensional neutral exospheres. The ExB drift produced by plasma motion picks up the ions while the effects of magnetic field draping, mass loading, wave particle scattering, and Coulomb collisions near a planetary body are ignored. Previously, one.dimensional exospheres were treated, resulting in closed form pickup ion distributions that explicitly depend on the ratio rg/H, where rg is the ion gyroradius and H is the neutral scale height at the exobase. In general, the pickup ion distributions, based on three.dimensional neutral exospheres, cannot be written in closed form, but can be computed numerically. They continue to reflect their source's spatial distributions in an implicit way. These ion distributions and their moments are applied to several bodies, including He(+) and Na(+) at the Moon, H(+2) and CH(+4) at Titan, and H+ at Venus. The best places to use these distributions are upstream of the Moon's surface, the ionopause of Titan, and the bow shock of Venus.

Hartle, R. E.↗

Benchmark of the KGMf with a coupled Boltzmann equation solver

The Kinetic Global Model framework (KGMf) is an open-source general-purpose global model (spatially averaged) simulation code developed to explore the reaction kinetics and pathways in plasma discharge systems. It contains species continuity and electron energy balance equations, with a time-dependent evaluated electron energy distribution function (EEDF) for electron impact reactions. The EEDF is utilized to determine the rate coefficients for electron impact reactions, which can have profound impact on the temporal evolutions of plasma parameters. Previously, the EEDF was commonly assumed as Maxwellian or an analytical function of the effective electron temperature. In this work, the KGMf is coupled with a Boltzmann equation (BE) solver to self-consistently compute the EEDF. The EEDF evolution frequency is determined based on relative changes of the reduced electric field. The KGMf is benchmarked with the ZDPlasKin code based on high-pressure low-temperature argon plasma discharge cases. Additionally, the temporal evolutions of reduced electric field, electron temperature, EEDF, reaction rates, and species densities, are obtained and compared under different discharge conditions, showing good agreement between the KGMf and the ZDPlasKin simulations. The application of the KGMf for predicting breakdown times in high power microwave discharges is also presented, which shows qualitative agreement with particle-in-cell simulations. The KGMf can be further applied for more complicated plasma discharge systems, where the reaction kinetics are more intricate (e.g., plasma-assisted combustion systems).

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

The effect of original and initial saturation on residual nonwetting phase capillary trapping efficiency

Injection of supercritical carbon dioxide (CO 2 ) into geological formations is a strategy for both atmospheric greenhouse gas reduction (climate change mitigation) and enhanced oil recovery. To understand CO 2 trapping efficiency, the capillary trapping behaviors that immobilize subsurface fluids were analyzed at the pore-scale using pairs of proxy fluids representing the range of in situ nonwetting and wetting fluid properties encountered in geologic storage reservoirs. The pairs of fluids were cycled through imbibition and drainage processes using a flow cell apparatus containing a sintered glass bead column. Computed x-ray microtomography (microCT) was used to identify immobilized nonwetting fluid after imbibition and drainage events. From microCT images, the trapped residual (post-secondary imbibition) nonwetting phase was spatially correlated to both the original (post-primary imbibition) and the initial (post-primary drainage) nonwetting phase; this relationship is referred to here as the original saturation dependence (SO-dependence) and initial saturation dependence (SI-dependence), respectively. Significant trends of decreasing SO- and SI-dependence with increasing wetting and nonwetting fluid phase viscosities were observed. This finding implies that the amount of CO 2 injected and ultimately trapped is dependent on the nonwetting phase (e.g. oil or gas) already present in the formation, as well as on the manner in which supercritical CO 2 is initially injected, which are factors not considered in current trapping models. To explore the potential effect that varying viscosity and interfacial tension (IFT) might have on this process, we report on a variety of fluid pairs with different viscosities and IFT.

42 ENGINEERING↗

Pro5 is not essential for the formation of ‘Ni-hook’ in nickel superoxide dismutase

The N-terminus of nickel-dependent superoxide dismutase (NiSOD) forms a structural motif known as the "Ni-hook," where the peptide wraps around the metal to bring cysteine-2 and cysteine-6 into spatial proximity, allowing these residues to coordinate in a cis-geometry. A highly conserved proline-5 residue in the Ni-hook adopts a cis-conformation that is widely considered important for its formation. Herein, we investigate this role by point mutation of Pro5 to alanine. The results obtained show that the variant exhibits wild-type-like redox catalysis and features a Ni(III) center very similar to that found in enzyme. Structural analysis using x-ray absorption spectroscopy of the nickel sites in as-isolated P5A-NiSOD reveals changes in the variant and are consistent with a six-coordinate Ni site with (N/O) 4 S 2 coordination. These changes are attributed to changes in the Ni(II) site structure. Nickel-binding studies using isothermal titration calorimetry reveal two binding events with K d = 25(20) nM, and 250(60) nM. These events are attributed to i) Ni(II) binding to a preformed Ni-hook containing cis-Pro5 and ii) the combination of trans- to cis- isomerization upon Ni(II) binding, respectively. Finally, the higher-affinity binding event is absent in P5A-NiSOD, an observation attributed to the low abundance of the cis-Ala5 isomer in the apo-protein.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

X-ray nano-imaging of defects in thin film catalysts via cluster analysis

Functional properties of transition-metal oxides strongly depend on crystallographic defects; crystallographic lattice deviations can affect ionic diffusion and adsorbate binding energies. Scanning x-ray nanodiffraction enables imaging of local structural distortions across an extended spatial region of thin samples. Yet, localized lattice distortions remain challenging to detect and localize using nanodiffraction, due to their weak diffuse scattering. Here, in this study, we apply an unsupervised machine learning clustering algorithm to isolate the low-intensity diffuse scattering in as-grown and alkaline-treated thin epitaxially strained SrIrO 3 films. We pinpoint the defect locations, find additional strain variation in the morphology of electrochemically cycled SrIrO 3 , and interpret the defect type by analyzing the diffraction profile through clustering. Our findings demonstrate the use of a machine learning clustering algorithm for identifying and characterizing hard-to-find crystallographic defects in thin films of electrocatalysts and highlight the potential to study electrochemical reactions at defect sites in operando experiments.

42 ENGINEERING↗

Decay of a laser beam under stimulated temperature scattering in toluene as a result of two-photon pump absorption and nonsteady-state laser pulse interaction with matter

The decay of the pump-beam spatial structure under conditions of stimulated temperature scattering of a pulsed laser beam in toluene as a result of two-photon absorption of pump radiation and time-dependent interaction of a laser pulse with a medium is considered. (paper)

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Velocity Dispersion σ aper Aperture Corrections as a Function of Galaxy Properties from Integral-field Stellar Kinematics of 10,000 MaNGA Galaxies

Abstract The second moment of the stellar velocity within the effective radius, denoted by σ e 2 , is a crucial quantity in galaxy studies, as it provides insight into galaxy properties and their mass distributions. However, large spectroscopic surveys typically do not measure σ e directly, instead providing σ aper , the second moment of the stellar velocity within a fixed fiber aperture. In this paper, we derive an empirical aperture correction formula, given by σ aper / σ e = ( R aper / R e ) α , using spatially resolved stellar kinematics extracted from approximately 10,000 Sloan Digital Sky Survey–Mapping Nearby Galaxies at Apache Point Observatory integral field unit observations. Our analysis reveals a strong dependence of α on the r -band absolute magnitude M r , g – i color, and Sérsic index n Ser , where α values are lower for brighter, redder galaxies with higher Sérsic indices. Our results demonstrate that the aperture correction derived from previous literature on early-type galaxies cannot be applied to predict the aperture corrections for galaxies with intermediate Sérsic indices. We provide a lookup table of α values for different galaxy types, with parameters in the ranges of −18 > M r > −24, 0.4 < g – i < 1.6, and 0 < n Ser < 8. A Python script is provided to obtain the correction factors from the lookup table.

Astronomy & Astrophysics↗

Role of electromagnetic gauge-field fluctuations in the selection between chiral and nematic superconductivity

Motivated by the observation of nematic superconductivity in several systems, we revisit the problem of the leading pairing instability of two-component unconventional superconductors on the triangular lattice—such as (p x , p y )-wave and (d x 2 –y 2 , d xy )-wave superconductors. Such a system has two possible superconducting states: the chiral state (e.g., p + ip or d + id), which breaks time-reversal symmetry, and the nematic state (e.g., p + p or d + d), which breaks the threefold rotational symmetry of the lattice. Weak-coupling calculations generally favor the chiral over the nematic superconducting state, raising the question of what mechanism can stabilize the latter. Here, we show that the electromagnetic field fluctuations can play a crucial role in selecting between these two states. Specifically, we derive and analyze the effective free energy for the two-component superconducting order parameter after integrating out the gauge-field fluctuations, which is formally justified if the spatial order parameter fluctuations can be neglected. A nonanalytic cubic term arises, as in the case of a conventional s -wave superconductor. However, unlike the latter, the cubic term depends on the relative phase and on the relative amplitudes between the two order parameter components, in such a way that it generally favors the nematic state. This result is a direct consequence of the fact that the stiffness of the superconducting order parameter is not isotropic. Competition with the quartic term, which favors the chiral state, leads to a renormalized phase diagram in which the nematic state displaces the chiral state over a wide region in the parameter space. Here, we analyze the stability of the fluctuation-induced nematic phase, generalize our results to tetragonal lattices, and discuss their applicability to candidate nematic superconductors, including twisted bilayer graphene.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Development and Analysis of Optimal Multilevel Solvers on Advanced Computers. Final Report

Constrained optimization in the context of time dependent, partial differential equations (PDE) leads to a symmetric, block-tridiagonal system of nonlinear equations that must be solved repeatedly in an iterative solution strategy. The blocks represent spatial discretization, while the connection between the blocks represents a forward and backward integration in time. The focus of this project is to apply a parallel-in-time (PiT) solution technique to the large block-triangular system.

97 MATHEMATICS AND COMPUTING↗

Verification of EvaluateFLux Utility Program

The EvaluateFlux program is a post-processing utility program for DIF3D, specifically DIF3D-VARIANT which handles Cartesian and hexagonal geometries. The EvaluateFlux program was developed to allow users to obtain flux and power traverses through the geometry domain, and its initial purpose was to facilitate foil analysis by evaluating the flux solution from DIF3D-VARIANT and combining it with foil cross section data. The EvaluateFlux program can calculate the neutron flux, as well as the reaction rates, at any user provided evaluation point. It does this by identifying the spatial mesh associated with the evaluation point and then evaluates the polynomial based neutron flux moments stored in the NHFLUX file at that point. The output of EvaluateFlux varies depending on the input setup. The maximum output includes the neutron flux and microscopic and macroscopic reaction rates at each evaluation point. The purpose of this work is to verify the outputs of EvaluateFlux. Simple models that have hand calculatable results are first defined and used to verify the EvaluateFlux outputs. More complex cases are then added where a duplicate program of EvaluateFlux that uses PrintTables outputs of the binary files is used to verify the EvaluateFlux outputs. In those complex cases, hand calculations of selected evaluation points were also displayed to confirm the software verification. For all the tests done, the hand calculations agreed well with those calculated by EvaluateFlux. For the larger complex problems, the duplicate program that can process hundreds of evaluation points was able to identify that zero points within some meshes have large errors. This aspect was attributed to the truncation error on the input provided to the duplicate program and is not a concern for the accuracy of the EvaluateFlux software.

97 MATHEMATICS AND COMPUTING↗

Verification of the EvaluateFlux Utility Program

The EvaluateFlux program is a post-processing utility program for DIF3D, specifically DIF3D-VARIANT which handles Cartesian and hexagonal geometries. The EvaluateFlux program was developed to allow users to obtain flux and power traverses through the geometry domain, and its initial purpose was to facilitate foil analysis by evaluating the flux solution from DIF3D-VARIANT and combining it with foil cross section data. The EvaluateFlux program can calculate the neutron flux, as well as the reaction rates, at any user provided evaluation point. It does this by identifying the spatial mesh associated with the evaluation point and then evaluates the polynomial based neutron flux moments stored in the NHFLUX file at that point. The output of EvaluateFlux varies depending on the input setup. The maximum output includes the neutron flux and microscopic and macroscopic reaction rates at each evaluation point. The purpose of this work is to verify the outputs of EvaluateFlux. Simple models that have hand calculatable results are first defined and used to verify the EvaluateFlux outputs. More complex cases are then added where a duplicate program of EvaluateFlux that uses PrintTables outputs of the binary files is used to verify the EvaluateFlux outputs. In those complex cases, hand calculations of selected evaluation points were also displayed to confirm the software verification. For all the tests done, the hand calculations agreed well with those calculated by EvaluateFlux. For the larger complex problems, the duplicate program that can process hundreds of evaluation points was able to identify that zero points within some meshes have large errors. This aspect was attributed to the truncation error on the input provided to the duplicate program and is not a concern for the accuracy of the EvaluateFlux software.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

The pointable multispectral linear array system and the large format camera for early 1980 shuttle payloads

Design options are described for the multispectral linear arrays which will provide data and facilities beyond those obtained by LANDSAT 1, 2, 3, and D (the thematic mapper) are considered. The features, uses, and risks of a 10M fine resolution pointable imager, a 15M/30M dual resolution pointable imager, a 30M selectable spectral resolution pointing imager, a 30M electronically pointable multi-temporal sampler, and a 40M/80M mixed resolution imager are described. A large format camera scheduled for 1981 shuttle launch which can obtain high spatial ground resolution (approximately 10M), a relative position accuracy of plus or minus 14M, and a relative heightening accuracy, by virtue of stereographic coverage, of between 7 and 28M depending on the base-to-height ratio used. A potential problem of the system is that it will orbit in the shuttle.

Slater, P. N.↗

Three-dimensional blade vortex interactions

A three-dimensional time dependent Navier-Stokes analysis was applied to the rotor blade vortex interaction problem. The numerical procedure is an iterative implicit procedure using three point central differences to represent spatial derivatives. A series of calculations were made to determine the time steps, pseudo-time steps, iterations, artificial dissipation level, etc. required to maintain a nondissipative vortex. Results show the chosen method to have excellent non-dissipative properties provided the correct parameters are chosen. This study was used to set parameters for both two- and three-dimensional blade vortex interaction studies. The case considered was the interaction between a vortex and a NACA0012 airfoil. The results showed the detailed physics during the interaction including the pressure pulse propagating from the blade. The simulated flow physics was qualitatively similar to that experimentally observed. The BVI phenomena is the result of the buildup and violent collapse of the shock waves and local supersonic pockets on the blade surfaces. The resulting pressure pulse build-up appears to be centered at the blade leading edge.

Davoudzadeh, Farhad↗