Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computer model calibration”

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 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
NASA NTRS1979

IGOS improvements for Seasat

Modifications to Battelle's Interactive Graphics Orbit Selection (IGOS) computer program to assist in the planning and evaluation of the Seasat-A Scatterometer System (SASS) flight program were studied. To meet the planning needs of the LaRC Seasat-A Scatterometer team, the following features/modifications were implemented in IGOS: (1) display and specification of time increments in orbital passes represented by the cross-hatching of ground swaths; (2) addition of pass number annotations on the horizontal axis of the STPLNG and STPTOD plots; (3) modification of the sensor model to include more than two swaths associated with a single sensor to approximate the SASS cell pattern; (4) inclusion of down range and cross-track swath geometry to display the characteristic skewed SASS pattern; (5) addition of a swath schedule to allow the display of the SASS mode changes and to calibrate gaps; and (6) development of a set of commands to generate the detailed swath data from sensor characteristics and orbit/earth motion.

Warmke, J. M.↗
DOE OSTI2026

Reformulating a Gurson-based dynamic damage model and demonstrating improved predictive power and numerical robustness

The original Tepla (TEnsile PLAsticity) ductile damage model, based on the Gurson yield surface, has long been used to model damage evolution and material failure under dynamic loading. Unfortunately, Tepla suffered from mesh sensitivity, numerical instability, and limited predictive capability. Here, we theoretically reformulate Tepla to address these issues. We especially focus on the prediction of porosity, which is the key state variable used for modeling ductile damage, as compared to more easily measured surface velocities, which are at best an indirect measure of damage. Key model changes include separating the viscosity during volumetric void growth from underlying shear strength behavior and switching to an iterative bisection solver. The new Tepla is then calibrated on incipient spall experiments on half-hard copper and tantalum, which demonstrate its ability to simultaneously fit the model to recovered porosity distributions and measured surface velocities, a stringent test. Improved numerical behavior, such as greatly reduced mesh sensitivity, is also shown in those simulations. Finally, the new Tepla model is applied to several high-explosive loaded, sweeping wave experiments, showing the ability of the model to predict behavior on tests with significantly different loading conditions and histories than the calibration data.

36 MATERIALS SCIENCE↗
NASA NTRS2007

White-light Interferometry using a Channeled Spectrum: II. Calibration Methods, Numerical and Experimental Results

In the companion paper, [Appl. Opt. 46, 5853 (2007)] a highly accurate white light interference model was developed from just a few key parameters characterized in terms of various moments of the source and instrument transmission function. We develop and implement the end-to-end process of calibrating these moment parameters together with the differential dispersion of the instrument and applying them to the algorithms developed in the companion paper. The calibration procedure developed herein is based on first obtaining the standard monochromatic parameters at the pixel level: wavenumber, phase, intensity, and visibility parameters via a nonlinear least-squares procedure that exploits the structure of the model. The pixel level parameters are then combined to obtain the required 'global' moment and dispersion parameters. The process is applied to both simulated scenarios of astrometric observations and to data from the microarcsecond metrology testbed (MAM), an interferometer testbed that has played a prominent role in the development of this technology.

phase correction↗
DOE OSTI2021

Cholesky-based experimental design for Gaussian process and kernel-based emulation and calibration.

Gaussian processes and other kernel-based methods are used extensively to construct approximations of multivariate data sets. The accuracy of these approximations is dependent on the data used. This paper presents a computationally efficient algorithm to greedily select training samples that minimize the weighted L p error of kernel-based approximations for a given number of data. The method successively generates nested samples, with the goal of minimizing the error in high probability regions of densities specified by users. The algorithm presented is extremely simple and can be implemented using existing pivoted Cholesky factorization methods. Training samples are generated in batches which allows training data to be evaluated (labeled) in parallel. For smooth kernels, the algorithm performs comparably with the greedy integrated variance design but has significantly lower complexity. Numerical experiments demonstrate the efficacy of the approach for bounded, unbounded, multi-modal and non-tensor product densities. We also show how to use the proposed algorithm to efficiently generate surrogates for inferring unknown model parameters from data using Bayesian inference.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Poromechanical cohesive interface element with combined Mode I-II cohesive zone elastoplasticity for simulating fracture in fluid-saturated porous media

