Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation analytics”

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 721 records · Page 40

Progress in the Modeling of NiAl-Based Alloys Using the BFS Method

The BFS method has been applied to the study of NiAl-based materials to assess the effect of alloying additions on structure. Ternary, quaternary and even pent-alloys based on Ni-rich NiAl with additions of Ti, Cr and Cu were studied. Two approaches were used, Monte Carlo simulations to determine ground state structures and analytical calculations of high symmetry configurations which give physical insight into preferred bonding. Site occupancy energetics for ternary and the more complicated case of quaternary additions were determined, and solubility limits and precipitate formation with corresponding information concerning structure and lattice parameter were also 'observed' computationally. The method was also applied to determine the composition of alloy surfaces and interfaces. Overall, the results demonstrate that the BFS method for alloys is a powerful tool for alloy design and with its simplicity and obvious advantages can be used to complement any experimental alloy design program.

Bozzolo, Guillermo↗

Mechanics of Lithosphere Deformation

The objective of this study was to perform theoretical modeling of mechanical deformation of the Earth's lithosphere making use of analytical solutions as well as numerical computer simulations.

Zuber, Maria T.↗

Design, Fabrication, and Testing of Composite Energy-Absorbing Keel Beams for General Aviation Type Aircraft

A lightweight energy-absorbing keel-beam concept was developed and retrofitted in a general aviation type aircraft to improve crashworthiness performance. The energy-absorbing beam consisted of a foam-filled cellular structure with glass fiber and hybrid glass/kevlar cell walls. Design, analysis, fabrication and testing of the keel beams prior to installation and subsequent full-scale crash testing of the aircraft are described. Factors such as material and fabrication constraints, damage tolerance, crush stress/strain response, seat-rail loading, and post crush integrity, which influenced the course of the design process are also presented. A theory similar to the one often used for ductile metal box structures was employed with appropriate modifications to estimate the sustained crush loads for the beams. This, analytical tool, coupled with dynamic finite element simulation using MSC.Dytran were the prime design and analysis tools. The validity of the theory as a reliable design tool was examined against test data from static crush tests of beam sections while the overall performance of the energy-absorbing subfloor was assessed through dynamic testing of 24 in long subfloor assemblies.

Kellas, Sotiris↗

Experimental Verification of the Theoretical Prediction of the Phase Structure of a Ni-Al-Ti-Cr-Cu Alloy

The Bozzolo-Ferrante-Smith (BFS) method for alloys was applied to the study of NiAl-based materials to assess the effect of alloying additions on structure. Ternary, quaternary and even pentalloys based on NiAl with additions of Ti, Cr and Cu were studied and experimental verification of the theoretical predictions including the phase structure of a Ni-Al-Ti-Cr-Cu alloy is presented. Two approaches were used, Monte Carlo simulations to determine low energy structures, and analytical calculations of the energy of high symmetry configurations which give physical insight into preferred structures. The energetics for site occupancy in ternary and quaternary systems were calculated leading to an indirect determination of solubility limits at 0 K. Precipitate formation with information concerning structure and lattice parameter were also 'observed' computationally and the general characteristics of a Ni-Al-Ti-Cr-Cu alloy were correctly predicted. The results indicate that the BFS method for alloys can be a useful tool for alloy design and can be used to complement experimental alloy design programs.

Wilson, A.↗

The Dissipation Mechanism in Collisionless Magnetic Reconnection

The dissipation mechanism of magnetic reconnection remains a subject of intense scientific interest. On one hand, one set of recent studies have shown that particle inertia-based processes, which include thermal and bulk inertial effects, provide the reconnection electric field in the diffusion region. On the other hand, a second set of studies emphasizes the role of wave-particle interactions in providing anomalous resistivity in the diffusion region. In this presentation, we present analytical theory results, as well as PIC simulations of guide-field magnetic reconnection. We will show that the thermal electron inertia-based dissipation mechanism, expressed through nongyrotropic electron pressure tensors, remains viable in three dimensions. We will demonstrate the thermal inertia effect through studies of electron distribution functions. Furthermore, we will show that the reconnection electric field provides a transient acceleration on particles traversing the inner reconnection region. This inertial effect can be described as a diffusion-like term of the current density, which matches key features of electron distribution functions.

