Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Monte Carlo Modelling”

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 991 records · Page 55
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Molecular dynamics simulation of effects of solutes on dislocation propagation in Ni-based superalloys

Ni-based superalloys are used in the hot sections of jet turbine engines because of their high strength, stability and resistance to oxidation. The properties of these alloys can be further optimized by adding alloying elements. Therefore, a fundamental knowledge on the effect of different elements on properties of Ni-based superalloys is required. Molecular dynamics simulation could shed light here but its application is hindered by the absence of reliable and computationally cheap semi-empirical potential of the interatomic interaction for 4 and more element alloys. We will present a new Ni-Al-Cr-Nb Finnis-Sinclair (FS) potential specially designed to simulate the dislocation propagation from to  phase. In order to construct this potential, we designed a special algorithm to incorporate the data on element partitioning in the potential development procedure. For example, it is known from experiment, that Cr is mostly present in the gamma phase. Figure 1 shows a snapshot obtained after equilibration of the model of the Ni68Al17Cr15 alloy at T=1000 K using the hybrid Monte-Carlo (MC)/molecular dynamics (MD) simulation with the developed semi-empirical potential. One can clearly see that the Cr partitioning is in agreement with the experimental data. We will discuss the developed algorithm to incorporate the solute partition data in details. Using the developed semi-empirical potential, we first investigated the effect of anti-site defects in the  phase on the single dislocation propagation. It was found that the dislocation velocity increases with the increasing of the anti-site defect concentration. This effect was attributed to smaller number of Al-Al pairs forming during the dislocation migration in the presence of the anti-site defects. Next, we investigated the effect of Nb on the dislocation pair propagation in the Kolbe mechanism. It was found that the addition of Nb leads to considerable decrease in the dislocation propagation rate. This is in agreement with the experimental data on the effect of Nb on the creep resistance of the Ni-based superalloys. We will discuss the origin of this effect.

Mikhail I. Mendelev↗
NASA NTRS2011

A New Aerodynamic Data Dispersion Method for Launch Vehicle Design

A novel method for implementing aerodynamic data dispersion analysis is herein introduced. A general mathematical approach combined with physical modeling tailored to the aerodynamic quantity of interest enables the generation of more realistically relevant dispersed data and, in turn, more reasonable flight simulation results. The method simultaneously allows for the aerodynamic quantities and their derivatives to be dispersed given a set of non-arbitrary constraints, which stresses the controls model in more ways than with the traditional bias up or down of the nominal data within the uncertainty bounds. The adoption and implementation of this new method within the NASA Ares I Crew Launch Vehicle Project has resulted in significant increases in predicted roll control authority, and lowered the induced risks for flight test operations. One direct impact on launch vehicles is a reduced size for auxiliary control systems, and the possibility of an increased payload. This technique has the potential of being applied to problems in multiple areas where nominal data together with uncertainties are used to produce simulations using Monte Carlo type random sampling methods. It is recommended that a tailored physics-based dispersion model be delivered with any aerodynamic product that includes nominal data and uncertainties, in order to make flight simulations more realistic and allow for leaner spacecraft designs.

Pinier, Jeremy T.↗
NASA NTRS2016

NASA's X-Plane Database and Parametric Cost Model v 2.0

The NASA Armstrong Cost Engineering Team with technical assistance from NASA HQ (SID)has gone through the full process in developing new CERs from Version #1 to Version #2 CERs. We took a step backward and reexamined all of the data collected, such as dependent and independent variables, cost, dry weight, length, wingspan, manned versus unmanned, altitude, Mach number, thrust, and skin. We used a well- known statistical analysis tool called CO$TAT instead of using "R" multiple linear or the "Regression" tool found in Microsoft Excel(TradeMark). We setup an "array of data" by adding 21" dummy variables;" we analyzed the standard error (SE) and then determined the "best fit." We have parametrically priced-out several future X-planes and compared our results to those of other resources. More work needs to be done in getting "accurate and traceable cost data" from historical X-plane records!

life cycle cost↗
NASA NTRSDate not supplied

Recent Updates to the Object Reentry Survival Analysis Tool (ORSAT) Version 7.1

