Search NASA⌕ Search

SEARCH · Search NASA

Results for “linear estimation”

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 253 records · Page 14
All sourcesNASA NTRSDOE OSTI
NASA NTRS1959

Analysis and Evaluation of Supersonic Underwing Heat Addition

The linearized theory for heat addition under a wing has been developed to optimize wing geometry, heat addition, and angle of attack. The optimum wing has all of the thickness on the underside of the airfoil, with maximum-thickness point well downstream, has a moderate thickness ratio, and operates at an optimum angle of attack. The heat addition is confined between the fore Mach waves from under the trailing surface of the wing. By linearized theory, a wing at optimum angle of attack may have a range efficiency about twice that of a wing at zero angle of attack. More rigorous calculations using the method of characteristics for particular flow models were made for heating under a flat-plate wing and for several wings with thickness, both with heat additions concentrated near the wing. The more rigorous calculations yield in practical cases efficiencies about half those estimated by linear theory. An analysis indicates that distributing the heat addition between the fore waves from the undertrailing portion of the wing is a way of improving the performance, and further calculations appear desirable. A comparison of the conventional ramjet-plus wing with underwing heat addition when the heat addition is concentrated near the wing shows the ramjet to be superior on a range basis up to Mach number of about B. The heat distribution under the wing and the assumed ramjet and airframe performance may have a marked effect on this conclusion. Underwing heat addition can be useful in providing high-altitude maneuver capability at high flight Mach numbers for an airplane powered by conventional ramjets during cruise.

wings↗
NASA NTRS2006

Predictions of Leukemia Risks to Astronauts from Solar Particle Events

