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 973 records · Page 54

Application of a laser scanner to three dimensional visual sensing tasks

The issues are described which are associated with using a laser scanner for visual sensing and the methods developed by the author to address them. A laser scanner is a device that controls the direction of a laser beam by deflecting it through a pair of orthogonal mirrors, the orientations of which are specified by a computer. If a calibrated laser scanner is combined with a calibrated camera, it is possible to perform three dimensional sensing by directing the laser at objects within the field of view of the camera. There are several issues associated with using a laser scanner for three dimensional visual sensing that must be addressed in order to use the laser scanner effectively. First, methods are needed to calibrate the laser scanner and estimate three dimensional points. Second, methods to estimate three dimensional points using a calibrated camera and laser scanner are required. Third, methods are required for locating the laser spot in a cluttered image. Fourth, mathematical models that predict the laser scanner's performance and provide structure for three dimensional data points are necessary. Several methods were developed to address each of these and has evaluated them to determine how and when they should be applied. The theoretical development, implementation, and results when used in a dual arm eighteen degree of freedom robotic system for space assembly is described.

Ryan, Arthur M.↗

Hyperspectral Imaging and Related Field Methods: Building the Science

The proposal requested funds for the computing power to bring hyperspectral image processing into undergraduate and graduate remote sensing courses. This upgrade made it possible to handle more students in these oversubscribed courses and to enhance CSES' summer short course entitled "Hyperspectral Imaging and Data Analysis" provided for government, industry, university and military. Funds were also requested to build field measurement capabilities through the purchase of spectroradiometers, canopy radiation sensors and a differential GPS system. These instruments provided systematic and complete sets of field data for the analysis of hyperspectral data with the appropriate radiometric and wavelength calibration as well as atmospheric data needed for application of radiative transfer models. The proposed field equipment made it possible to team-teach a new field methods course, unique in the country, that took advantage of the expertise of the investigators rostered in three different departments, Geology, Geography and Biology.

Goetz, Alexander F. H.↗

Detailed Chemical Abundances for a Benchmark Sample of M Dwarfs from the APOGEE Survey

Individual chemical abundances for 14 elements (C, O, Na, Mg, Al, Si, K, Ca, Ti, V, Cr, Mn, Fe, and Ni) are derived for a sample of M dwarfs using high-resolution, near-infrared H-band spectra from the Sloan Digital Sky Survey-IV/Apache Point Observatory Galactic Evolution Experiment (APOGEE) survey. The quantitative analysis included synthetic spectra computed with 1D LTE plane-parallel MARCS models using the APOGEE Data Release 17 line list to determine chemical abundances. The sample consists of 11 M dwarfs in binary systems with warmer FGK dwarf primaries and 10 measured interferometric angular diameters. To minimize atomic diffusion effects, [X/Fe] ratios are used to compare M dwarfs in binary systems and literature results for their warmer primary stars, indicating good agreement (<0.08 dex) for all studied elements. The mean abundance difference in primaries minus this work's M dwarfs is –0.05 ± 0.03 dex. It indicates that M dwarfs in binary systems are a reliable way to calibrate empirical relationships. A comparison with abundance, effective temperature, and surface gravity results from the APOGEE Stellar Parameter and Chemical Abundances Pipeline (ASPCAP) Data Release 16 finds a systematic offset of [M/H], T eff , log g = +0.21 dex, –50 K, and 0.30 dex, respectively, although ASPCAP [X/Fe] ratios are generally consistent with this study. The metallicities of the M dwarfs cover the range of [Fe/H] = –0.9 to +0.4 and are used to investigate Galactic chemical evolution via trends of [X/Fe] as a function of [Fe/H]. The behavior of the various elemental abundances [X/Fe] versus [Fe/H] agrees well with the corresponding trends derived from warmer FGK dwarfs, demonstrating that the APOGEE spectra can be used to examine Galactic chemical evolution using large samples of selected M dwarfs.

79 ASTRONOMY AND ASTROPHYSICS↗

Wind-tunnel preflight test program for aeroassist flight experiment