The Object Reentry Survival Analysis Tool (ORSAT) code is maintained and used by the NASA Orbital Debris Program Office (ODPO) and has been under continuous development and improvement since the mid-1990s. ORSAT is an object-oriented reentry simulation tool; it models a satellite as a collection of discrete components that follow independent trajectories upon the breakup of the parent object. Version 7.1 of the tool incorporates five years of new thermal and aerodynamic model development, multi-processor parametric study capability, codebase upgrades, and numerous bug-fixes. The thermal demise model was completely rewritten using a forward-time/central-space numerical stencil and incorporating a new pyrolysis model for fiber-reinforced plastic (FRP) materials. New aerodynamic and aeroheating models for hollow cylinders and hollow square prisms were developed using a combination of flow simulations in the direct simulation Monte Carlo (DSMC) Analysis Code (DAC) and Data Parallel Line Relaxation (DPLR) code and free-flight tests in the University of Texas at San Antonio’s Hypersonic Wind Tunnel. The latest version also incorporates a mechanical, strength-based demise model for FRP materials. Minor improvements include an update to the Fortran 2018 codebase; improved integration and speed with the Python-based, multi-core, parametric study tool, AutoORSAT; and fixes for many minor bugs. This new version of ORSAT will enable more accurate reentry risk assessments for modern satellites. This paper presents an overview of these changes and a summary of the verification and validation performed on the final code.

Benton R. Greene↗
NASA NTRSDate not supplied

Recent Updates to the Object Reentry Survival Analysis Tool (ORSAT) Version 7.1

The Object Reentry Survival Analysis Tool (ORSAT) code is maintained and used by the NASA Orbital Debris Program Office (ODPO) and has been under continuous development and improvement since the mid-1990s. ORSAT is an object-oriented reentry simulation tool; it models a satellite as a collection of discrete components that follow independent trajectories upon the breakup of the parent object. Version 7.1 of the tool incorporates five years of new thermal and aerodynamic model development, multi-processor parametric study capability, codebase upgrades, and numerous bug-fixes. The thermal demise model was completely rewritten using a forward-time/central-space numerical stencil and incorporating a new pyrolysis model for fiber-reinforced plastic (FRP) materials. New aerodynamic and aeroheating models for hollow cylinders and hollow square prisms were developed using a combination of flow simulations in the direct simulation Monte Carlo (DSMC) Analysis Code (DAC) and Data Parallel Line Relaxation (DPLR) code and free-flight tests in the University of Texas at San Antonio’s Hypersonic Wind Tunnel. The latest version also incorporates a mechanical, strength-based demise model for FRP materials. Minor improvements include an update to the Fortran 2018 codebase; improved integration and speed with the Python-based, multi-core, parametric study tool, AutoORSAT; and fixes for many minor bugs. This new version of ORSAT will enable more accurate reentry risk assessments for modern satellites. This paper presents an overview of these changes and a summary of the verification and validation performed on the final code.

Benton R. Greene↗
NASA NTRS1994

Space proton transport in one dimension

An approximate evaluation procedure is derived for a second-order theory of coupled nucleon transport in one dimension. An analytical solution with a simplified interaction model is used to determine quadrature parameters to minimize truncation error. Effects of the improved method on transport solutions with the BRYNTRN data base are evaluated. Comparisons with Monte Carlo benchmarks are given. Using different shield materials, the computational procedure is used to study the physics of space protons. A transition effect occurs in tissue near the shield interface and is most important in shields of high atomic number.

NASA Discipline Number 45-10↗
NASA NTRSDate not supplied

Assessment of Control Algorithms for Mars Entry Vehicles With Flap-Based Trajectory Control Under Uncertainty

Flap-based steering systems on blunt-body entry vehicles may improve flight performance relative to heritage bank-angle steering systems. Successful implementation of articulating aerodynamic flaps on a hypersonic entry vehicle requires an active control system to map angle of attack and sidelip angle commands to flap deflection commands. Here, a successive-linearization model prective control algorithm, as a linear quadratic regulator, are formulated, designed, and assessed to address this multiple input multiple output control problem. These two control algorithms are assessed under uncertainty in Monte Carlo simulations for various control profiles and vehicle configurations. Results indicate that while both control algorithms provide successful command tracking under uncertainty, the model predictive controller provides a lower mean and integral error, as well as greater robustness, relative to the linear-quadratic regulator.

Daniel L Engel↗
NASA NTRS2011

Unsupervised Detection of Planetary Craters by a Marked Point Process

