SEARCH · Search NASA

Results for “Discontinuity Factors”

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 73 records · Page 4

All sources NASA NTRS DOE OSTI

DOE OSTI2020

A quadratic programming flux correction method for high-order DG discretizations of S transport

In this work, we present a new flux-fixup approach for arbitrarily high-order discontinuous Galerkin (DG) discretizations of the S N transport equation, and we demonstrate the compatibility of this approach with the Variable Eddington Factor (VEF) method. The new fixup approach is sweep-compatible: during a transport sweep (block Gauss-Seidel iteration in which the scattering source is lagged), a local quadratic programming (QP) problem is solved in each spatial element to ensure that the solution satisfies certain physical constraints, including local particle balance. In this paper, we describe two choices of physical constraints, resulting in two variants of the method: QP Zero (QPZ) and QP Maximum Principle (QPMP). In QPZ, the finite element coefficients of the solution are constrained to be nonnegative. In QPMP, they are constrained to adhere to an approximate discrete maximum principle. There are two primary takeaways in this paper. First, when the positive Bernstein basis is used for DG discretization, the QPMP method eliminates negativities, preserves high-order accuracy for smooth problems, and significantly dampens unphysical oscillations in the solution. The latter feature – the dampening of unphysical oscillations – is an improvement upon standard, simpler fixup approaches such as the approach described in (denoted as the “zero and rescale” (ZR) method in this paper). This improvement comes at a moderate computational cost, but it is not prohibitive. Our results show that, even in an unrealistic worst-case scenario where 83% of the spatial elements require a fixup, the computational cost of performing a transport sweep with fixup is only ~31% greater than performing one without fixup. The second takeaway is that the VEF method can be used to accelerate the convergence of transport sweeps even when a fixup is applied. When optically thick regions are present, transport sweeps converge slowly, regardless of whether a fixup is applied, and acceleration is needed. However, attempting to apply standard diffusion synthetic acceleration (DSA) to fixed-up transport sweeps results in divergence for optically thick problems. Our results show that the same is not true for VEF. When VEF is combined with fixed-up transport sweeps, the result is a scheme that produces a nonnegative solution, converges independently of the mean free path, and, in the case of the QPMP fixup, adheres to an approximate discrete maximum principle.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1984

Crack problems for bonded nonhomogeneous materials under antiplane shear loading

The singular nature of the crack tip stress field in a nonhomogeneous medium with a shear modulus with a discontinuous derivative was investigated. The simplest possible loading and geometry, the antiplane shear loading of two bonded half spaces in which the crack is perpendicular to the interface is considered. It is shown that the square root singularity of the crack tip stress field is unaffected by the discontinuity in the derivative of the shear modulus. The problem is solved for a finite crack and results for the stress intensity factors are presented.

Erdogan, F.↗

NASA NTRS2002

A Method of Efficient Inclination Changes for Low-thrust Spacecraft

The evolution of low-thrust propulsion technologies has reached a point where such systems have become an economical option for many space missions. The development of efficient, low trip time control laws has received an increasing amount of attention in recent years, though few studies have examined the subject of inclination changing maneuvers in detail. A method for performing economical inclination changes through the use of an efficiency factor is derived front Lagrange's planetary equations. The efficiency factor can be used to regulate propellant expenditure at the expense of trip time. Such a method can be used for discontinuous-thrust transfers that offer reduced propellant masses and trip-times in comparison to continuous thrust transfers, while utilizing thrusters that operate at a lower specific impulse. Performance comparisons of transfers utilizing this approach with continuous-thrust transfers are generated through trajectory simulation and are presented in this paper.

Falck, Robert↗

DOE OSTI2022

Efficiency and Forward Voltage of Blue and Green Lateral LEDs with V-shaped Defects and Random Alloy Fluctuation in Quantum Wells

For nitride-based blue and green light-emitting diodes (LEDs), the forward voltage V for is larger than expected, especially for green LEDs. This is mainly due to the large barriers to vertical carrier transport caused by the total polarization discontinuity at multiple quantum well and quantum barrier interfaces. The natural random alloy fluctuation in quantum wells has proven to be an important factor reducing V for . However, this does not suffice in the case of green LEDs because of their larger polarization-induced barrier. V-shaped defects (V-defects) have been proposed as another key factor in reducing V for to allow lateral injection into multiple quantum wells, thus bypassing the multiple energy barriers incurred by vertical transport. In this paper, to model carrier transport in the whole LED, we consider both random-alloy and V-defect effects. A fully two-dimensional drift-diffusion charge-control solver is used to model both effects. The results indicate that the turn-on voltages for blue and green LEDs are both affected by random alloy fluctuations and the V-defect density. For green LEDs, V for decreases more due to V-defects, where the smaller polarization barrier at the V-defect sidewall is the major path for lateral carrier injection. Then, we discuss how the V-defect density and size affects the results.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

