SEARCH · Engineering Papers

Results for “model structural error”

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 1,783 records · Page 99

All sources NASA NTRS DOE OSTI

DOE OSTI2023

Temporal Subsampling Diminishes Small Spatial Scales in Recurrent Neural Network Emulators of Geophysical Turbulence

The immense computational cost of traditional numerical weather and climate models has sparked the development of machine learning (ML) based emulators. Because ML methods benefit from long records of training data, it is common to use data sets that are temporally subsampled relative to the time steps required for the numerical integration of differential equations. Here, we investigate how this often overlooked processing step affects the quality of an emulator's predictions. We implement two ML architectures from a class of methods called reservoir computing: (a) a form of Nonlinear Vector Autoregression (NVAR), and (b) an Echo State Network (ESN). Despite their simplicity, it is well documented that these architectures excel at predicting low dimensional chaotic dynamics. We are therefore motivated to test these architectures in an idealized setting of predicting high dimensional geophysical turbulence as represented by Surface Quasi-Geostrophic dynamics. In all cases, subsampling the training data consistently leads to an increased bias at small spatial scales that resembles numerical diffusion. Interestingly, the NVAR architecture becomes unstable when the temporal resolution is increased, indicating that the polynomial based interactions are insufficient at capturing the detailed nonlinearities of the turbulent flow. The ESN architecture is found to be more robust, suggesting a benefit to the more expensive but more general structure. Spectral errors are reduced by including a penalty on the kinetic energy density spectrum during training, although the subsampling related errors persist. Future work is warranted to understand how the temporal resolution of training data affects other ML architectures.

58 GEOSCIENCES↗

DOE OSTI2022

Soil-related developments of the Biome-BGCMuSo v6.2 terrestrial ecosystem model

Terrestrial biogeochemical models are essential tools to quantify climate–carbon cycle feedback and plant–soil relations from local to global scale. In this study, a theoretical basis is provided for the latest version of the Biome-BGCMuSo biogeochemical model (version 6.2). Biome-BGCMuSo is a branch of the original Biome-BGC model with a large number of developments and structural changes. Earlier model versions performed poorly in terms of soil water content (SWC) dynamics in different environments. Moreover, lack of detailed nitrogen cycle representation was a major limitation of the model. Since problems associated with these internal drivers might influence the final results and parameter estimation, additional structural improvements were necessary. In this paper the improved soil hydrology as well as the soil carbon and nitrogen cycle calculation methods are described in detail. Capabilities of the Biome-BGCMuSo v6.2 model are demonstrated via case studies focusing on soil hydrology, soil nitrogen cycle, and soil organic carbon content estimation. Soil-hydrology-related results are compared to observation data from an experimental lysimeter station. The results indicate improved performance for Biome-BGCMuSo v6.2 compared to v4.0 (explained variance increased from 0.121 to 0.8 for SWC and from 0.084 to 0.46 for soil evaporation; bias changed from -0.047 to -0.007 m 3 m -3 for SWC and from -0.68 to -0.2 mm d -1 for soil evaporation). Simulations related to nitrogen balance and soil CO 2 efflux were evaluated based on observations made in a long-term field experiment under crop rotation. The results indicated that the model is able to provide realistic nitrate content estimation for the topsoil. Soil nitrous oxide (N 2 O) efflux and soil respiration simulations were also realistic, with overall correspondence with the observations (for the N 2 O efflux simulation bias was between -0.13 and -0.1 mgNm -2 d -1 , and normalized root mean squared error (NRMSE) was 32.4 %–37.6 %; for CO 2 efflux simulations bias was 0.04–0.17 gCm -2 d -1 , while NRMSE was 34.1 %–40.1 %). Sensitivity analysis and optimization of the decomposition scheme are presented to support practical application of the model. The improved version of Biome-BGCMuSo has the ability to provide more realistic soil hydrology representation as well as nitrification and denitrification process estimation, which represents a major milestone.

58 GEOSCIENCES↗

NASA NTRS2003

First Year Wilkinson Microwave Anisotropy Probe (WMAP)Observations: Preliminary Maps and Basic Results

