Search NASA⌕ Search

SEARCH · Search NASA

Results for “model determination”

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 73 records · Page 4

VALIDATION, VERIFICATION, AND CALIBRATION THROUGH A CAUSAL LENS

This paper presents an alternative method based on causal inference to perform validation, verification, and calibration of simulation models. While classical validation and verification approaches focus on the identification of the associations between data elements using statistical and machine learning methods, the novel methods in this paper focus instead on the identification of causal relationships between data elements. Statistical and machine learning-based approaches are strictly data-driven, meaning that they provide quantitative comparison measures between datasets without explicitly considering the hypotheses behind them. This can lead to the erroneous conclusion that, if two data sets are close enough, then the models that generated them are similar. In addition, when experimental and simulated data differ to an extent that fails to meet the acceptance criteria, calibration techniques are used to tweak simulation model parameters to reduce the gap between the two types of data. This produces the false expectation that a simulation model will match reality. The methods presented in this paper move away from these strictly data-driven methods for validation and calibration toward more robust, model-driven methods based on causal inference. Causal inference aims to identify the possible mechanisms that might have generated data. Thus, this analysis targets the prediction of the effects when one (or more) of the identified mechanisms are altered. There are many approaches to identify, quantify and illustrate causal relationships. For the scope of this paper, directed graphs are employed as causal models. If the directed graph lacks cycles it is known as a directed acyclic graph (DAG). A node in such a graph represents an observed data element while a directed edge connecting two nodes represents a causal relationship between two variables. The developed causal methods are designed to extract causal models from simulation models and from experimental data. Causal models capture the causal relationships between data elements (e.g., simulated and experimental data). In this context, validation and verification are performed by comparing causal models. The proposed approach does not only inform system analysts on how a simulation model matches real-world data, but also identifies elements of the simulation model that should be revised when discrepancies between simulation and experimental data are observed. Through these causal methods, analysts have a means to identify the portion of the model equation(s) that are behind an edge connecting two variables. Hence, once the structural differences between causal models have been determined, model calibration can occur by changing only those model parameters that impact the identified causal relationships.

97 MATHEMATICS AND COMPUTING↗

Spacecraft Solar Particle Event (SPE) Shielding: Shielding Effectiveness as a Function of SPE model as Determined with the FLUKA Radiation Transport Code

Analysis of both satellite and surface neutron monitor data demonstrate that the widely utilized Exponential model of solar particle event (SPE) proton kinetic energy spectra can seriously underestimate SPE proton flux, especially at the highest kinetic energies. The more recently developed Band model produces better agreement with neutron monitor data ground level events (GLEs) and is believed to be considerably more accurate at high kinetic energies. Here, we report the results of modeling and simulation studies in which the radiation transport code FLUKA (FLUktuierende KAskade) is used to determine the changes in total ionizing dose (TID) and single-event environments (SEE) behind aluminum, polyethylene, carbon, and titanium shielding masses when the assumed form (i. e., Band or Exponential) of the solar particle event (SPE) kinetic energy spectra is changed. FLUKA simulations have fully three dimensions with an isotropic particle flux incident on a concentric spherical shell shielding mass and detector structure. The effects are reported for both energetic primary protons penetrating the shield mass and secondary particle showers caused by energetic primary protons colliding with shielding mass nuclei. Our results, in agreement with previous studies, show that use of the Exponential form of the event

Koontz, Steve↗

A comprehensive model to determine the effects of temperature and species fluctuations on reaction rates in turbulent reacting flows

A computationally-viable model describing the interaction between fluid-mechanical turbulence and finite-rate combustion reactions, principally in high-speed flows was developed. Chemical kinetic mechanisms, complete and global, were developed describing the finite rate reaction of fuels of interest to NASA. These fuels included principally hydrogen and silane, although a limited amount of work involved hydrocarbon fuels as well.

Chinitz, W.↗

A comprehensive model to determine the effects of temperature and species fluctuations on reaction rates in turbulent reacting flows

A principal element to be derived from modeling turbulent reacting flows is an expression for the reaction rates of the various species involved in any particular combustion process under consideration. A temperature-derived most-likely probability density function (pdf) was used to describe the effects of temperature fluctuations on the Arrhenius reaction rate constant. A most-likely bivariate pdf described the effects of temperature and species concentrations fluctuations on the reaction rate. A criterion is developed for the use of an "appropriate" temperature pdf. The formulation of models to calculate the mean turbulent Arrhenius reaction rate constant and the mean turbulent reaction rate is considered and the results of calculations using these models are presented.