With the launch of several planetary missions in the last decade, a large amount of planetary images is being acquired. Preferably, automatic and robust processing techniques need to be used for data analysis because of the huge amount of the acquired data. Here, the aim is to achieve a robust and general methodology for crater detection. A novel technique based on a marked point process is proposed. First, the contours in the image are extracted. The object boundaries are modeled as a configuration of an unknown number of random ellipses, i.e., the contour image is considered as a realization of a marked point process. Then, an energy function is defined, containing both an a priori energy and a likelihood term. The global minimum of this function is estimated by using reversible jump Monte-Carlo Markov chain dynamics and a simulated annealing scheme. The main idea behind marked point processes is to model objects within a stochastic framework: Marked point processes represent a very promising current approach in the stochastic image modeling and provide a powerful and methodologically rigorous framework to efficiently map and detect objects and structures in an image with an excellent robustness to noise. The proposed method for crater detection has several feasible applications. One such application area is image registration by matching the extracted features.

Troglio, G.↗
NASA NTRSDate not supplied

In-Depth Modeling and Simulation Analysis of Artemis Missions Using the Impact Probabilistic Risk Assessment Tool

BACKGROUND The Artemis campaign is a Moon exploration program with a series of six planned missions, five of which will be crewed. These five crewed missions will contain a single mission segment (space flight), or multiple mission segments involving space flight (Orion), lunar landing (LTV) and/or space habitat (Gateway). Each crewed segment faces the risk of unique medical conditions, necessitating medical sets/kits tailored to those specificities. To support and enable a data-driven and evidence-based decision-making process through out a mission’s life cycle, a software tool called IMPACT was developed. Using probabilistic risk assessment (PRA) methodologies, IMPACT (Informing Mission Planning via Analysis of Complex Tradespaces) is a novel tool built for analyzing the possibility of encountering complex medical risks during space flight, and for identifying the medical resources and capabilities needed to treat those potential at-risk medical conditions. IMPACT achieves this by performing hundreds of thousands of Monte Carlo simulations of missions to build aggregate pictures of medical risk. During an extended simulation modeling phase, IMPACT generated analytical results for medical risks, and the medical resources and capabilities to address those risks, for every segment of every crewed Artemis mission. This presentation will highlight the reliability, consistency and validity of IMPACT’s computational modeling techniques and will showcase the library of analytical outcomes generated for the Artemis missions. OVERVIEW During the early stages of IMPACT’s design, architecture and technical requirements collection, “scenarios” (use cases) - achievement goals required for acceptance testing, were identified by stakeholders. IMPACT successfully completed the scenario testing requirements and undertook an extensive operational run phase utilizing a wide range of input combinations with a goal of delivering a cohesive, trustworthy, reliable, vast, and diverse body of evidence. The intent of these modeling runs was to validate consistency in output, ensure solidity of executable operations and to streamline processes by identifying areas requiring efficiency improvements. Using the many missions of Artemis, IMPACT ran variations of operational runs to assess the output for acceptable, as well as unusual characteristics. This rigorous long-term “shakedown” analysis was implemented to help build a collective body of evidence to aid in securing a high level of confidence, reliability, and validity in the output, whether from the applicational components of IMPACT, or the entirety of the operational process. ANTICIPATED ANALYSIS AND CONCLUSION This presentation will discuss the various categories of input criteria; the comparisons in the application of these input criteria to various Artemis missions; the preparation and collection of the body of evidence, and reliability of the computational modeling techniques. This paper serves as an initial analytical overview of IMPACT’s probabilistic risk assessment (PRA) medical risk outputs covering Artemis missions and is not intended to be deemed the official medical response for the Artemis campaign.

Crew Composition↗
NASA NTRSDate not supplied

Evaluating Delta-V Dispersions Using Linear Covariance Techniques with Applications to Rendezvous and Docking

One of the attractive advantages of using linear covariance analysis (LinCov) is that it can accurately produce both navigation errors and trajectory dispersions in a fraction of the time comparable to the standard Monte Carlo analysis approach that is traditionally adopted. It also has the capability of generating delta-v dispersions which is perhaps the most influential performance metric since many aspects of a mission design revolve around the anticipated propellant usage. Accurately capturing this key performance criteria using linear covariance techniques has several subtle caveats raising reasonable doubts and a cautious sense of skepticism. This paper outlines multiple approaches of modeling delta-v usage in LinCov, highlights their limitations and advantages, and ultimately compares their corresponding results to the actual Monte Carlo performance. The theory is also extended to reliably account for an arbitrary thruster configuration and the resultant propellant consumptions due to both translational burns and attitude maneuvers. This provides near instant insight into the impacts of thruster layout design and different vehicle mass properties. These derived techniques are applied and verified using an NRHO rendezvous and docking scenario consistent with upcoming NASA Artemis missions for Orion.

Linear Covariance Analysis↗
NASA NTRS2002

Using Monte-Carlo Simulations to Study the Disk Structure in Cygnus X-1

