Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “empirical validation”

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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
NASA NTRS1997

Seasonal Phytoplankton Dynamics in the Eastern Tropical Atlantic

The coastal zone color scanner (CZCS) that operated aboard the Nimbus 7 satellite provided extensive coverage of phytoplankton pigment concentrations in the surface waters of the eastern tropical Atlantic (ETA) from March 1979 to February 1980 and coincided with four major research cruises to this region. Total primary production within the ETA (5 deg N-10 deg S, 25 deg W-10 deg E) was determined from CZCS pigment estimates and an empirical algorithm derived from concurrent in situ data taken along 4 deg W that relates near-surface chlorophyll concentration and integrated primary production. We estimated an average annual production for the ETA of 2.3 Gt C/yr with an associated 3.5-fold seasonal variation in the magnitude of this production. We describe the principal physical mechanisms controlling seasonal phytoplankton dynamics within the ETA and propose that in addition to seasonal change in the thermocline depth, one must also consider changes in the depth of the equatorial under current. An extensive validation effort indicates that the standard CZCS global products are a conservative estimate of pigment concentrations in ETA surface waters. Significant underestimates by the CZCS global products were observed in June and July which we attributed, in part, to aerosol correction errors and, more importantly, to errors caused by a significant reduction in the concentration of near-surface dissolved organic matter that resulted from strong equatorial upwelling.

Monger, Bruce↗
NASA NTRS2012

Observations of CMEs and Models of the Eruptive Corona

It is now realized that coronal mass ejections (CMEs) are the most energetic phenomenon in the heliosphere. Although early observations (in the 1970s and 19805) revealed most of the properties of CMEs, it is the extended and uniform data set from the Solar and Heliospheric Observatory (SOHO) mission that helped us consolidate our knowledge on CMEs. The Solar Terrestrial Relations Observatory (STEREO) mission has provided direct confirmation of the three-dimensional structure of CMEs. The broadside view provided by the STEREO coronagraphs helped us estimate the width of the halo CMEs and hence validate CME cone models. Current theoretical ideas on the internal structure of CMEs suggest that a flux rope is central to the CME structure, which has considerable observational support both from remote-sensing and in-situ observations. The flux-rope nature is also consistent with the post-eruption arcades with high-temperature plasma and the charge states observed within CMEs arriving at Earth. The quadrature observations also helped us understand the relation between the radial and expansion speeds of CMEs, which were only known from empirical relations in the past. This paper highlights some of these results obtained during solar cycle 23 and 24 and discusses implications for CME models.

Gopalswamy, Nat↗
DOE OSTI2025

Development of mechanistic, microstructure-informed BISON models for fission product-induced failure mechanisms in advanced nuclear fuels

The US continues to prioritize commercial deployment of advanced reactors—particularly those utilizing U-Zr metallic and TRISO particle fuels. Fission product-induced failure mechanisms for these systems include fuel–cladding chemical interaction in metallic fuel rods, which threatens cladding integrity, and Pd penetration in TRISO particles, which degrades SiC layer properties. Empirical models can be applied within the bounds of existing irradiation databases, but their utility is limited when considering new designs or investigating fundamental material behaviors. The Nuclear Energy Advanced Modeling and Simulation Program is therefore developing mechanistic models for these behaviors, which are expected to aid in fuel design, assist in development of failure mitigation strategies, and provide support for qualification and licensing. These efforts leverage multiscale capabilities to capture the microstructural features and processes that govern these behaviors. This talk summarizes past and ongoing efforts to develop and validate these models for the BISON fuel performance code.

FCCI↗
DOE OSTI2021

Empirical Evaluation of GPS Clock Accuracy for Isochronous Droop-Based Inverters: Preprint

Grid-forming inverters rely upon an internal reference to regulate the grid’s voltage and frequency. Without a reliable reference, issues such as loss of synchronism, frequency instability, and unequal power sharing can occur. One proposed internal phase reference for grid-forming inverters is generated by the controller’s clock, but this is affected by nonidealities in the crystal oscillator leading to time drift in the phase reference. To improve the controller’s clock accuracy, the Global Positioning System (GPS) can be used to calibrate the internal clock. Previously, GPS-based clock calibration was explored theoretically and various methods were developed for mitigating frequency/phase drift at the power system level. This paper experimentally validates the direct performance of GPS time-calibration at the controller level and justifies its use in maintaining an accurate phase reference for droop-based grid-forming inverters.

