Search NASA⌕ Search

SEARCH · Search NASA

Results for “Model initialization”

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,549 records · Page 86
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

Multi Resolution In-Situ Testing and Multiscale Simulation for Creep Fatigue Damage Analysis of Alloy 617

This report outlines the research activities that were carried out for the integrated experimental and simulation investigation of creep-fatigue damage mechanism and life prediction of Nickel-based alloy, Inconel 617 at high temperatures (950° and 850°). First, a novel experimental design using a hybrid control technique is proposed. The newly developed experimental technique can generate different combinations of creep and fatigue damage by changing the experimental design parameters. Next, detailed imaging analysis and statistical data analysis are performed to quantify the failure mechanisms of the creep fatigue of alloy 617 at high temperatures. It is observed that the creep damage is directly associated with the internal voids at the grain boundaries and the fatigue damage is directly related to the surface cracking. It is also observed that the classical time fraction approach does not has a good correlation with the experimental observed damage features. An effective time fraction parameter is seen to have an excellent correlation with the material microstructural damage. Thus, a new empirical damage interaction diagram is proposed based on the experimental observations. Following this, a macro level viscoplastic model coupled with damage is developed to simulate the stress/strain response under creep fatigue loadings. A damage rate function based on the hysteresis energy and creep energy is proposed to capture the softening behavior of the material and a good correlation with life prediction and material hysteresis behavior is observed. The simulation work is extended to include the microstructural heterogeneity. A crystal plasticity finite element model considering isothermal and large deformation conditions at the microstructural scale has been developed for fatigue, creep-fatigue as well as creep deformation and rupture at high temperature. The model considers collective dislocation glide and climb of the grains and progressive damage accumulation of the grain boundaries. The glide model incorporates a slip resistance evolution model that characterizes the solute-drag creep effects and can capture well the stress-strain and stress time response of fatigue and creep-fatigue tests at various strain ranges and hold times. In order to accurately capture the creep strains that accumulate particularly at relatively low stress levels, a dislocation climb model has been incorporated into the crystal plasticity modeling framework. The dislocation climb model parameters are calibrated and verified through experimental creep tests performed at 950°. In addition, a cohesive zone model has been fully implemented in the context of the crystal plasticity finite element model to capture the intergranular creep damage. The parameters of the cohesive zone model have been calibrated using available experimental data. The numerical simulations illustrate the capability of the proposed model in capturing damage initiation and growth under creep loads as compared to the experimental observations. The microscale analysis sheds light on the crack initiation sites and propagation patterns within the microstructure. The model is also utilized to investigate the hybrid-controlled creep-fatigue tests and has been found to capture reasonably well the stress-strain response with different hold times and hold stress magnitudes.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Utilizing physics-based input features within a machine learning model to predict wind speed forecasting error

Machine learning is quickly becoming a commonly used technique for wind speed and power forecasting. Many machine learning methods utilize exogenous variables as input features, but there remains the question of which atmospheric variables are most beneficial for forecasting, especially in handling non-linearities that lead to forecasting error. This question is addressed via creation of a hybrid model that utilizes an autoregressive integrated moving-average (ARIMA) model to make an initial wind speed forecast followed by a random forest model that attempts to predict the ARIMA forecasting error using knowledge of exogenous atmospheric variables. Variables conveying information about atmospheric stability and turbulence as well as inertial forcing are found to be useful in dealing with non-linear error prediction. Streamwise wind speed, time of day, turbulence intensity, turbulent heat flux, vertical velocity, and wind direction are found to be particularly useful when used in unison for hourly and 3 h timescales. The prediction accuracy of the developed ARIMA–random forest hybrid model is compared to that of the persistence and bias-corrected ARIMA models. The ARIMA–random forest model is shown to improve upon the latter commonly employed modeling methods, reducing hourly forecasting error by up to 5 % below that of the bias-corrected ARIMA model and achieving an R 2 value of 0.84 with true wind speed.

17 WIND ENERGY↗
NASA NTRS1983

Collisions and merging of disk galaxies

