Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Direct simulation Monte Carlo”

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 829 records · Page 46

Microreactor Optimization Using Simulation And Economics (mouse)

Microreactor Optimization Using Simulation and Economics (MOUSE) is a tool that integrates both nuclear microreactor design and reactor economics to provide comprehensive evaluations and optimizations. This tool enables stakeholders to explore the interplay between technical and economic variables, guiding them towards effective and competitive microreactor solutions. For the reactor core simulations, MOUSE leverages the OpenMC Monte Carlo Particle Transport Code to perform detailed core simulations for various microreactor designs. The included OpenMC models are 2D core designs of a Liquid Metal Thermal Microreactor (LMTR), a Gas-Cooled TRISO-Fueled Microreactor (GCMR), and a Heat Pipe Microreactor. Beyond core design, MOUSE includes simplified calculations for: - Calculating the masses of heat exchangers within the system. - Mechanical power of pumps. - Estimating the area occupied by various buildings within the nuclear plant. For the economic analysis, MOUSE provides detailed bottom-up cost estimates, encompassing a wide range of costs including preconstruction costs, direct costs, indirect costs, training costs, financial costs, operation & maintenance (O&M) costs, and fuel costs. These cost estimations are developed using data from the MARVEL project and additional literature sources, enabling the calculation of total capital costs and levelized cost of energy for both first-of-a-kind and nth-of-a-kind microreactors. MOUSE also enables analysis of the cost drivers and competitiveness in the electricity market. MOUSE allows users to modify a wide array of technical and economic parameters to evaluate different scenarios and their impacts. Examples of these parameters include: Fuels, coolants, or reflector materials Enrichment levels Control drum materials and geometry Fuel pin geometry and materials Moderator pin geometry and materials Reactor core and reflector dimensions Packing factor for the TRISO particles Nuclear reactor power and reactor burnup Number of sensors Shielding thickness Reactor vessel and guard vessel dimensions Operational staff requirements Number of emergency shutdowns Levelization period Interest rate Construction duration Since MOUSE is powered by the WATTS toolkit, it supports optimization studies, parametric analyses, and uncertainty calculations/propagation. The optimization techniques enable users to identify optimal design and economic configurations. The parametric analysis tools allow users to explore the sensitivity of various parameters, while uncertainty propagation helps quantify the impact of uncertainties on overall performance and cost. User Interface and Workflow: Currently, MOUSE is a command-line-based tool. Users can input various reactor design or economic parameters, modify the designs, run simulations, and visualize results through comprehensive data visualization and reporting capabilities. The typical workflow involves setting up the reactor model, defining economic parameters, running simulations, and analyzing the results to make informed decisions. By combining advanced design calculations with detailed economic modeling, MOUSE provides a robust framework for optimizing nuclear microreactor technologies, enhancing their competitiveness, and guiding stakeholders towards innovative and cost-effective solutions.

