SEARCH · Engineering Papers

Results for “consistent scheme”

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 919 records · Page 51

All sources NASA NTRS DOE OSTI

NASA NTRS1984

Newton like: Minimal residual methods applied to transonic flow calculations

A computational technique for the solution of the full potential equation is presented. The method consists of outer and inner iterations. The outer iterate is based on a Newton like algorithm, and a preconditioned Minimal Residual method is used to seek an approximate solution of the system of linear equations arising at each inner iterate. The present iterative scheme is formulated so that the uncertainties and difficulties associated with many iterative techniques, namely the requirements of acceleration parameters and the treatment of additional boundary conditions for the intermediate variables, are eliminated. Numerical experiments based on the new method for transonic potential flows around the NACA 0012 airfoil at different Mach numbers and different angles of attack are presented, and these results are compared with those obtained by the Approximate Factorization technique. Extention to three dimensional flow calculations and application in finite element methods for fluid dynamics problems by the present method are also discussed. The Inexact Newton like method produces a smoother reduction in the residual norm, and the number of supersonic points and circulations are rapidly established as the number of iterations is increased.

Wong, Y. S.↗

NASA NTRS1985

Newton-like minimal residual methods applied to transonic flow calculations

Wong, Y. S.↗

NASA NTRS1988

Computational technique for compressible vortex flows using the integral equation solution

The steady full-potential equation is written in the form of Poisson's equation, and the solution for the velocity field is expressed in terms of an integral equation. The integral solution consists of two surface integrals and one volume integral. The solution is obtained through successive iteration cycles. Each cycle of iteration consists of two sub-cycles, an inner cycle for wake relaxation and an out cycle for the strength of the source distribution integrals representing the flow compressibility. The density gradients in the source distribution is computed by using a type-differencing scheme of the Murman-Cole type. The method is applied to delta wings and the numerical examples show that a curved shock is captured on the wing suction side beneath the leading edge vortex sheet. Recently, a modified version of the scheme was applied to rectangular wings. In this modified scheme, the surface integral terms were computed by using a bilinear distribution of vorticity on triangular vortex panels which represent the wing and its wake. The results were compared with the available experimental data and they are in good agreement.

Kandil, Osama A.↗

NASA NTRS1999

1D Variational Assimilation of TOVS/ATOVS Level 1b Data

Recently, it has been shown that the use of observations from satellite-borne microwave and infrared radiometers in data assimilation systems consistently increases forecast skill. Considerable effort has been expended over the past two decades, particularly with the TIROS Operational Vertical Sounder (TOVS), to achieve this result. The positive impact on forecast skill has resulted from improvements in quality control algorithms, systematic error correction schemes, and data assimilation systems. Despite these recent advances, there are still many issues regarding the use of satellite data in data assimilation systems that remain unresolved. For example, at several centers, much of the TOVS data are not assimilated over land as a result of the difficulty in specifying the surface emissivity. In this study, we add the microwave surface emissivity to the state vector so that it is retrieved rather than specified by a model. The LDVAR temperature and humidity information is then assimilated interactively into the Goddard Earth Observing System - Data Assimilation System (GEOS-DAS). Results of 1DVAR assimilations with both TOVS and ATOVS will be shown. We will compare the results obtained with HIRS/MSU/SSU and HIRS/AMSU-A2.

Joiner, J.↗

NASA NTRS2001

A Strategy to Assess Aerosol Direct Radiative Forcing of Climate Using Satellite Radiation Measurements

