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 577 records · Page 32

Scaling and Benchmarking an Evolutionary Algorithm for Constructing Biophysical Neuronal Models

Single neuron models are fundamental for computational modeling of the brain's neuronal networks, and understanding how ion channel dynamics mediate neural function. A challenge in defining such models is determining biophysically realistic channel distributions. Here, we present an efficient, highly parallel evolutionary algorithm for developing such models, named NeuroGPU-EA. NeuroGPU-EA uses CPUs and GPUs concurrently to simulate and evaluate neuron membrane potentials with respect to multiple stimuli. We demonstrate a logarithmic cost for scaling the stimuli used in the fitting procedure. NeuroGPU-EA outperforms the typically used CPU based evolutionary algorithm by a factor of 10 on a series of scaling benchmarks. We report observed performance bottlenecks and propose mitigation strategies. Finally, we also discuss the potential of this method for efficient simulation and evaluation of electrophysiological waveforms.

59 BASIC BIOLOGICAL SCIENCES↗

The Snowmelt-Runoff Model (SRM) user's manual

A manual to provide a means by which a user may apply the snowmelt runoff model (SRM) unaided is presented. Model structure, conditions of application, and data requirements, including remote sensing, are described. Guidance is given for determining various model variables and parameters. Possible sources of error are discussed and conversion of snowmelt runoff model (SRM) from the simulation mode to the operational forecasting mode is explained. A computer program is presented for running SRM is easily adaptable to most systems used by water resources agencies.

Martinec, J.↗

O/S analysis of conceptual space vehicles

The application of recently developed computer models in determining operational capabilities and support requirements during the conceptual design of proposed space systems is discussed. The models used are the reliability and maintainability (R&M) model, the maintenance simulation model, and the operations and support (O&S) cost model. In the process of applying these models, the R&M and O&S cost models were updated. The more significant enhancements include (1) improved R&M equations for the tank subsystems, (2) the ability to allocate schedule maintenance by subsystem, (3) redefined spares calculations, (4) computing a weighted average of the working days and mission days per month, (5) the use of a position manning factor, and (6) the incorporation into the O&S model of new formulas for computing depot and organizational recurring and nonrecurring training costs and documentation costs, and depot support equipment costs. The case study used is based upon a winged, single-stage, vertical-takeoff vehicle (SSV) designed to deliver to the Space Station Freedom (SSF) a 25,000 lb payload including passengers without a crew.

Ebeling, Charles E.↗

Transient Finite Element Analyses Developed to Model Fan Containment Impact Events

Research is underway to establish an increased level of confidence in existing numerical techniques for predicting transient behavior when the fan of a jet engine is released and impacts the fan containment system. To evaluate the predictive accuracy that can currently be obtained, researchers at the NASA Lewis Research Center used the DYNA 3D computer code to simulate large-scale subcomponent impact tests that were conducted at the University of Dayton Research Institute (UDRI) Impact Physics Lab. In these tests, 20- by 40-in. flat metal panels, contoured to the shape of a typical fan case, were impacted by the root section of a fan blade. The panels were oriented at an angle to the path of the projectile that would simulate the conditions in an actual blade-out event. The metal panels were modeled in DYNA 3D using a kinematic hardening model with the strain rate dependence of the yield stress governed by the Cowper-Simons rule. Failure was governed by the effective plastic strain criterion. The model of the fan blade and case just after impact is shown. By varying the maximum effective plastic strain, we obtained good qualitative agreement between the model and the experiments. Both the velocity required to penetrate the case and the deflection during impact compared well. This indicates that the failure criterion and constitutive model may be appropriate, but for DYNA 3D to be useful as a predictive tool, methods to determine accurate model parameters must be established. Simple methods for measuring model parameters are currently being developed. In addition, alternative constitutive models and failure criteria are being investigated.

Pereira, J. Michael↗

Cluster electric current density measurements within a magnetic flux rope in the plasma sheet