Hanna, Botros [Idaho National Laboratory (INL), Id↗

Neutron Tagging From Neutrino Interactions in DUNE-ND 2x2 Prototype

The Deep Underground Neutrino Experiment (DUNE) is a long-baseline neutrino oscillation experiment that aims to measure whether CP is violated in the leptonic sector (if violated, how much) and unambiguously determine the neutrino mass ordering. DUNE consists of near and far detectors that rely on liquid argon time projection chamber (LArTPC) technology to observe neutrino interactions. The near detector (ND) will be placed in Fermilab, near the neutrino source, while the far detector (FD) will be deployed in Sanford Lab, 1.5 km deep underground, which is 1300 km away from the source. LArTPCs provide excellent particle identification and calorimetry; however, detecting neutrons is challenging, as they do not leave direct ionization signals in LArTPCs. Neutrons can carry away up to 25% of the neutrino energy, introducing a significant uncertainty in DUNE measurements. The DUNE near detector (ND) features a novel modular LArTPC with pixelated charge readout, which enhances event recons truction. The modular design enables precise correlation between ionization signals and light signals in a high-rate environment, improving the identification of delayed energy depositions from neutrons in neutrino interactions. We introduce a neutron tagging technique using the 2x2 Demonstrator, a small-scale prototype of the DUNE ND LArTPC. The analysis utilizes Monte Carlo simulations and deep-learning techniques to identify neutron-induced energy depositions and reconstruct low-energy activity.

Kufatty, Georgette [Florida State U.]↗

Reliability Analysis and Reliability-Based Design Optimization of Circular Composite Cylinders Under Axial Compression

This report describes the preliminary results of an investigation on component reliability analysis and reliability-based design optimization of thin-walled circular composite cylinders with average diameter and average length of 15 inches. Structural reliability is based on axial buckling strength of the cylinder. Both Monte Carlo simulation and First Order Reliability Method are considered for reliability analysis with the latter incorporated into the reliability-based structural optimization problem. To improve the efficiency of reliability sensitivity analysis and design optimization solution, the buckling strength of the cylinder is estimated using a second-order response surface model. The sensitivity of the reliability index with respect to the mean and standard deviation of each random variable is calculated and compared. The reliability index is found to be extremely sensitive to the applied load and elastic modulus of the material in the fiber direction. The cylinder diameter was found to have the third highest impact on the reliability index. Also the uncertainty in the applied load, captured by examining different values for its coefficient of variation, is found to have a large influence on cylinder reliability. The optimization problem for minimum weight is solved subject to a design constraint on element reliability index. The methodology, solution procedure and optimization results are included in this report.

Rais-Rohani, Masoud↗

Studying the effects of overlapping objects in dark energy

Observing the clustering of galaxies allows us to calculate cosmological parameters necessary for understanding dark energy. However, as the density of observed objects increases, the probability of these objects blending likewise increases, causing multiple galaxies to be observed as one. This affects the inferred values of parameters such as the galaxy bias (b) and the matter energy density (ΩM). To see whether the bias from incorrectly inferring the galaxy count is significant, we compare the correlation functions in simulated data for “true” and “observed” data sets with one-to-one and multiple-to-one correspondences, respectively. For each data set, we create two correlation functions: one “measured” function directly relying on the galaxies’ positions using the TreeCorr python library, and one “model” derived mathematically from the galaxies’ power spectrum using the Cosmological Core Library (CCL). By minimizing the residual between these two functions, we compute the ideal values for b and ΩM across the various possible redshifts that position the galaxies in three dimensional space. This minimization is done with an Markov chain Monte Carlo (MCMC) estimate that finds one value of ΩM and ten values for b corresponding to the ten redshift bins ranging from z = 0.2 to z = 1.2. We find that neither b nor ΩM is particularly affected by inclusion of blended galaxies. Though there is room for improvement, the data suggests that the fluctuations we found are a result of noise or limitations on the modeling rather than blending explicitly.

79 ASTRONOMY AND ASTROPHYSICS↗

3D Heisenberg universality in the van der Waals antiferromagnet NiPS 3

Van der Waals (vdW) magnetic materials are comprised of layers of atomically thin sheets, making them ideal platforms for studying magnetism at the two-dimensional (2D) limit. These materials are at the center of a host of novel types of experiments, however, there are notably few pathways to directly probe their magnetic structure. We confirm the magnetic order within a single crystal of NiPS 3 and show it can be accessed with resonant elastic X-ray diffraction along the edge of the vdW planes in a carefully grown crystal by detecting structurally forbidden resonant magnetic X-ray scattering. We find the magnetic order parameter has a critical exponent of β ~ 0.36, indicating that the magnetism of these vdW crystals is more adequately characterized by the three-dimensional (3D) Heisenberg universality class. We verify these findings with first-principles density functional theory, Monte-Carlo simulations, and density matrix renormalization group calculations.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Performance Considerations for the SIMPL Single Photon, Polarimetric, Two-Color Laser Altimeter as Applied to Measurements of Forest Canopy Structure and Composition

The Slope Imaging Multi-polarization Photon-counting Lidar (SIMPL) is a multi-beam, micropulse airborne laser altimeter that acquires active and passive polarimetric optical remote sensing measurements at visible and near-infrared wavelengths. SIMPL was developed to demonstrate advanced measurement approaches of potential benefit for improved, more efficient spaceflight laser altimeter missions. SIMPL data have been acquired for wide diversity of forest types in the summers of 2010 and 2011 in order to assess the potential of its novel capabilities for characterization of vegetation structure and composition. On each of its four beams SIMPL provides highly-resolved measurements of forest canopy structure by detecting single-photons with 15 cm ranging precision using a narrow-beam system operating at a laser repetition rate of 11 kHz. Associated with that ranging data SIMPL provides eight amplitude parameters per beam unlike the single amplitude provided by typical laser altimeters. Those eight parameters are received energy that is parallel and perpendicular to that of the plane-polarized transmit pulse at 532 nm (green) and 1064 nm (near IR), for both the active laser backscatter retro-reflectance and the passive solar bi-directional reflectance. This poster presentation will cover the instrument architecture and highlight the performance of the SIMPL instrument with examples taken from measurements for several sites with distinct canopy structures and compositions. Specific performance areas such as probability of detection, after pulsing, and dead time, will be highlighted and addressed, along with examples of their impact on the measurements and how they limit the ability to accurately model and recover the canopy properties. To assess the sensitivity of SIMPL's measurements to canopy properties an instrument model has been implemented in the FLIGHT radiative transfer code, based on Monte Carlo simulation of photon transport. SIMPL data collected in 2010 over the Smithsonian Environmental Research Center, MD are currently being modelled and compared to other remote sensing and in situ data sets. Results on the adaptation of FLIGHT to model micropulse, single'photon ranging measurements are presented elsewhere at this conference. NASA's ICESat-2 spaceflight mission, scheduled for launch in 2016, will utilize a multi-beam, micropulse, single-photon ranging measurement approach (although non-polarimetric and only at 532 nm). Insights gained from the analysis and modelling of SIMPL data will help guide preparations for that mission, including development of calibration/validation plans and algorithms for the estimation of forest biophysical parameters.

Dabney, Philip W.↗

Comparison of spherical harmonics method and discrete ordinates method for radiative transfer in a turbulent jet flame

Here, in this study, we systematically compared the accuracy and computational cost of two popular solution methods for the radiative transfer equation (RTE): the spherical harmonics method (P N ) and the discrete ordinates method (DOM). We first investigated convergence characteristics of different orders of P N and DOM in a series of 1D homogeneous configurations with varying optical thicknesses. Both solvers perform better for optically thicker cases. The accuracy of P N methods increases with its order, , but the gain in accuracy reduces with the increase in , i.e., improvement of P 7 over P 5 is less than that of P 3 over P 1 . This decreasing trend becomes more prominent as the optical thickness decreases. On the other hand, DOM’s accuracy increases almost linearly with the increase in the number of ordinates (or polar angles in this study) in all cases. While comparing the directional profile of radiative intensity, both solvers perform better when the radiative intensity is more isotropic. These solvers were then connected with a full spectrum k-distribution (FSK) spectral model and used to perform radiation-coupled simulations of a turbulent jet flame in an axi-symmetric cylindrical domain. Results obtained from P 1 to P 7 approximations for P N , and 2 x 4, 4 x 4, 4 x 8, 8 x 8 finite angles for DOM are compared with that from an optically thin model, and a reference solution from line-by-line (LBL) photon Monte Carlo (PMC) method. The choice of radiation solver shows a noticeable impact on the temperature distribution of the flame. The P N solvers lead to slightly higher radiant fractions and the DOM solvers lead to slightly lower radiant fractions than the PMC benchmark solution. Finally, the computational costs of each of these solvers are also reported and an intermittent evaluation / time blending scheme to improve the computational efficiency of radiation solvers in radiation-coupled simulations are also demonstrated.

42 ENGINEERING↗

Effects of low-frequency voltage on nonlinear standing wave excitation, plasma uniformity, and ion dynamics in dual-frequency asymmetric capacitive discharges

Abstract It is known that in very-high-frequency (VHF) capacitively coupled plasmas, the higher harmonics generated by nonlinear sheath motion can enhance the standing wave effect (SWE), which can lead to center-peaked plasma density profiles. In this work, an improved nonlinear electromagnetic model incorporating a transmission line model, an electron momentum balance model, a bulk plasma model, a collisionless nonlinear numerical sheath model, and an ion Monte-Carlo collision (MCC) model is developed to study the effects of low-frequency (LF) voltage V L on the nonlinear standing wave excitation, plasma uniformity, and ion energy and angular distribution functions (IEDFs and IADFs) in dual-frequency (DF) asymmetric capacitive argon discharges at relatively low pressure of 3 Pa. The plasma diffusion in the radial direction and ion dynamics within the LF oscillating sheath are self-consistently considered. The LF voltage V L at 2 MHz varies from 0 to 700 V while the HF voltage V H at 60 MHz is fixed at 100 V. Simulation results indicate that without the addition of an LF source (i.e. V L = 0 V), there are a considerable number of high-order harmonics with short wavelengths, leading to significant SWE and central peak in the radial plasma density profile. Nevertheless, the high-order harmonic excitations tend to be weakened and merely occur around the phase of the full LF sheath collapse due to a shorter characteristic damping time of the surface waves as V L increases. This, combined with increased surface wavelengths of both the driving frequency and the higher harmonics at a higher V L , leads to suppressed standing waves and improved plasma uniformity. Meanwhile, the simulations show that both the low and the high energy peaks of IEDF move towards higher energies, and the energy peak separation width ΔEbecomes wider with the increase of V L . The IEDF at the radial center of the powered electrode exhibits a broader ΔEthan that at the edge. For the IADF, an increased V L results in more ions incident on the electrode with a smaller deflection angle. Because of a thinner sheath and a higher sheath voltage at the electrode center, the peak value of IADF at the electrode center is greater than that at the edge.

Physics↗

An Exploratory Exercise in Taguchi Analysis of Design Parameters: Application to a Shuttle-to-space Station Automated Approach Control System

The chief goals of the summer project have been twofold - first, for my host group and myself to learn as much of the working details of Taguchi analysis as possible in the time allotted, and, secondly, to apply the methodology to a design problem with the intention of establishing a preliminary set of near-optimal (in the sense of producing a desired response) design parameter values from among a large number of candidate factor combinations. The selected problem is concerned with determining design factor settings for an automated approach program which is to have the capability of guiding the Shuttle into the docking port of the Space Station under controlled conditions so as to meet and/or optimize certain target criteria. The candidate design parameters under study were glide path (i.e., approach) angle, path intercept and approach gains, and minimum impulse bit mode (a parameter which defines how Shuttle jets shall be fired). Several performance criteria were of concern: terminal relative velocity at the instant the two spacecraft are mated; docking offset; number of Shuttle jet firings in certain specified directions (of interest due to possible plume impingement on the Station's solar arrays), and total RCS (a measure of the energy expended in performing the approach/docking maneuver). In the material discussed here, we have focused on single performance criteria - total RCS. An analysis of the possibility of employing a multiobjective function composed of a weighted sum of the various individual criteria has been undertaken, but is, at this writing, incomplete. Results from the Taguchi statistical analysis indicate that only three of the original four posited factors are significant in affecting RCS response. A comparison of model simulation output (via Monte Carlo) with predictions based on estimated factor effects inferred through the Taguchi experiment array data suggested acceptable or close agreement between the two except at the predicted optimum point, where a difference outside a rule-of-thumb bound was observed. We have concluded that there is most likely an interaction effect not provided for in the original orthogonal array selected as the basis for our experimental design. However, we feel that the data indicates that this interaction is a mild one and that inclusion of its effect will not alter the location of the optimum.

Deal, Don E.↗

Inflation of Unreefed and Reefed Extraction Parachutes

Data from the Orion and several other test programs have been used to reconstruct inflation parameters for 28 ft Do extraction parachutes as well as the parent aircraft pitch response during extraction. The inflation force generated by extraction parachutes is recorded directly during tow tests but is usually inferred from the payload accelerometer during Low Velocity Airdrop Delivery (LVAD) flight test extractions. Inflation parameters are dependent on the type of parent aircraft, number of canopies, and standard vs. high altitude extraction conditions. For standard altitudes, single canopy inflations are modeled as infinite mass, but the non-symmetric inflations in a cluster are modeled as finite mass. High altitude extractions have necessitated reefing the extraction parachutes, which are best modeled as infinite mass for those conditions. Distributions of aircraft pitch profiles and inflation parameters have been generated for use in Monte Carlo simulations of payload extractions.

Ray, Eric S.↗

Examination of synthetic gas puff imaging diagnostic data from a gyrokinetic turbulence code

A synthetic gas puff imaging (GPI) diagnostic has been developed for the purpose of validating the three-dimensional gyrokinetic turbulence code XGC. The synthetic diagnostic is described and applied to XGC simulations of two Alcator C-Mod discharges. The turbulence characteristics deduced from the resulting simulated GPI frames, using analysis techniques similar to those applied to experimental data, are compared with locally derived characteristics extracted directly from the XGC output. The comparison of the two is shown to be potentially impacted by misalignment between the GPI view and the magnetic field, the dependence of the light emission on the electron density and temperature, and spatial and temporal variations in the neutral gas cloud induced by the turbulent plasma fluctuations. We conclude that quantitative and, in some cases even qualitative, validation of turbulence simulations need to account for these effects. While we cannot directly compare our results with experimental data due to the absence of high quality GPI data from the shots and times simulated by XGC, we do relate the overall characteristics obtained from the synthetic GPI analysis to published Alcator C-Mod GPI data.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Proper Motions, Orbits, and Tidal Influences of Milky Way Dwarf Spheroidal Galaxies

We combine Gaia early data release 3 astrometry with accurate photometry and utilize a probabilistic mixture model to measure the systemic proper motion of 52 dwarf spheroidal (dSph) satellite galaxies of the Milky Way (MW). For the 46 dSphs with literature line-of-sight velocities we compute orbits in both a MW and a combined MW + Large Magellanic Cloud (LMC) potential and identify Car II, Car III, Hor I, Hyi I, Phx II, and Ret II as likely LMC satellites. 40% of our dSph sample has a >25% change in pericenter and/or apocenter with the MW + LMC potential. For these orbits, we use a Monte Carlo sample for the observational uncertainties for each dSph and the uncertainties in the MW and LMC potentials. We predict that Ant II, Boo III, Cra II, Gru II, and Tuc III should be tidally disrupting by comparing each dSph's average density relative to the MW density at its pericenter. dSphs with large ellipticity (CVn I, Her, Tuc V, UMa I, UMa II, UMi, Wil 1) show a preference for their orbital direction to align with their major axis even for dSphs with large pericenters. We compare the dSph radial orbital phase to subhalos in MW-like N-body simulations and infer that there is not an excess of satellites near their pericenter. With projections of future Gaia data releases, we find that dSph's orbital precision will be limited by uncertainties in the distance and/or MW potential rather than in proper motion precision. Finally, we provide our membership catalogs to enable community follow-up.

79 ASTRONOMY AND ASTROPHYSICS↗

LandScan mosaic enables high-resolution gridded population estimates with explicit uncertainty

Gridded population datasets represent high-resolution distributions of human occupancy, enabling informed decision-making across a broad range of fields. These data products are valuable for assessing environmental risk, urban development, disaster preparedness and resource allocation—areas where accurate population estimates directly enhance policy effectiveness and optimize resource distribution. Despite the importance of gridded population datasets, traditional population modeling approaches often overlook inherent uncertainties in the estimation process. This limitation can create a false sense of certainty in population estimates, potentially leading to flawed decisions by those who rely on the data. To address this methodological gap, we introduce a probabilistic machine learning modeling framework, LandScan Mosaic, that explicitly incorporates uncertainty into the population modeling process. Our approach systematically quantifies uncertainty in three key modeling parameters of the LandScan HD gridded population dataset: building use types, floor counts, and occupancy rates. By employing Monte Carlo simulations, we propagate these uncertainties through the modeling process, yielding probability distributions of population counts in place of deterministic point estimates. We demonstrate the practical application of this framework in Iloilo City, Philippines, using structured decision-making techniques and our probabilistic estimates to identify and prioritize areas most affected by projected flooding, supporting targeted interventions that address both economic and social risks. In doing so, we propose a population-specific approach for incorporating confidence into structured decision making processes. Through a comparative analysis with conventional deterministic approaches and point estimate approaches, including LandScan HD and WorldPop, we evaluate how the incorporation of machine learning and uncertainty influences decision rankings. This research advances population distribution modeling by offering a robust, quantitative approach that explicitly accounts for uncertainty in the underlying data, along with guidance for how users can apply uncertainty in their decision-making.

Environmental sciences↗

Analysis and modeling of tungsten emission and net erosion in the DIII-D divertor using updated atomic data

Tungsten (W) is one of the leading candidate materials for plasma-facing components. However, its main drawback is its high radiative efficiency; if W penetrates the plasma, it can lead to core degradation or even collapse. Since eroded tungsten tends to ionize in the sheath and redeposit promptly, the net erosion flux that escapes prompt redeposition can differ significantly from the gross erosion. This work presents a modeling framework to estimate net erosion and photon emission from W coatings exposed to the lower divertor of DIII-D using the DiMES material exposure probe. The approach couples RustBCA for sputtering yields with a Monte Carlo transport code (LPTMC) that models redeposition and W emission. Computation is carried out with new atomic data, based on R-matrix and Mons calculations, leading to lower ionization probabilities and a twofold increase in net erosion estimates compared to calculations done with OPEN-ADAS atomic data. The model results are benchmarked against experimental measurements, showing quantitative agreement for erosion, although the trends in W emission are reproduced only qualitatively. The model is also used to assess whether W II emission can serve as a direct measurement of the net erosion of W in the lower divertor of DIII-D. Simulations show that this is not valid if the electron pressure is above ~120 Pa or if the toroidal length of the eroded material is smaller than the parallel-to-B distance traveled by impurity ions before steady-state conditions are reached. Finally, simulations suggest that when W is sputtered by carbon ions with high impact energies (≳300 eV) in DIII-D, W net erosion scales with W gross erosion and can be numerically approximated using W I flux alone as input.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Analysis and Compensation of Tx and Rx IQ Imbalances in Uplink IM-OFDMA Systems

Index modulation-based orthogonal division multiple access (IM-OFDMA) has recently been proposed as a potential technique for future wireless communication systems due to its superior spectral efficiency and error performance over conventional multiple access schemes. However, its performance is still under investigation by researchers in light of a variety of scenarios and assumptions. Following this direction, in this Article, the individual and joint effects of the transmitter (Tx) and receiver (Rx) in-phase and quadrature imbalances (IQI) on the bit-error-rate (BER) performance of the uplink IM-OFDMA scheme are investigated by considering whether each user's RF front-end is identical or not. Moreover, to reduce the detrimental effect of IQI, a preamble-based estimation and compensation method is proposed for IM-OFDMA. Closed-form expressions for the average BER of IM-OFDMA are obtained by considering the physical effect of the Tx and Rx IQI. Also, using Monte Carlo simulations, the derived expressions are verified under different system configurations. Furthermore, analytical and simulation results reveal that Tx and Rx IQI cause an error floor in the BER performance of IM-OFDMA schemes. However, the proposed estimation and compensation method not only reduces the impact of IQI but also leads to a better BER performance compared to the case of no IQI, which mainly refers to the frequency diversity caused by the IQI.

42 ENGINEERING↗

Short-range magnetic correlations in quasicrystalline i -Tb-Cd

Here, we report on elastic and inelastic neutron scattering from single-grain isotopically enriched samples to elucidate the local magnetic correlations between Tb 3+ moments in quasicrystalline i -Tb-Cd. The inelastic neutron scattering measurements of the crystalline electric field excitations demonstrated that the Tb 3+ moments are directed primarily along the local fivefold axes of the Tsai-type cluster as was found for the TbCd 6 approximant phase. Based on the inelastic measurements, we consider a simple Ising-type model for the moment configurations on a single Tb 3+ icosahedron and enumerate the lowest energy moment configurations. We then calculate the diffuse scattering from these configurations and compare with the experimental magnetic diffuse scattering measurements to identify the most likely single cluster moment configurations and find reasonable agreement between the broad features observed in our scattering simulations. We further use a heuristic glass model as well as large-scale Monte-Carlo simulations of a multicluster model to consider the role of higher-order (longer-range) intercluster correlations for magnetic frustration and the magnetic scattering.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Uncertainty Modeling for Mars Ascent Vehicle’s Aerodatabase Development

The design of the Mars Ascent Vehicle - Mars Sample Return (MAV-MSR) trajectories requires an accurate assessment of flight performance. Typically, these trajectories are developed by flight mechanics analysis to meet complex mission requirements and then flight performance assessed through Monte Carlo simulations. Consequently, it is crucial to develop an aerodynamic aerodatabase as an input model for flight mechanics analysis to provide static and dynamic force and moment coefficients under specific flight conditions. The force and moment coefficients in the MAV aerodatabase are determined using the FUN3D computational fluid dynamics solver. Firstly, an overview of the aerodatabase is presented to demonstrate its applicability to trajectory-defined simulations. This effort is followed by an initial attempt to quantify uncertainties in the force and moment coefficients necessary for updating the current aerodatabase. The uncertainty model identifies uncertainty adders and multipliers for coefficient-based forces and moments through a direct comparison between FUN3D and wind tunnel test data from NASA Marshall Space Center's 14x14 inch Trisonic Wind Tunnel. These uncertainties aim to encompass various changes in Mach number, angle of attack, and aerodynamic roll angle.

Uncertainty analysis↗

Uncertainty Modeling for Mars Ascent Vehicle’s Aerodynamic Database Development

The design of the Mars Ascent Vehicle - Mars Sample Return (MAV-MSR) trajectories requires an accurate assessment of flight performance. Typically, these trajectories are developed by flight mechanics analysis to meet complex mission requirements and then flight performance assessed through Monte Carlo simulations. Consequently, it is crucial to develop an aerodynamic aerodatabase as an input model for flight mechanics analysis to provide static and dynamic force and moment coefficients under specific flight conditions. The force and moment coefficients in the MAV aerodatabase are determined using the FUN3D computational fluid dynamics solver. Firstly, an overview of the aerodatabase is presented to demonstrate its applicability to trajectory-defined simulations. This effort is followed by an initial attempt to quantify uncertainties in the force and moment coefficients necessary for updating the current aerodatabase. The uncertainty model identifies uncertainty adders and multipliers for coefficient-based forces and moments through a direct comparison between FUN3D and wind tunnel test data from NASA Marshall Space Center's 14x14 inch Trisonic Wind Tunnel. These uncertainties aim to encompass various changes in Mach number, angle of attack, and aerodynamic roll angle.

Uncertainty analysis↗