Leukemias consisting of acute and chronic myeloid leukemia and acute lymphatic lymphomas represent the earliest cancers that appear after radiation exposure, have a high lethality fraction, and make up a significant fraction of the overall fatal cancer risk from radiation for adults. Several considerations impact the recommendation of a preferred model for the estimation of leukemia risks from solar particle events (SPE's): The BEIR VII report recommends several changes to the method of calculation of leukemia risk compared to the methods recommended by the NCRP Report No. 132 including the preference of a mixture model with additive and multiplicative components in BEIR VII compared to the additive transfer model recommended by NCRP Report No. 132. Proton fluences and doses vary considerably across marrow regions because of the characteristic spectra of primary solar protons making the use of an average dose suspect. Previous estimates of bone marrow doses from SPE's have used an average body-shielding distribution for marrow based on the computerized anatomical man model (CAM). We have developed an 82-point body-shielding distribution that faithfully reproduces the mean and variance of SPE doses in the active marrow regions (head and neck, chest, abdomen, pelvis and thighs) allowing for more accurate estimation of linear- and quadratic-dose components of the marrow response. SPE's have differential dose-rates and a pseudo-quadratic dose response term is possible in the peak-flux period of an event. Also, the mechanistic basis for leukemia risk continues to improve allowing for improved strategies in choosing dose-rate modulation factors and radiation quality descriptors. We make comparisons of the various choices of the components in leukemia risk estimates in formulating our preferred model. A major finding is that leukemia could be the dominant risk to astronauts for a major solar particle event.

Cucinotta, F. A.↗
NASA NTRS1998

Estimating Subcanopy Soil Moisture with RADAR

The subcanopy soil moisture of a boreal old jack pine forest is estimated using polarimetric L- and P-band AIRSAR data. Model simulations have shown that for this stand, the principal scattering mechanism responsible for radar backscatter is the double-bounce mechanism between the tree trunks and the ground.

radar soil moisture vegetation non-linear estimati↗
NASA NTRS1986

On the design of classifiers for crop inventories

Crop proportion estimators that use classifications of satellite data to correct, in an additive way, a given estimate acquired from ground observations are discussed. A linear version of these estimators is optimal, in terms of minimum variance, when the regression of the ground observations onto the satellite observations in linear. When this regression is not linear, but the reverse regression (satellite observations onto ground observations) is linear, the estimator is suboptimal but still has certain appealing variance properties. In this paper expressions are derived for those regressions which relate the intercepts and slopes to conditional classification probabilities. These expressions are then used to discuss the question of classifier designs that can lead to low-variance crop proportion estimates. Variance expressions for these estimates in terms of classifier omission and commission errors are also derived.

Heydorn, R. P.↗
DOE OSTI2018

An Evaluation of Bang-Time Measurements from a Multichannel, Triaxial, nTOF Diagnostic for MagLIF Experiments at the Z facility

Neutron bang times for a series of MagLIF (Magnetic Liner Inertial Fusion) experiments with D 2 -filled targets have been measured at the Z facility. The emitted neutrons were detected as current-mode pulses in a multichannel, neutron time-of-flight (nTOF) diagnostic with conventional, scintillator-photomultiplier-tube (PMT) detectors. In these experiments, the detectors were fielded at known, fixed distances L (690-2510 cm) from the target, and on three, non-coplanar (but convergent) lines-of-sight (LOS). The primary goal of this diagnostic was to estimate a fiducial time (bang time) relative to an externally generated time-base for synchronizing all the diagnostics in an experiment. Recorded arrival times (A7) of the pulses were characterized experimentally by three numerical methods: a first-moment estimate (centroid) and two nodal measures — Savitzky-Golay (SG) smoothing and a single point peak estimate of the raw data. These times were corrected for internal detector time delays (transit and impulse-response function) — an adjustment that linked the recorded ATs to the corresponding arrival of uncollided neutrons at each detector. The bang time was then estimated by linearly regressing the arrival times against the associated distances to the source; t bang (on the system timescale) was taken as the temporal intercept of the regression equation at distance L = 0. This article reports the analysis for a representative shot #2584 for which (a) the recorded ATs — even without detector corrections — agreed by method in each channel to within 1-2 ns; (b) internal corrections were each ~3 — 5 ns; and (c) a 95% uncertainty (confidence) interval for t bang in this shot was estimated at ±3 ns with 4 degrees of freedom. A secondary goal for this diagnostic was to check that the bang time measurements corresponded to neutrons emitted by the D(d,n) 3 He reaction in a thermalized DD plasma. According to the theoretical studies by Brysk, such neutrons should be emitted with an isotropic Gaussian distribution of mean kinetic energy $ \overline{E}$ of 2.449 MeV; this energy translates to a mean neutron speed $ \overline{u}$ of 2.160 cm/ns [D. H. Munro, Nuclear Fusion, 56(3) 036001 (2016)]. In the MagLIF series of shots there was no evidence of spatial asymmetry in the time-distance regressions, and it was possible to extract the mean neutron speed from the slope of these fits. In shot 2584 $ \overline{u}$ was estimated at 2.152 cm/ns ± 0.010 cm/ns [95 % confidence, 4 dof] and the mean kinetic energy $ \overline{E}$ (with relativistic corrections) was 2.431 MeV ± 0.022 MeV [95 % confidence, 4 dof] — results supporting the assumption that D-D neutrons were, in fact, measured.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
NASA NTRS1976

Digital controllers for VTOL aircraft

Using linear-optimal estimation and control techniques, digital-adaptive control laws have been designed for a tandem-rotor helicopter which is equipped for fully automatic flight in terminal area operations. Two distinct discrete-time control laws are designed to interface with velocity-command and attitude-command guidance logic, and each incorporates proportional-integral compensation for non-zero-set-point regulation, as well as reduced-order Kalman filters for sensor blending and noise rejection. Adaptation to flight condition is achieved with a novel gain-scheduling method based on correlation and regression analysis. The linear-optimal design approach is found to be a valuable tool in the development of practical multivariable control laws for vehicles which evidence significant coupling and insufficient natural stability.

Stengel, R. F.↗
DOE OSTI2021

Injection data analysis using material balance time for CO 2 storage capacity estimation in deep closed saline aquifers

Estimating the ultimate storage capacity of deep saline aquifers is important to address the formation potential to store the envisioned large volumes of CO 2 . Injection data (i.e. injection rate, bottomhole pressure, and cumulative injected volume of CO 2 ) are routinely recorded during storage operations. These data contain valuable information on the subsurface (e.g. the reservoir pore volume and the formation storage capacity) that can be extracted. In this paper, we present a two-step graphical technique to infer the pore volume and the ultimate storage capacity of closed saline aquifers by analyzing the available injection data. First, the pore volume is inferred through adapting the concept of the material balance time. Material balance time is an approximate superposition time function developed to interpret production data from oil and gas wells operating at variable pressure/rate conditions during the boundary-dominated flow period. Using material balance techniques, the ultimate storage capacity is then estimated through linear extrapolation of the average pressure trend to the maximum allowable pressure the formation can withstand. The average pressure is not available in practice, but is can be obtained from the injection data. Two approaches are presented in this study to calculate the average pressure; namely the rigorous and the approximate approaches. Unlike the rigorous approach, the approximate approach does not require a prior knowledge of some reservoir properties (e.g. relative permeability, absolute permeability, formation porosity and thickness) to calculate the average pressure. To investigate its potential and reliability in analyzing CO 2 injection data, the proposed technique is applied to four synthetic cases representing different well operating conditions. Results indicate that the approximate approach consistently overestimates the actual (simulated) storage capacity as compared to the rigorous approach. The agreement - between the inferred and the simulated reservoir pore volume, and between the analytical and numerical estimates of storage capacity - validates the potential application of the technique to CO 2 storage in closed saline aquifers. The technique is further substantiated through application to a field data set utilized from a commercial-scale geological storage (CGS) project. Finally, field data interpretation shows that the proposed technique can be utilized to identify the degree of hydraulic continuity and reservoir compartmentalization within a target formation by interpreting the corresponding pressure and rate responses.

02 PETROLEUM↗
DOE OSTI2018

Convergence of the multigrid reduction in time algorithm for the linear elasticity equations

This paper presents some recent advances for parallel-in-time methods applied to linear elasticity. With recent computer architecture changes leading to stagnant clock speeds, but ever increasing numbers of cores, future speedups will be available through increased concurrency. Thus, sequential algorithms, such as time stepping, will suffer a bottleneck. This paper explores multigrid reduction in time (MGRIT) for an important application area, linear elasticity. Previously, efforts at parallel-in-time for elasticity have experienced difficulties, for example, the beating phenomenon. As a result, practical parallel-in-time algorithms for this application area currently do not exist. This paper proposes some solutions made possible by MGRIT (e.g., slow temporal coarsening and FCF-relaxation) and, more importantly, a different formulation of the problem that is more amenable to parallel-in-time methods. Using a recently developed convergence theory for MGRIT and Parareal, we show that the changed formulation of the problem avoids the instability issues and allows the reduction of the error using two temporal grids. We then extend our approach to the multilevel case, where we demonstrate how slow temporal coarsening improves convergence. Finally, the paper ends with supporting numerical results showing a practical algorithm enjoying speedup benefits over the sequential algorithm.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1982

A new solution to parameter adaptive estimation of random processes

This paper is concerned with the development of an adaptive state estimator that is capable of tracking switched linear plants that undergo rapid configuration changes. The particular adaptive estimator developed here is called the Sliding Window Detector/Estimator (SWDE) algorithm. Unlike previous algorithms, the SWDE algorithm is designed specifically for the switched-linear plant problem. It uses a joint detection/estimation approach to give a very close approximation to the unrealizable optimum switched-linear estimator. An extremely reliable and accurate estimator can be constructed by combining a modified Parameter Adaptive Estimation (PAE) algorithm with SWDE. The algorithm has been fully verified by extensive computer simulation, and the implementation advantages afforded by this method make it suitable for use in a wide variety of applications.

Zwicke, P. E.↗
DOE OSTI2020

An electrochemical study of hydrogen in molten 2LiF-BeF 2 (FLiBe) with addition of LiH

The chemical and transport behavior of hydrogen isotopes in molten 2LiF-BeF 2 (FLiBe) is of interest for the design of tritium management systems in nuclear fission and fusion reactors that use FLiBe. The chemical reaction of LiH with FLiBe is used to introduce hydrogen in the molten salt and electrochemical methods are used for in situ studies of hydrogen in FLiBe. LiH reacts with molten FLiBe to generate an electroactive species whose voltammetry peak is proportional to the added quantity of LiH. The cyclic voltammetry reaction potential of 2.009 ± 0.050 V vs Be/Be 2+ and the electron exchange of n = 0.8 ± 0.5 are consistent with the one-electron oxidation of dissolved H in the zero valence state, H o to H + . The concentration of H o is estimated by linear sweep voltammetry at 60–80% of the hydrogen introduced by the reaction of 0.42 mol% LiH with FLiBe and after eleven hours it remains above 40%. It is postulated that covalent BeH 2 is formed in FLiBe upon LiH addition, as a FLiBe-soluble quasi-stable intermediate product. Furthermore, the results provide an evaluation of LiH as a means of introducing dissolved hydrogen in FLiBe, enabling electrochemical methods as tools to advance the understanding of the chemistry of hydrogen isotopes in FLiBe.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2026

Modulation induced by very-large-scale motions on the inclination angle of wall-attached eddies: an atmospheric surface layer experiment

The forward leaning inclination angle, 𝛾, of coherent turbulent structures is a well-known feature of wall-bounded turbulent flows. Although invariant across friction Reynolds numbers within the range 𝑅𝑒 𝜏 =10 3 −10 6 , 𝛾 can vary significantly across turbulent scales within a high-Reynolds-number flow. Very-large-scale motions (VLSMs) are known to induce significant changes in the instantaneous shear profile, which is a conditioning event that could trigger variability in the inclination angle of smaller coherent turbulent structures. Although this aspect has been extensively studied via numerical and laboratory experiments, few studies have explored this feature for a very-high-Reynolds-number atmospheric flow. In this work, the inclination angle of turbulent structures within the atmospheric surface layer at a very high Reynolds number (𝑅𝑒 𝜏 =7.9 ×10 5 ) is investigated by deploying a scanning Doppler light detection and ranging and a super large particle image velocimetry (SLPIV) apparatus. The inclination angle of wall-attached eddies is inferred either from the two-point correlation of streamwise velocity (𝛾 =41.1∘) or with a scale-dependent approach through the spectral linear stochastic estimator (SLSE). The SLSE (and, thus, the scale-dependent inclination angle) is conditionally evaluated based on the high- and low-momentum events induced by VLSMs, both in the streamwise (𝑢′𝑉𝐿𝑆𝑀) and in the vertical (𝑤′ 𝑉𝐿𝑆𝑀) velocity components. As a result, lower inclination angles (𝛾 =30° −50°) are found for 𝑢$^{'}_{VLSM}$ >0 (𝑤$^{'}_{VLSM}$ <0), while higher values (50° −85°) are ascribed to 𝑢$^{'}_{VLSM}$ <0 (𝑤$^{'}_{VLSM}$ >0). This result emphasises the primary role that VLSMs play in shaping the wall-attached eddy geometry, which, in turn, is crucial to determine the Reynolds stress balance within the wall-attached eddy range.