On August 22, 2001 all 4 Cluster spacecraft nearly simultaneously penetrated a magnetic flux rope in the tail. The flux rope encounter took place in the central plasma sheet, Beta(sub i) approx. 1-2, near the leading edge of a bursty bulk flow. The "time-of-flight" of the flux rope across the 4 spacecraft yielded V(sub x) approx. 700 km/s and a diameter of approx.1 R(sub e). The speed at which the flux rope moved over the spacecraft is in close agreement with the Cluster plasma measurements. The magnetic field profiles measured at each spacecraft were first modeled separately using the Lepping-Burlaga force-free flux rope model. The results indicated that the center of the flux rope passed northward (above) s/c 3, but southward (below) of s/c 1, 2 and 4. The peak electric currents along the central axis of the flux rope predicted by these single-s/c models were approx.15-19 nA/sq m. The 4-spacecraft Cluster magnetic field measurements provide a second means to determine the electric current density without any assumption regarding flux rope structure. The current profile determined using the curlometer technique was qualitatively similar to those determined by modeling the individual spacecraft magnetic field observations and yielded a peak current density of 17 nA/m2 near the central axis of the rope. However, the curlometer results also showed that the flux rope was not force-free with the component of the current density perpendicular to the magnetic field exceeding the parallel component over the forward half of the rope, perhaps due to the pressure gradients generated by the collision of the BBF with the inner magnetosphere. Hence, while the single-spacecraft models are very successful in fitting flux rope magnetic field and current variations, they do not provide a stringent test of the force-free condition.

Slavin, J. A.↗

Preprocessing Inconsistent Linear System for a Meaningful Least Squares Solution

Mathematical models of many physical/statistical problems are systems of linear equations. Due to measurement and possible human errors/mistakes in modeling/data, as well as due to certain assumptions to reduce complexity, inconsistency (contradiction) is injected into the model, viz. the linear system. While any inconsistent system irrespective of the degree of inconsistency has always a least-squares solution, one needs to check whether an equation is too much inconsistent or, equivalently too much contradictory. Such an equation will affect/distort the least-squares solution to such an extent that renders it unacceptable/unfit to be used in a real-world application. We propose an algorithm which (i) prunes numerically redundant linear equations from the system as these do not add any new information to the model, (ii) detects contradictory linear equations along with their degree of contradiction (inconsistency index), (iii) removes those equations presumed to be too contradictory, and then (iv) obtain the minimum norm least-squares solution of the acceptably inconsistent reduced linear system. The algorithm presented in Matlab reduces the computational and storage complexities and also improves the accuracy of the solution. It also provides the necessary warning about the existence of too much contradiction in the model. In addition, we suggest a thorough relook into the mathematical modeling to determine the reason why unacceptable contradiction has occurred thus prompting us to make necessary corrections/modifications to the models - both mathematical and, if necessary, physical.

Sen, Syamal K.↗

Preprocessing in Matlab Inconsistent Linear System for a Meaningful Least Squares Solution

Mathematical models of many physical/statistical problems are systems of linear equations~ Due to measurement and possible human errors/mistakes in modeling/data, as well as due to certain assumptions to reduce complexity, inconsistency (contradiction) is injected into the model, viz. the linear system. While any inconsistent system irrespective of the degree of inconsistency has always a least-squares solution, one needs to check whether an equation is too much inconsistent or, equivalently too much contradictory. Such an equation will affect/distort the least-squares solution to such an extent that renders it unacceptable/unfit to be used in a real-world application. We propose an algorithm which (i) prunes numerically redundant linear equations from the system as these do not add any new information to the model, (ii) detects contradictory linear equations along with their degree of contradiction (inconsistency index), (iii) removes those equations presumed to be too contradictory, and then (iv) obtain the . minimum norm least-squares solution of the acceptably inconsistent reduced linear system. The algorithm presented in Matlab reduces the computational and storage complexities and also improves the accuracy of the solution. It also provides the necessary warning about the existence of too much contradiction in the model. In addition, we suggest a thorough relook into the mathematical modeling to determine the reason why unacceptable contradiction has occurred thus prompting us to make necessary corrections/modifications to the models - both mathematical and, if necessary, physical.

Sen, Symal K.↗

Foam-on-Tile Damage Model