Atmospheric aerosols have a complex internal chemical composition and optical properties. Therefore it is difficult to model their impact on redistribution and absorption of solar radiation, and the consequent impact on atmospheric dynamics and climate. The use in climate models of isolated aerosol parameters retrieved from satellite data (e.g. optical thickness) may result in inconsistent calculations, if the model assumptions differ from these of the satellite retrieval schemes. Here we suggest a strategy to assess the direct impact of aerosol on the radiation budget at the top and bottom of the atmosphere using satellite and ground based measurements of the spectral solar radiation scattered by the aerosol. This method ensures consistent use of the satellite data and increases its accuracy. For Kaufman and Tanre: Strategy for aerosol direct forcing anthropogenic aerosol in the fine mode (e.g. biomass burning smoke and urban pollution) consistent use of satellite derived optical thickness can yield the aerosol impact on the spectral solar flux with accuracy an order of magnitude better than the optical thickness itself. For example, a simulated monthly average smoke optical thickness of 0.5 at 0.55 microns (forcing of 40-50 W/sq m) derived with an error of 20%, while the forcing can be measured directly with an error of only 0-2 W/sq m. Another example, the effect of large dust particles on reflection of solar flux can be derived three times better than retrievals of optical thickness. Since aerosol impacts not only the top of the atmosphere but also the surface irradiation, a combination of satellite and ground based measurements of the spectral flux, can be the most direct mechanism to evaluate the aerosol effect on climate and assimilate it in climate models. The strategy is applied to measurements from SCAR-B and the Tarfox experiments. In SCAR-B aircraft spectral data are used to derive the 24 hour radiative forcing of smoke at the top of the atmosphere of (Delta)F(sub 24hr)/(Delta)tau = - 25 +/- 5 W/sq m. Ground based data give forcing at the surface of (Delta)F(sub 24hr)/(Delta)taur = -80 +/- 5 W/sq m. In TARFOX a mixture of maritime and regional pollution aerosol resulted in a varied forcing at the top of the atmosphere, (Delta)F(sub 24hr)/(Delta)tau, between -26 W/sq 2 and -50 W/sq m depending on mixture of coarse and accumulation modes, for Angstrom exponents of 1.0 and 0.2 respectively.

Kaufman, Yoram J.↗

NASA NTRS2021

Geomagnetic Secular Variation Forecast using the NASA GEMS Ensemble Kalman filter: A Candidate SV Model for IGRF‑13

We have produced a 5-year mean secular variation (SV ) of the geomagnetic field for the period 2020–2025. We use the NASA Geomagnetic Ensemble Modeling System (GEMS), which consists of the NASA Goddard geodynamo model and ensemble Kalman filter (EnKF) with 400 ensemble members. Geomagnetic field models are used as observations for the assimilation, including gufm1 (1590–1960), CM4 (1961–2000) and CM6 (2001–2019). The forecast involves a bias correction scheme that assumes that the model bias changes on timescales much longer than the forecast period, so that they can be removed by successive forecast series. The algorithm was validated on the time period 2010-2015 by comparing with CM6 before being applied to the 2020–2025 time period. This forecast has been sub-mitted as a candidate predictive model of IGRF-13 for the period 2020–2025.

Andrew Tangborn↗

NASA NTRS2000

Does The Earth Have an Adaptive Infrared Iris?

Observations and analyses of water vapor and clouds in the tropics over the past decade suggest a different approach to radiative climate feedbacks: namely, that high clouds and high free-tropospheric relative humidity are largely tied to each other, and that the main feedback consists in changing the relative areas of cloudy/moist regions vis a vis clear/dry regions in response to the surface temperature of the cloudy/moist regions - as opposed to altering the humidity in either of the regions. This is an intrinsically 2-dimensional (horizontal and vertical) effect which does not readily enter simple 1-dimensional (vertical) radiative-convective schemes which emphasize average humidity, etc. Preliminary analyses of cloud data for the eastern part of the Western Pacific from the Japanese GMS-5(Geostationary Meteorological Satellite), are supportive of this suggestion - pointing to a 15% reduction in cloudy/moist area for a 1C increase of the sea surface temperature as measured by the cloud-weighted SST (sea surface temperature). The implication of this result is examined using a simple 2-dimensional radiative-convective model. The calculations show that such a change in the tropics would lead to a strong negative feedback in the global climate, with a feedback factor of about -1.7, which, if correct, would easily dominate the positive water vapor feedback found in current models. This new feedback mechanism, in effect, constitutes an adaptive infrared iris that opens and closes in order to control the OLR (outgoing longwave radiation) in response to changes in surface temperature in a manner similar to the way in which an eye's iris opens and closes in response to changing light levels. The climate sensitivity resulting from this thermostatic mechanism is consistent with the independent determination by Lindzen and Giannitisis (1998). Preliminary attempts to replicate observations with GCMs (General Circulation Models) suggest that models lack such a negative cloud/moist areal feedback.