As the first dynamically determined black hole X-ray binary system, Cygnus X-1 has been studied extensively. However, its broad-band spectra in hard state with BeppoSAX is still not well understood. Besides the soft excess described by the multi-color disk model (MCD), the power- law component and a broad excess feature above 10 keV (disk reflection component), there is also an additional soft component around 1 keV, whose origin is not known currently.We propose that the additional soft component is due to the thermal Comptonization process between the s oft disk photon and the warm plasma cloud just above the disk.i.e., a warm layer. We use Monte-Carlo technique t o simulate this Compton scattering process and build several table models based on our simulation results.

Yao, Y.↗
NASA NTRSDate not supplied

Validating Drag and Heating Coefficients for Hollow Reentry Objects in Continuum Flow Using a Mach 7 Ludwieg Tube

Drag and heating coefficient databases and models are crucial to destructive reentry simulation. The NASA Orbital Debris Program Office (ODPO) develops, maintains, and performs analysis with the Object Reentry Survival Analysis Tool (ORSAT), which comprises drag and heating models for free molecular, transitional, and continuum flow regimes. These models have, in the past, only included solid, convex, blunt shapes (such as boxes, spheres, and cylinders). Previous work led by ODPO includes the extension of these models to hollow cylinders and square boxes in free molecular and transitional flow using the Direct Simulation Monte Carlo (DSMC) method. Since 2019, the ODPO has continued its program of DSMC simulations and extended the project to include analyses with the NASA Data Parallel Line Relaxation (DPLR) program on hollow cylinders and boxes (with varying wall thickness-diameter ratio). In fall 2022, the ODPO began a collaboration with the University of Texas San Antonio (UTSA) to use the Mach 7 Ludwieg Tube facility to validate the model built using numerical simulations. This facility can replicate (at a scale of approximately 100:1) the conditions seen by reentering objects near typical demise altitudes. We present here the drag and heating coefficients derived from the continued DSMC simulations, the new DPLR simulations, and the 26-test series at UTSA.

Chris Ostrom↗
NASA NTRS2023

Validating Drag and Heating Coefficients for Hollow Reentry Objects in Continuum Flow Using a Mach 7 Ludwieg Tube

Drag and heating coefficient databases and models are crucial to destructive reentry simulation. The NASA Orbital Debris Program Office (ODPO) develops, maintains, and performs analysis with the Object Reentry Survival Analysis Tool (ORSAT), which comprises drag and heating models for free molecular, transitional, and continuum flow regimes. These models have, in the past, only included solid, convex, blunt shapes (such as boxes, spheres, and cylinders). Previous work led by ODPO includes the extension of these models to hollow cylinders and square boxes in free molecular and transitional flow using the Direct Simulation Monte Carlo (DSMC) method. Since 2019, the ODPO has continued its program of DSMC simulations and extended the project to include analyses with the NASA Data Parallel Line Relaxation (DPLR) program on hollow cylinders and boxes (with varying wall thickness-diameter ratio). In fall 2022, the ODPO began a collaboration with the University of Texas San Antonio (UTSA) to use the Mach 7 Ludwieg Tube facility to validate the model built using numerical simulations. This facility can replicate (at a scale of approximately 100:1) the conditions seen by reentering objects near typical demise altitudes. We present here the drag and heating coefficients derived from the continued DSMC simulations, the new DPLR simulations, and the 26-test series at UTSA.

Chris Ostrom↗
NASA NTRS1998

Direct Acceleration of Pickup Ions at The Solar Wind Termination Shock: The Production of Anomalous Cosmic Rays

We have modeled the injection and acceleration of pickup ions at the solar wind termination shock and investigated the parameters needed to produce the observed Anomalous Cosmic Ray (ACR) fluxes. A non-linear Monte Carlo technique was employed, which in effect solves the Boltzmann equation and is not restricted to near-isotropic particle distribution functions. This technique models the injection of thermal and pickup ions, the acceleration of these ions, and the determination of the shock structure under the influence of the accelerated ions. The essential effects of injection are treated in a mostly self-consistent manner, including effects from shock obliquity, cross- field diffusion, and pitch-angle scattering. Using recent determinations of pickup ion densities, we are able to match the absolute flux of hydrogen in the ACRs by assuming that pickup ion scattering mean free paths, at the termination shock, are much less than an AU and that modestly strong cross-field diffusion occurs. Simultaneously, we match the flux ratios He(+)/H(+) or O(+)/H(+) to within a factor approx. 5. If the conditions of strong scattering apply, no pre-termination-shock injection phase is required and the injection and acceleration of pickup ions at the termination shock is totally analogous to the injection and acceleration of ions at highly oblique interplanetary shocks recently observed by the Ulysses spacecraft. The fact that ACR fluxes can be modeled with standard shock assumptions suggests that the much-discussed "injection problem" for highly oblique shocks stems from incomplete (either mathematical or computer) modeling of these shocks rather than from any actual difficulty shocks may have in injecting and accelerating thermal or quasi-thermal particles.