NASA NTRS1999

Numerical System Solver Developed for the National Cycle Program

As part of the National Cycle Program (NCP), a powerful new numerical solver has been developed to support the simulation of aeropropulsion systems. This software uses a hierarchical object-oriented design. It can provide steady-state and time-dependent solutions to nonlinear and even discontinuous problems typically encountered when aircraft and spacecraft propulsion systems are simulated. It also can handle constrained solutions, in which one or more factors may limit the behavior of the engine system. Timedependent simulation capabilities include adaptive time-stepping and synchronization with digital control elements. The NCP solver is playing an important role in making the NCP a flexible, powerful, and reliable simulation package.

Binder, Michael P.↗

NASA NTRS2017

Heat Loads Due to Small Penetrations in Multilayer Insulation Blankets

The main penetrations (supports and piping) through multilayer insulation systems for cryogenic tanks have been previously addressed by heat flow measurements. Smaller penetrations due to fasteners and attachments are now experimentally investigated. The use of small pins or plastic garment tag fasteners to each the handling and construction of multilayer insulation (MLI) blankets goes back many years. While it has long been understood that penetrations and other discontinuities degrade the performance of the MLI blanket, quantification of this degradation has generally been lumped into gross performance multipliers (often called degradation factors or scale factors). Small penetrations contribute both solid conduction and radiation heat transfer paths through the blanket. The conduction is down the stem of the structural element itself while the radiation is through the hole formed during installation of the pin or fastener. Analytical models were developed in conjunction with MLI perforation theory and Fouriers Law. Results of the analytical models are compared to experimental testing performed on a 10 layer MLI blanket with approximately 50 small plastic pins penetrating the test specimen. The pins were installed at 76-mm spacing inches in both directions to minimize the compounding of thermal effects due to localized compression or lateral heat transfer. The testing was performed using a liquid nitrogen boil-off calorimeter (Cryostat-100) with the standard boundary temperatures of 293 K and 78 K. Results show that the added radiation through the holes is much more significant than the conduction down the fastener. The results are shown to be in agreement with radiation theory for perforated films.

Structural attachments↗

NASA NTRSDate not supplied

A Deep Neural Network for Achieving Spectrally Consistent and Seamless Infrared Radiance Measurements Across Geostationary Satellite Domains

The NASA Clouds and the Earth's Radiant Energy System (CERES) project provides the scientific community with observed top-of-atmosphere (TOA) shortwave and longwave fluxes for climate monitoring and climate model validation. To achieve this goal, CERES relies on TOA broadband fluxes derived from geostationary satellite (GEO) imagery to account for the diurnal flux variations between the CERES observation intervals. Consistent global flux derivation depends on accurate and consistent cloud retrievals. Scene-dependent spectral measurement inconsistency of the instruments that make up the contiguous ring of GEO observations (GEO-Ring), as well as limb darkening effects, can cause discontinuities in derived cloud properties and radiative fluxes at the boundaries of adjacent imager domains. Although the algorithms utilize radiative transfer models to account for instrument-band-dependent atmospheric correction and viewing zenith angle (VZA) dependency, small discontinuities may persist due to uncertainties inherent to the multiple imager-specific algorithms. Furthermore, while hyperspectral-instrument-based spectral band adjustment factors may effectively account for spectrally induced bias, they are less effective at reducing variance owed to the specific composition of the viewed scene, which is challenging to robustly characterize. As such, this article highlights the use of a deep neural network (DNN) to resolve spectral-and VZA-induced biases between GEO-Ring imagers. The DNN uses available infrared (IR) channels from the GEO instruments, along with viewing and solar illumination geometry, to estimate homogenized, VIIRS-like IR radiances for use in the GEO cloud algorithm. This approach is effective at mitigating scene-dependent spectral variance and VZA dependency, resulting in consistent radiance measurements across the GEO-Ring, thereby leading toward a more seamless global cloud assessment.

deep learning↗

NASA NTRS1995

Correction Factor for Determining the London Penetration Depth from Strip Resonators

