Search NASA⌕ Search

SEARCH · Search NASA

Results for “Approximation methods”

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 775 records · Page 43

Cloud Edge Properties Measured by the ARM Shortwave Spectrometer over Ocean and Land

We use the spectrally-invariant method to study the variability of cloud optical thickness () and droplet effective radius (reff) in transition zones between the cloudy and clear sky columns observed by Shortwave Array Spectroradiometer-Zenith (SASZe) at the Southern Great Plains Central Facility site (SGP C1) and during the Marine ARM GPCI Investigation of Clouds (MAGIC) field campaign. The spectrally-invariant method approximates the spectra in the transition zone as a linear combination of definitely clear and definitely cloudy spectra. The slope and intercept of the linear relations characterize and reff in the transition region, respectively. The radiative transfer model simulations show that that (i) the slope of the visible (VIS) band is positively correlated with while (ii) the intercept of the near-infrared (NIR) band has a high negative correlation with reff. We have analyzed 22 cloud edge cases from the SGP and MAGIC and found that from cloud to clear in the transitions (a) the slopes of the VIS band decrease, indicating the decrease of towards cloud edges; (b) the intercepts of the NIR band show a much more significant increase at the SGP than from the MAGIC. The results from observed cases suggest that while decreases for all cases, the decrease in reff is much more significant for cloud over land at the SGP site compared to the ocean counterpart during the MAGIC campaign.

Yang, Weidong↗

Plasma-environment effects on K lines of astrophysical interest III. IPs, K thresholds, radiative rates, and Auger widths in Fe ix – Fe xvi

Aims. In the context of black-hole accretion disks, we aim to compute the plasma-environment e ects on the atomic parameters used to model the decay of K-vacancy states in moderately charged iron ions, namely Fe ix – Fe xvi. Methods. We used the fully relativistic multiconfiguration Dirac–Fock method approximating the plasma electron–nucleus and electron–electron screenings with a time-averaged Debye–Hückel potential. Results. We report modified ionization potentials, K-threshold energies, wavelengths, radiative emission rates, and Auger widths for plasmas characterized by electron temperatures and densities in the ranges 10^(5)-10^(7) K and 10^(18)-10^(22)/cu. cm. Conclusions. This study confirms that the high-resolution X-ray spectrometers onboard the future XRISM and Athena space missions will be capable of detecting the lowering of the K edges of these ions due to the extreme plasma conditions occurring in accretion disks around compact objects.

J. Deprince↗

Plasma environment effects on K lines of astrophysical interest IV. IPs, K thresholds, radiative rates, and Auger widths in Fe ii–Fe viii

Aims. Within the framework of compact-object accretion disks, we calculate plasma environment e ects on the atomic structure and decay parameters used in the modeling of K lines in lowly charged iron ions, namely Fe ii–Fe viii. Methods. For this study, we used the fully relativistic multiconfiguration Dirac–Fock method approximating the plasma electron–nucleus and electron-electron screenings with a time-averaged Debye-Hückel potential. Results. We report modified ionization potentials, K-threshold energies, wavelengths, radiative emission rates, and Auger widths for plasmas characterized by electron temperatures and densities in the ranges 10(exp 5)-10(exp 7) K and 10(exp 18)-10(exp 22) per cu.cm. In addition, we propose two universal fitting formulae to predict the IP and K-threshold lowerings in any elemental ion. Conclusions. We conclude that the high-resolution X-ray spectrometers onboard the future XRISM and ATHENA space missions will be able to detect the lowering of the K edges of these Fe ions due to the extreme plasma conditions occurring in the accretion disks around compact objects.

J. Deprince↗

Progress in Time-Series Soil Moisture Retrieval Using L- and S-Band Radar Backscatter