42 ENGINEERING↗
DOE OSTI2025

Decomposing a Compound Flood Event in an Urban Pacific Northwest Estuary: Primary Drivers and Projections for the Future

The Duwamish River Estuary (DRE) of Washington is prone to compound flooding during atmospheric river (AR) events. The processes contributing to such flooding (coastal and fluvial) have remained opaque to municipalities that are increasingly impacted. Here, we conduct a suite of coupled atmosphere-hydrology-ocean model simulations with varying forcing combinations (tide, surge, and/or river discharge) to identify the primary drivers of compound flooding during a recent AR event. We also test year 2,100 climate forcing to project how flooding and drivers may change in the future for the same event. We identify a clear distinction between dynamics in the downstream, engineered portion of the DRE compared to the upstream, “natural” river. Downstream, tides dominate water levels but contributions from storm surge and nonlinear tide-surge interaction elevate tide-only high waters from no flooding to major flooding. Upstream, total water levels during the event are ∼6 cm higher than downstream due to an increasing influence of river discharge over surge and tides. Notably, nonlinear surge-river and tide-river interactions act to reduce upstream water levels up to 50% compared to estimates which linearly sum tides, surge, and river, likely reducing flood vulnerability. Under two future climate scenarios: one with only sea level rise (SLR) and another with SLR plus atmospheric warming, we find little change in mechanism contributions to water levels. Expanded flooding in both cases is largely due to SLR, as a ∼50% increase to river discharge under the warming scenario has no impact downstream and marginally increases (∼3 cm) water level upstream.