Lindzen, Richard S.↗

NASA NTRS2000

Multiresolution With Super-Compact Wavelets

The solution data computed from large scale simulations are sometimes too big for main memory, for local disks, and possibly even for a remote storage disk, creating tremendous processing time as well as technical difficulties in analyzing the data. The excessive storage demands a corresponding huge penalty in I/O time, rendering time and transmission time between different computer systems. In this paper, a multiresolution scheme is proposed to compress field simulation or experimental data without much loss of important information in the representation. Originally, the wavelet based multiresolution scheme was introduced in image processing, for the purposes of data compression and feature extraction. Unlike photographic image data which has rather simple settings, computational field simulation data needs more careful treatment in applying the multiresolution technique. While the image data sits on a regular spaced grid, the simulation data usually resides on a structured curvilinear grid or unstructured grid. In addition to the irregularity in grid spacing, the other difficulty is that the solutions consist of vectors instead of scalar values. The data characteristics demand more restrictive conditions. In general, the photographic images have very little inherent smoothness with discontinuities almost everywhere. On the other hand, the numerical solutions have smoothness almost everywhere and discontinuities in local areas (shock, vortices, and shear layers). The wavelet bases should be amenable to the solution of the problem at hand and applicable to constraints such as numerical accuracy and boundary conditions. In choosing a suitable wavelet basis for simulation data among a variety of wavelet families, the supercompact wavelets designed by Beam and Warming provide one of the most effective multiresolution schemes. Supercompact multi-wavelets retain the compactness of Haar wavelets, are piecewise polynomial and orthogonal, and can have arbitrary order of approximation. The advantages of the multiresolution algorithm are that no special treatment is required at the boundaries of the interval, and that the application to functions which are only piecewise continuous (internal boundaries) can be efficiently implemented. In this presentation, Beam's supercompact wavelets are generalized to higher dimensions using multidimensional scaling and wavelet functions rather than alternating the directions as in the 1D version. As a demonstration of actual 3D data compression, supercompact wavelet transforms are applied to a 3D data set for wing tip vortex flow solutions (2.5 million grid points). It is shown that high data compression ratio can be achieved (around 50:1 ratio) in both vector and scalar data set.

Lee, Dohyung↗

DOE OSTI2025

Studying baryon acoustic oscillations using photometric redshifts from the DESI Legacy Imaging survey DR9

Context. The Dark Energy Spectroscopic Instrument (DESI) Legacy Imaging Survey DR9 (DR9 hereafter), with its extensive dataset of galaxy locations and photometric redshifts, presents an opportunity to study baryon acoustic oscillations (BAOs) in the region covered by the ongoing spectroscopic survey with DESI. Aims. We aim to investigate differences between different parts of the DR9 footprint. Furthermore, we want to measure the BAO scale for luminous red galaxies within them. Our selected redshift range of 0.6–0.8 corresponds to the bin in which a tension between DESI Y1 and eBOSS was found. Methods. We calculated the anisotropic two-point correlation function in a modified binning scheme to detect the BAOs in DR9 data. We then used template fits based on simulations to measure the BAO scale in the imaging data. Results. Our analysis reveals the expected correlation function shape in most of the footprint areas, showing a BAO scale consistent with Planck’s observations. Aside from identified mask-related data issues in the southern region of the South Galactic Cap, we find a notable variance between the different footprints. Conclusions. We find that this variance is consistent with the difference between the DESI Y1 and eBOSS data, and it supports the argument that that tension is caused by sample variance. Additionally, we also uncovered systematic biases not previously accounted for in photometric BAO studies. We emphasize the necessity of adjusting for the systematic shift in the BAO scale associated with typical photometric redshift uncertainties to ensure accurate measurements.

