Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximate inverses”

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 271 records · Page 15

Microstructural characterization of a star-linear polymer blend under shear flow by using rheo-SANS

We present an investigation into the dynamic relaxation mechanisms of a polybutadiene blend composed of a four-arm star (10 wt. %) and a linear polymer matrix in the presence of an applied shear flow. Our focus was the response of the star polymer, which cannot be unambiguously assessed via linear viscoelastic measurements since the signature of the star polymer can barely be detected due to the dominant contribution of the linear matrix. By utilizing small-angle neutron scattering (SANS) coupled with a Couette shear device and a deuterated matrix polymer, we investigated the dynamics of the minority star component of the blend. Our results confirm that the stars deform anisotropically with increasing shear rate. We have compared the SANS data with predictions from the well-established scattering adaptation of the state-of-the-art tube model for entangled linear polymer melts undergoing shear, i.e., Graham, Likhtman, Milner, and McLeish (GLaMM) approach, appropriately modified following earlier studies in order to apply to the star. This modified model, GLaMM-R, includes the physics necessary to understand stress relaxation in both the linear and nonlinear flow regimes, i.e., contour length fluctuations, constraint release, convective constraint release, and chain retraction. The full scattering signal is due to the minority star component and, although the contribution of the linear chains is hidden from the neutron scattering, they still influence the star polymer molecular dynamics, with the applied shear rate ranging from approximately 8 to 24 s –1 , below the inverse relaxation time of the linear component. This study provides another confirmation that the combination of rheology and neutron scattering is an indispensable tool for investigating the nonlinear dynamics of complex polymeric systems.

36 MATERIALS SCIENCE↗

Tracking performance of unbalanced QPSK demodulators. II - Biphase Costas loop with active arm filters

In a Costas loop study for biphase modulation conducted by Simon and Lindsey (1977), it was demonstrated that considerable improvement in tracking performance could be obtained by employing active arm filters of the integrate-and-dump type as opposed to passive arm filters. An investigation is conducted concerning the possibility to obtain a similar performance improvement for an unbalanced quadriphase-shift-keying (QPSK) modulation. It is found that the biphase Costas loop can be used as an efficient demodulator of QPSK in cases in which the ratio of data rates is of the same order of magnitude as the inverse of the power ratio. These cases involve approximately equal signal energies in the two channels.

Simon, M. K.↗

The continuous similarity model of bulk soil-water evaporation

The continuous similarity model of evaporation is described. In it, evaporation is conceptualized as a two stage process. For an initially moist soil, evaporation is first climate limited, but later it becomes soil limited. During the latter stage, the evaporation rate is termed evaporability, and mathematically it is inversely proportional to the evaporation deficit. A functional approximation of the moisture distribution within the soil column is also included in the model. The model was tested using data from four experiments conducted near Phoenix, Arizona; and there was excellent agreement between the simulated and observed evaporation. The model also predicted the time of transition to the soil limited stage reasonably well. For one of the experiments, a third stage of evaporation, when vapor diffusion predominates, was observed. The occurrence of this stage was related to the decrease in moisture at the surface of the soil. The continuous similarity model does not account for vapor flow. The results show that climate, through the potential evaporation rate, has a strong influence on the time of transition to the soil limited stage. After this transition, however, bulk evaporation is independent of climate until the effects of vapor flow within the soil predominate.

Clapp, R. B.↗

Recent improvements in efficiency, accuracy, and convergence for implicit approximate factorization algorithms

In 1977 and 1978, general purpose centrally space differenced implicit finite difference codes in two and three dimensions have been introduced. These codes, now called ARC2D and ARC3D, can run either in inviscid or viscous mode for steady or unsteady flow. Since the introduction of the ARC2D and ARC3D codes, overall computational efficiency could be improved by making use of a number of algorithmic changes. These changes are related to the use of a spatially varying time step, the use of a sequence of mesh refinements to establish approximate solutions, implementation of various ways to reduce inversion work, improved numerical dissipation terms, and more implicit treatment of terms. The present investigation has the objective to describe the considered improvements and to quantify advantages and disadvantages. It is found that using established and simple procedures, a computer code can be maintained which is competitive with specialized codes.