An impact model was developed to predict how three specific foam types would damage the Space Shuttle Orbiter insulating tiles. The inputs needed for the model are the foam type, the foam mass, the foam impact velocity, the foam impact incident angle, the type being impacted, and whether the tile is new or aged (has flown at least one mission). The model will determine if the foam impact will cause damage to the tile. If it can cause damage, the model will output the damage cavity dimensions (length, depth, entry angle, exit angle, and sidewall angles). It makes the calculations as soon as the inputs are entered (less than 1 second). The model allows for the rapid calculation of numerous scenarios in a short time. The model was developed from engineering principles coupled with significant impact testing (over 800 foam impact tests). This model is applicable to masses ranging from 0.0002 up to 0.4 pound (0.09 up to 181 g). A prior tool performed a similar function, but was limited to the assessment of a small range of masses and did not have the large test database for verification. In addition, the prior model did not provide outputs of the cavity damage length, entry angle, exit angle, or sidewall angles.

Koharchik, Michael↗

Development and Calibration of a System-Integrated Rotorcraft Finite Element Model for Impact Scenarios

Two full-scale crash tests of an MD-500 helicopter were conducted in 2009 and 2010 at NASA Langley's Landing and Impact Research Facility in support of NASA s Subsonic Rotary Wing Crashworthiness Project. The first crash test was conducted to evaluate the performance of an externally mounted composite deployable energy absorber (DEA) under combined impact conditions. In the second crash test, the energy absorber was removed to establish baseline loads that are regarded as severe but survivable. The presence of this energy absorbing device reduced the peak impact acceleration levels by a factor of three. Accelerations and kinematic data collected from the crash tests were compared to a system-integrated finite element model of the test article developed in parallel with the test program. In preparation for the full-scale crash test, a series of sub-scale and MD-500 mass simulator tests were conducted to evaluate the impact performances of various components and subsystems, including new crush tubes and the DEA blocks. Parameters defined for the system-integrated finite element model were determined from these tests. Results from 19 accelerometers placed throughout the airframe were compared to finite element model responses. The model developed for the purposes of predicting acceleration responses from the first crash test was inadequate when evaluating more severe conditions seen in the second crash test. A newly developed model calibration approach that includes uncertainty estimation, parameter sensitivity, impact shape orthogonality, and numerical optimization was used to calibrate model results for the full-scale crash test without the DEA. This combination of heuristic and quantitative methods identified modeling deficiencies, evaluated parameter importance, and proposed required model changes. The multidimensional calibration techniques presented here are particularly effective in identifying model adequacy. Acceleration results for the calibrated model were compared to test results and the original model results. There was a noticeable improvement in the pilot and copilot region, a slight improvement in the occupant model response, and an over-stiffening effect in the passenger region. One lesson learned was that this approach should be adopted early on, in combination with the building-block approaches that are customarily used, for model development and pretest predictions. Complete crash simulations with validated finite element models can be used to satisfy crash certification requirements, potentially reducing overall development costs.

Annett, Martin S.↗

Nonlinear Unsteady Aerodynamic Modeling Using Empirical Orthogonal Functions

Empirical orthogonal function modeling is explained and applied to identify compact discrete-time nonlinear unsteady aerodynamic models from data generated by an unsteady three-dimensional compressible Navier-Stokes flow solver for an airfoil undergoing various pitching motions. Model structures, model parameter estimates, and model parameter uncertainty estimates for nondimensional lift, drag, and pitching moment coefficient models were determined autonomously and directly from the data. Prediction tests using data that were not used in the modeling process showed that the identified models exhibited excellent prediction capability, which is a strong indicator of an accurate model.

empirical↗

Fairfax County Urban Development: Identifying Urban Heat Mitigation Strategies for Climate Adaptation Planning in Fairfax County, Virginia