L- and S-band observations from NASA's Passive/Active L/S band (PALS) sensor from the SMEX02 campaign were used to estimate soil moisture. The retrieval process is based on the “alpha approximation” method. This method utilizes a time-series of normalized radar backscatter measurements as well as ancillary information to estimate soil moisture over the Walnut Creek watershed. The resulting retrieved soil moistures are compared to in-situ soil moisture measurements at multiple test sites within the watershed. The calculations show reasonable results for both L- and S-band and provide further insight into the use of L- and S-bands for the upcoming NASA/ISRO mission.

Soil moisture↗

4pPA1 - Turbulence Effects on Shaped Booms: Central Composite Design of Modeled Atmospheric Turbulence Parameters for Sonic Boom Propagation

Propagation of sonic booms through turbulence reduces mean sonic boom perception metric levels and also causes considerable variability. NASA’s PCBoom suite of sonic boom acoustic propagation modules includes an approximate method for accounting for the effects of turbulence on traditional N-wave sonic booms. The current implementation is ineffective for shaped sonic booms or low-booms, and it also has limited values for turbulence and ambient input parameters. NASA’s future X-59 low-boom community noise surveys require an accurate estimate of the effects of turbulence in regions across the USA, so the module must be improved. This work presents the methods of selecting which ambient and turbulence parameters should be included in an improved PCBoom turbulence module. Turbulence and ambient data were collected from two atmospheric model databases, the Climate Forecast System Version 2 and European Centre for Medium-Range Weather Forecast Reanalysis Version 5 (ERA5), hourly from 7 AM to 7 PM local time for 10 years at 19 locations across the USA. A fully-factorial propagation analysis using these parameters would be exceedingly computationally expensive. Instead, a central composite design was chosen resulting in 45 combinations of ambient and turbulence parameters. These 45 cases effectively sample the space balancing computational burden.

turbulence↗

Real-time Ridesharing for Transportation Hubs with Demand and Supply Uncertainty

Transportation hubs in major cities generate a significant amount of trips by taxis and for-hail vehicles (FHV), with many of the trips sharing similar destinations. This suggests promising opportunities to leverage the collective travel needs with dedicated ridesharing solutions to reduce the externalities of excessive traffic at transportation hubs. In this study, we develop a novel dynamic ridesharing approach to serve trips from the transportation hub by considering (1) demand (new passengers) and supply (newly available vehicles) in the near future and (2) the uncertainty of future predictions. Our approach consists of two stages. In the first stage, we develop a data structure called hub mobility tree to generate potential combinations of shareable trips as candidate schedules efficiently. Then the generated schedules are used in the second stage to formulate the stochastic hub-based ridesharing problem (SHRP), which is a stochastic integer programming problem with the objective to maximize the total expected ridesharing profit over time. Due to the prohibitive number of shareable trips, we then approximately solve SHRP by the sample average approximate method (SAA), and a dual Lagrangian technique is implemented to further improve the scalability of the solution approach. We demonstrate the performances of the proposed method by simulating the ridesharing service at JFK airport using NYC taxi and FHV data. The results indicate that the proposed method outperforms the myopic ridesharing (maximize profit for a single time step) and the rolling horizon method with point estimation of future demand and supply.

dynamic ridesharing↗

Description of spin–orbit coupling in excited states with two-component methods based on approximate coupled-cluster theory

A generalization of the approximated coupled-cluster singles and doubles method and the algebraic diagrammatic construction scheme up to second order to two-component spinors obtained from a relativistic Hartree–Fock calculation is reported. Computational results for zero-field splittings of atoms and monoatomic cations, triplet lifetimes of two organic molecules, and the spin-forbidden part of the UV/Vis absorption spectrum of tris(ethylenediamine)cobalt(III) are presented.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Unsteady Stokes flow near boundaries: the point-particle approximation and the method of reflections