A significant disagreement is often seen between the theoretical temperature dependent magnetic penetration depth profile and experimentally derived calculations based on stripline type resonators. This short paper shows that the disagreement can be attributed to the susceptance coupled into the resonator from the gap discontinuity as well as the feed line. When the effect is taken into account, the natural resonant frequency of the resonator is increased, and the frequency shift due to kinetic inductance can be calculated much more accurately. While it is necessary to include this effect to determine the penetration depth, it is shown that the impact on unloaded quality factor is generally negligible. The situation when the strip characteristic impedance is not matched to the generator is included.

Romanofsky, Robert R.↗

DOE OSTI2024

Langmuir turbulence in suspended kelp farms

This study investigates the influence of suspended kelp farms on ocean mixed layer hydrodynamics in the presence of currents and waves. We use the large eddy simulation method, where the wave effect is incorporated by solving the wave-averaged equations. Distinct Langmuir circulation patterns are generated within various suspended farm configurations, including horizontally uniform kelp blocks and spaced kelp rows. Intensified turbulence arises from the farm-generated Langmuir circulation, as opposed to the standard Langmuir turbulence observed without a farm. The creation of Langmuir circulation within the farm is attributed to two primary factors depending on farm configuration: (i) enhanced vertical shear due to kelp frond area density variability, and (ii) enhanced lateral shear due to canopy discontinuity at lateral edges of spaced rows. Both enhanced vertical and lateral shear of streamwise velocity, representing the lateral and vertical vorticity components, respectively, can be tilted into downstream vorticity to create Langmuir circulation. This vorticity tilting is driven by the Craik–Leibovich vortex force associated with the Stokes drift of surface gravity waves. In addition to the farm-generated Langmuir turbulence, canopy shear layer turbulence is created at the farm bottom edge due to drag discontinuity. The intensity of different types of turbulence depends on both kelp frond area density and the geometric configuration of the farm. The farm-generated turbulence has substantial consequences for nutrient supply and kelp growth. These findings also underscore the significance of the presence of obstacle structures in modifying ocean mixed layer characteristics.

Bo, Tong (ORCID:0000000260300561)↗

NASA NTRS1988

On governing equations for crack layer propagation

Results of analysis on damage distribution of a crack layer, in a model material, supported the self-similarity hypothesis of damage evolution which has been adopted by the crack layer theory. On the basis of measurements of discontinuity density and the double layer potential technique, a solution to the crack damage interaction problem has been developed. Evaluation of the stress intensity factor illustrated the methodology. Analysis of experimental results showed that Arrhenius type constitutive relationship described very well the expansion of the active zone of a crack layer.

Chudnovsky, A.↗

DOE OSTI2021

Linear and Nonlinear Solvers for Simulating Multiphase Flow within Large-Scale Engineered Subsurface Systems

Simulation of multiphase flow in the subsurface is well-known to be computationally challenging. While there have been many studies that have explored approaches to overcoming these challenges, they often utilize relatively simple case studies. In this paper, we focus on the unique numerical challenges posed by modeling large-scale engineered subsurface systems, characterized by discrete features embedded in a heterogeneous natural subsurface setting. The man-made features such as shafts, tunnels, and barriers often cause multiple challenges in modeling the domain for multiphase porous media flow. This flow scenario can have a wide range of applications such as nuclear waste repositories, enhanced recovery of a petroleum reservoir, geothermal engineering, and carbon sequestration. An example of these severe numerical challenges is the case of performance assessment (PA) for Waste Isolation Pilot Plant (WIPP), the only operating deep geological repository in the US, which simulates extreme material properties of bedded salt rock formation and extreme contrast due to open excavation next to the formation. The models have extremes not only of permeability and porosity but also of the constitutive models needed for multiphase flow; additionally, they have process models like salt creep closure reducing porosity over time, fracturing in clay and anhydrite interbeds of the bedded salt, gas generation from the waste materials, and unintentional human borehole intrusions in some scenarios. Numerical simulations require the solution of coupled systems of nonlinear PDEs; in our work, we use the open-source simulator PFLOTRAN which is based on Finite Volume discretization. The solution of the nonlinear equations requires use of the Newton-Raphson iteration at each time step, which entails the solution of the linearized Jacobian system at each iteration. The effects of all the processes (i.e., large number of unknowns, highly nonlinear constitutive relations, large contrasts in material properties in short distances) lead to an ill-conditioned Jacobian matrix that severely challenges traditional linear solver, i.e., stabilized biconjugate gradient with block Jacobi incomplete LU preconditioner (BCGS-ILU) leading to non-convergence for traditional Newton-Raphson nonlinear solver causing unacceptably long computation time for each model. This paper presents linear solvers such as constrained pressure residual (CPR) two-stage preconditioner with alternate-block-factorization (ABF) and quasi- implicit pressure and explicit saturation (QIMPES) decouplers and flexible generalized residual solver (FGMRES). The new general-purpose nonlinear solver, Newton trust-region dogleg Cauchy (NTRDC), is also introduced to resolve extreme nonlinearities in the models. We demonstrate the effectiveness of each method relative to the default BCGS-Newton solver. The two best cases had nearly 50 times speed-up and achieved completion of a simulation in 14 hours that never completed due to non-convergence with the default solver. We also investigate the strong scalability of each method and discuss some of the deficiencies found for Block Jacobi preconditioner using parallel domain decomposition, and node packing effects of modern processor architecture.