atmospheric river↗
DOE OSTI2024

The cluster decomposition of the configurational energy of multicomponent alloys

Abstract The cluster expansion method (CEM) is a widely used lattice-based technique in the study of multicomponent alloys. Despite its prevalent use, a clear understanding of expansion terms is lacking. We present a modern mathematical formalism of the CEM and introduce thecluster decomposition—a unique and basis-independent decomposition for functions of the atomic configuration in a crystal. We identify the cluster decomposition as an invariant ANOVA decomposition; and demonstrate how functional analysis of variance and sensitivity analysis can be used to interpret interactions among species. Furthermore, we show how the mathematical structure of the cluster decomposition enables numerical evaluation that scales with the number of clusters and is independent of the number of species. Overall, our work enables rigorous interpretations of interactions among species, provides opportunities to explore parameter estimation beyond linear regression, introduces a numerical efficient implementation, and enables analysis of cluster expansions based on established mathematical and statistical principles.

Chemistry↗
DOE OSTI2023

Spectra and characteristics of somatic mutations induced by ionizing radiation in hematopoietic stem cells

Spectra and frequencies of spontaneous and X-ray-induced somatic mutations were revealed with mouse long-term hematopoietic stem cells (LT-HSCs) by whole-genome sequencing of clonal cell populations propagated in vitro from single isolated LT-HSCs. SNVs and small indels were the most common types of somatic mutations, and increased up to twofold to threefold by whole-body X-irradiation. Base substitution patterns in the SNVs suggested a role of reactive oxygen species in radiation mutagenesis, and signature analysis of single base substitutions (SBS) revealed a dose-dependent increase of SBS40. Most of spontaneous small deletions were shrinkage of tandem repeats, and X-irradiation specifically induced small deletions out of tandem repeats (non-repeat deletions). Presence of microhomology sequences in non-repeat deletions suggested involvement of microhomology mediated end-joining repair mechanisms as well as nonhomologous end-joining in radiation-induced DNA damages. We also identified multisite mutations and structural variants (SV), i.e., large indels, inversions, reciprocal translocations, and complex variants. The radiation-specificity of each mutation type was evaluated from the spontaneous mutation rate and the per-Gy mutation rate estimated by linear regression, and was highest with non-repeat deletions without microhomology, followed by those with microhomology, SV except retroelement insertions, and multisite mutations; these types were thus revealed as mutational signatures of ionizing radiation. Further analysis of somatic mutations in multiple LT-HSCs indicated that large fractions of postirradiation LT-HSCs originated from single LT-HSCs that survived the irradiation and then expanded in vivo to confer marked clonality to the entire hematopoietic system, with varying clonal expansion and dynamics depending on radiation dose and fractionation.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗
DOE OSTI2020