Problems of particle dynamics involving unsteady Stokes flows in confined geometries are typically harder to solve than their steady counterparts. Approximation techniques are often the only resort. Felderhof (see e.g. J. Phys. Chem. B, vol. 109 (45), 2005, pp. 21406–21412; J. Fluid Mech. , vol. 637, 2009, pp. 285–303) has developed a point-particle approximation framework to solve such problems, especially in the context of Brownian motion. Despite excellent agreement with past experiments, this framework produces unsteady drag coefficients that depend on particle density. This is inconsistent, since the problem can be formulated mathematically without any reference to the particle’s density. We address this inconsistency in our work. Upon implementing our modifications, the framework passes consistency checks that it previously failed. Further, it is not obvious that such an approximation should work for short-time-scale motion. We investigate its validity by deriving it from a general formalism based on integral equations through a series of systematic approximations. We also compare results from the point-particle framework against a calculation performed using the method of reflections, for the specific case of a sphere near a full-slip plane boundary. We find from our analysis that the reasons for the success of the point-particle approximation are subtle and have to do with the nature of the unsteady Oseen tensor. Finally, we provide numerical predictions for Brownian motion near a full-slip and a no-slip plane wall based on the point-particle approximation as used by Felderhof, our modified point-particle approximation and the method of reflections. We show that our modifications to Felderhof’s framework would become significant for systems of metallic nanoparticles in liquids.

Mechanics↗

SU-F-T-192: Study of Robustness Analysis Method of Multiple Field Optimized IMPT Plans for Head & Neck Patients

Purpose: Proton therapy is more sensitive to uncertainties than photon treatments due to protons’ finite range depending on the tissue density. Worst case scenario (WCS) method originally proposed by Lomax has been adopted in our institute for robustness analysis of IMPT plans. This work demonstrates that WCS method is sufficient enough to take into account of the uncertainties which could be encountered during daily clinical treatment. Methods: A fast and approximate dose calculation method is developed to calculate the dose for the IMPT plan under different setup and range uncertainties. Effects of two factors, inversed square factor and range uncertainty, are explored. WCS robustness analysis method was evaluated using this fast dose calculation method. The worst-case dose distribution was generated by shifting isocenter by 3 mm along x,y and z directions and modifying stopping power ratios by ±3.5%. 1000 randomly perturbed cases in proton range and x, yz directions were created and the corresponding dose distributions were calculated using this approximated method. DVH and dosimetric indexes of all 1000 perturbed cases were calculated and compared with the result using worst case scenario method. Results: The distributions of dosimetric indexes of 1000 perturbed cases were generated and compared with the results using worst case scenario. For D95 of CTVs, at least 97% of 1000 perturbed cases show higher values than the one of worst case scenario. For D5 of CTVs, at least 98% of perturbed cases have lower values than worst case scenario. Conclusion: By extensively calculating the dose distributions under random uncertainties, WCS method was verified to be reliable in evaluating the robustness level of MFO IMPT plans of H&N patients. The extensively sampling approach using fast approximated method could be used in evaluating the effects of different factors on the robustness level of IMPT plans in the future.

60 APPLIED LIFE SCIENCES↗

Numerical methods and hypoexponential approximations for gamma distributed delay differential equations

Abstract Gamma distributed delay differential equations (DDEs) arise naturally in many modelling applications. However, appropriate numerical methods for generic gamma distributed DDEs have not previously been implemented. Modellers have therefore resorted to approximating the gamma distribution with an Erlang distribution and using the linear chain technique to derive an equivalent system of ordinary differential equations (ODEs). In this work, we address the lack of appropriate numerical tools for gamma distributed DDEs in two ways. First, we develop a functional continuous Runge–Kutta (FCRK) method to numerically integrate the gamma distributed DDE without resorting to Erlang approximation. We prove the fourth-order convergence of the FCRK method and perform numerical tests to demonstrate the accuracy of the new numerical method. Nevertheless, FCRK methods for infinite delay DDEs are not widely available in existing scientific software packages. As an alternative approach to solving gamma distributed DDEs, we also derive a hypoexponential approximation of the gamma distributed DDE. This hypoexponential approach is a more accurate approximation of the true gamma distributed DDE than the common Erlang approximation but, like the Erlang approximation, can be formulated as a system of ODEs and solved numerically using standard ODE software. Using our FCRK method to provide reference solutions, we show that the common Erlang approximation may produce solutions that are qualitatively different from the underlying gamma distributed DDE. However, the proposed hypoexponential approximations do not have this limitation. Finally, we apply our hypoexponential approximations to perform statistical inference on synthetic epidemiological data to illustrate the utility of the hypoexponential approximation.