Hesse, Michael↗

Joint Bayesian Component Separation and CMB Power Spectrum Estimation

We describe and implement an exact, flexible, and computationally efficient algorithm for joint component separation and CMB power spectrum estimation, building on a Gibbs sampling framework. Two essential new features are (1) conditional sampling of foreground spectral parameters and (2) joint sampling of all amplitude-type degrees of freedom (e.g., CMB, foreground pixel amplitudes, and global template amplitudes) given spectral parameters. Given a parametric model of the foreground signals, we estimate efficiently and accurately the exact joint foreground- CMB posterior distribution and, therefore, all marginal distributions such as the CMB power spectrum or foreground spectral index posteriors. The main limitation of the current implementation is the requirement of identical beam responses at all frequencies, which restricts the analysis to the lowest resolution of a given experiment. We outline a future generalization to multiresolution observations. To verify the method, we analyze simple models and compare the results to analytical predictions. We then analyze a realistic simulation with properties similar to the 3 yr WMAP data, downgraded to a common resolution of 3 deg FWHM. The results from the actual 3 yr WMAP temperature analysis are presented in a companion Letter.

numerical↗

The Electron Diffusion Region: Forces and Currents

The dissipation mechanism of magnetic reconnection remains a subject of intense scientific interest. On one hand, one set of recent studies have shown that particle inertia-based processes, which include thermal and bulk inertial effects, provide the reconnection electric field in the diffusion region. On the other hand, a second set of studies emphasizes the role of wave-particle interactions in providing anomalous resistivity in the diffusion region. In this presentation, we present analytical theory results, as well as PIC simulations of guide-field magnetic reconnection. We will show that the thermal electron inertia-based dissipation mechanism, expressed through nongyrotropic electron pressure tensors, remains viable in three dimensions. We will demonstrate the thermal inertia effect through studies of electron distribution functions. Furthermore, we will show that the reconnection electric field provides a transient acceleration on particles traversing the inner reconnection region. This inertia1 effect can be described as a diffusion-like term of the current density, which matches key features of electron distribution functions.

Hesse, Michael↗

What Turns Reconnection On and Off

The dissipation mechanism of magnetic reconnection remains a subject of intense scientific interest. On one hand, one set of recent studies have shown that particle inertia-based processes, which include thermal and bulk inertial effects, provide the reconnection electric field in the diffusion region. On the other hand, a second set of studies emphasizes the role of wave-particle interactions in providing anomalous resistivity in the diffusion region. In this presentation, we present analytical theory results, as well as PIC simulations of guide-field magnetic reconnection. We will show that the thermal electron inertia-based dissipation mechanism, expressed through nongyrotropic electron pressure tensors, remains viable in three dimensions. We will demonstrate the thermal inertia effect through studies of electron distribution functions. Furthermore, we will show that the reconnection electric field provides a transient acceleration on particles traversing the inner reconnection region. This inertial effect can be described as a diffusion-like term of the current density. which matches key features of electron distribution functions.

Hesse, Michael↗

Theory and Simulations of Solar System Plasmas

"Theory and simulations of solar system plasmas" aims to highlight results from microscopic to global scales, achieved by theoretical investigations and numerical simulations of the plasma dynamics in the solar system. The theoretical approach must allow evidencing the universality of the phenomena being considered, whatever the region is where their role is studied; at the Sun, in the solar corona, in the interplanetary space or in planetary magnetospheres. All possible theoretical issues concerning plasma dynamics are welcome, especially those using numerical models and simulations, since these tools are mandatory whenever analytical treatments fail, in particular when complex nonlinear phenomena are at work. Comparative studies for ongoing missions like Cassini, Cluster, Demeter, Stereo, Wind, SDO, Hinode, as well as those preparing future missions and proposals, like, e.g., MMS and Solar Orbiter, are especially encouraged.