A proposed Aeroassist Flight Experiment (AFE) will simulate return from geosynchronous orbit to low-earth orbit to provide fundamental information about the hypersonic flow conditions surrounding aeroassisted orbital transfer vehicles (AOTV's). Future AOTV design must rely heavily on computational fluid dynamic (CFD) computer codes now under development, since the entry conditions are beyond the simulation capabilities of ground-based test facilities. However, existing codes and best available ground-based test facilities must provide for design of the AFE. The AFE ground-based test program provides data for calibration of existing and developing pertinent CFD codes. Confidence in code prediction capability increases with ability to predict forces and moments, pressure and heat-transfer distributions, shock shapes, and surface streamline directions. The ground-based test program provides these data over a wide range of hypersonic test parameters. This paper describes the test program, models, facilities, and some representative results.

Wells, William L.↗

Reduction and Analysis of Phosphor Thermography Data With the IHEAT Software Package

Detailed aeroheating information is critical to the successful design of a thermal protection system (TPS) for an aerospace vehicle. This report describes NASA Langley Research Center's (LaRC) two-color relative-intensity phosphor thermography method and the IHEAT software package which is used for the efficient data reduction and analysis of the phosphor image data. Development of theory is provided for a new weighted two-color relative-intensity fluorescence theory for quantitatively determining surface temperatures on hypersonic wind tunnel models; an improved application of the one-dimensional conduction theory for use in determining global heating mappings; and extrapolation of wind tunnel data to flight surface temperatures. The phosphor methodology at LaRC is presented including descriptions of phosphor model fabrication, test facilities and phosphor video acquisition systems. A discussion of the calibration procedures, data reduction and data analysis is given. Estimates of the total uncertainties (with a 95% confidence level) associated with the phosphor technique are shown to be approximately 8 to 10 percent in the Langley's 31-Inch Mach 10 Tunnel and 7 to 10 percent in the 20-Inch Mach 6 Tunnel. A comparison with thin-film measurements using two-inch radius hemispheres shows the phosphor data to be within 7 percent of thin-film measurements and to agree even better with predictions via a LATCH computational fluid dynamics solution (CFD). Good agreement between phosphor data and LAURA CFD computations on the forebody of a vertical takeoff/vertical lander configuration at four angles of attack is also shown. In addition, a comparison is given between Mach 6 phosphor data and laminar and turbulent solutions generated using the LAURA, GASP and LATCH CFD codes. Finally, the extrapolation method developed in this report is applied to the X-34 configuration with good agreement between the phosphor extrapolation and LAURA flight surface temperature predictions. The phosphor process outlined in the paper is believed to provide the aerothermodynamic community with a valuable capability for rapidly obtaining (4 to 5 weeks) detailed heating information needed in TPS design.

Merski, N. Ronald↗

Determination of Ballistic Limit for IM7/8552 Using LS-DYNA MAT 261

The goal of the NASA ACC High Energy Dynamic Impact Project is to determine the state of the art of dynamic fracture simulations for high velocity impact for composite fuselage shielding applications. Using a building block approach, several computational models considered under NASA ACC are being validated against test data, starting at unconfigured panels and progressing to configured panels under combined out-of-plane and in-plane loading due to ballistic impact. The computational models being evaluated in this project include MAT 162, MAT 213, MAT 261, SPG, and Peridynamics. In this paper, the simulation results using LS-Dyna Material MAT 261 are presented. In particular, a series of blind predictions for unconfigured panels were performed to determine the ballistic limit or V50 velocity. MAT 261 employs failure approach that is generally physically-based using fracture toughness criteria. The overall material model relies on typical ply-level stiffness properties, similar to MAT 162 and other composite continuum damage material models. The fracture toughness values are based on standard tests, and thus are not subject to extensive calibration. This approach is more efficient than performing extensive optimization studies for calibration of parameters. Also, this approach of relying on physical properties reduces the uncertainty of results, as questions concerning the quality and extent of the calibration studies is no longer relevant. However, it was found that carefully controlled coupon-level tests are needed to accurately obtain the required fracture toughness values. Additionally, it should be noted that there is one significant parameter in MAT 261 that does appear to require calibration, and that is the overall failure strain. This is the strain at which the element is deleted, and is not the same as the strain at which damage begins to accumulate. This failure strain is termed EFS (Effective Failure Strain), and is the maximum effective stain for element failure. Simulations have shown that this value will significantly affect impact response and failure. The paper presents the effect of this element failure strain parameter, along with possible uncertainties in fracture toughness values. With an adjusted appropriate value for EFS, it is seen that simulation results compare well with impact test data for predicted penetration velocity.

Byar, Alan↗

ICESat-2 Pointing Calibration and Geolocation Performance

ICESat-2 science requirements are dependent on the accurate real-time pointing control (i.e., geolocation control) and postprocessed geolocation knowledge of the laser altimeter surface returns. Prelaunch pointing alignment errors and postlaunch pointing alignment variation result in large geolocation errors that must be calibrated on orbit. In addition, the changing sun-orbit geometry causes thermal-mechanical forced laser frame alignment variations at the orbit period and trends from days, weeks, and months. Early mission analysis computed precise postlaunch laser beam alignment calibration. The alignment calibration was uploaded to the spacecraft and enabled the pointing control performance to achieve 4.4 ± 6.0 m, a significant improvement over the 45 m (1 σ) mission requirement. Laser frame alignment calibrations are used to reduce the alignment bias and time variation, as well as the orbital variation contributions to geolocation knowledge error from 6 to 1.7 m (1 σ). Relative beam alignment of the six beams is calibrated and shown to contribute between 0.5 ± 0.1 m and 2.4 ± 0.2 m of remaining geolocation knowledge error. Independent geolocation assessment based on comparison to high-resolution digital elevation models agrees well with the calibration geolocation error estimates. The analysis demonstrates the ICESat-2 mission is performing far better than its geolocation knowledge requirement of 6.5 m (1 σ) after the laser frame alignment bias variation and orbital variation calibrations have been applied. Remaining geolocation error is beam dependent and ranges from 2.5 m for beam 6 to 4.4 m for beam 2 (mean + 1 σ).

S B Luthcke↗

Computational general relativistic force-free electrodynamics: II. Characterization of numerical diffusivity

Scientific codes are an indispensable link between theory and experiment; in (astro-)plasma physics, such numerical tools are one window into the universe’s most extreme flows of energy. The discretization of Maxwell’s equations – needed to make highly magnetized (astro)physical plasma amenable to its numerical modeling – introduces numerical diffusion. It acts as a source of dissipation independent of the system’s physical constituents. Understanding the numerical diffusion of scientific codes is the key to classifying their reliability. It gives specific limits in which the results of numerical experiments are physical. We aim at quantifying and characterizing the numerical diffusion properties of our recently developed numerical tool for the simulation of general relativistic force-free electrodynamics by calibrating and comparing it with other strategies found in the literature. Our code correctly models smooth waves of highly magnetized plasma. We evaluate the limits of general relativistic force-free electrodynamics in the context of current sheets and tearing mode instabilities. We identify that the current parallel to the magnetic field (j ∥ ), in combination with the breakdown of general relativistic force-free electrodynamics across current sheets, impairs the physical modeling of resistive instabilities. We find that at least eight numerical cells per characteristic size of interest (e.g., the wavelength in plasma waves or the transverse width of a current sheet) are needed to find consistency between resistivity of numerical and of physical origins. High-order discretization of the force-free current allows us to provide almost ideal orders of convergence for (smooth) plasma wave dynamics. The physical modeling of resistive layers requires suitable current prescriptions or a sub-grid modeling for the evolution of j ∥ .

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling of the mechanical alloying process

Two programs have been developed to compute the dimensional and property changes that occur with repetitive impacts during the mechanical alloying process. The more sophisticated of the programs also maintains a running count of the fractions of particles present and from this calculates a population distribution. The programs predict powder particle size and shape changes in accord with the accepted stages of powder development during mechanical alloying of ductile species. They also predict hardness and lamellar thickness changes with processing, again with reasonable agreement with experimental results. These predictions offer support of the model (and thereby give insight into the possible 'actual' happenings of mechanical alloying) and hence allow refinement and calibration of the myriad aspects of the model. They also provide a vehicle for establishing control over the dimensions and properties of the output powders used for consolidation, thereby facilitating optimization of the consolidation process.

Maurice, D.↗

Deterministic Calibration of MFiX-PIC, Part 1: Settling Bed

The Particle-in-cell (PIC) numerical approach for modeling granular solids in fluid flow has gained significant interest in recent years. Valued for its often shorter time-to-solution, the PIC formulation relies on modeling statistical groupings of particles called parcels in cooperation with a solids stress model to affect local solids velocity. This is in contrast to the discrete element model (DEM) where every particle in a system is modelled individually and directly coupled to local solids velocity through Newtonian mechanics. The U.S. Department of Energy (DOE), National Energy Technology Laboratory (NETL) develops and maintains Multiphase Flow with Interphase eXchanges (MFiX), a collection of open-source computational fluid dynamics (CFD) solvers. Included in the MFiX suite are traditional two-fluid model (TFM) and DEM solvers, and a recently added PIC solver (NETL, 2021). In general, PIC methodologies offer an accuracy trade-off in lieu of computational speed; and therefore, it is important to assess the credibility of MFiX-PIC simulations. For this purpose, a systematic verification, validation and uncertainty quantification (VVUQ) effort was initiated at NETL to assess the new PIC solver

42 ENGINEERING↗

Hydrological analysis in watersheds with a variable-resolution global climate model (VR-CESM)

Traditionally, watershed-scale hydrology is simulated by distributed hydrological models with offline meteorological forcing data, or by regional regional climate models that link atmospheric and land hydrology interactions. Global climate models (GCMs) are rarely used to study watershed-scale hydrology due to the relatively coarse grid resolution, computationally expensive downscaling, and simplified physical processes. Recently, however, watershed-scale hydrology analysis has become possible in GCMs due to the development of variable-resolution GCMs that dynamically couple the hydrological processes between atmosphere and land systems at fine resolutions in selected regions and coarse resolution elsewhere. In this study, we used the variable-resolution Community Earth System Model (VR-CESM) with refined-resolution (14 km) in the western U.S. and eastern China to evaluate watershed-scale hydrology. Furthermore, we compared the historical VR-CESM outputs with gauge measurements and other hydrological models (e.g., National Water Model in the U.S.) and calibrated the subsurface runoff capacities in four mountainous watersheds. An RCP8.5 projection from 2007 to 2050 is used to estimate the impact of changing precipitation and snow climatology on watershed hydrology. We also analyzed the long-term runoff variability and the possibility of extreme runoff events as simulated by the VR-CESM. Although calibration is not possible in larger-scale watersheds, VR-CESM simulates the long-term annual variability of watersheds and provides insights on climate change impacts on hydrology. We conclude that refined-resolution VR-CESM can be used for watershed-scale hydrology analysis to understand water resources and water balance, in addition to traditional watershed-scale hydrological models. It enables hydrological analysis at multiple watersheds in one simulation and can help understand the two-way dynamics between land surface hydrology and atmospheric processes, and is especially practical for projecting climate change impacts. However, it is challenging to apply VR-CESM for hydrologic analysis in regulated watersheds as human factors (e.g., pumping, irrigation, water diversion) have not been fully addressed in VR-CESM.

54 ENVIRONMENTAL SCIENCES↗

Geopotential models of the Earth from satellite tracking, altimeter and surface gravity observations: GEM-T3 and GEM-T3S

Improved models of the Earth's gravitational field have been developed from conventional tracking data and from a combination of satellite tracking, satellite altimeter and surface gravimetric data. This combination model represents a significant improvement in the modeling of the gravity field at half-wavelengths of 300 km and longer. Both models are complete to degree and order 50. The Goddard Earth Model-T3 (GEM-T3) provides more accurate computation of satellite orbital effects as well as giving superior geoidal representation from that achieved in any previous GEM. A description of the models, their development and an assessment of their accuracy is presented. The GEM-T3 model used altimeter data from previous satellite missions in estimating the orbits, geoid, and dynamic height fields. Other satellite tracking data are largely the same as was used to develop GEM-T2, but contain certain important improvements in data treatment and expanded laser tracking coverage. Over 1300 arcs of tracking data from 31 different satellites have been used in the solution. Reliable estimates of the model uncertainties via error calibration and optimal data weighting techniques are discussed.

Lerch, F. J.↗

A Geometric Approach to Modeling Microstructurally Small Fatigue Crack Formation: Simulation and Prediction of Crack Nucleation in AA 7075-T651 - 2

The objective of this paper is to develop further a framework for computationally modeling microstructurally small fatigue crack growth in AA 7075-T651 [1]. The focus is on the nucleation event, when a crack extends from within a second-phase particle into a surrounding grain, since this has been observed to be an initiating mechanism for fatigue crack growth in this alloy. It is hypothesized that nucleation can be predicted by computing a non-local nucleation metric near the crack front. The hypothesis is tested by employing a combination of experimentation and nite element modeling in which various slip-based and energy-based nucleation metrics are tested for validity, where each metric is derived from a continuum crystal plasticity formulation. To investigate each metric, a non-local procedure is developed for the calculation of nucleation metrics in the neighborhood of a crack front. Initially, an idealized baseline model consisting of a single grain containing a semi-ellipsoidal surface particle is studied to investigate the dependence of each nucleation metric on lattice orientation, number of load cycles, and non-local regularization method. This is followed by a comparison of experimental observations and computational results for microstructural models constructed by replicating the observed microstructural geometry near second-phase particles in fatigue specimens. It is found that orientation strongly influences the direction of slip localization and, as a result, in uences the nucleation mechanism. Also, the baseline models, replication models, and past experimental observation consistently suggest that a set of particular grain orientations is most likely to nucleate fatigue cracks. It is found that a continuum crystal plasticity model and a non-local nucleation metric can be used to predict the nucleation event in AA 7075-T651. However, nucleation metric threshold values that correspond to various nucleation governing mechanisms must be calibrated.

Hochhalter, Jake D.↗

Unraveling the Mechanics: Computational Modeling of Residual Stresses in U10Mo Fuel with Aluminum Cladding

LEU (Low Enriched Uranium) plate-type fuel elements consist of a high-density, low-enrichment U–Mo alloy-based fuel foil encapsulated in an aluminum cladding and is fabricated through the Hot Isostatic Pressing (HIP) technique. Understanding the mechanism of possible failure modes of nuclear fuel is critical to mitigate potential consequences. One of the major contributing factors in various failure modes is stress and it greatly affects the integrity under temperature, pressure, and irradiation. These stresses could originate from various sources, including manufacturing, and operation, and can lead to deformation, cracking, or even complete failure of fuel elements. This study focuses on the development of a computational model that accurately predicts the residual stresses generated in the U10Mo fuel during fabrication. It has been observed that cladding creep has a substantial impact on the residual stresses in the U10Mo fuel post the HIP fabrication process. Furthermore, during the HIP bonding process the fuel plate system is heated to a temperature of 560 oC and as a result, the aluminum cladding transitions from (Al 6061-T6 to -O). This presents a challenge in capturing the change in material properties accurately. Therefore, a temperature dependent creep model such as hyperbolic sine creep model is considered to estimate the creep properties of the Aluminum cladding. The proposed calibrated creep model accurately predicts the residual stresses in the U10Mo fuel foil and agrees well with the experiments.

42 ENGINEERING↗

Tundra SGS – Critical Importance of Site-Specific Injectivity Testing for Carbon Dioxide Injection Well Simulation Calibration

Paper presented at 16th International Conference on Greenhouse Gas Control Technologies (GHGT-16), October 23–27, 2022, Lyon, France. Brine injectivity tests were conducted on both the Broom Creek (BC) and Deadwood (DW) storage formations in one of the cored stratigraphic test wells prior to finalized simulations. The injectivity tests entailed injection of brine into the formations, followed by shutting in the wells. The shut-in duration represents the fall-off period which was used to carry out pressure transient analyses of the tests. The analysis of injectivity tests indicated that core-indicated and modeled permeability distributions should be conservatively multiplied by 2.5× and 5× in the BC and DW geologic models, respectively, prior to numerical simulations of CO 2 injection using Computer Modelling Group’s (CMG’s) GEM equation-of-state reservoir simulator. The scaler permeability multipliers allowed the fluid flow simulations to fully represent well test-measured injection well rates in the respective formations.

20 FOSSIL-FUELED POWER PLANTS↗

Continuum-Level Modeling of Li-Ion Battery SEI by Upscaling Atomistically Informed Reaction Mechanisms

Understanding and controlling solid-electrolyte interphase (SEI) formation to stabilize cell performance is a significant challenge for next-generation Li-ion battery technologies. In recent years, computational modeling has become an essential tool in providing fundamental insights into SEI properties and dynamics. However, neither atomistic nor continuum-level approaches alone can capture the complexities of SEI chemistry across all relevant length and time scales. In this work, a continuum-level model is developed that is informed by reaction mechanisms obtained from first-principle calculations. The atomistically informed continuum-level model is used to understand electrolyte degradation, including the decomposition of ethylene carbonate (EC), ethyl methyl carbonate (EMC), and fluoroethylene carbonate (FEC). The model presented here is the most chemically complex continuum-level SEI model in the literature to date. The SEI model is calibrated against experimental irreversible leakage currents and shows qualitative agreement with expected SEI growth trends. The model framework is expected to accelerate fundamental understanding of SEI formation, facilitate mechanism development feedback, and dynamically interact with experimental insights.

continuum-level model↗

ICESat-2 Pointing Calibration and Geolocation Performance

ICESat-2 science requirements are dependent on the accurate real-time pointing control (i.e. geolocation control) and post-processed geolocation knowledge of the laser altimeter surface returns. Pre-launch pointing alignment errors and post-launch pointing alignment variation result in large geolocation errors that must be calibrated on orbit. In addition, the changing sun-orbit geometry causes thermal-mechanical forced laser frame alignment variations at the orbit period and trends from days, weeks and months. Early mission analysis computed precise post-launch laser beam alignment calibration. The alignment calibration was uploaded to the spacecraft and enabled the pointing control performance to achieve 4.4 ± 6.0 m, a significant improvement over the 45 m (1 σ) mission requirement. Laser frame alignment calibrations are used to reduce the alignment bias and time variation, as well as the orbital variation contributions to geolocation knowledge error from 6 m to 1.7 m (1 σ). Relative beam alignment of the six beams is calibrated and shown to contribute between 0.5 ± 0.1 m and 2.4 ± 0.2 m of remaining geolocation knowledge error. Independent geolocation assessment based on comparison to high-resolution digital elevation models agrees well with the calibration geolocation error estimates. The analysis demonstrates the ICESat-2 mission is performing far better than its geolocation knowledge requirement of 6.5 m (1 σ) after the laser frame alignment bias variation and orbital variation calibrations have been applied. Remaining geolocation error is beam dependent and ranges from 2.5 m for beam 6 to 4.4 m for beam 2 (mean + 1 σ).

S. B. Luthcke↗

Joint estimation of biogeochemical model parameters from multiple experiments: A bayesian approach applied to mercury methylation

Here, to characterize complex biogeochemical systems, results from multiple experiments, where each targets a specific subprocess, are commonly combined. The resulting datasets are interpreted through the calibration of biogeochemical models for process inference and predictions. Commonly used calibration approaches of fitting datasets from individual experiments to subprocess models one at a time is prone to missing information shared between datasets and incomplete uncertainty propagation. We propose a Bayesian joint-fitting scheme addressing the above-mentioned concerns by jointly fitting all the available datasets, thus calibrating the entire biogeochemical model in one go using Markov Chain Monte Carlo (MCMC). The identification of null spaces in the parameter distributions from MCMC guided the simplification of certain subprocess models. For example, fast kinetic sorption was replaced by equilibrium sorption, and Monod demethylation was replaced by first-order demethylation. Joint fitting of datasets resulted in complete uncertainty propagation with parameter estimates informed by all available data.

54 ENVIRONMENTAL SCIENCES↗