Pulliam, T. H.↗

Rapid delivery of small payloads to Mars

To sustain a continued human presence on the surface of Mars or to achieve critical mission objectives the rapid delivery of small payloads may be necessary. Consideration is given to the physics of delivering payloads of about 10 kg over a nominal distance of 1 AU (149.5 x 10 to the 6th km) in 10 days. It is proposed that the most effective method of delivery is the use of a laser lightsail vehicle. The lightsail would be accelerated by a high-power laser from the vicinity of earth to velocities of approximately 174 km/s. Coasting at this velocity for most of the trip the vehicle would be decelerated by a similar laser upon arrival at Mars. The continuous laser power required is calculated to be about 47-billion W and the total duty time of both lasers, determined by optimization of the total energy, is calculated to be 3.9 hours. The laser power scales approximately linearly with payload mass and as the inverse square of the trip time.

Meyer, T. R.↗

On the development of efficient algorithms for three dimensional fluid flow

The difficulties of constructing efficient algorithms for three-dimensional flow are discussed. Reasonable candidates are analyzed and tested, and most are found to have obvious shortcomings. Yet, there is promise that an efficient class of algorithms exist between the severely time-step sized-limited explicit or approximately factored algorithms and the computationally intensive direct inversion of large sparse matrices by Gaussian elimination.

Maccormack, R. W.↗

Estimation on nonlinear damping in second order distributed parameter systems

An approximation and convergence theory for the identification of nonlinear damping in abstract wave equations is developed. It is assumed that the unknown dissipation mechanism to be identified can be described by a maximal monotone operator acting on the generalized velocity. The stiffness is assumed to be linear and symmetric. Functional analytic techniques are used to establish that solutions to a sequence of finite dimensional (Galerkin) approximating identification problems in some sense approximate a solution to the original infinite dimensional inverse problem.

Banks, H. T.↗

Kelvin-Helmholtz instability in an Alfven resonant layer of a solar coronal loop

A Kelvin-Helmholtz instability has been identified numerically on an azimuthally symmetric Alfven resonant layer in an axially bounded, straight cylindrical coronal loop. The set of equations is solved numerically as an initial value problem. The linear growth rate of this instability is shown to be approximately proportional to the Alfven driving amplitude and inversely proportional to the width of the Alfven resonant layer. It is also shown that the linear growth rate increases linearly with m - 1 up to a certain m, reaches its maximum value for the mode whose half wavelength is comparable to the Alfven resonant layer width, and decreases at higher azimuthal mode number.

Uchimoto, E.↗

An analysis for high Reynolds number inviscid/viscid interactions in cascades

An efficient steady analysis for predicting strong inviscid/viscid interaction phenomena such as viscous-layer separation, shock/boundary-layer interaction, and trailing-edge/near-wake interaction in turbomachinery blade passages is needed as part of a comprehensive analytical blade design prediction system. Such an analysis is described. It uses an inviscid/viscid interaction approach, in which the flow in the outer inviscid region is assumed to be potential, and that in the inner or viscous-layer region is governed by Prandtl's equations. The inviscid solution is determined using an implicit, least-squares, finite-difference approximation, the viscous-layer solution using an inverse, finite-difference, space-marching method which is applied along the blade surfaces and wake streamlines. The inviscid and viscid solutions are coupled using a semi-inverse global iteration procedure, which permits the prediction of boundary-layer separation and other strong-interaction phenomena. Results are presented for three cascades, with a range of inlet flow conditions considered for one of them, including conditions leading to large-scale flow separations. Comparisons with Navier-Stokes solutions and experimental data are also given.

Barnett, Mark↗

Retrieval of Surface Directional Reflectance Properties Using Ground Level Multiangle Measurements