Goldstein, Melvyn L.↗

The Electron Diffusion Region: Forces and Currents

The dissipation mechanism of magnetic reconnection remains a subject of intense scientific interest. On one hand, one set of recent studies have shown that particle inertia-based processes, which include thermal and bulk inertial effects, provide the reconnection electric field in the diffusion region. On the other hand, a second set of studies emphasizes the role of wave-particle interactions in providing anomalous resistivity in the diffusion region. In this presentation, we present analytical theory results, as well as PIC simulations of guide-field magnetic reconnection. We will show that the thermal electron inertia-based dissipation mechanism, expressed through nongyrotropic electron pressure tensors, remains viable in three dimensions. We will demonstrate the thermal inertia effect through studies of electron distribution functions. Furthermore, we will show that the reconnection electric field provides a transient acceleration on particles traversing the inner reconnection region. This inertial effect can be described as a diffusion-like term of the current density, which matches key features of electron distribution functions.

Hesse, Michael↗

Understanding the Impact of Return-Current Losses on the X-Ray Emission from Solar Flares

I obtain and examine the implications of one-dimensional analytic solutions for return-current losses on an initially power-law distribution of energetic electrons with a sharp low-energy cutoff in flare plasma with classical (collisional) resistivity. These solutions show, for example, that return-current losses are not sensitive to plasma density, but are sensitive to plasma temperature and the low energy cutoff of the injected nonthermal electron distribution. A characteristic distance from the electron injection site, x(sub rc), is derived. At distances less than x(sub rc) the electron flux density is not reduced by return-current losses, but plasma heating can be substantial in this region, in the upper, coronal part of the flare loop. Before the electrons reach the collisional thick-target region of the flare loop, an injected power-law electron distribution with a low-energy cutoff maintains that structure, but with a flat energy distribution below the cutoff energy, which is now determined by the total potential drop experienced by the electrons. Modifications due to the presence of collisional losses are discussed. I compare these results with earlier analytical results and with more recent numerical simulations. Emslie's 1980 conjecture that there is a maximum integrated X-ray source brightness on the order of 10(exp -15) photons per square centimeter per second per square centimeter is examined. I find that this is not actually a maximum brightness and its value is parameter dependent, but it is nevertheless a valuable benchmark for identifying return-current losses in hard X-ray spectra. I discuss an observational approach to identifying return-current losses in flare data, including identification of a return-current "bump" in X-ray light curves at low photon energies.

Holman, Gordon D.↗

Macro Scale Independently Homogenized Subcells for Modeling Braided Composites

An analytical method has been developed to analyze the impact response of triaxially braided carbon fiber composites, including the penetration velocity and impact damage patterns. In the analytical model, the triaxial braid architecture is simulated by using four parallel shell elements, each of which is modeled as a laminated composite. Currently, each shell element is considered to be a smeared homogeneous material. The commercial transient dynamic finite element code LS-DYNA is used to conduct the simulations, and a continuum damage mechanics model internal to LS-DYNA is used as the material constitutive model. To determine the stiffness and strength properties required for the constitutive model, a top-down approach for determining the strength properties is merged with a bottom-up approach for determining the stiffness properties. The top-down portion uses global strengths obtained from macro-scale coupon level testing to characterize the material strengths for each subcell. The bottom-up portion uses micro-scale fiber and matrix stiffness properties to characterize the material stiffness for each subcell. Simulations of quasi-static coupon level tests for several representative composites are conducted along with impact simulations.

Blinzler, Brina J.↗

On the Electron Diffusion Region in Asymmetric Reconnection with a Guide Magnetic Field