We present full sky microwave maps in five frequency bands (23 to 94 GHz) from the WMAP first year sky survey. Calibration errors are less than 0.5% and the low systematic error level is well specified. The cosmic microwave background (CMB) is separated from the foregrounds using multifrequency data. The sky maps are consistent with the 7 in. full-width at half-maximum (FWHM) Cosmic Background Explorer (COBE) maps. We report more precise, but consistent, dipole and quadrupole values. The CMB anisotropy obeys Gaussian statistics with -58 less than f(sub NL) less than 134 (95% CL). The 2 less than or = l less than or = 900 anisotropy power spectrum is cosmic variance limited for l less than 354 with a signal-to-noise ratio greater than 1 per mode to l = 658. The temperature-polarization cross-power spectrum reveals both acoustic features and a large angle correlation from reionization. The optical depth of reionization is tau = 0.17 +/- 0.04, which implies a reionization epoch of t(sub r) = 180(sup +220, sub -80) Myr (95% CL) after the Big Bang at a redshift of z(sub r) = 20(sup +10, sub -9) (95% CL) for a range of ionization scenarios. This early reionization is incompatible with the presence of a significant warm dark matter density. A best-fit cosmological model to the CMB and other measures of large scale structure works remarkably well with only a few parameters. The age of the best-fit universe is t(sub 0) = 13.7 +/- 0.2 Gyr old. Decoupling was t(sub dec) = 379(sup +8, sub -7)kyr after the Big Bang at a redshift of z(sub dec) = 1089 +/- 1. The thickness of the decoupling surface was Delta(sub z(sub dec)) = 195 +/- 2. The matter density of the universe is Omega(sub m)h(sup 2) = 0.135(sup +0.008, sub -0.009) the baryon density is Omega(sub b)h(sup 2) = 0.0224 +/- 0.0009, and the total mass-energy of the universe is Omega(sub tot) = 1.02 +/- 0.02. There is progressively less fluctuation power on smaller scales, from WMAP to fine scale CMB measurements to galaxies and finally to the Ly-alpha forest. This is accounted for with a running spectral index, significant at the approx. 2(sigma) level. The spectral index of scalar fluctuations is fit as n(sub s) = 0.93 +/-0.03 at wavenumber k(sub o) = 0.05/Mpc ((sub eff) approx. = 700), with a slope of dn(sub s)/d I(sub nk) = -0.031(sup + 0.016, sub -0.018) in the best-fit model.

Bennett, C. L.↗

NASA NTRSDate not supplied

Tropical Convection Through the Lens of the INCUS Mission

The overarching goal of the recently selected NASA INvestigation of Convective UpdraftS (INCUS) mission is to enhance our understanding of why, when and where tropical convective storms form, and why only some of these storms produce extreme weather. Convective storms provide an important pathway for the transport of air and water between Earth’s surface and the upper troposphere. This vertical transport of air and water, often referred to as convective mass flux (CMF), plays a critical role in Earth’s weather and climate system through its impacts on microphysical and precipitation rates, detrainment and upper tropospheric moistening, high cloud feedbacks, and the large-scale circulation. Potential changes to the CMF as a function of the local environment or with changing climates may also have significant implications for severe weather such as flood-producing rainfall, damaging hail and lightning. In spite of the critical role of this vertical transport of water and air, representation of CMF remains a major source of error in weather and climate models, thereby limiting our ability to accurately predict convective storms and their impacts. The tropics-wide observations from INCUS will enhance our understanding of tropical convective storm processes and will provide guidance for representing these processes in numerical models across scales. INCUS is comprised of three SmallSat platforms each carrying a RainCube-heritage Ka-band 7-beam scanning radar. The satellite platforms will be 30 and 90 seconds apart, thus providing three time intervals (30, 90 and 120 seconds) over which observations will be made. INCUS will investigate CMF using a novel time-differenced radar reflectivity profile approach. In addition to the Ka-band radar, a single TEMPEST-D-heritage cross-track-scanning passive microwave radiometer will be housed on the middle SmallSat. The radiometer will provide extensive storm context for the radar observations, as well as observations of the convective anvils. The combination of the radars and radiometer on INCUS will deliver unprecedented three-dimensional views of tropical convective storms. INCUS is the first systematic investigation of the rapidly evolving CMF within tropical convective storms, the observations of which are expected to significantly enhance both our understanding and prediction of storm structure, their dynamics and microphysical processes, and the ways in which these evolve over storm lifetimes. This presentation will highlight the observational capabilities and scientific approach of the INCUS mission.