79 ASTRONOMY AND ASTROPHYSICS↗

NASA NTRS1972

Low G accelerometer testing

Eight different types of low-g accelerometer tests are covered on the Bell miniature electrostatically suspended accelerometer (MESA) which is known to be sensitive to less than 10 to the minus 7th power earth's gravity. These tests include a mass attracting scheme, Leitz dividing head, Wild theodolite, precision gage blocks, precision tiltmeters, Hilger Watts autocollimator, Razdow Mark 2 autocollimator, and laser interferometer measuring system. Each test is described and a comparison of the results is presented. The output of the MESA was as linear and consistent as any of the available devices were capable of measuring. Although the extent of agreement varied with the test equipment used, it can only be concluded that the indicated errors were attributable to the test equipment coupled with the environmental conditions.

Vaughan, M. S.↗

NASA NTRS1986

A refined mixed shear flexible finite element for the nonlinear analysis of laminated plates

The present study is concerned with the development of a mixed shear flexible finite element with relaxed continuity for the geometrically linear and nonlinear analysis of laminated anisotropic plates. The formulation of the element is based on a refined higher-order theory. This theory satisfies the zero transverse shear stress boundary conditions on the top and bottom faces of the plate. Shear correction coefficients are not needed. The developed element consists of 11 degrees-of-freedom per node, taking into account three displacements, two rotations, and six moment resultants. An evaluation of the element is conducted with respect to the accuracy obtained in the bending of laminated anistropic rectangular plates with different lamination schemes, loadings, and boundary conditions.

Putcha, N. S.↗

NASA NTRS1990

A predictive theoretical model for electron tunneling pathways in proteins

A practical method is presented for calculating the dependence of electron transfer rates on details of the protein medium intervening between donor and acceptor. The method takes proper account of the relative energetics and mutual interactions of the donor, acceptor, and peptide groups. It also provides a quantitative search scheme for determining the important tunneling pathways (specific sequences of localized bonding and antibonding orbitals of the protein which dominate the donor-acceptor electronic coupling) in native and tailored proteins, a tool for designing new proteins with prescribed electron transfer rates, and a consistent description of observed electron transfer rates in existing redox labeled metalloproteins and small molecule model compounds.

Onuchic, Jose Nelson↗

NASA NTRS2001

Assimilation of Cloud- and Land-Affected TOVS/ATOVS Level 1B Radiances in DAO's Next Generation Finite-Volume Data Assimilation System

The Physical-space/Finite-volume Data Assimilation System (fvDAS) is the next generation global atmospheric data assimilation system in development at the Data Assimilation Office (DAO) at NASA's Goddard Space Flight Center. It is based on a new finite-volume general circulation model jointly developed by NASA and NCAR, and on the Physical-Space Statistical Analysis System (PSAS) developed at the DAO. In this talk we will focus on the assimilation of data from the (Advanced) TIROS Operational Vertical Sounder (ATOVS), with emphasis on the impact of cloud- and land-affected level 1B radiances. Recently, it has been shown that the use of observations from satellite-borne microwave and infrared radiometers in data assimilation systems consistently increases forecast skill. Considerable effort has been expended over the past two decades, particularly with the (Advanced) TIROS Operational Vertical Sounder (ATOVS), to achieve this result. The positive impact on forecast skill has resulted from improvements in quality control algorithms, systematic error correction schemes, and more sophisticated data assimilation algorithms. Despite these advances, there are still many issues regarding the use of satellite data in data assimilation systems that remain unresolved. In particular, most operational centers still do not assimilate cloud- and land-affected TOVS data. In this study, we evaluate the impact of assimilating cloud-and land-affected TOVS/ATOVS level 1B data in DAO's next generation fvDAS, using a 1D variational scheme. We will discuss the impact of these data on both tropospheric and stratospheric forecasts, as well as on the general aspects of the earth climate system.

