Engineering PapersSearch

SEARCH · Engineering Papers

Results for “Reduced Order Modeling”

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 1,045 records · Page 58

Uncertainty and Sensitivity Analysis of Afterbody Radiative Heating Predictions for Earth Entry

The objective of this work was to perform sensitivity analysis and uncertainty quantification for afterbody radiative heating predictions of Stardust capsule during Earth entry at peak afterbody radiation conditions. The radiation environment in the afterbody region poses significant challenges for accurate uncertainty quantification and sensitivity analysis due to the complexity of the flow physics, computational cost, and large number of un-certain variables. In this study, first a sparse collocation non-intrusive polynomial chaos approach along with global non-linear sensitivity analysis was used to identify the most significant uncertain variables and reduce the dimensions of the stochastic problem. Then, a total order stochastic expansion was constructed over only the important parameters for an efficient and accurate estimate of the uncertainty in radiation. Based on previous work, 388 uncertain parameters were considered in the radiation model, which came from the thermodynamics, flow field chemistry, and radiation modeling. The sensitivity analysis showed that only four of these variables contributed significantly to afterbody radiation uncertainty, accounting for almost 95% of the uncertainty. These included the electronic- impact excitation rate for N between level 2 and level 5 and rates of three chemical reactions in uencing N, N(+), O, and O(+) number densities in the flow field.

West, Thomas K., IV

Modeling Studies of PVT Growth of ZnSe: Current Status and Future Course

Bulk growth of wide band gap II-VI semiconductors by physical vapor transport (PVT) has been developed and refined over the past several years at NASA Marshall Space Flight Center. Results from a modeling study of PVT crystal growth of ZnSe are reported in this paper. The PVT process is numerically investigated using a two-dimensional formulation of the governing equations and associated boundary conditions. Both the incompressible Boussinesq approximation and a compressible model are tested to determine the influence of gravity on the process and to discern the differences between the two approaches. The influence of a residual gas is included in the models. The results show that both the incompressible and compressible approximations provide comparable results and the presence of a residual gas tends to measurably reduce the mass flux in the system. Detailed flow, thermal and concentration profiles are provided. The simulations show that the Stefan flux dominates the system flow field and the subtle gravitational effects can be gauged by subtracting this flux from the calculated profiles. Shear flows, due to solutal buoyancy, of the order of 50 microns/s for the liorizont,-d growth orientation and 10 microns/s for the vertical orientation are predicted. Whether these flows can fully account for the observed gravity related growth morphological effects and inhomogeneous solute and dopant distributions is a matter of conjecture. A template for future modeling efforts in this area is suggested which incorporates a mathematical approach to the tracking of the growth front based on energy of formation concepts.

Ramachandran, N.

Bubble Generation in a Flowing Liquid Medium and Resulting Two-Phase Flow in Microgravity