97 MATHEMATICS AND COMPUTING↗

Fourier Analysis of the Multilevel QD Method with {alpha}-Approximation For Time-Dependent Radiative Transfer Problems

Thermal radiation transfer (TRT) problems describe interaction of photon radiation with matter. They are defined by the time-dependent radiative transfer (RT) equation for the specific intensity I coupled with the energy balance (EB) equation. This class of problems is characterised by high dimensionality, multiple scales and strong nonlinearity. In this paper we analyze the multilevel quasi-diffusion (QD) method for solving TRT problems. This is a nonlinear method of moments with exact closures. The QD method enables one to reduce dimensionality of the original radiative transfer problem. It is formulated by a nonlinear system of equations consisting of (i) the high-order RT (HORT) equation for the intensity of photon radiation and (ii) low-order equations for the angular and energy moments of the intensity. The nonlinear closures are formulated by means of stable factors. The low-order problem brings the radiative transfer model to the level of multi-physics equations. The multilevel QD (MLQD) method possesses attractive features for efficient multi-physics coupling. To calculate the high-order RT solution we apply the {alpha}-approximation which assumes that the intensity changes exponentially in time over a time step. The {alpha}-approximation leads to the steady-state RT equation with a modified opacity and enables one to avoid storing the high-dimensional solution from the previous time step. The rate of change is evaluated by the solution of the time-dependent low-order QD (LOQD) equations for moments. We note that the HORT equation in the {alpha}-approximation contains an extra nonlinear term. It affects the iterative behaviour of the MLQD method. The obtained computational results indicated that the rate of convergence is slower compared to the rate when the time-dependent HORT is used to calculate the intensity. In this paper we perform an iterative stability analysis of the MLQD method in a discretized form and analyze convergence of outer (transport) iterations on a time step. We consider multigroup TRT problems in 1D slab geometry. The fully implicit scheme is used to discretize equations in time. The HORT equation is approximated with the method of step characteristics. A finite volume method is applied for spatial discretization of multigroup LOQD equations. The grey (one-group) LOQD equations are algebraically consistent with the multigroup LOQD equations. To perform Fourier analysis a model TRT problem in an infinite homogeneous medium is considered. The system of nonlinear equations of the MLQD method is linearized near an equilibrium solution at some given temperature. (authors)

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Compressible laminar boundary layers with large acceleration and cooling

The paper extends an approximate prediction method involving the integral form of the momentum equation to deduce the flow quantities of interest when compressibility effects become important and heat transfer may occur. The approximate method is applicable to a two-dimensional laminar boundary layer on an impermeable surface of negligible curvature for a specified integral momentum equation. The proposed approximate method involves the use of similar solutions in conjunction with the integral momentum equation for isentropic freestream flow of a perfect gas by assuming that viscosity is proportional to stagnation condition and Prandtl number is unity. The study pertains to larger values of acceleration (acceleration parameter to 20 rather than 2) than previously considered to account for rapidly accelerating flows such as in supersonic nozzles. A discussion of some applications confirms the advantages of the method in establishing better confidence in the predictions.

Back, L. H.↗

Less-Complex Method of Classifying MPSK