Susan C. van den Heever↗

DOE OSTI2020

Examination of the inversion of isobaric analogue states in mirror nuclei

A recent study by Hoff and collaborators [Nature (London) 580, 52 (2020)] presented evidence that the ground-state spin of [Math Processing Error] is different from that of its mirror, [Math Processing Error], likely due to an inversion of the ground- and first-excited states. Here, we assess the likelihood of such an inversion arising from normal isospin-symmetry breaking and, more broadly, whether such phenomena challenge our understanding of charge-symmetry-breaking forces in atomic nuclei. By placing the result within the context of previous experimental and theoretical work we demonstrate that this inversion lies entirely within the bounds of normal isospin-symmetry-breaking behavior. We further note that, in the context of isospin, neither level inversions nor the nuclear ground state hold any special significance.

59 ≤ A ≤ 89↗

DOE OSTI2023

A General Spatiotemporal Imputation Framework for Missing Sensor Data

Many applications from precision agriculture, environmental monitoring and transportation networks rely on data collected across space and time over a large geographic area. Missing data poses a significant challenge for any data-driven inference and control tasks. Data imputation or the estimation of missing data can help fill these gaps by utilizing inherent spatial relationships and temporal patterns. A variety of spatiotemporal imputation models have been developed to address missing data in spatiotemporal datasets. However, these classical methods rely on the assumption that the underlying data follows a smooth trend and fail to provide accurate estimates when there is a large number of missing points in the data. Even though there are machine learning driven tensor completion approaches such as convolutional neural network based tensor completion (CoSTCo) that capture the non-linear relationships in the dataset, the transductive nature makes the algorithm less scalable. Thus, existing approaches for estimating the missing information do not effectively capture all dimensions of the spatiotemporal data structure, resulting in erroneous predictions and poor performance. The main contributions of this paper are: (1) We propose a novel inductive framework (G-LSTM) for missing data imputation that integrates a graph neural network with LSTMs to effectively capture both spatial and temporal dependencies. (2) Experimental results on a traffic dataset demonstrate that the proposed GNN integrated with an LSTM framework achieves improved imputation and maintains steady performance even when there are extreme missing conditions in comparison with the state-of-the-art imputation framework (i.e, CoSTCo). (3) The simulation results on a traffic network show up to 69% reduction in mean absolute error and 61% reduction in root mean square error when compared to CoSTCo.

Tharzeen, Aabila↗

NASA NTRS2000

Role of muscle pulleys in producing eye position-dependence in the angular vestibuloocular reflex: a model-based study

It is well established that the head and eye velocity axes do not always align during compensatory vestibular slow phases. It has been shown that the eye velocity axis systematically tilts away from the head velocity axis in a manner that is dependent on eye-in-head position. The mechanisms responsible for producing these axis tilts are unclear. In this model-based study, we aimed to determine whether muscle pulleys could be involved in bringing about these phenomena. The model presented incorporates semicircular canals, central vestibular pathways, and an ocular motor plant with pulleys. The pulleys were modeled so that they brought about a rotation of the torque axes of the extraocular muscles that was a fraction of the angle of eye deviation from primary position. The degree to which the pulleys rotated the torque axes was altered by means of a pulley coefficient. Model input was head velocity and initial eye position data from passive and active yaw head impulses with fixation at 0 degrees, 20 degrees up and 20 degrees down, obtained from a previous experiment. The optimal pulley coefficient required to fit the data was determined by calculating the mean square error between data and model predictions of torsional eye velocity. For active head impulses, the optimal pulley coefficient varied considerably between subjects. The median optimal pulley coefficient was found to be 0.5, the pulley coefficient required for producing saccades that perfectly obey Listing's law when using a two-dimensional saccadic pulse signal. The model predicted the direction of the axis tilts observed in response to passive head impulses from 50 ms after onset. During passive head impulses, the median optimal pulley coefficient was found to be 0.21, when roll gain was fixed at 0.7. The model did not accurately predict the alignment of the eye and head velocity axes that was observed early in the response to passive head impulses. We found that this alignment could be well predicted if the roll gain of the angular vestibuloocular reflex was modified during the initial period of the response, while pulley coefficient was maintained at 0.5. Hence a roll gain modification allows stabilization of the retinal image without requiring a change in the pulley effect. Our results therefore indicate that the eye position-dependent velocity axis tilts could arise due to the effects of the pulleys and that a roll gain modification in the central vestibular structures may be responsible for countering the pulley effect.