Extreme high temperatures lead to increased instances of cardiovascular disease, pulmonary disease, and even death, as well as increased energy consumption and infrastructure costs. People in urbanized areas experience higher temperatures than rural areas due to diminished vegetation and increased impervious surfaces which absorb and radiate heat. Fairfax County, Virginia has embarked on Resilient Fairfax, a program aimed at addressing climate adaptation and resilience. The DEVELOP team partnered with the Fairfax County Office of Environmental and Energy Coordination (OEEC) to assess the extent of the urban heat island effect on the county and its most vulnerable populations. The team used data from Landsat 8 Operational Land Imager (OLI) and Thermal Infrared Sensor (TIRS), as well as the ECOsystem Spaceborne Thermal Radiometer Experiment on Space Station (ECOSTRESS) for the years 2013 to 2021 and found that the hottest spots were in densely urbanized areas, with temperatures as much as 47°F above that of undeveloped reference areas. The team used the Integrated Valuation of Ecosystem Services and Tradeoffs (InVEST) urban cooling model and determined that areas with higher tree canopy cover had greater heat mitigation capacity. Estimates from the InVEST model showed that a 4.5% increase in canopy cover across the county could result in a temperature reduction of up to 2.4°F in some areas. The results will allow partners to assess heat distribution across Fairfax County and implement effective mitigation strategies, including locating prime locations for cooling centers and increasing canopy cover.

W. Pierce Holloway↗

Reliability models applicable to space telescope solar array assembly system