Foy, E.↗

A comprehensive model to determine the effects of temperature and species fluctuations on reactions in turbulent reacting flows

The joint probability distribution function (pdf), which is a modification of the bivariate Gaussian pdf, is discussed and results are presented for a global reaction model using the joint pdf. An alternative joint pdf is discussed. A criterion which permits the selection of temperature pdf's in different regions of turbulent, reacting flow fields is developed. Two principal approaches to the determination of reaction rates in computer programs containing detailed chemical kinetics are outlined. These models represent a practical solution to the modeling of species reaction rates in turbulent, reacting flows.

Antaki, P. J.↗

A technique using a nonlinear helicopter model for determining trims and derivatives

A technique is described for determining the trims and quasi-static derivatives of a flight vehicle for use in a linear perturbation model; both the coupled and uncoupled forms of the linear perturbation model are included. Since this technique requires a nonlinear vehicle model, detailed equations with constants and nonlinear functions for the CH-47B tandem rotor helicopter are presented. Tables of trims and derivatives are included for airspeeds between -40 and 160 knots and rates of descent between + or - 10.16 m/sec (+ or - 200 ft/min). As a verification, the calculated and referenced values of comparable trims, derivatives, and linear model poles are shown to have acceptable agreement.

Ostroff, A. J.↗

Application of two creep fatigue life models for the prediction of elevated temperature crack initiation of a nickel base alloy

Cyclic Damage Accumulation (CDA) and Total Strain-Strain Range Partitioning (TS-SRP) models for predicting the creep-fatigue crack initiation life of high temperature alloys are presented. The models differ in their fundamental assumptions regarding the controlling parameters for fatigue crack initiation and in the amount of data required to determine model constants. The CDA model represents a ductility exhaustion approach and uses stress quantities to calculate the cyclic fatigue damage. The TS-SRP model is based on the use of total mechanical strain and earlier concepts of the Strain Range Partitioning Method. Both models were applied to a well controlled fatigue data set at a high temperature nickel base alloy, B1900 + Hf, tested at 1600 F and 1800 F. The tests were divided into a baseline data set required to determine model constants and a verification data set for evaluation of the predictive capability of the models. Both models correlated the baseline data set to within factors of two in life, and predicted the verification data set to within a factor of three or better. In addition, sample calculations to demonstrate the application of each model and discusions of the predictive capabilities and areas requiring further development are presented.

Moreno, V.↗

Computer Aided Modeling to Determine the Effectiveness of Resistive Exercises as Countermeasures for Bone Mineral Density Loss

Due to the loss of gravitational loading, astronauts have a tendency to lose bone mineral density in their lumbar spine and lower extremities on orbit. NASA requires astronauts to perform exercises during space flight to help reduce the amount of demineralization. To test these exercises on earth, 17 week bed rest studies are conducted that consist of specific diet and exercise regimes. Developing a finite element model of these exercises will help to quantify the stress distribution imposed by of each of these exercises. To help develop this model, MRI images are acquired from individuals participating in the bed rest studies. The MRIs can be used to create a subject specific model of each individual for testing. The MRIs are processed in the Magnetic Resonance Imaging Data Transfer System program to develop a three-dimensional finite element model of the femur for evaluation. Modifications were made to the MRIDTS that simplified the model creation process. These modifications made it possible to construct two separate models of different portions of a bone simultaneously and then later connect them manually. This helped alleviate the warping problem associated with the drastic changes in geometry found in some body parts, such as the joints. The code was also modified to incorporate material properties of various bone components into the model. Interior meshing was also incorporated into the program to allow for both the cortical shell and the entire bone to be modeled. A prototype model of the right femur of an adult female is being constructed and tested to determine the feasibility of finite element analysis as a tool for evaluating exercise effectiveness. The model is being run through the ANSYS finite element program on the Alabama Super Computer Network. After the model is validated, models of bedrest subjects can be generated to investigate exercise countermeasures.

Murphy, Benjamin M.↗

A comprehensive model to determine the effects of temperature and species fluctuations on reaction rates in turbulent reacting flows