clock drift↗
NASA NTRS2012

Development of Modeling Capabilities for Launch Pad Acoustics and Ignition Transient Environment Prediction

This paper presents development efforts to establish modeling capabilities for launch vehicle liftoff acoustics and ignition transient environment predictions. Peak acoustic loads experienced by the launch vehicle occur during liftoff with strong interaction between the vehicle and the launch facility. Acoustic prediction engineering tools based on empirical models are of limited value in efforts to proactively design and optimize launch vehicles and launch facility configurations for liftoff acoustics. Modeling approaches are needed that capture the important details of the plume flow environment including the ignition transient, identify the noise generation sources, and allow assessment of the effects of launch pad geometric details and acoustic mitigation measures such as water injection. This paper presents a status of the CFD tools developed by the MSFC Fluid Dynamics Branch featuring advanced multi-physics modeling capabilities developed towards this goal. Validation and application examples are presented along with an overview of application in the prediction of liftoff environments and the design of targeted mitigation measures such as launch pad configuration and sound suppression water placement.

West, Jeff↗
DOE OSTI2023

ORNL_AISD-Ex: Quantum chemical prediction of UV/Vis absorption spectra for over 10 million organic molecules

We performed calculations of electronic excitation energies and associated oscillator strengths based on the time-dependent density-functional tight-binding (TD-DFTB) method [1]. The SMILES (Simplified molecular-input line-entry system) strings of the molecules from the AISD HOMO-LUMO database [2] were converted to a 3D atomistic structure and stored in a PDB file after preliminary geometry optimization using the Merck Molecular Force Field (MMFF94) in RDKit [3,4]. The primary information stored in the PDB file archive consists of Cartesian coordinates for each atom of the molecule in their 3D location in space, along with summary information about the structure, sequence, and experiment. We then performed molecular geometry optimization using the density-functional tight-binding (DFTB) method [5] in the electronic ground state, followed by single-point excited states calculations, as described below. We note that, since RDKit employs a random choice for the generation of molecular conformers, the molecular geometries obtained in this dataset could be different from the ones that were generated when the AISD HOMO-LUMO dataset was generated. The computed excitation energies and associated oscillator strengths can be converted to predict UV/Vis absorption spectra, where excitation energies correspond to absorption peak positions, and oscillator strengths are a good measure of the probability of absorption of visible or UV light in transitions between electronic ground and excited states. The conversion of SMILES strings to 3D Cartesian coordinates of fully DFTB-optimized molecules was successful for 10,502,904 out of 10,502,917 molecules. For these molecules, both geometry optimizations and excited states calculations were successful. The DFTB calculations did not complete for 13 molecules of the original AISD HOMO-LUMO dataset. We still provide information about the geometry of these molecules. The molecules are diverse for chemical compositions (which span 5 non-hydrogen elements: oxygen, carbon, nitrogen, fluorine, sulfur) and molecular size (the smallest molecule contains 5 non-hydrogen atoms, and the largest molecule contains 71 non-hydrogen atoms). The DFTB method [5] is an approximation to density functional theory (DFT), utilizing a minimal basis set in conjunction with a two-center approximation to the electronic Hamiltonian and overlap matrix elements. The DFTB total energy is the sum of an electronic and a repulsive energy contribution, and their calculation requires optimized electronic parameters and diatomic repulsive potential energy functions. All DFTB calculations were performed using the DFTB+ code [6] (version 21.2) and the wrapper for DFTB+ in the Atomic Simulation Environment (ASE) (version 3.22.1) [7], which performed an internal conversion of Cartesian coordinates from PDB to the .gen file format. For the geometry optimizations on the electronic ground state potential energy surface of the molecules, we have chosen the third-order DFTB (DFTB3) method [5c] and employed the matching 3ob set of electronic parameters and repulsive potentials [8]. The empirical γ-damping for hydrogen bond correction, and Grimme's D3 empirical dispersion correction with Becke-Johnson damping (D3(BJ)) [9] dispersion correction was included to improve the description of non-covalent interactions. For excited states single-point energy calculations, we employed the TD-DFTB method in conjunction with the DFTB2 method [5b] and the matching mio [5b,10] and halorg [11] parameter sets. We opted to request the simultaneous calculation of 50 excited states for singlet transition to investigate sufficient number of excited states, based on linear response theory using the Casida equation [Ref: T. A. Niehaus, S. Suhai, F. Della Sala, P Lugli, M. Elstner, G. Seifert, and Th. Frauenheim. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 63:085108, 2001] and the ARPACK diagonalizer [R. B. Lehoucq, D. C. Sorensen, and C. Yang. Arpack users guide: Solution of large-scale eigenvalue problems by implicitly restarted arnoldi methods, 1997. 46, 51]. The dataset contains 1001 tar.gz files. Tar files are named as “ornl_aisd_ex_1.tar.gz†through “ornl_aisd_ex_1000.tar.gzâ€. Additionally, the 13 failed molecules are in “ornl_aisd_ex_unprocessed.tar.gzâ€. Except for the tar files listed below, each tar file contains 10,500 molecules. Tar files numbered 34, 121, 128, 352, 360, 429, 495, 509, 518, 627, 676, 668, and 862 contain 10,499 molecules each. The last tar file numbered 1000 contains 13,417 molecules. The total size of the uncompressed dataset is over 283 Gigabytes. The code for calculating the electronic excitation energies and statistical analysis of the dataset is provided at the following GitLab repository: https://github.com/ORNL/Analysis-of-Large-Scale-Molecular-Datasets-with-Python Calculating the UV spectrum of a molecule requires performing 3 main operations: 1. Converting the smiles string representation of a molecule into a geometric structure where each atom is assigned XYZ coordinates. The geometric structure is written to the file smiles.pdb. 2. Using smiles.pdb to compute the relaxed geometry of the molecule, which corresponds with the position of the atoms at the position of equilibrium at the ground state. This generates the files band.out, detailed.out, and geo_end.gen. 3. Using geo_end.gen to calculate the UV spectrum of the molecule which is written into the file EXC.DAT. Every molecule in the dataset has its own directory. The files contained in each molecule directory are as follows: 1. geo_end.gen 2. detailed.out 3. band.out 4. EXC.DAT 5. smiles.pdb REFERENCES [1] Niehaus, T. A.; Suhai, S.; Della Salla, F.; Lugli, P.; Elstner, M.; Seifert, G.; Frauenheim, Th. Tight-binding approach to time-dependent density-functional response theory. Phys. Rev. B, 2001, 63, 085108/1-9. [2] Blanchard, A.; Gounley, J.; Metha, K.; Yoo, P.; Irle, S. AISD HOMO-LUMO. DOI: 10.13139/ORNLNCCS/1869409 [3] RDKit: Cheminformatics and Machine Learning Software. 2013, [http://www.rdkit.org] [4] Tosco, P.; Stiefl, N. and Landrum, G. Bringing the MMFF force field to the RDKit: implementation and validation. J Cheminform. 2014, 6, 1–4. [5] a) Porezag, D.; Frauenheim, T.; Kohler, T.; Seifert, G.; Kaschner, Construction of tight-binding-like potentials on the basis of density-functional theory: Application to carbon, R. Phys. Rev. B 1995, 51, 12947-12957; b) Elstner, M.; Porezag, D.; Jungnickel, G.; Elsner, J.; Haugk, M.; Frauenheim, Th.; Suhai, S.; Seifert, G.; Phys. Rev. B 1998, 58, 7260-7268; c) Gaus, M.; Cui, Q.; Elstner, M. DFTB3: Extension of the Self-Consistent-Charge Density-Functional Tight-Binding Method (SCC-DFTB), J. Chem. Theory Comput. 2011, 7, 931-948; d) Cui, Q.; Elstner, M. Density functional tight binding: values of semi-empirical methods in an ab initio era, Phys. Chem. Chem. Phys. 2014, 16, 14368-14377. [6] Hourahine, B. et al. DFTB+, a software package for efficient approximate density functional theory based atomistic simulations, J. Chem. Phys. 2020, 152, 124101/1-19. [7] Larsen, A. H. et al. The atomic simulation environment—a Python library for working with atoms. J. Phys.: Cond. Matter 2017, 29, 273002. [8] Kubillus, M.; Kubar, T.; Gaus, M.; Rezac, J.; Elstner, M. Parameterization of the DFTB3 Method for Br, Ca, Cl, F, I, K, and Na in Organic and Biological Systems, J. Chem. Theory Comput. 2015, 11, 332-342. [9] Brandenburg, J. G.; Grimme, S. Accurate Modeling of Organic Molecular Crystals by Dispersion-Corrected Density Functional Tight Binding (DFTB), J. Phys. Chem. Lett. 2014, 5, 1785−1789. [10] a) Niehaus, T. A.; Elstner, M.; Frauenheim, Th.; Suhai, S. Application of an approximate density-functional method to sulfur containing compounds. J. Mol. Struct.: THEOCHEM 2001, 541, 185-94; b) Elstner, M.; Hobza, P.; Frauenheim, Th.; Suhai, S.; Kaxiras, E. Hydrogen bonding and stacking interactions of nucleic acid base pairs: A density-functional-theory based treatment. J. Chem. Phys. 2001, 114, 5149-55. [11] Kubar, T.; Bodrog, Z.; Gaus, M.; Köhler, C.; Aradi, B.; Frauenheim, Th.; Elstner, M. Parametrization of the SCC-DFTB Method for Halogens. J. Chem. Theory Comput. 2013, 9, 2939-49.