A combined Mode I-II cohesive zone (CZ) elasto-plastic constitutive model, and a two-dimensional (2D) cohesive interface element (CIE) are formulated and implemented at small strain within an ABAQUS User Element (UEL) for simulating 2D crack nucleation and propagation in fluid-saturated porous media. Here, the CZ model mitigates problems of convergence for the global Newton-Raphson solver within ABAQUS, which when combined with a viscous stabilization procedure allows for simulation of post-peak response under load control for coupled poromechanical finite element analysis, such as concrete gravity dam stability analysis. Verification examples are presented, along with a more complex ambient limestone-concrete wedge fracture experiment, water-pressurized concrete wedge experiment, and concrete gravity dam stability analyses. A calibration procedure for estimating the CZ parameters is demonstrated with the limestone-concrete wedge fracture process. For the water-pressurized concrete wedge fracture experiment it is shown that the inherent time-dependence of the poromechanical CIE analysis provides a good match with experimental force versus displacement results at various crack mouth opening rates, yet misses the pore water pressure evolution ahead of the crack tip propagation. This is likely a result of the concrete being partially-saturated in the experiment, whereas the finite element analysis assumes fully water saturated concrete. For the concrete gravity dam analysis, it is shown that base crack opening and associated water uplift pressure leads to a reduced Factor of Safety, which is confirmed by separate analytical calculations.

97 MATHEMATICS AND COMPUTING↗
NASA NTRS1979

The double seismic zone in downgoing slabs and the viscosity of the mesosphere

The seismic zone beneath several island arcs between about 100 and 200 km depth consists of an upper zone having down-dip compression and a lower zone having down-dip tension. Several numerical models of the Aleutian arc were computed to test the hypothesis that these double seismic zones are due to sagging of the slab under its own weight. This sagging occurs because the asthenosphere (between about 100 and 200 km) provides little support or resistance to the slab, which is supported from below by the more viscous mesosphere and from above by the lithosphere. The viscosity of the mesosphere was constrained to the interval between 0.25 x 10 to the 22nd and 0.5 x 10 to the 22nd P by noting that the slab would have mainly down-dip compression at higher viscosities and mainly down-dip tension at lower viscosities. The deviatoric stress in the slab and the fault plane between the slab and the island arc is about 200-300 bars (expressed as shear stress). The models were calibrated to the observed depth and gravity anomalies in the trench.

Sleep, N. H.↗
NASA NTRSDate not supplied

Evaluation of Spectral Band Adjustment Factors for Cross-Calibration of Visible Imagers

The NASA CERES EBAF dataset provides TOA SW and LW fluxes for long-term monitoring of the Earth’s energy balance and to validate climate models. The EBAF products, based on the Terra, Aqua, and NOAA20 CERES instrument observed radiances, rely on coincident measurements from the onboard MODIS or VIIRS imagers to determine cloud properties used for angular distribution model scene selection, which is required to convert the CERES observed radiances into flux values. Furthermore, CERES utilizes geostationary imager (GEO) hourly fluxes and clouds to infer the regional averaged daily flux in the SYN1deg product. A seamless transition of fluxes and clouds can only occur if the analogous MODIS, VIIRS, and GEO channels are properly inter-calibrated. The analogous band SRFs differ noticeably and require scene dependent spectral band adjustment factors (SBAF) for proper radiometric scaling between them. Given their disparate but overlapping SRFs, the coincident VIIRS I1 and M5 band TOA reflectance measurements provide the optimal opportunity to validate SBAFs over many surface and cloud conditions. The CERES project maintains SCIAMACHY, GOME-2, and Hyperion scene-stratified hyper-spectral reflectance measurements that can be convolved with sensor pair SRFs to compute the corresponding SBAF. This study highlights the SCIAMACHY, GOME-2, and Hyperion based SBAFs for different MODIS, VIIRS, and GEO inter-calibration targets, including all-sky tropical ocean (ATO), Libya-4, Dome-C, and deep convective clouds. To mitigate the observed sensor radiance fluctuations due to atmospheric absorption and aerosol variations, the Earth target observed radiances are correlated with multiple atmospheric parameters, such as precipitable water and ozone. The remaining M5 and I1 mean spectral band radiance difference should be resolved by the SBAF correction. The formulation of hyper-spectral sensor based SBAFs and validation methods can be verified in future utilizing high-accuracy, SI-traceable hyper-spectral measurements from CLARREO Pathfinder.

David R Doelling↗
DOE OSTI2023

Integrating LEO and GEO Observations: Toward Optimal Summertime Satellite Precipitation Retrieval