The present investigation reports a study of bubble generation under reduced gravity conditions, using both a co-flow and a cross-flow configuration. This study may be used in the conceptual design of a space-based thermal management system. Ensuing two-phase flow void fraction can be accurately monitored using a single nozzle gas injection system within a continuous liquid flow conduit, as utilized in the present investigation. Accurate monitoring of void fraction leads to precise control of heat and mass transfer coefficients related to a thermal management system; hence providing an efficient and highly effective means of removing heat aboard spacecraft or space stations. Our experiments are performed in parabolic flight aboard the modified DC-9 Reduced Gravity Research Aircraft at NASA Lewis Research Center, using an air-water system. For the purpose of bubble dispersion in a flowing liquid, we use both a co-flow and a cross-flow configuration. In the co-flow geometry, air is introduced through a nozzle in the same direction with the liquid flow. On the other hand, in the cross-flow configuration, air is injected perpendicular to the direction of water flow, via a nozzle protruding inside the two-phase flow conduit. Three different flow conduit (pipe) diameters are used, namely, 1.27 cm, 1.9 cm and 2.54 cm. Two different ratios of nozzle to pipe diameter (D(sub N))sup * are considered, namely (D(sub N))sup * = 0.1 and 0.2, while superficial liquid velocities are varied from 8 to 70 cm/s depending on flow conduit diameter. It is experimentally observed that by holding all other flow conditions and geometry constant, generated bubbles decrease in size with increase in superficial liquid velocity. Detached bubble diameter is shown to increase with air injection nozzle diameter. Likewise, generated bubbles grow in size with increasing pipe diameter. Along the same lines, it is shown that bubble frequency of formation increases and hence the time to detachment of a forming bubble decreases, as the superficial liquid velocity is in-creased. Furthermore, it is shown that the void fraction of the resulting two-phase flow increases with volumetric gas flow rate Q(sub d), pipe diameter and gas injection nozzle diameter, while they decrease with surrounding liquid flow. The important role played by flowing liquid in detaching bubbles in a reduced gravity environment is thus emphasized. We observe that the void fraction can be accurately controlled by using single nozzle gas injection, rather than by employing multiple port injection, since the later system gives rise to unpredictable coalescence of adjacent bubbles. It is of interest to note that empirical bubble size and corresponding void fraction are somewhat smaller for the co-flow geometry than the cross-flow configuration at similar flow conditions with similar pipe and nozzle diameters. In order to supplement the empirical data, a theoretical model is employed to study single bubble generation in the dynamic (Q(sub d) = 1 - 1000 cu cm/s) and bubbly flow regime within the framework of the co-flow configuration. This theoretical model is based on an overall force balance acting on the bubble during the two stages of generation, namely the expansion and the detachment stage. Two sets of forces, one aiding and the other inhibiting bubble detachment are identified. Under conditions of reduced gravity, gas momentum flux enhances, while the surface tension force at the air injection nozzle tip inhibits bubble detachment. In parallel, liquid drag and inertia can act as both attaching and detaching forces, depending on the relative velocity of the bubble with respect to the surrounding liquid. Predictions of the theoretical model compare well with our experimental results. However, at higher superficial liquid velocities, as the bubble loses its spherical form, empirical bubble size no longer matches the theoretical predictions. In summary, we have developed a combined experimental and theoretical work, which describes the complex process of bubble generation and resulting two-phase flow in a microgravity environment. Results of the present study can be used in a wide range of space-based applications, such as thermal energy and power generation, propulsion, cryogenic storage and long duration life support systems, necessary for programs such as NASA's Human Exploration for the Development of Space (HEDS).

Pais, S. C.

Benchmarking MCNP models to support gamma irradiation experiments

Seismic isolators and damping devices have been used for 30+ years to reduce seismic forces and deployed in 10,000+ structures, including bridges, buildings, and mission-critical infrastructure. Studies have shown that seismic isolation of nuclear power plants can reduce the seismic risk by several orders of magnitude and substantially reduce the overnight capital cost. Advanced reactor developers are considering seismic isolation as an integral design feature in their reactor design. It can either be deployed at the building level, isolating the entire reactor building, or at the component level, isolating individual pieces of safety-class equipment, such as a reactor vessel or a steam generator. For equipment isolation, located close to the reactor vessel, the seismic isolators and/or damping devices may be exposed to gamma and neutron radiation, with a possible effect on their mechanical properties which would affect the response of the isolated equipment under earthquake shaking. The Department of Energy, as part of its Nuclear Energy University Program (NEUP), in collaboration with the Idaho National Laboratory (INL), is funding a project at the University at Buffalo (UB) to investigate the effect of gamma radiation on the mechanical properties of various seismic protective devices. (When isolating equipment, the seismic protective devices will be located outside the reactor vessel where the neutron exposure will be insignificant.

22 - GENERAL STUDIES OF NUCLEAR REACTORS

An In Situ , Automated High-Explosives Aging Method Utilizing Two-Dimensional Gas Chromatography–Mass Spectrometry

Understanding chemical changes that occur in high explosives as they age is of great importance to the safe employment and storage of these compounds. Traditional methods of aging high explosives even under accelerated aging conditions are time intensive with durations on the order of months to years. The nature of traditional aging analyses reduces each sample to a snapshot data point often separated widely in time, requiring many assumptions as to how the degradation products develop. Further complicating matters, several analytical techniques are typically employed for each sample analysis in order to ascertain an entire picture of the decomposition pathways. To address these shortcomings with existing methods, a new method of accelerated aging of high explosives utilizing comprehensive two-dimensional gas chromatography coupled to high-resolution mass spectrometry (GC × GC-HRMS) was developed using 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) as a model compound for method development. This in situ automated method reduces the time scale of aging to a matter of hours using the inlet of the GC × GC as the aging vessel. GC × GC in combination with HRMS allowed for the collection of both evolved gases and other decomposition products produced during the entire aging process in real time with HRMS providing far greater certainty in identification of explosives aging products. Additionally, this method allowed for a higher throughput of samples with greatly simplified sample preparation. Chemometric analysis of the GC × GC-HRMS data set via the alteration analysis (ALA) enabled discovery of statistically significant chemical changes providing insight into the variation of decomposition pathways with varying aging temperatures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH

Iterative Mechanism Solutions with Scenario and ADAMS

This slide presentation reviews the use of iterative solutions using Scenario for Motion (UG NX 2 Motion) to assist in designing the Mars Science Laboratory (MSL). The MSL will have very unique design requirements, and in order to meet these requirements the system must have the ability to design for static stability, simulate mechanism kinematics, simulate dynamic behaviour and be capable of reconfiguration, and iterations as designed. The legacy process used on the Mars Exploration rovers worked, but it was cumbersome using multiple tools, limited configuration control, with manual process and communication, and multiple steps. The aim is to develop a mechanism that would reduce turn around time, and make more reiterations possible, to improve the quality and quantity of data, and to enhance configuration control. Currently for NX Scenario for Motion uses are in the articulation studies, the simulations of traverse motions,and subsystem simulations. The design of the Rover landing model requires accurate results, flexible elements, such as beams, and the use of the full ADAMS solver has been used. In order to achieve this, when required, there has been a direct translation from Scenario to ADAMS, with additional data in ascii format. The process that has been designed to move from Scenario to ADAMS is reviewed.

Mars Science Laboratory (MSL)

Characteristics of Non-Premixed Turbulent Flames in Microgravity

The momentum of the fuel (and/or air) jet is important in classifying gas-jet diffusion flame behavior. Normal-gravity data on gas-jet flames show that the flame height (non-dimensionalized with respect to an effective diameter) can be correlated to a density weighted Froude number in the buoyancy-dominated limit. In the momentum-dominated limit this non-dimensional flame height asymptotes to a constant value. The momentum-dominated limit under normal gravity conditions is usually obtained for very high injection velocities which in turn results in high values of the injection Reynolds number. This results in a complicated flame structure because of the large number of turbulence scales involved. In order to gain better insight into the structure of these flames it would be useful to reduce the injection Reynolds number while still maintaining turbulent conditions. This can be done in microgravity where momentum-dominated turbulent flames are obtained at much smaller velocities than in normal gravity. In this paper, experimental results on the effects of nozzle diameter and fuel dilution on flame height are discussed. The experimental values are compared with predictions from a numerical procedure utilizing the standard k-epsilon turbulence model. Flame height scaling with nozzle size and dilution is established. Differences between model predictions and measurements are presented. In order to explain these differences, evolutions of turbulent spectra and Taylor microscale along the flame axis are considered.

Hegde, U.

Modeling cross-beam energy transfer with sector ray tracing

Ray-based cross-beam energy transfer (CBET) models are an essential feature of the radiation-hydrodynamic codes used to simulate inertial confinement fusion implosions, but full 3D ray-based CBET calculations can have a prohibitively high computational cost. Sector ray tracing can be used to reduce the cost by orders of magnitude in cases where the coronal plasma and laser drive can be approximated as spherically symmetric. An extension of sector ray tracing (section ray tracing) can be used to relax the assumption of a spherically symmetric drive while still retaining most of the computational savings of a sector ray trace. We discuss the foundations of sector and section ray tracing and compare them to full ray tracing.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY

Structure of a fluid interface near the critical point

The structure of the liquid-liquid interface of three very different mixtures (carbon disulfide + methanol, methanol + cyclohexane + deuterated cyclohexane, and nitrobenzene + n-decane) has been studied using ellipsometry in the reduced temperature range between 0.0009 and 0.042. Although the elliptical thickness varies by an order of magnitude between these mixtures, the data from all three mixtures can be scaled to the same universal constant by a combined mean-field plus capillary-wave model of the interface. The universal constant determined experimentally is significantly less than the theoretical value.

Schmidt, James W.

Guided waves in a multilayered composite and ultrasonic nondestructive evaluation

Analysis of wave propagation in a laminated composite plate with many layers is made difficult by the anisotropy of the laminae, different lay-ups, many layers and defects in the fabricated plate. A good understanding of the wave field in such a plate without defects is a prerequisite for use of ultrasonic techniques to evaluate the various defects that can substantially reduce the strength and service life of this structure. Results of modeling studies of wave propagation and scattering in a multilayered composite plate are presented. The analysis technique combines the finite element discretization through the thickness with wave form representation along the plate in order to calculate the dispersive modes of propagation in the plate. To study scattering of these waves by a matrix crack that grows into delaminations, a hybrid finite element representation of the field near the crack and the modal representation of the scattered exterior field are used. Results are presented showing the dependence of the reflection and transmission coefficients on the size of the crack with delamination which can be used to size the defect using ultrasonic techniques.

Datta, S. K.