NASA Discipline Neuroscience↗

DOE OSTI2022

Reduced order models for thermal radiative transfer problems based on moment equations and data-driven approximations of the Eddington tensor

Here a new group of structure and asymptotic preserving reduced-order models (ROMs) for multidimensional nonlinear thermal radiative transfer (TRT) problems is presented. They are formulated by means of the nonlinear projective approach and data compression techniques. The nonlinear projection is applied to the Boltzmann transport equation (BTE) to derive a hierarchy of low-order moment equations. Approximation of the Eddington tensor that provides exact closure for the system of moment equations is found with projection-based data-driven methodologies. These include the (i) proper orthogonal decomposition (POD), (ii) dynamic mode decomposition (DMD) and (iii) a variant of the DMD. A parameterization is derived for this ROM for the temperature of radiation incoming to the problem domain (the radiation drive temperature). This parameterization is informed from results of a dimensionless study of the TRT problem. Analysis of the ROMs is performed on the classical Fleck-Cummings TRT multigroup test problem in 2D geometry with a radiation-driven Marshak wave. Numerical results are presented to demonstrate the performance of these ROMs for the simulation of evolving radiation and heat waves. Results show these models to be sufficiently accurate for practical computations with rather low-rank representations of the Eddington tensor. As the rank of the approximation is increased, the errors of solutions generated by the ROMs gradually decreases.

42 ENGINEERING↗

NASA NTRS1998

An Investigation of Reliability Models for Ceramic Matrix Composites and their Implementation into Finite Element Codes

The development of modeling approaches for the failure analysis of ceramic-based material systems used in high temperature environments was the primary objective of this research effort. These materials have the potential to support many key engineering technologies related to the design of aeropropulsion systems. Monolithic ceramics exhibit a number of useful properties such as retention of strength at high temperatures, chemical inertness, and low density. However, the use of monolithic ceramics has been limited by their inherent brittleness and a large variation in strength. This behavior has motivated material scientists to reinforce the monolithic material with a ceramic fiber. The addition of a second ceramic phase with an optimized interface increases toughness and marginally increases strength. The primary purpose of the fiber is to arrest crack growth, not to increase strength. The material systems of interest in this research effort were laminated ceramic matrix composites, as well as two- and three- dimensional fabric reinforced ceramic composites. These emerging composite systems can compete with metals in many demanding applications. However, the ongoing metamorphosis of ceramic composite material systems, and the lack of standardized design data has in the past tended to minimize research efforts related to structural analysis. Many structural components fabricated from ceramic matrix composites (CMC) have been designed by "trial and error." The justification for this approach lies in the fact that during the initial developmental phases for a material system fabrication issues are paramount. Emphasis is placed on demonstrating feasibility rather than fully understanding the processes controlling mechanical behavior. This is understandable during periods of rapid improvements in material properties for any composite system. But to avoid the ad hoc approach, the analytical methods developed under this effort can be used to develop rational structural design protocols.

Duffy, Stephen F.↗

NASA NTRS2004

The Precision Formation Flying Integrated Analysis Tool (PFFIAT)