A complex system may consist of a number of subsystems with several components in series, parallel, or combination of both series and parallel. In order to predict how well the system will perform, it is necessary to know the reliabilities of the subsystems and the reliability of the whole system. The objective of the present study is to develop mathematical models of the reliability which are applicable to complex systems. The models are determined by assuming k failures out of n components in a subsystem. By taking k = 1 and k = n, these models reduce to parallel and series models; hence, the models can be specialized to parallel, series combination systems. The models are developed by assuming the failure rates of the components as functions of time and as such, can be applied to processes with or without aging effects. The reliability models are further specialized to Space Telescope Solar Arrray (STSA) System. The STSA consists of 20 identical solar panel assemblies (SPA's). The reliabilities of the SPA's are determined by the reliabilities of solar cell strings, interconnects, and diodes. The estimates of the reliability of the system for one to five years are calculated by using the reliability estimates of solar cells and interconnects given n ESA documents. Aging effects in relation to breaks in interconnects are discussed.

Patil, S. A.↗

Low order equivalent models of highly augmented aircraft determined from flight data using maximum likelihood estimation

This paper presents the results of a study of the feasibility of using low order equivalent mathematical models of a highly augmented aircraft, the F-8 digital fly-by-wire (DFBW), for flying qualities research. Increasingly complex models were formulated and evaluated using flight data and maximum likelihood estimation techniques. The aircraft actuator was modeled alone first. Next the equivalent derivatives were used to model the longitudinal unaugmented F-8 DFBW aircraft dynamics. The most complex model incorporated a pure time shift of the pilot input, a first order lag, and the basic longitudinal airframe model. This same model was then implemented for the F-8 DFBW aircraft in a highly augmented mode. Excellent matching of the dynamics resulted for this model, indicating that low order equivalent models which are good representations of the highly augmented F-8 DFBW aircraft can be formulated with these methods.

Shafer, M. F.↗

Multi-Scale Computational Modeling of Two-Phased Metal Using GMC Method

A multi-scale computational model for determining plastic behavior in two-phased CMSX-4 Ni-based superalloys is developed on a finite element analysis (FEA) framework employing crystal plasticity constitutive model that can capture the microstructural scale stress field. The generalized method of cells (GMC) micromechanics model is used for homogenizing the local field quantities. At first, GMC as stand-alone is validated by analyzing a repeating unit cell (RUC) as a two-phased sample with 72.9% volume fraction of gamma'-precipitate in the gamma-matrix phase and comparing the results with those predicted by finite element analysis (FEA) models incorporating the same crystal plasticity constitutive model. The global stress-strain behavior and the local field quantity distributions predicted by GMC demonstrated good agreement with FEA. High computational saving, at the expense of some accuracy in the components of local tensor field quantities, was obtained with GMC. Finally, the capability of the developed multi-scale model linking FEA and GMC to solve real life sized structures is demonstrated by analyzing an engine disc component and determining the microstructural scale details of the field quantities.

nickel alloys↗

Model uncertainty in accelerator application simulations

Monte-Carlo nuclear reaction and transport codes are widely used to devise accelerator-based nuclear physics experiments; at the same time, many experiments are performed to validate the Monte-Carlo codes, which can be used for the design of full-scale nuclear power applications or the design of new benchmark experiments. Dedicated model benchmark studies investigate a broad range of nuclear reactions and quantities. Examples of these include isotope formation or secondary particle fluxes that result from the interactions of GeV-range hadrons with monoisotopic targets, which can be used to assess the respective systematic uncertainty of models. Such benchmark studies, as well as many nuclear application experiments and simulations carried out by various groups over the last few decades, enable us to draw methodological lessons. In this work, model uncertainty determined based on available experimental data allow us to identify the effects of practitioner expertise as well as the design of codes (user access to micro-scale parameters) on the range of uncertainties. We found that in cases when simulations are performed by code developers or users that are very experienced in performing simulations, the model to experiment quantity ratios generally agree with the limits determined by dedicated benchmark studies. In other cases, the ratios generally tend to be either smaller (underestimation of model error) or larger (overestimation of model error). A plausible explanation of the aforementioned effects is suggested.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Determination of Extrapolation Distance with Measured Pressure Signatures from Two Low-Boom Models

A study to determine a limiting distance to span ratio for the extrapolation of near-field pressure signatures is described and discussed. This study was to be done in two wind-tunnel facilities with two wind-tunnel models. At this time, only the first half had been completed, so the scope of this report is limited to the design of the models, and to an analysis of the first set of measured pressure signatures. The results from this analysis showed that the pressure signatures measured at separation distances of 2 to 5 span lengths did not show the desired low-boom shapes. However, there were indications that the pressure signature shapes were becoming 'flat-topped'. This trend toward a 'flat-top' pressure signatures shape was seen to be a gradual one at the distance ratios employed in this first series of wind-tunnel tests.

Mack, Robert J.↗

Direct uranium enrichment assay in gaseous uranium hexafluoride with laser induced breakdown spectroscopy

A set of eleven U atomic emission lines, between 421.099 and 421.460 nm, were characterized for direct enrichment assay of gaseous UF 6 samples with laser induced breakdown spectroscopy (LIBS). Several iterative, multivariate nonlinear spectral-fitting algorithms were evaluated for their efficacies to extract the enrichment information from the measured LIBS spectra. Furthermore, wavelength-dependent weight factors, which take into consideration of the spectral-line position in the fitting model, the determined isotopic ratio, and the isotopic shift of the emission line, are essential for the spectral-fitting model to function adequately. The analytical accuracies and precision were typically within 0.5% in absolute [ 235 U/( 235 U + 238 U)] ratios.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Migration energy barriers and diffusion anisotropy of point defects on tungsten surfaces

We report tungsten is one of the most promising candidates for plasma-facing materials in future fusion devices, owing to its high performance under extreme irradiation conditions. However, irradiation-induced surface morphology varies significantly, depending on the irradiation type, fluence, and flux. Therefore, it is critical to examine the dynamics of point defects on tungsten (W) surfaces to understand how irradiation affects surface morphology. In this study, we employ the Self-Evolving Atomistic Kinetic Monte Carlo (SEAKMC) method to search for potential migration paths of point defects on W (1 0 0), (1 1 0) and (1 1 1) surfaces. The first-principles calculations are then used to accurately determine the migration energy barriers. The obtained paths and barriers are incorporated into a kinetic Monte Carlo (KMC) model to determine trajectories and diffusivities, which are described by diffusion tensors to demonstrate their anisotropic features. Multiple diffusion mechanisms have been identified on different surfaces, with various anisotropy factors. Particularly, point defects on the W (1 1 0) surface have the highest diffusivities, with anisotropy factors independent of temperatures. In comparison, the anisotropy of the W (1 0 0) surface decreases as temperature increases, while the W (1 1 1) surface is isotropic for point defects. This study provides insights into defect transport properties on different surfaces, which are essential for understanding the early stages of irradiation-induced microstructural evolutions and surface morphology of tungsten.

36 MATERIALS SCIENCE↗