Joiner, J.↗

NASA NTRS2022

Unique Observational Constraints on the Seasonal and Longitudinal Variability of the Earth’s Planetary Albedo and Cloud Distribution Inferred From EPIC Measurements

Thorough comparison to observations is key to developing a credible climate model forecasting capability. Deep Space Climate Observatory (DSCOVR) measurements of Earth’s reflected solar and emitted thermal radiation provide a unique observational perspective that permits a more reliable model/data comparison than is possible with the otherwise available satellite data. The uniqueness is in the DSCOVR satellite’s viewing geometry, which enables continuous viewing of the Earth’s sunlit hemisphere from its Lissajous orbit around the Lagrangian L1 point. The key instrument is the Earth Polychromatic Imaging Camera (EPIC), which views the Earth’s sunlit hemisphere with 1024-by-1024-pixel imagery in 10 narrow spectral bands from 317 to 780 nm, acquiring up to 22 high spatial resolution images per day. The additional feature is that the frequency of EPIC image acquisition is nearly identical to that of the climate GCM data generation scheme where climate data for the entire globe are ‘instantaneously’ calculated at 1-h radiation time-step intervals. Implementation of the SHS (Sunlit Hemisphere Sampling) EPIC-view geometry for the in-line GCM output data sampling establishes a precise self-consistency in the space-time data sampling between EPIC observational and GCM output data generation and sampling. The remaining problem is that the GCM generated data are radiative fluxes, while the EPIC measurements are backscatter-dependent radiances. Radiance to flux conversion is a complex problem with no simple way to convert GCM radiative fluxes into spectral radiances. The more expedient approach is to convert the EPIC spectral radiances into broadband radiances by MODIS/CERES-based regression relationships and then into solar radiative fluxes using the CERES angular distribution models. Averaging over the sunlit hemisphere suppresses the meteorological weather noise, but preserves the intra-seasonal larger scale variability. Longitudinal slicing by the Earth’s rotation permits a self-consistent model/data comparison of the longitudinal model/data differences in the variability of the reflected solar radiation. Ancillary EPIC Composite data provide additional cloud property information for climate model diagnostics. Comparison of EPIC-derived seasonal and longitudinal variability of the Earth’s planetary albedo with the GISS ModelE2 results shows systematic overestimate of cloud reflectivity over the Pacific Ocean with corresponding underestimates over continental land areas.

DSCOVR↗

DOE OSTI2020

Entanglement properties of disordered quantum spin chains with long-range antiferromagnetic interactions

Entanglement measures are useful tools in characterizing otherwise unknown quantum phases and indicating transitions between them. Here we examine the concurrence and entanglement entropy in quantum spin chains with random long-range couplings, spatially decaying with a power-law exponent $\textit{α}$. Using the strong disorder renormalization group (SDRG) technique, we find by analytical solution of the master equation a strong disorder fixed point, characterized by a fixed point distribution of the couplings with a finite dynamical exponent, which describes the system consistently in the regime $α > \frac{1}{2}$ . A numerical implementation of the SDRG method yields a power-law spatial decay of the average concurrence, which is also confirmed by exact numerical diagonalization. However, we find that the lowest-order SDRG approach is not sufficient to obtain the typical value of the concurrence. We therefore implement a correction scheme which allows us to obtain the leading-order corrections to the random singlet state. This approach yields a power-law spatial decay of the typical value of the concurrence, which we derive both by a numerical implementation of the corrections and by analytics. Next, using numerical SDRG, the entanglement entropy (EE) is found to be logarithmically enhanced for all $\textit{α}$, corresponding to a critical behavior with an effective central charge $\textit{c}$ = ln(2), independent of $\textit{α}$. This is confirmed by an analytical derivation. Using numerical exact diagonalization (ED), we confirm the logarithmic enhancement of the EE and a weak dependence on $\textit{α}$. For a wide range of partition size $\textit{l}$, the EE fits a critical behavior with a central charge close to $\textit{c}$ = 1, which is the same as for the clean Haldane-Shastry model with a power-law-decaying interaction with $\textit{α}$ = 2. Only for small $l \ll L$, in a range which increases with the number of spins N, we find deviations which are rather consistent with the strong disorder fixed point central charge $\textit{c}$ = ln(2). Furthermore, we find using ED that the concurrence shows power-law decay, albeit with smaller power exponents than obtained by SDRG. Finally, we also present results obtained with DMRG and find agreement with ED for sufficiently small $\textit{α}$ < 2, whereas for larger $\textit{α}$ DMRG tends to underestimate the entanglement entropy and finds a faster decaying concurrence.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2024