Several space missions presently in the concept phase (e.g. Stellar Imager, Submillimeter Probe of Evolutionary Cosmic Structure, Terrestrial Planet Finder) plan to use multiple spacecraft flying in precise formation to synthesize unprecedently large aperture optical systems. These architectures present challenges to the attitude and position determination and control system; optical performance is directly coupled to spacecraft pointing with typical control requirements being on the scale of milliarcseconds and nanometers. To investigate control strategies, rejection of environmental disturbances, and sensor and actuator requirements, a capability is needed to model both the dynamical and optical behavior of such a distributed telescope system. This paper describes work ongoing at NASA Goddard Space Flight Center toward the integration of a set of optical analysis tools (Optical System Characterization and Analysis Research software, or OSCAR) with the Formation Flying Test Bed (FFTB). The resulting system is called the Precision Formation Flying Integrated Analysis Tool (PFFIAT), and it provides the capability to simulate closed-loop control of optical systems composed of elements mounted on multiple spacecraft. The attitude and translation spacecraft dynamics are simulated in the FFTB, including effects of the space environment (e.g. solar radiation pressure, differential orbital motion). The resulting optical configuration is then processed by OSCAR to determine an optical image. From this image, wavefront sensing (e.g. phase retrieval) techniques are being developed to derive attitude and position errors. These error signals will be fed back to the spacecraft control systems, completing the control loop. A simple case study is presented to demonstrate the present capabilities of the tool.

Stoneking, Eric↗

NASA NTRS2004

The Precision Formation Flying Integrated Analysis Tool (PFFIAT)

Several space missions presently in the concept phase (e.g. Stellar Imager, Sub- millimeter Probe of Evolutionary Cosmic Structure, Terrestrial Planet Finder) plan to use multiple spacecraft flying in precise formation to synthesize unprecedently large aperture optical systems. These architectures present challenges to the attitude and position determination and control system; optical performance is directly coupled to spacecraft pointing with typical control requirements being on the scale of milliarcseconds and nanometers. To investigate control strategies, rejection of environmental disturbances, and sensor and actuator requirements, a capability is needed to model both the dynamical and optical behavior of such a distributed telescope system. This paper describes work ongoing at NASA Goddard Space Flight Center toward the integration of a set of optical analysis tools (Optical System Characterization and Analysis Research software, or OSCAR) with the Formation J?lying Test Bed (FFTB). The resulting system is called the Precision Formation Flying Integrated Analysis Tool (PFFIAT), and it provides the capability to simulate closed-loop control of optical systems composed of elements mounted on multiple spacecraft. The attitude and translation spacecraft dynamics are simulated in the FFTB, including effects of the space environment (e.g. solar radiation pressure, differential orbital motion). The resulting optical configuration is then processed by OSCAR to determine an optical image. From this image, wavefront sensing (e.g. phase retrieval) techniques are being developed to derive attitude and position errors. These error signals will be fed back to the spacecraft control systems, completing the control loop. A simple case study is presented to demonstrate the present capabilities of the tool.

Stoneking, Eric↗

NASA NTRS2016

Exercise Effects on the Brain and Sensorimotor Function in Bed Rest