Preconditioner, Nonlinear, Porous media, Multiphas↗

NASA NTRS1985

Crack problems for bonded nonhomogeneous materials under antiplane shear loading

The singular nature of the crack tip stress field in a nonhomogeneous medium having a shear modulus with a discontinuous derivative was investigated. The problem is considered for the simplest possible loading and geometry, namely the antiplane shear loading of two bonded half spaces in which the crack is perpendicular to the interface. It is shown that the square-root singularity of the crack tip stress field is unaffected by the discontinuity in the derivative of the shear modulus. The problem is solved for a finite crack and extensive results are given for the stress intensity factors.

Erdogan, F.↗

NASA NTRS1985

The crack problem for bonded nonhomogeneous materials under antiplane shear loading

Erdogan, F.↗

NASA NTRS1983

Magnetic Bearings for Inertial Energy Storage

The selection of a noncontacting bearing technique with no wear out phenomena and which is vacuum compatible which is the decisive factor in selecting magnetic bearings for kinetic energy storage was investigated. Unlimited cycle life without degradation is a primary goal. Storage efficiency is a key parameter which is defined as the ratio of the energy remaining to energy stored after a fixed time interval at no load conditions. Magnetic bearings, although noncontacting, are not perfectly frictionless in that magnetic losses due to eddy currents and hysteresis can occur. Practical magnetic bearings, however, deviate from perfect symmetry and have discontinuities and asymmetric flux paths either by design or when controlled in the presence of disturbances, which cause losses. These losses can be kept smaller in the bearings than in a high power motor/generator, however, are a significant factor in selecting the magnetic bearing type.

Studer, P. A.↗

DOE OSTI2022

Factors Controlling a Synthetic Aperture Radar (SAR) Derived Root-Zone Soil Moisture Product over The Seward Peninsula of Alaska

Root-zone soil moisture exerts a fundamental control on vegetation, energy balance, and the carbon cycle in Arctic ecosystems, but it is still not well understood in vast, remote, and understudied regions of discontinuous permafrost. The root-zone soil moisture product (30 m resolution) used in this analysis was retrieved from a time-series P-Band (420–440 MHz) synthetic aperture radar (SAR) backscatter observations (August 2017 & October 2017). While similar approaches have been taken to retrieve surface (0 cm to 5 cm) soil moisture from L-Band (1.2 GHz) SAR backscatter, this is one of the first known attempts at reaching the root-zone in permafrost regions. Here, we analyze secondary factors (excluding primary factors, such as precipitation) controlling summer (August) soil moisture at depths of 6 cm, 12 cm, and 20 cm over a 4500 km2 area on the Seward Peninsula of Alaska. Using a random forest model, we quantify the impact of topography, vegetation, and meteorological factors on soil moisture distributions. In developing the random forest model, we explore a variety of feature scales (30 m, 60 m, 90 m, 120 m, 180 m, and 240 m), tune hyperparameters (the structure of individual decision trees making up the ensemble including the number and depth of trees), and perform the final feature selection using cross-validated recursive feature elimination. Results suggest that root-zone soil moisture on the Seward Peninsula is primarily controlled by vegetation at 6 cm, but deeper in the soil column topography and meteorological factors, such as predominant winter wind direction and summer insolation, play a larger role. The random forest model accounts for 40% to 60% of the variation observed (R2 = 0.44 at 6 cm, R2 = 0.52 at 12 cm, R2 = 0.58 at 20 cm). These results indicate that vegetation is the dominant control on soil moisture shallow in the soil column, but the impact of vegetation does not extend to deeper layers retrieved from P-Band SAR backscatter.