Knowledge of the directional reflectance properties of natural surfaces such as soils and vegetation canopies is essential for classification studies and canopy model inversion. Atmospheric correction schemes, using various levels of approximation, are described to retrieve surface bidirectional reflectance factors (BRFs) and directional-hemispherical reftectances (albedos)from multiangle radiance measurements taken at ground level. The retrieval schemes are tested on simulated data incorporating realistic surface BRFs and atmospheric models containing aerosols. Sensitivity of the atmospherically corrected BRFs and associated directional-hemispherical reflectances to various aerosol properties and the sun-view geometry is illustrated. A measurement strategy for obtaining highly accurate surface reflectance properties also is examined in the context of instrument radiametric calibration, knowledge of the atmospheric properties, and sun-view angular coverage.

Martonchik, John V.↗

An Inverse Interpolation Method Utilizing In-Flight Strain Measurements for Determining Loads and Structural Response of Aerospace Vehicles

An important and challenging technology aimed at the next generation of aerospace vehicles is that of structural health monitoring. The key problem is to determine accurately, reliably, and in real time the applied loads, stresses, and displacements experienced in flight, with such data establishing an information database for structural health monitoring. The present effort is aimed at developing a finite element-based methodology involving an inverse formulation that employs measured surface strains to recover the applied loads, stresses, and displacements in an aerospace vehicle in real time. The computational procedure uses a standard finite element model (i.e., "direct analysis") of a given airframe, with the subsequent application of the inverse interpolation approach. The inverse interpolation formulation is based on a parametric approximation of the loading and is further constructed through a least-squares minimization of calculated and measured strains. This procedure results in the governing system of linear algebraic equations, providing the unknown coefficients that accurately define the load approximation. Numerical simulations are carried out for problems involving various levels of structural approximation. These include plate-loading examples and an aircraft wing box. Accuracy and computational efficiency of the proposed method are discussed in detail. The experimental validation of the methodology by way of structural testing of an aircraft wing is also discussed.

Shkarayev, S.↗

Retrieval of Surface Directional Reflectance and Hemispherical Albedo Using Multi-angle Measurements, I. Ground Level Observations

Knowledge of the directional reflectance properties of natural surfaces such as soils and vegetation canopies is essential for classification studies and canopy model inversions. Atmospheric correction schemes, using various levels of approximation, are described to retrieve surface bidirectional reflectance factors (BRFs) and directional-hemispherical reflectances (albedos) from multi-angle radiance measurements taken at ground level.

factors (BRFs) radiometric calibration sun-view an↗

A Mesoscale Analysis of Column-Integrated Aerosol Properties in Northern India During the TIGERZ 2008 Pre-Monsoon Period and a Comparison to MODIS Retrievals