Long duration spaceflight microgravity results in cephalad fluid shifts and deficits in posture control and locomotion. Effects of microgravity on sensorimotor function have been investigated on Earth using head down tilt bed rest (HDBR). HDBR serves as a spaceflight analogue because it mimics microgravity in body unloading and bodily fluid shifts. Preliminary results from our prior 70 days HDBR studies showed that HDBR is associated with focal gray matter (GM) changes and gait and balance deficits, as well as changes in brain functional connectivity. In consideration of the health and performance of crewmembers we investigated whether exercise reduces the effects of HDBR on GM, functional connectivity, and motor performance. Numerous studies have shown beneficial effects of exercise on brain health. We therefore hypothesized that an exercise intervention during HDBR could potentially mitigate the effects of HDBR on the central nervous system. Eighteen subjects were assessed before (12 and 7 days), during (7, 30, and ~70 days) and after (8 and 12 days) 70 days of 6-degrees HDBR at the NASA HDBR facility in UTMB, Galveston, TX, US. Each subject was randomly assigned to a control group or one of two exercise groups. Exercise consisted of daily supine exercise which started 20 days before the start of HDBR. The exercise subjects participated either in regular aerobic and resistance exercise (e.g. squat, heel raise, leg press, cycling and treadmill running), or aerobic and resistance exercise using a flywheel apparatus (rowing). Aerobic and resistance exercise intensity in both groups was similar, which is why we collapsed the two exercise groups for the current experiment. During each time point T1-weighted MRI scans and resting state functional connectivity scans were obtained using a 3T Siemens scanner. Focal changes over time in GM density were assessed using voxel based morphometry (VBM8) under SPM. Changes in resting state functional connectivity was assessed using both a region of interest (ROI, or seed-to-voxel) approach as well as a whole brain intrinsic connectivity (i.e., voxel-to-voxel) analysis. For the ROI analysis we selected 11 ROIs of brain regions that are involved in sensorimotor function (i.e., L. Insular C., L. Putamen, R. Premotor C., L.+R. Primary Motor C., R. Vestibular C., L. Posterior Cingulate G., R. Cerebellum Lobule V + VIIIb + Crus I, and the R. Superior Parietal G.) and correlated their time course of brain activation during rest with all other voxels in the brain. The whole brain connectivity analysis tests changes in the strength of the global connectivity pattern between each voxel and the rest of the brain. Functional mobility was assessed using an obstacle course. Vestibular contribution to balance was measured using Neurocom Sensory Organization Test 5. Behavioral measures were assessed pre-HDBR, and 0, 8 and 12 days post-HDBR. Linear mixed models were used to test for effects of time, group, and group-by-time interactions. Family-wise error corrected VBM revealed significantly larger increases in GM volume in the right primary motor cortex in bed rest control subjects than in bed rest exercise subjects. No other significant group by time interactions in gray matter changes with bed rest were observed. Functional connectivity MRI revealed that the increase in connectivity during bed rest of the left putamen with the bilateral midsagittal precunes and the right cingulate gyrus was larger in bed rest control subjects than in bed rest exercise subjects. Furthermore, the increase in functional connectivity with bed rest of the right premotor cortex with the right inferior frontal gyrus and the right primary motor cortex with the bilateral premotor cortex was smaller in bed rest control subjects than in bed rest exercise subjects. Functional mobility performance was less affected by HDBR in exercise subjects than in control subjects and post HDBR exercise subjects recovered faster than control subjects. The group performance differences and GM changes were not related. Exercise in HDBR partially mitigates the adverse effect of HDBR on functional mobility, particularly during the post-bed rest recovery phase. In addition, exercise appears to result in differential brain structural and functional changes in motor regions such as the primary motor cortex, the premotor cortex and the putamen. Whether these central nervous system changes are related to motor behavioral changes including gait and balance warrants further research.

Koppelmans, V.↗

DOE OSTI2016

Investigating the age distribution of fracture discharge using multiple environmental tracers, Bedrichov Tunnel, Czech Republic

We investigated the transit time distribution (TTD) of discharge collected from fractures in the Bedrichov Tunnel, Czech Republic, using lumped parameter models and multiple environmental tracers. We then utilize time series of δ 18 O, δ 2 H and 3 H along with CFC measurements from individual fractures in the Bedrichov Tunnel of the Czech Republic to investigate the TTD, and the uncertainty in estimated mean travel time in several fracture networks of varying length and discharge. We also compare several TTDs, including the dispersion distribution, the exponential distribution, and a developed TTD which includes the effects of matrix diffusion. The effect of seasonal recharge is explored by comparing several seasonal weighting functions to derive the historical recharge concentration. We identify best fit mean ages for each TTD by minimizing the error-weighted, multi-tracer χ2 residual for each seasonal weighting function. We use this methodology to test the ability of each TTD and seasonal input function to fit the observed tracer concentrations, and the effect of choosing different TTD and seasonal recharge functions on the mean age estimation. We find that the estimated mean transit time is a function of both the assumed TTD and seasonal weighting function. Best fits as measured by the χ2 value were achieved for the dispersion model using the seasonal input function developed here for two of the three modeled sites, while at the third site, equally good fits were achieved with the exponential model and the dispersion model and our seasonal input function. The average mean transit time for all TTDs and seasonal input functions converged to similar values at each location. The sensitivity of the estimated mean transit time to the seasonal weighting function was equal to that of the TTD. These results indicated that understanding seasonality of recharge is at least as important as the uncertainty in the flow path distribution in fracture networks and that unique identification of the TTD and mean transit time is difficult given the uncertainty in the recharge function. But, the mean transit time appears to be relatively robust to the structural model uncertainty. The results presented here should be applicable to other studies using environmental tracers to constrain flow and transport properties in fractured rock systems.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2020