Reaction rates in turbulent, reacting flows are reviewed. Assumed probability density functions (pdf) modeling of reaction rates is being investigated in relation to a three variable pdf employing a 'most likely pdf' model. Chemical kinetic mechanisms treating hydrogen air combustion is studied. Perfectly stirred reactor modeling of flame stabilizing recirculation regions was used to investigate the stable flame regions for silane, hydrogen, methane, and propane, and for certain mixtures thereof. It is concluded that in general, silane can be counted upon to stabilize flames only when the overall fuel air ratio is close to or greater than unity. For lean flames, silane may tend to destabilize the flame. Other factors favoring stable flames are high initial reactant temperatures and system pressure.

Goldstein, D.↗

Using Particle Image Velocimetry to Determine Turbulence Model Parameters

The primary parameter of a standard k-ϵ model, C μ , was calculated from stereoscopic particle image velocimetry (PIV) data for a supersonic jet exhausting into a transonic crossflow. This required the determination of turbulent kinetic energy, turbulent eddy viscosity, and turbulent energy dissipation rate. Image interrogation was optimized, with different procedures used for mean strain rates and Reynolds stresses, to produce useful turbulent eddy viscosity fields. The eddy viscosity was calculated by a least-squares fit to all components of the three-dimensional strain-rate tensor that were available from the PIV data. This eliminated artifacts and noise observed when using a single strain component. Local dissipation rates were determined via Kolmogorov’s similarity hypotheses and the second-order structure function. The eddy viscosity and dissipation rates were then combined to determine C μ . Considerable spatial variation was observed in C μ , with the highest values found in regions where turbulent kinetic energy was relatively low but where turbulent mixing was important, e.g., along the high-strain jet edges and in the wake region. This suggests that use of a constant C μ in modeling may lead to poor Reynolds stress predictions at mixing interfaces. A data-driven modeling approach that can predict this spatial variation of C μ based on known state variables may lead to improved simulation results without the need for calibration.

42 ENGINEERING↗

A model to determine open or closed cellular convection

The horizontal asymmetry of cellular convection is demonstrated to be detectable by consideration of the vertical asymmetry of the driving force. A two-dimensional numerical model for convection in an internally heated and cooled fluid is presented, based on equations for the conservation of temperature and vorticity. Attention is focused on the steady-state finite-amplitude solutions for fixed Rayleigh number, Prandtl number, and aspect ratio. Temperature in the model corresponds to potential temperature under dry conditions and the equivalent potential temperature under saturated conditions. Slowly varying convection driven by asymmetric boundary fluxes is expressed in terms of steady-state convection driven by an asymmetric internal heat source. It is shown that heating near the ground produces open convection patterns where most of the fluid is descending, while cooling near the top of the flow leads to closed cellular patterns with a preponderance of ascending fluid. Extension of the model to three dimensions is indicated.

Helfand, H. M.↗

An Improved Model For Determining Salinity Recharge Time for Deep Borehole Disposal - 20405

In a previous publication (K. P. Travis, D. Burley and F. G. F. Gibb, WM2017 Conference, Phoenix Arizona paper no. 17480) we introduced a numerical model to estimate the time it would take for a pressure perturbation to subside following its creation from the construction of a deep geological borehole in water-saturated rock. The model was based on the notion of a point source of momentum - a solution of the time-dependent pressure-diffusion equation. Using the model, we estimated that it would take on the order of 10 k years for physical equilibrium to become re-established following the sinking of a 5 km borehole in granite. The significance of such a model is that it places an upper bound on the lifetime required of a borehole sealing system - a necessary input to a borehole post-closure safety case assessment. No engineered sealing system has ever been devised which is capable of retaining its sealing properties for 10 half-lives of long-lived radioisotopes in spent fuel or high-level waste. One of the key advantages of disposing of nuclear waste via Deep Borehole Disposal (DBD) is the natural sealing provided by density stratified groundwater. Upon sinking of a borehole, and subsequent filling with fresh water, brine or drilling mud, this natural barrier may be temporarily damaged. Over time, fresh brine from the far-field will flow towards or away from the hole (driven by a pressure gradient) and will eventually re-establish the original salinity gradient. It follows that the engineered seals need last only as long as the time required for this salinity gradient to reset itself. One of the limitations of our previous model was the use of a static boundary condition on the borehole wall. The model used a boundary pressure which varied quadratically with depth (arising from differences between the pressure of a column of fresh water in the borehole and that of a column of brine in the host rock). However as brine replaces fresh water in the borehole (driven by a pressure gradient), the boundary function must change with time. We now introduce an improved model which takes this time dependent boundary condition into account. The model also takes into account the time taken for a mixture of brine and fresh water to re-establish chemical equilibrium through the process of diffusion. The paper contains the mathematical details of our new iterative model as well as results showing the time taken to reach steady state, and concentration profiles for the components of the brine-filled borehole as a function of time together with a discussion on the implications of the results for developing a post-closure safety assessment for DBD. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Preliminary Deposition Modeling: For Determining the Deposition of Corrosive Contaminants on SNF Canisters