36 MATERIALS SCIENCE↗
DOE OSTI2024

A three-dimensional model of terrain-induced updrafts for movement ecology studies

Abstract Background Spatially explicit simulation models of animal movements through the atmosphere necessarily require a representation of the spatial and temporal variation of atmospheric conditions. In particular, for movements of soaring birds that rely extensively on vertical updrafts to avoid flapping flight, accurate and reliable estimation of the vertical component of wind is critical. The interaction between wind and complex terrain shapes both the horizontal and vertical wind fields, highlighting the need to model the coupling between local terrain features and atmospheric conditions at scales relevant to animal movement. Methods In this work, we propose a new empirical model for estimating the orographic updraft field. The model is developed using computational fluid dynamics simulations of canonical atmospheric conditions over moderately complex terrain. To isolate buoyancy and thermal effects, and focus on terrain-induced effects, we use only simulations of a neutrally stratified atmosphere to develop the model. The model, which we name Engineering Vertical Velocity Estimator (EVVE), is simple to implement and is a function of the underlying terrain elevation map, the desired height above ground level (AGL), and wind conditions at a reference height (80 m). We validate the model with data from the Alaiz mountain (Spain) field campaign. Results Compared to observations, the proposed improved model estimates the updrafts at 120 m AGL with a mean error of 0.11 m/s ( $$\sigma =0.28$$ σ = 0.28 m/s), compared to 0.85 m/s ( $$\sigma =0.58$$ σ = 0.58 m/s) for its baseline. For typical land-based wind turbine hub heights of 80 m AGL, the proposed model has a mean error of 0.04 m/s ( $$\sigma =0.25$$ σ = 0.25 m/s), compared to baseline 0.54 m/s ( $$\sigma =0.45$$ σ = 0.45 m/s) estimations. We illustrate an application of the model in movement ecology by comparing simulated tracks and presence maps of golden eagles ( Aquila chrysaetos ) moving across two distinct landscapes. The tracks and presence maps are obtained using a simple heuristic-based movement model, with the updraft field given by the proposed model and a wind vector-based estimation approach that is currently in wide use in movement ecology studies of raptors and other soaring birds. Conclusions We highlight that movement model results can be sensitive to the underlying orographic updraft model, especially in studies of fine-scale movements in regions of complex topography. We suggest adopting the proposed model rather than the wind vector estimation method for studies of soaring bird movements.