The Indo-Gangetic Plain (IGP) of the northern Indian subcontinent produces anthropogenic pollution from urban, industrial and rural combustion sources nearly continuously and is affected by convection-induced winds driving desert and alluvial dust into the atmosphere during the premonsoon period. Within the IGP, the NASA Aerosol Robotic Network (AERONET) project initiated the TIGERZ measurement campaign in May 2008 with an intensive operational period from May 1 to June 23, 2008. Mesoscale spatial variability of aerosol optical depth (AOD, tau) measurements at 500mn was assessed at sites around Kanpur, India, with averages ranging from 0.31 to 0.89 for spatial variability study (SVS) deployments. Sites located downwind from the city of Kanpur indicated slightly higher average aerosol optical depth (delta Tau(sub 500)=0.03-0.09). In addition, SVS AOD area-averages were compared to the long-tenn Kanpur AERONET site data: Four SVS area-averages were within +/- 1 cr of the climatological mean of the Kanpur site, while one SVS was within 2sigma below climatology. For a SVS case using AERONET inversions, the 440-870mn Angstrom exponent of approximately 0.38, the 440-870mn absorption Angstrom exponent (AAE) of 1.15-1.53, and the sphericity parameter near zero suggested the occurrence of large, strongly absorbing, non-spherical aerosols over Kanpur (e.g., mixed black carbon and dust) as well as stronger absorption downwind of Kanpur. Furthermore, the 3km and lOkm Terra and Aqua MODIS C005 aerosol retrieval algorithms at tau(sub 550) were compared to the TIGERZ data set. Although MODIS retrievals at higher quality levels were comparable to the MODIS retrieval uncertainty, the total number of MODIS matchups (N) were reduced with subsequent quality levels (N=25, QA>=0; N=9,QA>=l; N=6, QA>=2; N=1, QA=3) over Kanpur during the premonsoon primarily due to the semi-bright surface, complex aerosol mixture and cloud-contaminated pixels. The TIGERZ 2008 data set provided a unique opportunity to measure the spatial and temporal variations of aerosol loading in the IGP. The strong aerosol absorption derived from ground-based sun/sky radiometer measurements suggested the presence of a predominately black carbon and dust mixture during the pre-monsoon period. Consistent with the elevated heat-pump hypothesis, these absorbing aerosols found across Kanpur and the greater IGP region during the pre-monsoon period likely induced regional atmospheric warming, which lead to a more rapid advance of the southwest Asian monsoon and above normal precipitation over northern India in June 2008.

Giles, D. M.↗

On the Solution of ℓ 0 -Constrained Sparse Inverse Covariance Estimation Problems

The sparse inverse covariance matrix is used to model conditional dependencies between variables in a graphical model to fit a multivariate Gaussian distribution. Estimating the matrix from data are well known to be computationally expensive for large-scale problems. Sparsity is employed to handle noise in the data and to promote interpretability of a learning model. Although the use of a convex ℓ 1 regularizer to encourage sparsity is common practice, the combinatorial ℓ 0 penalty often has more favorable statistical properties. In this paper, we directly constrain sparsity by specifying a maximally allowable number of nonzeros, in other words, by imposing an ℓ 0 constraint. Here, we introduce an efficient approximate Newton algorithm using warm starts for solving the nonconvex ℓ 0 -constrained inverse covariance learning problem. Numerical experiments on standard data sets show that the performance of the proposed algorithm is competitive with state-of-the-art methods.

$\ell_0$-Constrained↗

Improved Representation of the Global Dust Cycle Using Observational Constraints on Dust Properties and Abundance

Even though desert dust is the most abundant aerosol by mass in Earth’s atmosphere, atmospheric models 25 struggle to accurately represent its spatial and temporal distribution. These model errors are partially caused byfundamental difficulties in simulating dust emission in coarse-resolution models and in accurately representing dust microphysical properties. Here we mitigate these problems by developing a new methodology that yields an improved representation of the global dust cycle. We present an analytical framework that uses inverse modeling to integrate an ensemble of global model simulations with observational constraints on the dust size distribution, 30 extinction efficiency, and regional dust aerosol optical depth. We then compare the inverse model results against independent measurements of dust surface concentration and deposition flux and find that errors are reduced by approximately a factor of two relative to current model simulations of the Northern Hemisphere dust cycle. The inverse model results show smaller improvements in the less dusty Southern Hemisphere, most likely because both the model simulations and the observational constraints used in the inverse model are less accurate. On a global 35 basis, we find that the emission flux of dust with geometric diameter up to 20 μm (PM20) is approximately 5,000 Tg/year, which is greater than most models account for. This larger PM20 dust flux is needed to match observational constraints showing a large atmospheric loading of coarse dust. We obtain gridded data sets of dust emission, vertically integrated loading, dust aerosol optical depth, (surface) concentration, and wet and dry deposition fluxes that are resolved by season and particle size. As our results indicate that this data set is more accurate than current 40 model simulations and the MERRA-2 dust reanalysis product, it can be used to improve quantifications of dust impacts on the Earth system.

aerosols, dust↗

Improved representation of the global dust cycle using observational constraints on dust properties and abundance