Reduced Stress Tensor and Dissipation and the Transport of Lamb Vector

We develop a methodology to ensure that the stress tensor, regardless of its number of independent components, can be reduced to an exactly equivalent one which has the same number of independent components as the surface force. It is applicable to the momentum balance if the shear viscosity is constant. A direct application of this method to the energy balance also leads to a reduction of the dissipation rate of kinetic energy. Following this procedure, significant saving in analysis and computation may be achieved. For turbulent flows, this strategy immediately implies that a given Reynolds stress model can always be replaced by a reduced one before putting it into computation. Furthermore, we show how the modeling of Reynolds stress tensor can be reduced to that of the mean turbulent Lamb vector alone, which is much simpler. As a first step of this alternative modeling development, we derive the governing equations for the Lamb vector and its square. These equations form a basis of new second-order closure schemes and, we believe, should be favorably compared to that of traditional Reynolds stress transport equation.

Wu, Jie-Zhi

Assessment of the Mars 2020 Entry, Descent, and Landing Simulation

On February 18, 2021, the Mars 2020 Perseverance rover successfully landed inside Jezero Crater at 18.44463 deg North latitude and 77.45088 deg East longitude. At 1026 kg, Perseverance is the largest, most sophisticated rover ever delivered to another planet. This event marked the ninth successful landing and fifth rover to be delivered at Mars. The Program to Optimize Simulated Trajectories II (POST2) is a trajectory simulation tool maintained by the NASA Langley Research Center. POST2 was the prime EDL performance simulation for Mars 2020. This tool has significant heritage on the previous Pathfinder, MER, Phoenix, InSight, and MSL landings. This paper presents a few initial comparisons between EDL flight telemetry and POST2 simulation predictions. These comparisons are important in order to understand how each of the individual models performed, as well as the integrated simulation as a whole. This information is fed forward to future missions, which benefit from the knowing where additional resources or study are needed and where uncertainties may be reduced to enable improved performance.

David W Way

Assessment of the Mars 2020 Entry, Descent, and Landing Simulation

On February 18, 2021, the Mars 2020 Perseverance rover successfully landed inside Jezero Crater at 18.44463 deg North latitude and 77.45088 deg East longitude. At 1026 kg, Perseverance is the largest, most sophisticated rover ever delivered to another planet. This event marked the ninth successful landing and fifth rover to be delivered at Mars. The Program to Optimize Simulated Trajectories II (POST2) is a trajectory simulation tool maintained by the NASA Langley Research Center. POST2 was the prime EDL performance simulation for Mars 2020. This tool has significant heritage on the previous Pathfinder, MER, Phoenix, InSight, and MSL landings. This paper presents a few initial comparisons between EDL flight telemetry and POST2 simulation predictions. These comparisons are important in order to understand how each of the individual models performed, as well as the integrated simulation as a whole. This information is fed forward to future missions, which benefit from the knowing where additional resources or study are needed and where uncertainties may be reduced to enable improved performance.

David Way

Adapting ASPEN for Orbital Express

By studying the Orbital Express mission, modeling the spacecraft and scenarios, and testing the system, a technique has been developed that uses recursive decomposition to represent procedural actions declaratively, schema-level uncertainty reasoning to make uncertainty reasoning tractable, and lightweight, natural language processing to automatically parse procedures to produce declarative models. Schema-level uncertainty reasoning has, at its core, the basic assumption that certain variables are uncertain, but not independent. Once any are known, then the others become known. This is important where a variable is uncertain for an action and many actions of the same type exist in the plan. For example, if the number of retries to purge pump lines was unknown (but bounded), and each attempt required a sub-plan, then, once the correct number of attempts required for a purge was known, it would likely be the same for all subsequent purges. This greatly reduces the space of plans that needs to be searched to ensure that all executions are feasible. To accommodate changing scenario procedures, each is ingested into a tabular format in temporal order, and a simple natural-language parser is used to read each step and to derive the impact of that step on memory, power, and communications. Then an ASPEN (Activity Scheduling and Planning Environment) model is produced based on this analysis. The model is tested and further changed by hand, if necessary, to reflect the actual procedure. This results in a great savings of time used for modeling procedures. Many processes that need to be modeled in ASPEN (a declarative system) are, in fact, procedural. ASPEN includes the ability to model activities in a hierarchical fashion, but this representation breaks down if there is a practically unbounded number of sub-activities and decomposition topologies. However, if recursive decomposition is allowed, HTN-like encodings are enabled to represent most procedural phenomena. For example, if a switch requires a variable (but known at the time of the attempt) number of attempts to switch on, one can recurse on the number of remaining switch attempts and decompose into either the same switching activity with one less required attempt, or not decompose at all (or decompose into a dummy task), resulting in the end of the decomposition. In fact, any bounded procedural behavior can be modeled using recursive decompositions assuming that the variables impinging the disjunctive decomposition decision are computable at the time that the decision is made. This enables one to represent tasks that are controlled outside of the scheduler, but that the scheduler must accommodate, without requiring one to give a declarative model of the procedural behavior.