17 WIND ENERGY↗
NASA NTRS1984

Models of auroral-zone conductances

The magnetosphere-ionosphere system is strongly coupled, with magnetospheric Birkeland currents feeding ionospheric Pedersen and Hall currents. Central to any computer simulation of this system is a detailed, valid conductivity model. An accurate conductivity model is also vital in order to infer Birkeland currents and electric field patterns from inversions of magnetometer chain data. Several recent attempts at constructing conductivity models are presented and their strengths and weaknesses discussed. Incoherent scatter radar measurements can determine height profiles of electron content, from which Pedersen and Hall conductances may be calculated. These yield excellent spatial and good temporal resolution; however, they are limited in field of view. A global pattern requires either 24 hours of data or a chain of stations. Synoptic empirical models (quantized by indices such as Kp or AE) typically are limited by their large bin size (1 deg invariant latitude x 1 hour MLT), and cannot reproduce arcs. Estimating conductivity globally from Dynamics Explorer auroral images is promising, and can yield reasonable time scales (of about 10 minutes); however, this procedure is still only now being tested.

Reiff, P. H.↗
NASA NTRS2008

Systematic Independent Validation of Inner Heliospheric Models

This presentation is the first in a series which will provide independent validation of community models of the outer corona and inner heliosphere. In this work we establish a set of measures to be used in validating this group of models. We use these procedures to generate a comprehensive set of results from the Wang- Sheeley-Arge (WSA) model which will be used as a baseline, or reference, against which to compare all other models. We also run a test of the validation procedures by applying them to a small set of results produced by the ENLIL Magnetohydrodynamic (MHD) model. In future presentations we will validate other models currently hosted by the Community Coordinated Modeling Center(CCMC), including a comprehensive validation of the ENLIL model. The Wang-Sheeley-Arge (WSA) model is widely used to model the Solar wind, and is used by a number of agencies to predict Solar wind conditions at Earth as much as four days into the future. Because it is so important to both the research and space weather forecasting communities, it is essential that its performance be measured systematically, and independently. In this paper we offer just such an independent and systematic validation. We report skill scores for the model's predictions of wind speed and IMF polarity for a large set of Carrington rotations. The model was run in all its routinely used configurations. It ingests line of sight magnetograms. For this study we generated model results for monthly magnetograms from the National Solar Observatory (SOLIS), Mount Wilson Observatory and the GONG network, spanning the Carrington rotation range from 1650 to 2068. We compare the influence of the different magnetogram sources, performance at quiet and active times, and estimate the effect of different empirical wind speed tunings. We also consider the ability of the WSA model to identify sharp transitions in wind speed from slow to fast wind. These results will serve as a baseline against which to compare future versions of the model