The stable disks required as initial states for an experiment are produced by embedding the disk (which is visible) in a halo (invisible). The disk is a luminous tracer that represents only 1 percent of the total mass of the system. The stability of the initial state of the model is confirmed experimentally. Collisions with various combinations of initial orbital energy and angular momentum and of disk orientations are investigated. Two initially parabolic cases with different impact parameters are followed in time until the interactions lead to a merger. The interpenetration of the two galaxies gives rise to a contraction and subsequent disruption of the entire system. A variety of responses is found for the visible disks, ranging from stretched-out nearly linear features to rapidly propagating ringlike patterns. It is noted that these forms depend critically on initial disk orientations but do not depend strongly on the sense of disk rotation. The ringlike patterns are found even for collisions with significant impact parameters.

Smith, B. F.↗
NASA NTRS2002

2-D Magnetohydrodynamic Modeling of A Pulsed Plasma Thruster

Experiments are being performed on the NASA Marshall Space Flight Center (MSFC) MK-1 pulsed plasma thruster. Data produced from the experiments provide an opportunity to further understand the plasma dynamics in these thrusters via detailed computational modeling. The detailed and accurate understanding of the plasma dynamics in these devices holds the key towards extending their capabilities in a number of applications, including their applications as high power (greater than 1 MW) thrusters, and their use for producing high-velocity, uniform plasma jets for experimental purposes. For this study, the 2-D MHD modeling code, MACH2, is used to provide detailed interpretation of the experimental data. At the same time, a 0-D physics model of the plasma initial phase is developed to guide our 2-D modeling studies.

Thio, Y. C. Francis↗
NASA NTRS2002

An Experimental Study of a Pulsed Electromagnetic Plasma Accelerator

Experiments are being performed on the NASA Marshall Space Flight Center (MSFC) pulsed electromagnetic plasma accelerator (PEPA-0). Data produced from the experiments provide an opportunity to further understand the plasma dynamics in these thrusters via detailed computational modeling. The detailed and accurate understanding of the plasma dynamics in these devices holds the key towards extending their capabilities in a number of applications, including their applications as high power (greater than 1 MW) thrusters, and their use for producing high-velocity, uniform plasma jets for experimental purposes. For this study, the 2-D MHD modeling code, MACH2, is used to provide detailed interpretation of the experimental data. At the same time, a 0-D physics model of the plasma initial phase is developed to guide our 2-D modeling studies.

Thio, Y. C. Francis↗
NASA NTRS2004

Stressed Oxidation Life Prediction for C/SiC Composites

The residual strength and life of C/SiC is dominated by carbon interface and fiber oxidation if seal coat and matrix cracks are open to allow oxygen ingress. Crack opening is determined by the combination of thermal, mechanical and thermal expansion mismatch induced stresses. When cracks are open, life can be predicted by simple oxidation based models with reaction controlled kinetics at low temperature, and by gas phase diffusion controlled kinetics at high temperatures. Key life governing variables in these models include temperature, stress, initial strength, oxygen partial pressure, and total pressure. These models are described in this paper.

Levine, Stanley R.↗
NASA NTRSDate not supplied

Application of Advanced Earth Observations and Model Simulations to Improve Air Quality Monitoring in the Hindu-Kush-Himalayan Region

Air pollution is a serious environmental health concern in the Hindu Kush Himalayan (HKH) region of south-central Asia, as rapid industrialization and population growth have led to increased anthropogenic emissions from transportation, residential, industrial, energy, and biomass burning sources. Natural emissions from dust and forest fires are additional sources of air pollutants that can exacerbate air quality in the region. The combination of the complex pollutant mixtures and atmospherically stable weather conditions during the winter monsoon can visibility reductions and hazardous air quality from persistent haze episodes. The Kathmandu Valley is especially vulnerable to extreme haze issues due to the surrounding mountains that restrict air movement and retains pollutants in the atmosphere. This study uses state-of-the-art satellite observations and modeling capabilities in conjunction with ground-based networks to provide a comprehensive data toolkit for advancing air quality monitoring and forecasting decisions in the HKH region. The toolkit includes new generation satellite observations from the TROPOspheric Monitoring Instrument (TROPOMI), Geostationary Environment Monitoring Spectrometer (GEMS), and Advanced Meteorological Imager, which provide high spatiotemporal information on NO2, HCHO, SO2, O3, and aerosol optical depth (AOD). Particulate matter with diameters less than 2.5 micrometers (PM2.5) are derived from the satellite-retrieved AOD using ground-based observations and forecast model data. The satellite observations are also used to initialize and constrain forecast model systems designed for the HKH region. This talk will highlight the performance of the air quality toolkit for enhancing decision-making processes during exceptional air quality events in the region. Note: Presentation includes additional attachment of full presentation with sound and animation (best when viewed as slide show) with runtime of 15 min 32 secs