An alternative to an optimal method of automated classification of signals modulated with M-ary phase-shift-keying (M-ary PSK or MPSK) has been derived. The alternative method is approximate, but it offers nearly optimal performance and entails much less complexity, which translates to much less computation time. Modulation classification is becoming increasingly important in radio-communication systems that utilize multiple data modulation schemes and include software-defined or software-controlled receivers. Such a receiver may "know" little a priori about an incoming signal but may be required to correctly classify its data rate, modulation type, and forward error-correction code before properly configuring itself to acquire and track the symbol timing, carrier frequency, and phase, and ultimately produce decoded bits. Modulation classification has long been an important component of military interception of initially unknown radio signals transmitted by adversaries. Modulation classification may also be useful for enabling cellular telephones to automatically recognize different signal types and configure themselves accordingly. The concept of modulation classification as outlined in the preceding paragraph is quite general. However, at the present early stage of development, and for the purpose of describing the present alternative method, the term "modulation classification" or simply "classification" signifies, more specifically, a distinction between M-ary and M'-ary PSK, where M and M' represent two different integer multiples of 2. Both the prior optimal method and the present alternative method require the acquisition of magnitude and phase values of a number (N) of consecutive baseband samples of the incoming signal + noise. The prior optimal method is based on a maximum- likelihood (ML) classification rule that requires a calculation of likelihood functions for the M and M' hypotheses: Each likelihood function is an integral, over a full cycle of carrier phase, of a complicated sum of functions of the baseband sample values, the carrier phase, the carrier-signal and noise magnitudes, and M or M'. Then the likelihood ratio, defined as the ratio between the likelihood functions, is computed, leading to the choice of whichever hypothesis - M or M'- is more likely. In the alternative method, the integral in each likelihood function is approximated by a sum over values of the integrand sampled at a number, 1, of equally spaced values of carrier phase. Used in this way, 1 is a parameter that can be adjusted to trade computational complexity against the probability of misclassification. In the limit as 1 approaches infinity, one obtains the integral form of the likelihood function and thus recovers the ML classification. The present approximate method has been tested in comparison with the ML method by means of computational simulations. The results of the simulations have shown that the performance (as quantified by probability of misclassification) of the approximate method is nearly indistinguishable from that of the ML method (see figure).

Hamkins, Jon↗

Approximate two layer (inviscid/viscous) methods to model aerothermodynamic environments

Aeroassisted orbital transfer vehicles (AOTVs) and HSTs fly hypersonically at altitudes higher than those of Apollo and Space Shuttle type, whose flowfield could be considered a continuum with insignificant chemical nonequilibrium effects; AOTVs and HSTs will by contrast fly in the transitional regime between continuum and free molecule flow, where chemical nonequilibrium within the shock layer can be significant. Attention is presently given to the numerical and approximate methods of hypersonic flowfield calculation, emphasizing the role of aerodynamic heating. The development status of boundary layer methods, fully viscous methods, and gas models is evaluated.

De Jarnette, Fred R.↗

Generalized quantum master equations in and out of equilibrium: When can one win?

Generalized quantum master equations (GQMEs) are an important tool in modeling chemical and physical processes. For a large number of problems, it has been shown that exact and approximate quantum dynamics methods can be made dramatically more efficient, and in the latter case more accurate, by proceeding via the GQME formalism. However, there are many situations where utilizing the GQME approach with an approximate method has been observed to return the same dynamics as using that method directly. Here, for systems both in and out of equilibrium, we provide a more detailed understanding of the conditions under which using an approximate method can yield benefits when combined with the GQME formalism. In particular, we demonstrate the necessary manipulations, which are satisfied by exact quantum dynamics, that are required to recast the memory kernel in a form that can be analytically shown to yield the same result as a direct application of the dynamics regardless of the approximation used. By considering the connections between these forms of the kernel, we derive the conditions that approximate methods must satisfy if they are to offer different results when used in conjunction with the GQME formalism. These analytical results thus provide new insights as to when proceeding via the GQME approach can be used to improve the accuracy of simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Accumulated approximation: A new method for structural optimization by iterative improvement

A new method for the solution of non-linear mathematical programming problems in the field of structural optimization is presented. It is an iterative scheme which for each iteration refines the approximation of objective and constraint functions by accumulating the function values of previously visited design points. The method has proven to be competitive for a number of well-known examples of which one is presented here. Furthermore because of the accumulation strategy, the method produces convergence even when the sensitivity analysis is inaccurate.

Rasmussen, John↗