A molecular dynamics simulation study of thermal conductivity anisotropy in β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (β-HMX)

Molecular dynamics (MD) simulations were used to predict the thermal conductivity of β-octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (β-HMX) along directions normal to the (011), (110), and (010) crystal planes. These directions were selected based on the measured morphological importance of the corresponding crystal surfaces. A reverse non-equilibrium MD approach was used wherein a constant heat flux is imposed along a prescribed direction and the resulting steady-state temperature gradient determined. The coefficient of thermal conductivity λ is the quotient of heat flux and temperature gradient (i.e. Fourier's law). Finite-size effects and sensitivity to imposed heat flux were investigated. The results reveal a modest dependence of the conductivity on crystal orientation, significant finite-size effects, and low sensitivity to imposed flux so long as the Fourier's law analysis is limited to the spatial interval in the simulation cell for which the temperature gradient is constant. Infinite-length thermal conductivities were estimated by linear regression of λ -1 (L) versus reciprocal cell length L -1 for each direction. The predicted values are systematically larger, but within a factor of two, than most published experimental determinations, the latter of which were obtained for pressed-powder or composite samples rather than oriented single crystals.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2020

Estimation of elastic properties of an additively manufactured lattice using resonant ultrasound spectroscopy

Development of improved approaches in the characterization of additively manufactured structures continues to be a topic of interest for the advanced manufacturing community. This article will investigate an approach using resonant ultrasound spectroscopy (RUS) to determine the effective elastic constants of an orthotropic lattice structure. The evaluation is performed on a cube shaped 316 L stainless steel test specimen, constructed using selective laser melting techniques. The approach uses RUS techniques in conjunction with the assumption that in the frequency regime of interest, the wavelength of the diagnostic ultrasound is greater than the discrete structural features of the unit cell of the lattice; thus, the AM structure can be treated as an anisotropic continuum with effective material properties and symmetry inherited from the unit cell. The RUS analysis estimates the nine elastic coefficients associated with orthotropic sample symmetry, which, in turn, are used to determine the elastic moduli and Poisson ratios. Current results show good agreement between experiments and modeled data. Comparisons to published results are also in good agreement, indicating the potential applicability of this characterization technique for estimating the linear elastic properties of innovative additive manufactured metal lattice structures.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Kinematics reconstruction in solenoidal spectrometers operated in active target mode

Here, we discuss the reconstruction of low-energy nuclear reaction kinematics from charged-particle tracks in solenoidal spectrometers working in Active Target Time Projection Chamber mode. In this operation mode, reaction products are tracked within the active gas medium of the Active Target with a three dimensional space point cloud. We have inferred the reaction kinematics from the point cloud using an algorithm based on a linear quadratic estimator (Kalman filter). The performance of this algorithm has been evaluated using experimental data from nuclear reactions measured with the Active Target Time Projection Chamber (AT-TPC) detector.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