Chouinard, Caroline

Rare Earth Element Partitioning in Lunar Minerals: An Experimental Study

The partitioning behavior of rare earth elements (REE) between minerals and melts is widely used to interpret the petrogenesis and geologic context of terrestrial and extra-terrestrial samples. REE are important tools for modelling the evolution of the lunar interior. The ubiquitous negative Eu anomaly in lunar basalts is one of the main lines of evidence to support the lunar magma ocean (LMO) hypothesis, by which the plagioclase-rich lunar highlands were formed as a flotation crust during differentiation of a global-scale magma ocean. The separation of plagioclase from the mafic cumulates is thought to be the source of the Eu depletion, as Eu is very compatible in plagioclase. Lunar basalts and volcanic glasses are commonly depleted in light REEs (LREE), and more enriched in heavy REEs (HREE). However, there is very little experimental data available on REE partitioning between lunar minerals and melts. In order to interpret the source of these distinctive REE patterns, and to model lunar petrogenetic processes, REE partition coefficients (D) between lunar minerals and melts are needed at conditions relevant to lunar processes. New data on D(sub REE) for plagioclase, and pyroxenes are now available, but there is limited available data for olivine/melt D(sub REE), particularly at pressures higher than 1 bar, and in Fe-rich and reduced compositions - all conditions relevant to the lunar mantle. Based on terrestrial data, REE are highly incompatible in olivine (i.e. D much less than 1), however olivine is the predominant mineral in the lunar interior, so it is important to understand whether it is capable of storing even small amounts of REE, and how the REEs might be fractionatied, in order to understand the trace element budget of the lunar interior. This abstract presents results from high-pressure and temperature experiments investigating REE partitioning between olivine and melt in a composition relevant to lunar magmatism.

McIntosh, E. C.

PSCF+: An Extended and Improved Open-Source Software Package for Polymer Self-Consistent Field Calculations

This work introduces PSCF+, a freely available open‑source software tool for modeling how block copolymers self‑assemble into ordered nanoscale structures. The software allows researchers to study a wide range of polymer models and interaction types, making it possible to directly compare theoretical predictions with molecular and simulation‑based studies in a consistent way. PSCF+ incorporates several recently developed numerical methods that significantly reduce computer memory requirements and speed up calculations, enabling efficient use of modern GPU hardware. Here, it also provides automated tools to calculate free‑energy landscapes and phase boundaries, which are essential for understanding and predicting material behavior. Together, these advances make PSCF+ a powerful and flexible platform for studying polymer self‑assembly and support the design of new polymeric materials for scientific and technological applications.

Chemical calculations

Economic Storage Size Optimization for Electric Vehicle Extreme-Fast Charging Stations

En-route charging infrastructure for electric vehicles is critical to support transportation needs. These charging stations are likely to have high loads and especially sharp peak loads given fast charging capabilities needed to meet transportation schedules. In order to reduce both strain on distribution grid infrastructure and charging station operational costs, many stations are likely to employ behind the meter storage. This paper demonstrates a behind the meter storage sizing optimization that employs an open-source agent-based vehicle behavior model (BEAM) to determine the best sizing across many scenarios. This optimization and analysis is novel in that it examines how storage size impacts not only charging station cost and peak load, but also vehicle queue times. The optimization is also applied across a wide analysis region with sufficient diversity and numbers to provide novel statistical analysis of optimal sizes.

Aka, Julius

Computer assisted thermal-vacuum testing

In testing complex systems and components under dynamic thermal-vacuum environments, it is desirable to optimize the environment control sequence in order to reduce test duration and cost. This paper describes an approach where a computer is utilized as part of the test control operation. Real time test data is made available to the computer through time-sharing terminals at appropriate time intervals. A mathematical model of the test article and environmental control equipment is then operated on using the real time data to yield current thermal status, temperature analysis, trend prediction and recommended thermal control setting changes to arrive at the required thermal condition. The data acquisition interface and the time-sharing hook-up to an IBM-370 computer is described along with a typical control program and data demonstrating its use.

Petrie, W.