MacNeice, P. J.↗
NASA NTRS2019

Empirical Absolute Calibration Model for Multiple Pseudo-Invariant Calibration Sites

This work extends an empirical absolute calibration model initially developed for the Libya 4 Pseudo-Invariant Calibration Site (PICS) to five additional Saharan Desert PICS (Egypt 1, Libya 1, Niger 1, Niger 2, and Sudan 1), and demonstrates the efficacy of the resulting models at predicting sensor top-of-atmosphere (TOA) reflectance. It attempts to generate absolute calibration models for these PICS that have an accuracy and precision comparable to or better than the current Libya 4 model, with the intent of providing additional opportunities for sensor calibration. In addition, this work attempts to validate the general applicability of the model to other sites. The method uses Terra Moderate Resolution Imaging Spectroradiometer (MODIS) as the reference radiometer and Earth Observing-1 (EO-1) Hyperion image data to provide a representative hyperspectral reflectance profile of the PICS. Data from a region of interest (ROI) in an “optimal region” of 3% temporal, spatial, and spectral stability within the PICS are used for developing the model. The developed models were used to simulate observations of the Landsat 7 (L7) Enhanced Thematic Mapper Plus (ETM+), Landsat 8 (L8) Operational Land Imager (OLI), Sentinel 2A (S2A) MultiSpectral Instrument (MSI) and Sentinel 2B (S2B) MultiSpectral Instrument (MSI) from their respective launch date through 2018. The models developed for the Egypt 1, Libya 1 and Sudan 1 PICS have an estimated accuracy of approximately 3% and precision of approximately 2% for the sensors used in the study, comparable to the current Libya 4 model. The models developed for the Niger 1 and Niger 2 sites are significantly less accurate with similar precision.

Raut, Bipin↗
NASA NTRS2021

Comparison of soil dielectric mixing models for Soil Moisture Retrieval using SMAP Brightness Temperature over croplands in India