Abstract Reliable quantitative precipitation estimation with a rich spatiotemporal resolution is vital for understanding the Earth’s hydrological cycle. Precipitation estimation over land and coastal regions is necessary for addressing the high degree of spatial heterogeneity of water availability and demand, and for resolving the extremes that modulate and amplify hazards such as flooding and landslides. Advancements in computation power along with unique high spatiotemporal and spectral resolution data streams from passive meteorological sensors aboard geosynchronous Earth-orbiting (GEO) and low Earth-orbiting (LEO) satellites offer exciting opportunities to retrieve information about surface precipitation phenomena using data-driven machine learning techniques. In this study, the capabilities of U-Net–like architecture are investigated to map instantaneous, summertime surface precipitation intensity at the spatial resolution of 2 km. The calibrated brightness temperature products from the Global Precipitation Measurement (GPM) Microwave Imager (GMI) radiometer are combined with multispectral images (visible, near-infrared, and infrared bands) from the Advanced Baseline Imager (ABI) aboard the GOES-R satellites as main inputs to the U-Net–like precipitation algorithm. Total precipitable water and 2-m temperature from the Global Forecast System (GFS) model are also used as auxiliary inputs to the model. The results show that the U-Net–like algorithm can capture fine-scale patterns and intensity of surface precipitation at high spatial resolution over stratiform and convective precipitation regimes. The evaluations reveal the potential of extracting relevant, high spatial features over complex surface types such as mountainous regions and coastlines. The algorithm allows users to interpret the inputs’ importance and can serve as a starting point for further exploration of precipitation systems within the field of hydrometeorology.

Meteorology & Atmospheric Sciences↗
DOE OSTI2021

Building Intelligent Cyberinfrastructure to Learn Iteratively from both Observations and Models for Understanding Watershed Dynamics

Focal Area(s): Predictive modeling through the use of AI techniques and AI-derived model components; the use of AI and other tools to design a prediction system comprising of a hierarchy of models (e.g., AI-driven model/component/parameterization selection). Science Challenge: Watershed processes, such as the fate and transport of sediment, carbon and nutrients across landscapes and their fluxes to water bodies (e.g., streams, rivers and lakes), have important implications for global and regional carbon and nutrient dynamics, biogeochemical functioning of terrestrial ecosystems, and soil functions. The magnitude of lateral surface/subsurface transport and fluxes of sediment, carbon and nutrients are key factors controlling the vulnerability of watersheds to climate extremes such as droughts, wildfires, and floods. Recent field observations and other scientific evidence suggest that the magnitudes of lateral transport and fluxes of sediment, carbon and nutrients are governed primarily by the spatial and vertical heterogeneity of landscape and soil properties and by pedogenic processes. However, the current generation of land surface and watershed models do not mechanistically couple the terrestrial and hydrologic systems, nor do they represent sufficiently the spatial and vertical heterogeneity of land surface and subsurface properties. On the other hand, increasing complexity of coupled watershed and land surface models requires more data to parameterize, calibrate and validate. Remote sensing (RS) provides a means to acquire spatial data and characterize their heterogeneity at the watershed scale, overcoming a major limitation associated with conventional point measurements. To improve the representation of land-surface and surface/subsurface process coupling and sub-grid heterogeneity in watershed models, it is essential to build our predictive understanding by learning from both the multi-scale multi-process modeling and diverse multi-scale data while leveraging powerful artificial intelligence (AI) techniques.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

The Role of Subsurface Flow on Evapotranspiration: A Global Sensitivity Analysis

Abstract Water resources are impacted by water‐energy balance fluxes at the land surface, most notably evapotranspiration ( ET ), the largest component of surface energy balance. While integrated hydrologic models show promise in quantifying the nonlinear dynamics at this interface, the model results are plagued by parametric uncertainties. Given the high computational demand of running multiple parameter spaces of these models, little is known about how these uncertainties propagate into land surface processes. In this work, we perform a global sensitivity analysis by computing Sobol and AMAE indices using a surrogate model constructed with a polynomial chaos expansion to assess the impacts of the subsurface physical properties on ET . We do so by modeling a semisynthetic test case. Our results show that the effects of vertical hydraulic conductivity, porosity, and the water retention curve parameter Van Genuchten α mainly control ET . However, we note that while evaporation ( E ) shows behavior similar to the ET with high sensitivity to the parameters controlling the flow in the unsaturated zone, transpiration ( T ) is very sensitive to the saturated zone parameters and groundwater flow, especially during periods without rain. Our results also show that heterogeneities in land cover impact the flow in the saturated zone. Furthermore, our work demonstrates that reduced‐order models can be developed to make integrated models more accessible for rigorous sensitivity analyses and calibration purposes.

Maina, Fadji Zaouna↗
NASA NTRSDate not supplied

Effect of Roughness on Transition on the BOLT-1a Geometry in Supersonic Flow