Aaron Naeger↗
DOE OSTI2017

Type Ib and IIb Supernova Progenitors in Interacting Binary Systems

We explore properties of Type Ib and IIb SN progenitors that are produced by stable mass transfer in binary systems using a new grid of stellar evolution models from an initial primary mass in the range of 10--18 M{sub ⊙} at solar and Large Magellanic Cloud metallicities. We find that blue and yellow supergiant SN IIb progenitors (e.g., of SN 2008ax, SN 2011dh, and SN 2016gkg) have a hydrogen envelope mass less than about 0.15 M{sub ⊙}, mostly resulting from early Case B mass transfer with relatively low initial masses and/or low metallicity. Red supergiant (RSG) SN IIb progenitors (e.g., of SN 1993J, SN 2013df) are produced via late Case B mass transfers and have a more massive hydrogen envelope (M{sub H,env}>0.15 M{sub ⊙}). SN Ib progenitors are predominantly produced by early Case B mass transfers. Our models predict that SN IIb progenitors are systematically more luminous in the optical (−8.0≲M{sub V}≲−5.0) than the majority of SN Ib progenitors (M{sub V}≳−5.0) for our considered initial mass range. However, the optically bright progenitor of SN Ib iPTF13bvn (i.e., M{sub V}≃−6.5) can be well explained by a relatively low-mass progenitor with a final mass of ∼3.0 M{sub ⊙}. The event rate of blue and yellow SN IIb progenitors would increase as metallicity decreases, while the event rate of SN Ib progenitors would decrease instead. By contrast, the population of RSG SN IIb progenitors would not be significantly affected by metallicity.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Spectra and flow of light nuclei in relativistic heavy ion collisions at energies available at the BNL Relativistic Heavy Ion Collider and at the CERN Large Hadron Collider

Within the framework of the coalescence model based on the phase-space distributions of protons and neutrons generated from the iebe-vishnu hybrid model with ampt initial conditions, we study the spectra and elliptic flow of deuterons and helium-3 in relativistic heavy ion collisions at the BNL Relativistic Heavy Ion Collider (RHIC) and the CERN Larger Hadron Collider (LHC). Results from our model calculations for Au+Au collisions at √ s NN = 200 GeV at RHIC and Pb+Pb collisions at √ s NN = 2.76 TeV at the LHC are compared with available experimental data. Good agreements are generally seen between theoretical results and experimental data, except that the calculated yield of helium-3 in Pb+Pb collisions at √ s NN = 2.76 TeV underestimates the data by about a factor of 2. Possible reasons for these discrepancies are discussed. Here, we also make predictions on the spectra and elliptic flow of deuterons and helium-3 in Pb+Pb collisions at √ s NN = 5.02 TeV that are being studied at LHC.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2010

The Impact of Ocean Observations in Seasonal Climate Prediction

The ocean provides the most significant memory for the climate system. Hence, a critical element in climate forecasting with coupled models is the initialization of the ocean with states from an ocean data assimilation system. Remotely-sensed ocean surface fields (e.g., sea surface topography, SST, winds) are now available for extensive periods and have been used to constrain ocean models to provide a record of climate variations. Since the ocean is virtually opaque to electromagnetic radiation, the assimilation of these satellite data is essential to extracting the maximum information content. More recently, the Argo drifters have provided unprecedented sampling of the subsurface temperature and salinity. Although the duration of this observation set has been too short to provide solid statistical evidence of its impact, there are indications that Argo improves the forecast skill of coupled systems. This presentation will address the impact these different observations have had on seasonal climate predictions with the GMAO's coupled model.

Rienecker, Michele↗
NASA NTRS1988

The effects of rotational flow, viscosity, thickness, and shape on transonic flutter dip phenomena