The accurate estimation of soil moisture (SM) using microwave remote sensing depends mostly on careful selection of retrieval parameters among which the soil dielectric mixing model is the important one. These models are often categorized into empirical, semi-empirical or volumetric based on their methodologies and input data requirements. To study in detail, the comparative performance of four dielectric mixing models -- Wang & Schmugge model, Hallikainen model, Dobson model and Mironov model were used with Soil Moisture Active Passive (SMAP) L-band brightness temperature and Single Channel Algorithm for SM retrieval over agricultural landscapes in India. The highest performance statistics combination in terms of Root Mean Square Error (RMSE), correlation coefficient (R^(2)) and percentage bias (PBIAS) against the concurrent in-situ SM measurements were calculated at the selected validation sites. The overall results indicate that the best performance was given by the Mironov model (RMSE = 0.07 cu. m/cu. m), followed by Wang & Schmugge model (RMSE = 0.08 cu. m3/cu. m), Hallikainen model (RMSE = 0.09 cu. m/cu. m), Dobson model (RMSE = 0.10 cu. m/cu. m) and original SMAP radiometer SM (RMSE = 0.12 cu. m/cu. m). Findings of this study provides important insights into application and performance of dielectric mixing models in mapping surface SM variations. This study also underlines the pivotal role of local conditions for SM retrieval which should be carefully included in the algorithms.

Swati Suman↗
DOE OSTI2023

Verification and benchmarking relativistic electron beam transport through a background gas

It is necessary to establish confidence in high-consequence codes containing an extensive suite of physics algorithms in the regimes of interest. Verification problems allow code developers to assess numerical accuracy and increase confidence that specific sets of model physics were implemented correctly in the code. The two main verification techniques are code verification and solution verification. In this work, we present verification problems that can be used in other codes to increase confidence in simulations of relativistic beam transport. Specifically, we use the general plasma code EMPIRE to model and compare with the analytical solution to the evolution of the outer radial envelope of a relativistic charged particle beam. Additionally, we also outline a benchmark test of a relativistic beam propagating through a vacuum and pressurized gas cell, and present the results between EMPIRE and the hybrid code GAZEL. Further, we discuss the subtle errors that were caught with these problems and detail lessons learned.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Uncertainty quantification for Multiphase-CFD simulations of bubbly flows: a machine learning-based Bayesian approach supported by high-resolution experiments

In this paper, we developed a machine learning-based Bayesian approach to inversely quantify and reduce the uncertainties of multiphase computational fluid dynamics (MCFD) simulations for bubbly flows. The proposed approach is supported by high-resolution two-phase flow measurements, including those by double-sensor conductivity probes, high-speed imaging, and particle image velocimetry. Local distributions of key physical quantities of interest (QoIs), including the void fraction and phasic velocities, are obtained to support the Bayesian inference. In the process, the epistemic uncertainties of the closure relations are inversely quantified while the aleatory uncertainties from stochastic fluctuations of the system are evaluated based on experimental uncertainty analysis. The combined uncertainties are then propagated through the MCFD solver to obtain uncertainties of the QoIs, based on which probability-boxes are constructed for validation. The proposed approach relies on three machine learning methods: feedforward neural networks and principal component analysis for surrogate modeling, and Gaussian processes for model form uncertainty modeling. The whole process is implemented within the framework of an open-source deep learning library PyTorch with graphics processing unit (GPU) acceleration, thus ensuring the efficiency of the computation. The results demonstrate that with the support of high-resolution data, the uncertainties of MCFD simulations can be significantly reduced. The proposed approach has the potential for other applications that involve numerical models with empirical parameters.

42 ENGINEERING↗
DOE OSTI2021

Analysis of Element Size and Element Type on SS304L Material Model Performance

Finite element solutions are only valid if the information inputted into the models is a valid assumption. This is why it is important to determine if the model behaves as expected especially in the material response. Material models that consider deformation in the elastic and plastic regions are desirable. As a material deforms plastically, hardening will occur until the ultimate strength is reached due to the reduction of the cross sectional area in a uniaxial tensile test environment. A phenomenological approximation of this behavior which is known as power law hardening (plh) can be expressed as σ = σ y + K $ϵ^{n}_{p}$ where K and n are empirically determined constants, and σ y is the proportional limit of the material. The first objective of this study is to verify the behavior for a calibrated plh material model for selective deviatoric hexahedral elements when applied to a different material stock. The Sierra SM input for this calibrated model can be found in the appendix. The second objective of this study is to determine the differences in failure characteristics for considerably large shell elements versus that of fine hexahedral elements.