X-ray spectroscopic investigation of crystal fields in Ce 2 Rh 1 − x Ir x In 8 heavy fermions

The higher dimensionality in the crystal fields of the Ce 2 M In 8 ( M = Rh , Ir ) compounds and its interplay with hybridization and disorder are key ingredients to understand the complex phase diagrams by this family, which have been explored extensively by macroscopic techniques. Here, we present an investigation of the crystal-electric field schemes of Ce 2 Rh 1 − x Ir x In 8 using x-ray absorption spectroscopy. Our full multiplet calculations for the 4 f 1 configuration of Ce 3 + to describe the temperature-dependent linear dichroism in Ce 2 M In 8 are consistent with a Γ 7 1 = 1 − α 2 · | ∓ 3 2 〉 − | α | · | ± 5 2 〉 ground state containing a predominant | ± 3 / 2 〉 contribution that increases further with x . This enhancement is believed to favor superconductivity in Ce-based heavy fermion materials, observed in previous results in the Ce M In 5 family. Our recent observations shed light on the unexpected emergence of the ambient-pressure superconducting dome in the center of the composition phase diagram and its subsequent suppression on the Ir-rich side due to the early onset of fluctuations associated with the structurally more disordered state, inferred from previous neutron magnetic diffraction experiments. Published by the American Physical Society 2024

Christovam, D. S. (ORCID:0000000250698107)↗

NASA NTRS1977

Automated array assembly task, phase 1

An assessment of state-of-the-art technologies that are applicable to silicon solar cell and solar cell module fabrication is provided. The assessment consists of a technical feasibility evaluation and a cost projection for high-volume production of silicon solar cell modules. The cost projection was approached from two directions; a design-to-cost analysis assigned cost goals to each major process element in the fabrication scheme, and a cost analysis built up projected costs for alternate technologies for each process element. A technical evaluation was used in combination with the cost analysis to identify a baseline low cost process. A novel approach to metal pattern design based on minimum power loss was developed. These design equations were used as a tool in the evaluation of metallization technologies.

Carbajal, B. G.↗

NASA NTRS1980

Development and validation of a hybrid-computer simulator for a transonic cryogenic wind tunnel

A study was undertaken to model the cryogenic wind tunnel process, to validate the model by the use of experimental data from the Langley 0.3 Meter Transonic Cryogenic Tunnel, and to construct an interactive simulator of the cryogenic tunnel using the validated model. Additionally, this model was used for designing closed loop feedback control laws for regulation of temperature and pressure in the 0.3 meter cryogenic tunnel. The global mathematical model of the cryogenic tunnel that were developed consists of coupled, nonlinear differential governing equations based on an energy state concept of the physical cryogenic phenomena. Process equations and comparisons between actual tunnel responses and computer simulation predictions were examined. Also included are the control laws and simulator responses obtained using the feedback schemes for closed loop control of temperature and pressure were also included.

Thibodeaux, J. J.↗