The transonic flutter dip phenomena on thin airfoils, which are employed for propfan blades, is investigated using an integrated Euler/Navier-Stokes code and a two degrees of freedom typical section structural model. As a part of the code validation, the flutter characteristics of the NACA 64A010 airfoil are also investigated. In addition, the effects of artificial dissipation models, rotational flow, initial conditions, mean angle of attack, viscosity, airfoil thickness and shape on flutter are investigated. The results obtained with a Euler code for the NACA 64A010 airfoil are in reasonable agreement with published results obtained by using transonic small disturbance and Euler codes. The two artificial dissipation models, one based on the local pressure gradient scaled by a common factor and the other based on the local pressure gradient scaled by a spectral radius, predicted the same flutter speeds except in the recovery region for the case studied. The effects of rotational flow, initial conditions, mean angle of attack, and viscosity for the Reynold's number studied seem to be negligible or small on the minima of the flutter dip.

Reddy, T. S. R.↗
NASA NTRS1988

The effects of rotational flow, viscosity, thickness, and shape on transonic flutter dip phenomena

The transonic flutter dip phenomena on thin airfoils, which are employed for propfan blades, is investigated using an integrated Euler/Navier-Stokes code and a two degrees of freedom typical section structural model. As a part of the code validation, the flutter characteristics of the NACA 64A010 airfoil are also investigated. In addition, the effects of artificial dissipation models, rotational flow, initial conditions, mean angle of attack, viscosity, airfoil thickness, and shape on flutter are investigated. The results obtained with an Euler code for the NACA 64A010 airfoil are in reasonable agreement with published results obtained by using transonic small disturbance and Euler codes. The two artifical dissipation models, one based on the local pressure gradient scaled by a common factor and the other based on the local presure gradient scaled by a spectral radius, predicted the same flutter speeds except in the recovery region for the case studied. The effects of rotation flow, initial conditions, mean angle of attack, and viscosity for the Reynolds number studied seem to be negligible or small on the minima of the flutter dip.

Reddy, T. S. R.↗
NASA NTRS1993

A nonsingular model of the open magnetosphere

We present a modified version of the Toffoletto and Hill (1989) open magnetosphere model that incorporates a tail-like interconection field with a discontinuity 10 represent the slow-mode expansion fan that defines the high-latitude tail magnetopause. (The interconnection field is defined as the perturbation on an initially closed magnetosphere model to make it open.) The expansion fan controls the open field line region in the tail, and the intersection of the fan with the tail current sheet is, by design, the x line. The new interconnection field allows greater control of the tail field structure; in particular, it enables us to eliminate the nightside mapping singularity that occurs in previous models when the interplanetary magnetic field is nonsouthward. Also, in contrast to earlier models, the far tail x line extends farther downstream on the flanks than in the center of the tail, consistent with observations.

Toffoletto, F. R.↗
DOE OSTI2021

Modeling of High-Temperature Corrosion of Zirconium Alloys Using the eXtended Finite Element Method (X-FEM)

Oxidation modeling in modern nuclear fuel performance codes is currently limited by the lack of coupling with mechanics, thus preventing proper description of how high-temperature oxidation impacts mechanical properties. This is mostly due to the fact that the finite difference formalism adopted in corrosion models is incompatible with the direct coupling with mechanics in the finite element modeling employed in modern nuclear fuel performance codes. In this study, a physically based zirconium alloy corrosion model called the Coupled-Current Charge Compensation (C4) model, which was initially developed for operating temperature conditions, has been updated to include high-temperature corrosion in order to provide additional critical information (e.g., oxygen concentration profile) under loss-of-coolant accident (LOCA) conditions—information lacking in existing empirical models. The C4 model was implemented in the MOOSE finite-element framework developed at Idaho National Laboratory, enabling it to be used in the BISON nuclear fuel performance code based on the MOOSE framework. To precisely track the different interfaces at a relatively low computational cost, the eXtended Finite Element Method (X-FEM) was applied in MOOSE. The model’s results were compared to those of existing empirical models as well as metallographic analysis of high-temperature oxidized Zircaloy-4 coupons. Oxygen diffusivities in the a and ß phases resulting from this comparison closely agree with those found in the literature. The C4 model implemented with X-FEM in MOOSE now has the capability to accurately predict oxide, oxygen-stabilized a, and prior ß phase layer growth kinetics under isothermal exposure at high temperature (1000–1500°C). Furthermore, in contrast with the empirical models, the C4 model accounts for the finite thickness of the fuel cladding. It can predict the oxygen concentration profile evolution through the whole cladding, enabling evaluation of the remaining ductile thickness—a crucial variable for modeling the mechanical behavior of the fuel cladding under LOCA. Furthermore, this implementation allows direct coupling with mechanics, at a low computing cost, using finite-element-based nuclear fuel performance codes such as BISON.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Experiments and simulations of isochorically heated warm dense carbon foam at the Texas Petawatt Laser