This report fulfills the M3 milestone M3SF-20PN010207045. During fiscal year (FY) 2019, Pacific Northwest National Laboratory (PNNL) worked to develop initial deposition and particle tracking models which could be used to better understand how corrosive contaminants deposit on the surfaces of spent nuclear fuel (SNF) canisters during dry storage. These models may be useful in understanding the likelihood of chloride induced stress corrosion cracking (CISCC) on SNF canisters

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Analysis of geologic terrain models for determination of optimum SAR sensor configuration and optimum information extraction for exploration of global non-renewable resources. Pilot study: Arkansas Remote Sensing Laboratory, part 1, part 2, and part 3

Computer-generated radar simulations and mathematical geologic terrain models were used to establish the optimum radar sensor operating parameters for geologic research. An initial set of mathematical geologic terrain models was created for three basic landforms and families of simulated radar images were prepared from these models for numerous interacting sensor, platform, and terrain variables. The tradeoffs between the various sensor parameters and the quantity and quality of the extractable geologic data were investigated as well as the development of automated techniques of digital SAR image analysis. Initial work on a texture analysis of SEASAT SAR imagery is reported. Computer-generated radar simulations are shown for combinations of two geologic models and three SAR angles of incidence.

Kaupp, V. H.↗

Model-independent determination of the dipole response of 66 Zn using quasimonoenergetic and linearly polarized photon beams

Background: Photon strength functions are an important ingredient in calculations relevant for the nucleosynthesis of heavy elements. The relation to the photoabsorption cross section allows to experimentally constrain photon strength functions by investigating the photoresponse of atomic nuclei. Purpose: Here, we determine the photoresponse of 66 Zn in the energy region of 5.6 MeV to 9.9 MeV and analyze the contribution of the 'elastic' decay channel back to the ground state. In addition, for the elastic channel electric and magnetic dipole transitions were separated. Methods: Nuclear resonance fluorescence experiments were performed using a linearly polarized quasi-monoenergetic photon beam at the High Intensity γ-ray Source. Photon beam energies from 5.6 to 9.9 MeV with an energy spread of about 3% were selected in steps of 200–300 keV. Two high purity germanium detectors were used for the subsequent γ-ray spectroscopy. Results: Full photoabsorption cross sections are extracted from the data making use of the monoenergetic character of the photon beam. For the ground-state decay channel, the average contribution of electric and magnetic dipole strengths is disentangled. The average branching ratio back to the ground state is determined as well. Conclusions: The new results indicate lower cross sections when compared to the values extracted from a former experiment using bremsstrahlung on 66 Zn. In the latter, the average branching ratio to the ground state is estimated from statistical-model calculations in order to analyze the data. Corresponding estimates from statistical-model calculations underestimate this branching ratio compared to the values extracted from the present analysis, which would partly explain the high cross sections determined from the bremsstrahlung data.

59 ≤ A ≤ 89↗

A comprehensive model to determine the effects of temperature and species fluctuations on reaction rates in turbulent reacting flows

The use of probability theory to determine the effects of turbulent fluctuations on reaction rates in turbulent combustion systems is briefly reviewed. Results are presented for the effect of species fluctuations in particular. It is found that turbulent fluctuations of species act to reduce the reaction rates, in contrast with the temperature fluctuations previously determined to increase Arrhenius reaction rate constants. For the temperature fluctuations, a criterion is set forth for determining if, in a given region of a turbulent flow field, the temperature can be expected to exhibit ramp like fluctuations. Using the above results, along with results previously obtained, a model is described for testing the effects of turbulent fluctuations of temperature and species on reaction rates in computer programs dealing with turbulent reacting flows. An alternative model which employs three variable probability density functions (temperature and two species) and is currently being formulated is discussed as well.

Chinitz, W.↗