36 MATERIALS SCIENCE↗
NASA NTRS1972

Relativistic stellar stability: An empirical approach

The PPN formalism which encompasses the post-Newtonian limit of nearly every metric theory of gravity is used to analyze stellar stability. This analysis enables one to infer, for any given gravitation theory, the extent to which post-Newtonian effects induce instabilities in white dwarfs, in neutron stars, and in supermassive stars. It also reveals the extent to which our current empirical knowledge of post-Newtonian gravity (based on solar-system experiments) actually guarantees that relativistic instabilities exist. In particular, it shows that: (1) for conservative theories of gravity, current solar-system experiments guarantee that the critical adiabatic index, for the stability of stars against radial pulsations exceeds the Newtonian value of 4/3 and (2) for nonconservative theories, current experiments do not permit any firm conclusion about the sign of the critical adiabatic index, and (3) in the PPN approximation to every metric theory, the standard Schwarzschild criterion for convection is valid.

Ni, W.↗
NASA NTRS1976

A mathematical model of an active control landing gear for load control during impact and roll-out

A mathematical model of an active control landing gear (ACOLAG) was developed and programmed for operation on a digital computer. The mathematical model includes theoretical subsonic aerodynamics; first-mode wing bending and torsional characteristics; oleo-pneumatic shock strut with fit and binding friction; closed-loop, series-hydraulic control; empirical tire force-deflection characteristics; antiskid braking; and sinusoidal or random runway roughness. The mathematical model was used to compute the loads and motions for a simulated vertical drop test and a simulated landing impact of a conventional (passive) main landing gear designed for a 2268-kg (5000-lbm) class airplane. Computations were also made for a simply modified version of the passive gear including a series-hydraulic active control system. Comparison of computed results for the passive gear with experimental data shows that the active control landing gear analysis is valid for predicting the loads and motions of an airplane during a symmetrical landing. Computed results for the series-hydraulic active control in conjunction with the simply modified passive gear show that 20- to 30-percent reductions in wing force, relative to those occurring with the modified passive gear, can be obtained during the impact phase of the landing. These reductions in wing force could result in substantial increases in fatigue life of the structure.

Mcgehee, J. R.↗
NASA NTRS2006

A New Empirical Thermospheric Density Model JB2006 using New Solar Indices

A new empirical atmospheric density model is developed using the CIRA72 (Jacchia 71) model as the basis for the diffusion equations. New solar indices based on orbit based sensor data are used for the solar irradiances in the extreme and far ultraviolet wavelengths. New exospheric temperature and semiannual density equations are employed to represent the major thermospheric density variations. Temperature correction equations are also developed for diurnal and latitudinal effects, and finally density correction factors are used for model corrections required at high altitude (1500-4000 km). The new model, Jacchia-Bowman 2006, is validated through comparisons of accurate daily density drag data previously computed for numerous satellite. For 400 km altitude the standard deviation of 16% for the standard Jacchia model is reduced to 10% for the new JB2006 model for periods of low geomagnetic storm activity.

Bowman, Bruce R.↗
NASA NTRS2007

Estimating the Soil Temperature Profile from a Single Depth Observation: A Simple Empirical Heatflow Solution

Two data sets of experimental field observations with a range of meteorological conditions are used to investigate the possibility of modeling near-surface soil temperature profiles in a bare soil. It is shown that commonly used heat flow methods that assume a constant ground heat flux can not be used to model the extreme variations in temperature that occur near the surface. This paper proposes a simple approach for modeling the surface soil temperature profiles from a single depth observation. This approach consists of two parts: 1) modeling an instantaneous ground flux profile based on net radiation and the ground heat flux at 5cm depth; 2) using this ground heat flux profile to extrapolate a single temperature observation to a continuous near surface temperature profile. The new model is validated with an independent data set from a different soil and under a range of meteorological conditions.

Holmes, Thomas↗