Infrared thermography measurements on an 8.8% scale BOLT-1a model were made in the NASA Langley Probe Calibration Tunnel for freestream unit Reynolds numbers from 7.7 × 10^ -6 /m to 30.7 × 10^ -6 /m, a freestream Mach number of 3.5, and nominal angle of attack of 0 degrees. The model was printed from polycarbonate and tested before and after fine-grit sanding to evaluate the effect of surface roughness on transition. Global surface roughness measurements were obtained using a 3D optical profilometer. The temperature distributions on the upper and lower experimental surfaces of the test article were simultaneously captured using two infrared cameras. Prior to testing, the facility’s freestream environment was characterized. At a freestream unit Reynolds number of approximately 17 × 10^ -6 /m, the nozzle-wall boundary layer began to transition non-uniformly, inducing transition on the model. The boundary-layer state was determined by analyzing Stanton number distributions computed from the surface temperature measurements. The present supersonic conditions that were examined on BOLT-1a were found to display both similarities and disparities in comparison to previous hypersonic experiments.

boundary layer transition↗
DOE OSTI2021

A NICER View of the Massive Pulsar PSR J0740+6620 Informed by Radio Timing and XMM-Newton Spectroscopy

We report on Bayesian estimation of the radius, mass, and hot surface regions of the massive millisecond pulsar PSR J0740+6620, conditional on pulse-profile modeling of Neutron Star Interior Composition Explorer X-ray Timing Instrument event data. We condition on informative pulsar mass, distance, and orbital inclination priors derived from the joint North American Nanohertz Observatory for Gravitational Waves and Canadian Hydrogen Intensity Mapping Experiment/Pulsar wideband radio timing measurements of Fonseca et al. We use XMM-Newton European Photon Imaging Camera spectroscopic event data to inform our X-ray likelihood function. The prior support of the pulsar radius is truncated at 16 km to ensure coverage of current dense matter models. We assume conservative priors on instrument calibration uncertainty. We constrain the equatorial radius and mass of PSR J0740+6620 to be ${12.39}_{-0.98}^{+1.30}$ km and ${2.072}_{-0.066}^{+0.067}$ M ⊙ respectively, each reported as the posterior credible interval bounded by the 16% and 84% quantiles, conditional on surface hot regions that are non-overlapping spherical caps of fully ionized hydrogen atmosphere with uniform effective temperature; a posteriori, the temperature is ${\mathrm{log}}_{10}(T\,[{\rm{K}}])={5.99}_{-0.06}^{+0.05}$ for each hot region. All software for the X-ray modeling framework is open-source and all data, model, and sample information is publicly available, including analysis notebooks and model modules in the Python language. Furthermore, our marginal likelihood function of mass and equatorial radius is proportional to the marginal joint posterior density of those parameters (within the prior support) and can thus be computed from the posterior samples.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1991

Computation of a circular-to-rectangular transition duct flow field

This paper presents the results of a Computational Fluid Dynamics (CFD) calibration study of flow through a circular-to-rectangular transition duct. The design of these ducts is critical to the optimum performance of aircraft with rectangular exhaust nozzles, since these ducts transfer the flow from the gas turbine engine to the exhaust nozzle. Two duct inflow conditions are considered, the first with straight inflow and the second with swirling inflow. Both flows contain realistic wall boundary layers. The first case permits examination of the effects of the geometric transition on the flowfield, while the second case adds in the rotational flow effect that can be present from the gas turbine engine. The flowfields associated with the two cases have been shown, both with CFD and experiment, to have significant differences. The Navier-Stokes CFD code PARC with the Baldwin and Lomax turbulence model was used in this study. The turbulence model had to be modified for both cases in order to achieve accurate determination of the edge of the wall bounded vorticity layers and thus turbulent viscosity. The results of this calibration study will be valuable to aircraft designers who rely on CFD to assist in the design and evaluation of propulsion systems.

Sirbaugh, J. R.↗
DOE OSTI2025

Adaptive Client Selection in Federated Learning: A Network Anomaly Detection Use Case

Federated Learning (FL) has become a ubiquitous approach for training machine learning models on decentralized data, addressing the myriad privacy concerns inherent in traditional centralized methods. However, the efficiency of FL depends on effective client selection and robust privacy preservation mechanisms. Inadequate client selection may lead to suboptimal model performance, while insufficient privacy measures risk exposing sensitive data. This paper proposes a client selection framework for FL that integrates differential privacy and fault tolerance. Our adaptive approach dynamically adjusts the number of selected clients based on model performance and system constraints, ensuring privacy through calibrated noise addition. We evaluate our method on a network anomaly detection use case using the UNSW-NB15 and ROAD datasets. Results show up to a 7% increase in accuracy and a 25% reduction in training time compared to FedL2P. Moreover, we highlight the trade-offs between privacy budgets and model performance, with higher privacy budgets reducing noise and improving accuracy. Our fault tolerance mechanism, while causing a slight performance drop, enhances robustness to client failures. Statistical validation using Mann-Whitney U tests confirms the significance of these improvements (p < 0.05).