An extended radar relative calibration adjustment (eRCA) technique for higher-frequency radars and range–height indicator (RHI) scans

Abstract. This study extends the relative calibration adjustment technique for calibration of weather radars to higher-frequency radars as well as range–height indicator (RHI) scans. The calibration of weather radars represents one of the most dominant sources of error for their use in a variety of fields including quantitative precipitation estimation and model comparisons. While most weather radars are routinely calibrated, the frequency of calibration is often less than required, resulting in miscalibrated time periods. While full absolute calibration techniques often require the radar to be taken offline for a period of time, there have been online calibration techniques discussed in the literature. The relative calibration adjustment (RCA) technique uses the statistics of the ground clutter surrounding the radar as a monitoring source for the stability of calibration but has only been demonstrated to work at S- and C-band for plan-position indicator (PPI) scans at a constant elevation. In this work the RCA technique is modified to work with higher-frequency radars, including Ka-band cloud radars. At higher frequencies the properties of clutter can be much more variable. This work introduces an extended clutter selection procedure that incorporates the temporal stability of clutter and helps to improve the operational stability of RCA for relatively higher-frequency radars. The technique is also extended to utilize range–height scans from radars where the elevation is varied rather than the azimuth. These types of scans are often utilized with research radars to examine the vertical structure of clouds. The newly extended technique (eRCA) is applied to four Department of Energy Atmospheric Radiation Measurement (DOE ARM) weather radars ranging in frequency from C- to Ka-band. Cross comparisons of three co-located radars with frequencies C, X, and Ka at the ARM Cloud, Aerosol, and Complex Terrain Interactions (CACTI) site show that the technique can determine changes in calibration. Using an X-band radar at the ARM Eastern North Atlantic (ENA) site, we show how the technique can be modified to be more resilient to clutter fields that show increased variability, in this case due to sea clutter. The results show that this technique is promising for a posteriori data calibration and monitoring.

47 OTHER INSTRUMENTATION↗

DOE OSTI2018

Parton Dynamics Inferred from High-Mass Drell-Yan Dimuons Induced by 120 GeV p+D Interactions

Fermilab Experiment 906/SeaQuest (E906/SeaQuest) is the latest in a well established tradition of studying leptoproduction from the annihilation of a quark and anti-quark, known as the Drell-Yan process. The broad goal of E906/SeaQuest is measuring various properties of nucleon structure in order to learn more about quarks and Quantum Chromodynamics (QCD), the mathematical description of the strong force. The present work investigated violations of the Lam-Tung relation between virtual photon polarization and quark and lepton angular momentum. The violation of Lam-Tung can be explained as the signature of quark-nucleon spin-orbit coupling through the use of the Transverse-Momentum-Dependent (TMD) framework, which assumes that the initial transverse momentum of quarks is smaller than the hard scattering scale, but also non-negligible. An analysis of the angular moments in Drell-Yan collected by E906/SeaQuest was performed with four different configurations in order to estimate the systematic errors attributed to each correction. After correction for background and error propagation, the final extraction of the azimuthal moment excluding contributions from the trigger was ν = 0.151 ± 0.88(stat.) ± 0.346(syst.) at an average transverse momentum of 0.87 ± 0.50 GeV/c and an average dimuon mass of 5.48 ± 0.70 GeV. In the future, the magnitude of the systematic errors on the extraction could potentially be reduced by improving the quality of the trigger efficiency calculation, improving the intensity dependent event reconstruction efficiency, considering the changes in acceptance due to a beam shift relative to the E906/SeaQuest spectrometer, and improving the modeling of background.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1977