An experimental and simulation study of warm dense carbon foams at ambient density (n e ~10 21 cm –3 ) is presented. This study of isochorically heated foams is motivated by their potential application in carbon-atmosphere white-dwarf envelopes, where there are modeling uncertainties due to the equation of state. The foams are heated on an approximately picosecond time scale with a laser-accelerated proton beam. The cooling and expansion of the heated foams can be modeled with appropriately initialized radiation-hydrodynamics codes; xRAGE code is used in this work. The primary experimental diagnostic is the streaked optical pyrometer, which images a narrow band of radiation from the rear surface of the heated material. Presented are xRAGE modeling results for both solid aluminum targets and carbonized resorcinol-formaldehyde foam targets, showing that the foam appears to cool slowly on the pyrometer because of partial transparency. So that simulations of cooling foam are processed properly, it is necessary to account for finite optical depth in the photosphere calculation, and the methods for performing that calculation are presented in depth.

36 MATERIALS SCIENCE↗
DOE OSTI2022

1D Modeling for the Development of a Digital Twin for Hydropower Systems with Applications for a Testing Rig

Modeling and simulation constitute two important parts in constructing a sensible digital twin to mimic the dynamics of hydropower systems. Because of the physical nature of a hydropower system, the basic modeling should cover the water flow systems from the reservoir to the penstock (inlet water pipes), penstock to hydro turbine, and hydro turbine to generator and from the linkage of the hydropower systems to the grid. This report describes initial attempts to model these dynamic components, including the formulation of the linearized state space model for hydro turbine systems using the well-known six-coefficients linearization method and the formulation of the voltage and power dynamical models of a synchronous generator. To validate the accuracy of the proposed model, data collected from Norwegian University of Science and Technology were used to obtain relevant parameters for the model, and the desired simulation results were obtained. Additionally, neural network modeling and learning were also developed and applied to model the generation torque and water flow rate subsystems to demonstrate a potential learning scheme, which can be used in the development of a digital twin.

13 HYDRO ENERGY↗
DOE OSTI2023

Predicting encapsulant delamination in photovoltaic modules bridging photochemical reaction kinetics and fracture mechanics

Abstract Photovoltaic (PV) modules are subjected to environmental stressors (UV exposure, temperature, and humidity) that cause degradation within the encapsulant and its interfaces with adjacent glass and cell substrates. To save experimental time and to enable long‐term assessment with intensive degradation only taking place after many years, the development of predictive models is indispensable. Previous works have modeled the delamination of the ethylene vinyl acetate (EVA) encapsulant/glass and encapsulant/cell interfaces under field aging conditions with fundamental photochemical degradation reactions that lead to molecular scission and loss of interfacial adhesion, characterized by the fracture resistance, G c . However, these models were fundamentally limited in that the following aspects were not incorporated: (i) molecular crosslinking in the field, (ii) synergistic autocatalytic interactions of degradation mechanisms, (iii) connection between degraded encapsulant structure and its mechanical properties, and (iv) rigorous treatment of the plasticity contribution to G c with finite element models. Here, we present a time‐dependent multiscale model that addresses these limitations and is applicable to a wide range of encapsulants and interfaces. For the reference EVA encapsulant and its interfaces with the glass and cell, the presented model predicts an initial rise in G c in the first 3 years of field aging from crosslinking, then a subsequent sharp decline from degradation mechanisms. We used nanoindentation to measure the changes in EVA mechanical properties over exposure time to tune the model parameters. The model predictions of G c and mechanical properties match with experimental data and show an improvement compared to previous models. The model can even predict switches in failure interfaces, such as the observed EVA/cell to EVA/glass transition. We also conducted a sensitivity analysis study by varying the degradation and crosslinking kinetic parameters to demonstrate their effects on G c . Model extensions to polyolefin elastomer‐ and silicone‐encapsulants and their interfaces are also demonstrated.

Liu, Kuan↗