Marfo, William [University of Texas at El Paso,Dep↗
NASA NTRS1992

Computer simulations of comet- and asteroidlike bodies passing through the Venusian atmosphere: Preliminary results on atmospheric and ground shock effects

We have completed computer simulations that model shock effects in the venusian atmosphere caused during the passage of two cometlike bodies 100 m and 1000 m in diameter and an asteroidlike body 10 km in diameter. Our objective is to examine hypervelocity-generated shock effects in the venusian atmosphere for bodies of different types and sizes in order to understand the following: (1) their deceleration and depth of penetration through the atmosphere; and (2) the onset of possible ground-surface shock effects such as splotches, craters, and ejecta formations. The three bodies were chosen to include both a range of general conditions applicable to Venus as well as three specific cases of current interest. These calculations use a new multiphase computer code (DICE-MAZ) designed by California Research & Technology for shock-dynamics simulations in complex environments. The code was tested and calibrated in large-scale explosion, cratering, and ejecta research. It treats a wide range of different multiphase conditions, including material types (vapor, melt, solid), particle-size distributions, and shock-induced dynamic changes in velocities, pressures, temperatures (internal energies), densities, and other related parameters, all of which were recorded in our calculations.

Roddy, D.↗
NASA NTRS1983

Rye Canyon X-ray noise test: One-third octave-band data

Acoustic data were obtained for the 25 ft. diameter X-wing rotor model during performance testing of the rotor system in hover. Data collected at the outdoor whirl tower test facility with a twelve microphone array were taken for approximately 150 test conditions comprised of various combinations of RPM, blade pressure ratio (BPR), and blade angle of attack (collective). The three test parameters had four values of RPM from 404 to 497, twelve values of BPR from 1.0 to 2.1, and six values of collective from 0.0 deg to 8.5 deg. Fifteen to twenty seconds of acoustic data were reduced to obtain an average 1/3 octave band spectrum for each microphone for each test condition. The complete, as measured, 1/3 octave band results for all the acoustic data are listed. Another part of the X-wing noise test was the acoustic calibration of the Rye Canyon whirl tower bowl. Corrections were computed which, when applied to as measured data, yield estimates of the free field X-wing noise. The free field estimates provide a more realistic measure of the rotor system noise levels. Trend analysis of the three test parameters on noise level were performed.

Willshire, W. L., Jr.↗
DOE OSTI2021

Columnar data analysis with ATLAS analysis formats

Future analysis of ATLAS data will involve new small-sized analysis formats to cope with the increased storage needs. The smallest of these, named DAOD_PHYSLITE, has calibrations already applied to allow fast downstream analysis and avoid the need for further analysis-specific intermediate formats. This allows for application of the "columnar analysis" paradigm where operations are applied on a per-array instead of a per-event basis. We will present methods to read the data into memory, using Uproot, and also discuss I/O aspects of columnar data and alternatives to the ROOT data format. Furthermore, we will show a representation of the event data model using the Awkward Array package and present proof of concept for a simple analysis application.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Estimation of Arrivals on Green at Signalized Intersections Using Stop-Bar Video Detection

Across the world, traffic congestion is increasing with alarming rapidity. Traffic signal control effectiveness, in coordinated networks, is often investigated in relation to the type of vehicle arrivals at the signalized intersections. Recently, several transportation agencies have switched from traditional loop detectors to video detection. When video cameras are accompanied by computer vision, one can extract more information about traffic “dynamics” than by using traditional inductive loop detectors. Collecting arrival times of multiple vehicles after the first arrival at the stop-bar detector might be challenging when using inductive loop detectors (since after the first arrival, detector status is always occupied). However, emerging video detection systems allow tracking of each vehicle’s entrance time in the detection zone, departure time from the detection zone, and the type of vehicle. This information can be used to estimate vehicular arrival and departure times, which then can be fed into machine learning algorithms to estimate arrivals on green (AOG). However, such research ideas have not been documented so far. Thus, this paper presents an estimation model for AOG, which was developed using multigene genetic programming. A robust experimental dataset was collected from a highly calibrated and validated microsimulation model of an 11-intersection corridor in Chattanooga, TN. The results of the model’s performance analysis showed the high accuracy of the training-, testing-, and validation datasets. The practical benefit of this model is that it can be applied to estimate arrival types at intersections where only stop-bar video detection exists.

Engineering↗