Particle-in-cell simulations in a 2.5-D geometry and analytical theory are employed to study the electron diffusion region in asymmetric reconnection with a guide magnetic field. The analysis presented here demonstrates that similar to the case without guide field, in-plane flow stagnation and null of the in-plane magnetic field are well separated. In addition, it is shown that the electric field at the local magnetic X point is again dominated by inertial effects, whereas it remains dominated by nongyrotropic pressure effects at the in-plane flow stagnation point. A comparison between local electron Larmor radii and the magnetic gradient scale lengths predicts that distribution should become nongyrotropic in a region enveloping both field reversal and flow stagnation points. This prediction is verified by an analysis of modeled electron distributions, which show clear evidence of mixing in the critical region.

reconnection↗

Analytic Spacecraft Attitude and Rate Estimation Performance During Attitude Sensor Outages

Analytic expressions for spacecraft attitude and rate estimation performance of an attitude estimation filter in terms of sensor specifications are useful tools for spacecraft design. Farrenkopf (1978) famously found analytic expressions for steady-state pre-update and post-update attitude and gyro bias estimate error variances for an attitude estimation filter for a single-axis spacecraft with a Rate Output Gyro (ROG). Markley and Reynolds (2000) extended the analysis for a Rate-Integrating Gyro (RIG) with angle white noise. These expressions allow for the rapid evaluation of system performance during preliminary mission design phases. One contribution of this paper is the analytic calculation of the steady-state pre-update and post-update angular rate estimate uncertainty for both the ROG and RIG cases. The primary contribution of this paper is the extension of the results for both the ROG and the RIG cases to the situation of an attitude sensor outage. This situation arises frequently in practice; for example when a star sensor’s field of view is occluded, when a star sensor’s readings are unreliable during a thruster burn that vibrates the spacecraft, or during star sensor outages due to radiation upsets. Analytic expressions for the attitude estimate uncertainty, gyro bias estimate uncertainty, and angular rate estimate uncertainty are given in terms of the attitude sensor outage interval, the star tracker measurement noise, and gyro noise parameters. Validity of the analytic results is demonstrated via Monte Carlo simulation.

Galante, Joseph M.↗

Analytic Spacecraft Attitude and Rate Estimation Performance During Attitude Sensor Outages

Analytic expressions for spacecraft attitude and rate estimation performance of an attitude estimation filter in terms of sensor specifications are useful tools for spacecraft design. Farrenkopf (1978) famously found analytic expressions for steady-state pre-update and post-update attitude and gyro bias estimate error variances for an attitude estimation filter for a single-axis spacecraft with a Rate Output Gyro (ROG). Markley and Reynolds (2000) extended the analysis for a Rate-Integrating Gyro (RIG) with angle white noise. These expressions allow for the rapid evaluation of system performance during preliminary mission design phases. One contribution of this paper is the analytic calculation of the steady-state pre-update and post-update angular rate estimate uncertainty for both the ROG and RIG cases. The primary contribution of this paper is the extension of the results for both the ROG and the RIG cases to the situation of an attitude sensor outage. This situation arises frequently in practice; for example when a star sensor’s field of view is occluded, when a star sensor’s readings are unreliable during a thruster burn that vibrates the spacecraft, or during star sensor outages due to radiation upsets. Analytic expressions for the attitude estimate uncertainty, gyro bias estimate uncertainty, and angular rate estimate uncertainty are given in terms of the attitude sensor outage interval, the star tracker measurement noise, and gyro noise parameters. Validity of the analytic results is demonstrated via Monte Carlo simulation.

Galante, Joseph M.↗

Experimental and Computational Investigation of Flow Boiling in Microgravity