Dann, Julian↗

DOE OSTI2020

Efficient Low-Order Refined Preconditioners for High-Order Matrix-Free Continuous and Discontinuous Galerkin Methods

In this paper, we design preconditioners for the matrix-free solution of high-order continuous and discontinuous Galerkin discretizations of elliptic problems based on finite element method--spectral element method (FEM-SEM) equivalence and additive Schwarz methods. The high-order operators are applied without forming the system matrix, making use of sum factorization for efficient evaluation. The system is preconditioned using a spectrally equivalent low-order (p = 1) finite element operator discretization on a refined mesh. The low-order refined mesh is anisotropic and not shape regular in the polynomial degree of the high-order operator, requiring specialized solvers to treat the anisotropy. We make use of an element-structured, geometric multigrid V-cycle with ordered ILU(0) smoothing. The preconditioner is parallelized through an overlapping additive Schwarz method that is robust in h and p. The method is extended to interior penalty and Bassi and Rebay (BR2) discontinuous Galerkin discretizations, for which it is also robust in the size of the penalty parameter. Finally, numerical results are presented on a variety of examples, verifying the uniformity of the preconditioner.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1993

Resolution of the COBE Earth sensor anomaly

Since its launch on November 18, 1989, the Earth sensors on the Cosmic Background Explorer (COBE) have shown much greater noise than expected. The problem was traced to an error in Earth horizon acquisition-of-signal (AOS) times. Due to this error, the AOS timing correction was ignored, causing Earth sensor split-to-index (SI) angles to be incorrectly time-tagged to minor frame synchronization times. Resulting Earth sensor residuals, based on gyro-propagated fine attitude solutions, were as large as plus or minus 0.45 deg (much greater than plus or minus 0.10 deg from scanner specifications (Reference 1)). Also, discontinuities in single-frame coarse attitude pitch and roll angles (as large as 0.80 and 0.30 deg, respectively) were noted several times during each orbit. However, over the course of the mission, each Earth sensor was observed to independently and unexpectedly reset and then reactivate into a new configuration. Although the telemetered AOS timing corrections are still in error, a procedure has been developed to approximate and apply these corrections. This paper describes the approach, analysis, and results of approximating and applying AOS timing adjustments to correct Earth scanner data. Furthermore, due to the continuing degradation of COBE's gyroscopes, gyro-propagated fine attitude solutions may soon become unavailable, requiring an alternative method for attitude determination. By correcting Earth scanner AOS telemetry, as described in this paper, more accurate single-frame attitude solutions are obtained. All aforementioned pitch and roll discontinuities are removed. When proper AOS corrections are applied, the standard deviation of pitch residuals between coarse attitude and gyro-propagated fine attitude solutions decrease by a factor of 3. Also, the overall standard deviation of SI residuals from fine attitude solutions decrease by a factor of 4 (meeting sensor specifications) when AOS corrections are applied.

Sedler, J.↗

NASA NTRS2009

Health Monitoring of a Rotating Disk Using a Combined Analytical-Experimental Approach

Rotating disks undergo rigorous mechanical loading conditions that make them subject to a variety of failure mechanisms leading to structural deformities and cracking. During operation, periodic loading fluctuations and other related factors cause fractures and hidden internal cracks that can only be detected via noninvasive types of health monitoring and/or nondestructive evaluation. These evaluations go further to inspect material discontinuities and other irregularities that have grown to become critical defects that can lead to failure. Hence, the objectives of this work is to conduct a collective analytical and experimental study to present a well-rounded structural assessment of a rotating disk by means of a health monitoring approach and to appraise the capabilities of an in-house rotor spin system. The analyses utilized the finite element method to analyze the disk with and without an induced crack at different loading levels, such as rotational speeds starting at 3000 up to 10 000 rpm. A parallel experiment was conducted to spin the disk at the desired speeds in an attempt to correlate the experimental findings with the analytical results. The testing involved conducting spin experiments which, covered the rotor in both damaged and undamaged (i.e., notched and unnotched) states. Damaged disks had artificially induced through-thickness flaws represented in the web region ranging from 2.54 to 5.08 cm (1 to 2 in.) in length. This study aims to identify defects that are greater than 1.27 cm (0.5 in.), applying available means of structural health monitoring and nondestructive evaluation, and documenting failure mechanisms experienced by the rotor system under typical turbine engine operating conditions.

Abdul-Aziz, Ali↗