Industrial Productivity

NASTRAN is an offshoot of the computer-design technique used in construction of airplanes and spacecraft. [n this technique engineers create a mathematical model of the aeronautical or space vehicle and "fly" it on the ground by means of computer simulation. The technique enables them to study performance and structural behavior of a number of different designs before settling on the final configuration and proceeding with construction. From this base of aerospace experience, NASA-Goddard developed the NASTRAN general purpose computer program, which offers an exceptionally wide range of analytic capability with regard to structures. NASTRAN has been applied to autos, trucks, railroad cars, ships, nuclear power reactors, steam turbines, bridges, and office buildings. NASA-Langley provides program maintenance services regarded as vital by many NASTRAN users. NASTRAN is essentially a predictive tool. It takes an electronic look at a computerire$.dedgn and reports how the structure will react under a great many different conditions. It can, for example, note areas where high stress levels will occur-potential failure points that need strengthening. Conversely, it can identify over-designed areas where weight and material might be saved safely. NASTRAN can tell how pipes stand up under strong fluid flow, how metals are affected by high temperatures, how a building will fare in an earthquake or how powerful winds will cause a bridge to oscillate. NASTRAN analysis is quick and inexpensive. It minimizes trial-and-error in the design process and makes possible better, safe, lighter structures affording large-scale savings in development time and materials. Some examples of the broad utility NASTRAN is finding among industrial firms are shown on these pages.

Source record↗

DOE OSTI2022

Real-time equation-of-motion CC cumulant and CC Green’s function simulations of photoemission spectra of water and water dimer

Newly developed coupled-cluster (CC) methods enable simulations of ionization potentials and spectral functions of molecular systems in a wide range of energy scales ranging from core-binding to valence. This paper discusses the results obtained with the real-time equation-of-motion CC cumulant (RT-EOM-CC) approach and CC Green’s function (CCGF) approaches in applications to the water and water dimer molecules. Here, we compare the ionization potentials obtained with these methods for the valence region with the results obtained with the coupled-cluster with singles, doubles, and perturbative triples formulation as a difference of energies for N and N – 1 electron systems. All methods show good agreement with each other. They also agree well with the experiment with errors usually below 0.1 eV for the ionization potentials. We also analyze unique features of the spectral functions, associated with the position of satellite peaks, obtained with the RT-EOM-CC and CCGF methods employing single and double excitations, as a function of the monomer OH bond length and the proton transfer coordinate in the dimer. Finally, we analyze the impact of the basis set effects on the quality of calculated ionization potentials and find that the basis set effects are less pronounced for the augmented-type sets.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2020

Fast and Universal Kohn-Sham Density Functional Theory Algorithm for Warm Dense Matter to Hot Dense Plasma

Understanding many processes, e.g., fusion experiments, planetary interiors, and dwarf stars, depends strongly on microscopic physics modeling of warm dense matter and hot dense plasma. This complex state of matter consists of a transient mixture of degenerate and nearly free electrons, molecules, and ions. This regime challenges both experiment and analytical modeling, necessitating predictive ab initio atomistic computation, typically based on quantum mechanical Kohn-Sham density functional theory (KS-DFT). However, cubic computational scaling with temperature and system size prohibits the use of DFT through much of the warm dense matter regime. A recently developed stochastic approach to KS-DFT can be used at high temperatures, with the exact same accuracy as the deterministic approach, but the stochastic error can converge slowly and it remains expensive for intermediate temperatures (< 50 eV). Here we have developed a universal mixed stochastic-deterministic algorithm for DFT at any temperature. This approach leverages the physics of KS-DFT to seamlessly integrate the best aspects of these different approaches. We demonstrate that this method significantly accelerated self-consistent field calculations for temperatures from 3 to 50 eV, while producing stable molecular dynamics and accurate diffusion coefficients.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