This study explores use of Computational Fluid Dynamics (CFD) to predict near-saturated flow boiling of FC-72 in microgravity. The computational method employs transient analysis to predict detailed interfacial behavior and heat transfer characteristics along a rectangular channel heated along two opposite walls. Predicted results are validated against experimental temperature measurements and high-speed video images captured during a series of parabolic aircraft maneuvers for three sets of operating conditions which include variations of both mass velocity and wall heat flux. The computational method is based on the multi-phase volume of fluid (VOF) model, which is combined with appropriate phase change and turbulence models, and accounts for both shear-lift force on bubbles and conjugate heat transfer along the heating walls. A key advantage of the CFD method is ability to capture details that are very difficult to measure experimentally, including detailed spatial variations of bubble shape, void fraction, mixture fluid temperature, liquid velocity, and vapor velocity, results for which are presented for each of the three test cases. Different flow regimes predicted along the heated length exhibit a number of dominant mechanisms including bubble nucleation, bubble growth, coalescence, vapor blankets, interfacial waviness, and residual liquid sub-layer, all of which agree well with experiment. Vapor velocity is shown to increase appreciably along the heated length because of increased void fraction, while liquid velocity experiences large fluctuations. Non-equilibrium effects are accentuated with increasing mass velocity, contributing minor deviations of fluid temperature from simulations compared to those predicted by the analytical method. Predicted wall temperature is fairly uniform in the middle of the heated length but increases in the entrance region, due to sensible heat transfer in the subcooled liquid, and decreases toward the exit, mostly because of flow acceleration resulting from increased void fraction.

two-phase cooling; flow boiling; CFD; microgravity↗

LUVOIR-ECLIPS Closed-Loop Adaptive Optics Performance and Contrast Predictions

One of the primary science goals of the Large UV/Optical/Infrared Surveyor (LUVOIR) mission concept is to detect and characterize Earth-like exoplanets orbiting nearby stars with direct imaging. The success of its coronagraph instrument ECLIPS (Extreme Coronagraph for Living Planetary Systems) depends on the ability to stabilize the wavefront from a large segmented mirror such that optical path differences are limited to tens of picometers RMS during an exposure time of a few hours. In order to relax the constraints on the mechanical stability, ECLIPS will be equipped with a wavefront sensing and control (WS\&C) architecture to correct wavefront errors up to temporal frequencies $\gtrsim$1 Hz. These errors may be dominated by spacecraft structural dynamics exciting vibrations at the segmented primary mirror. In this work, we present detailed simulations of the WS\&C system within the ECLIPS instrument and the resulting contrast performance. This study assumes wavefront aberrations based on a finite element model of a simulated telescope with spacecraft structural dynamics. Wavefront residuals are then computed according to a model of the adaptive optics system that includes numerical propagation to simulate a realistic wavefront sensor and an analytical model of the temporal performance. An end-to-end numerical propagation model of ECLIPS is then used to estimate the residual starlight intensity distribution at the science detector. We show that the contrast performance depends strongly on the target star magnitude and the spatio-temporal distribution of wavefront errors from the telescope. In cases with significant vibration, we advocate for the use of laser metrology to mitigate high temporal frequency wavefront errors and increase the mission yield.

Sirbu, Dan↗

Application of Long-Path Gas Cell Assisted Infrared Spectroscopy for Detection of Lunar Volatiles and Advanced Curation Efforts

The motivation behind analyzing gaseous samples from the Permanently Shadowed Regions (PSRs) on the lunar surface is built upon the volatiles observed during the LCROSS mission. Advanced curation techniques under cold conditions are essential to the critical scientific analysis of these returned samples. The detection of various volatiles in returned samples by gas-phase infrared (IR) spectroscopy is known to be a powerful tool to understand the resources present in planetary bodies and Astro-materials. The identification and quantification of volatile compounds in returned planetary samples and regolith simulants by advanced and custom-built analytical devices are areas of primary research interest in our group and the Planetary Exploration & Astromaterials Research Laboratory (PEARL) facility. In this context, we discuss recent advancements and experimental results for quick analysis of the major volatiles detected by the LCROSS (water, CO 2 , CH 3 OH, and NH 3 ) by a long-path gas cell/FTIR (Fourier Transform Infrared) spectrometer system. The applicability of this system is demonstrated by qualitative analysis of a series of volatiles measured individually. In addition, the linear dependence of absorption intensity on headspace sample volumes was checked by measuring gas-phase IR absorption spectra of increasing volumes of pure CO 2 gas aliquots. Finally, the instrument performance was evaluated by measuring the relative concentration of NH 3 (g) in methanol solution.

Curation↗