Ellison, Donald C.↗
NASA NTRS2013

Dynamic Modeling Accuracy Dependence on Errors in Sensor Measurements, Mass Properties, and Aircraft Geometry

A nonlinear simulation of the NASA Generic Transport Model was used to investigate the effects of errors in sensor measurements, mass properties, and aircraft geometry on the accuracy of dynamic models identified from flight data. Measurements from a typical system identification maneuver were systematically and progressively deteriorated and then used to estimate stability and control derivatives within a Monte Carlo analysis. Based on the results, recommendations were provided for maximum allowable errors in sensor measurements, mass properties, and aircraft geometry to achieve desired levels of dynamic modeling accuracy. Results using other flight conditions, parameter estimation methods, and a full-scale F-16 nonlinear aircraft simulation were compared with these recommendations.

Grauer, Jared A.↗
NASA NTRS2006

Adapting Covariance Propagation to Account for the Presence of Modeled and Unmodeled Maneuvers

This paper explores techniques that can be used to adapt the standard linearized propagation of an orbital covariance matrix to the case where there is a maneuver and an associated execution uncertainty. A Monte Carlo technique is used to construct a final orbital covariance matrix for a 'prop-burn-prop' process that takes into account initial state uncertainty and execution uncertainties in the maneuver magnitude. This final orbital covariance matrix is regarded as 'truth' and comparisons are made with three methods using modified linearized covariance propagation. The first method accounts for the maneuver by modeling its nominal effect within the state transition matrix but excludes the execution uncertainty by omitting a process noise matrix from the computation. The second method does not model the maneuver but includes a process noise matrix to account for the uncertainty in its magnitude. The third method, which is essentially a hybrid of the first two, includes the nominal portion of the maneuver via the state transition matrix and uses a process noise matrix to account for the magnitude uncertainty. The first method is unable to produce the final orbit covariance except in the case of zero maneuver uncertainty. The second method yields good accuracy for the final covariance matrix but fails to model the final orbital state accurately. Agreement between the simulated covariance data produced by this method and the Monte Carlo truth data fell within 0.5-2.5 percent over a range of maneuver sizes that span two orders of magnitude (0.1-20 m/s). The third method, which yields a combination of good accuracy in the computation of the final covariance matrix and correct accounting for the presence of the maneuver in the nominal orbit, is the best method for applications involving the computation of times of closest approach and the corresponding probability of collision, PC. However, applications for the two other methods exist and are briefly discussed. Although the process model ("prop-burn-prop") that was studied is very simple - point-mass gravitational effects due to the Earth combined with an impulsive delta-V in the velocity direction for the maneuver - generalizations to more complex scenarios, including high fidelity force models, finite duration maneuvers, and maneuver pointing errors, are straightforward and are discussed in the conclusion.

Schiff, Conrad↗
NASA NTRS1988

Development of a satellite SAR image spectra and altimeter wave height data assimilation system for ERS-1

The applicability of ERS-1 wind and wave data for wave models was studied using the WAM third generation wave model and SEASAT altimeter, scatterometer and SAR data. A series of global wave hindcasts is made for the surface stress and surface wind fields by assimilation of scatterometer data for the full 96-day SEASAT and also for two wind field analyses for shorter periods by assimilation with the higher resolution ECMWF T63 model and by subjective analysis methods. It is found that wave models respond very sensitively to inconsistencies in wind field analyses and therefore provide a valuable data validation tool. Comparisons between SEASAT SAR image spectra and theoretical SAR spectra derived from the hindcast wave spectra by Monte Carlo simulations yield good overall agreement for 32 cases representing a wide variety of wave conditions. It is concluded that SAR wave imaging is sufficiently well understood to apply SAR image spectra with confidence for wave studies if supported by realistic wave models and theoretical computations of the strongly nonlinear mapping of the wave spectrum into the SAR image spectrum. A closed nonlinear integral expression for this spectral mapping relation is derived which avoids the inherent statistical errors of Monte Carlo computations and may prove to be more efficient numerically.

Hasselmann, Klaus↗