Even though desert dust is the most abundant aerosol by mass in Earth's atmosphere, atmospheric models struggle to accurately represent its spatial and temporal distribution. These model errors are partially caused by fundamental difficulties in simulating dust emission in coarse-resolution models and in accurately representing dust microphysical properties. Here we mitigate these problems by developing a new methodology that yields an improved representation of the global dust cycle. We present an analytical framework that uses inverse modeling to integrate an ensemble of global model simulations with observational constraints on the dust size distribution, extinction efficiency, and regional dust aerosol optical depth. We then compare the inverse model results against independent measurements of dust surface concentration and deposition flux and find that errors are reduced by approximately a factor of 2 relative to current model simulations of the Northern Hemisphere dust cycle. The inverse model results show smaller improvements in the less dusty Southern Hemisphere, most likely because both the model simulations and the observational constraints used in the inverse model are less accurate. On a global basis, we find that the emission flux of dust with a geometric diameter up to 20 µm (PM20) is approximately 5000 Tg yr−1, which is greater than most models account for. This larger PM20 dust flux is needed to match observational constraints showing a large atmospheric loading of coarse dust. We obtain gridded datasets of dust emission, vertically integrated loading, dust aerosol optical depth, (surface) concentration, and wet and dry deposition fluxes that are resolved by season and particle size. As our results indicate that this dataset is more accurate than current model simulations and the MERRA-2 dust reanalysis product, it can be used to improve quantifications of dust impacts on the Earth system.

Dust cycle↗

What Controls the Entrainment Rate of Dry Buoyant Thermals with Varying Initial Aspect Ratio?

Here, this study uses theory and numerical simulations to analyze the nondimensional spreading rate a (change in radius with height) of buoyant thermals as they rise and entrain surrounding environmental fluid. A focus is on how a varies with initial thermal aspect ratio A r , defined as height divided by width of the initial buoyancy perturbation. An analytic equation for thermal ascent rate w t that depends on a is derived from the thermal-volume-averaged momentum bud get equation. The thermal top height when w t is maximum, defining a critical height z c , is inversely proportional to a. The height z c also corresponds to the thermal top height when buoyant fluid along the thermal’s vertical axis is fully replaced by entrained nonbuoyant environmental fluid rising from below the thermal. The time scale for this process is controlled by the vertical velocity of parcels rising upward through the thermal’s core. This parcel vertical velocity is approximated from Hill’s analytic spherical vortex, yielding an analytic inverse relation between a and A r . Physically, this a–A r relation is connected to changes in circulation as A r is modified. Numerical simulations of thermals with A r varied from 0.5 to 2 give a values close to the analytic theoretical relation, with a factor of ~3 decrease in a as A r is increased from 0.5 to 2. The theory also explains why a of initially spherical thermals from past laboratory and modeling studies is about 0.15. Overall, this study provides a theoretical underpinning for understanding the entrainment behavior of thermals, relevant to buoyantly driven atmospheric flows.

54 ENVIRONMENTAL SCIENCES↗

An approximation theory for the identification of nonlinear distributed parameter systems

An abstract approximation framework for the identification of nonlinear distributed parameter systems is developed. Inverse problems for nonlinear systems governed by strongly maximal monotone operators (satisfying a mild continuous dependence condition with respect to the unknown parameters to be identified) are treated. Convergence of Galerkin approximations and the corresponding solutions of finite dimensional approximating identification problems to a solution of the original finite dimensional identification problem is demonstrated using the theory of nonlinear evolution systems and a nonlinear analog of the Trotter-Kato approximation result for semigroups of bounded linear operators. The nonlinear theory developed here is shown to subsume an existing linear theory as a special case. It is also shown to be applicable to a broad class of nonlinear elliptic operators and the corresponding nonlinear parabolic partial differential equations to which they lead. An application of the theory to a quasilinear model for heat conduction or mass transfer is discussed.

Banks, H. T.↗