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 793 records · Page 44

Optimal attitude maneuver execution for the Advanced Composition Explorer (ACE) mission

The Advanced Composition Explorer (ACE) spacecraft will require frequent attitude reorientations in order to maintain the spacecraft high gain antenna (HGA) within 3 deg of earth-pointing. These attitude maneuvers will be accomplished by employing a series of ground-commanded thruster pulses, computed by ground operations personnel, to achieve the desired change in the spacecraft angular momentum vector. With each maneuver, attitude nutation will be excited. Large nutation angles are undesirable from a science standpoint. It is important that the thruster firings be phased properly in order to minimize the nutation angle at the end of the maneuver so that science collection time is maximized. The analysis presented derives a simple approximation for the nutation contribution resulting from a series of short thruster burns. Analytic equations are derived which give the induced nutation angle as a function of the number of small thruster burns used to execute the attitude maneuver and the phasing of the burns. The results show that by properly subdividing the attitude burns, the induced nutation can be kept low. The analytic equations are also verified through attitude dynamics simulation and simulation results are presented. Finally, techniques for quantifying the post-maneuver nutation are discussed.

Woodard, Mark A.↗

Automotive Design

Analytical Design Service Corporation, Ann Arbor, MI, used NASTRAN (a NASA Structural Analysis program that analyzes a design and predicts how parts will perform) in tests of transmissions, engine cooling systems, internal engine parts, and body components. They also use it to design future automobiles. Analytical software can save millions by allowing computer simulated analysis of performance even before prototypes are built.

Source record↗

Dynamics of Quantum Adiabatic Evolution Algorithm for Number Partitioning

We have developed a general technique to study the dynamics of the quantum adiabatic evolution algorithm applied to random combinatorial optimization problems in the asymptotic limit of large problem size n. We use as an example the NP-complete Number Partitioning problem and map the algorithm dynamics to that of an auxiliary quantum spin glass system with the slowly varying Hamiltonian. We use a Green function method to obtain the adiabatic eigenstates and the minimum exitation gap, gmin = O(n2(sup -n/2)), corresponding to the exponential complexity of the algorithm for Number Partitioning. The key element of the analysis is the conditional energy distribution computed for the set of all spin configurations generated from a given (ancestor) configuration by simultaneous flipping of a fixed number of spins. For the problem in question this distribution is shown to depend on the ancestor spin configuration only via a certain parameter related to the energy of the configuration. As the result, the algorithm dynamics can be described in terms of one-dimensional quantum diffusion in the energy space. This effect provides a general limitation of a quantum adiabatic computation in random optimization problems. Analytical results are in agreement with the numerical simulation of the algorithm.

Smelyanskiy, Vadius↗

Dynamics of Quantum Adiabatic Evolution Algorithm for Number Partitioning

We have developed a general technique to study the dynamics of the quantum adiabatic evolution algorithm applied to random combinatorial optimization problems in the asymptotic limit of large problem size n. We use as an example the NP-complete Number Partitioning problem and map the algorithm dynamics to that of an auxiliary quantum spin glass system with the slowly varying Hamiltonian. We use a Green function method to obtain the adiabatic eigenstates and the minimum excitation gap. g min, = O(n 2(exp -n/2), corresponding to the exponential complexity of the algorithm for Number Partitioning. The key element of the analysis is the conditional energy distribution computed for the set of all spin configurations generated from a given (ancestor) configuration by simultaneous flipping of a fixed number of spins. For the problem in question this distribution is shown to depend on the ancestor spin configuration only via a certain parameter related to 'the energy of the configuration. As the result, the algorithm dynamics can be described in terms of one-dimensional quantum diffusion in the energy space. This effect provides a general limitation of a quantum adiabatic computation in random optimization problems. Analytical results are in agreement with the numerical simulation of the algorithm.

Smelyanskiy, V. N.↗

Performance Modeling of Network-Attached Storage Device Based Hierarchical Mass Storage Systems

Network attached storage devices improve I/O performance by separating control and data paths and eliminating host intervention during the data transfer phase. Devices are attached to both a high speed network for data transfer and to a slower network for control messages. Hierarchical mass storage systems use disks to cache the most recently used files and a combination of robotic and manually mounted tapes to store the bulk of the files in the file system. This paper shows how queuing network models can be used to assess the performance of hierarchical mass storage systems that use network attached storage devices as opposed to host attached storage devices. Simulation was used to validate the model. The analytic model presented here can be used, among other things, to evaluate the protocols involved in 1/0 over network attached devices.

Menasce, Daniel A.↗

Maximum Likelihood Time-of-Arrival Estimation of Optical Pulses via Photon-Counting Photodetectors

Many optical imaging, ranging, and communications systems rely on the estimation of the arrival time of an optical pulse. Recently, such systems have been increasingly employing photon-counting photodetector technology, which changes the statistics of the observed photocurrent. This requires time-of-arrival estimators to be developed and their performances characterized. The statistics of the output of an ideal photodetector, which are well modeled as a Poisson point process, were considered. An analytical model was developed for the mean-square error of the maximum likelihood (ML) estimator, demonstrating two phenomena that cause deviations from the minimum achievable error at low signal power. An approximation was derived to the threshold at which the ML estimator essentially fails to provide better than a random guess of the pulse arrival time. Comparing the analytic model performance predictions to those obtained via simulations, it was verified that the model accurately predicts the ML performance over all regimes considered. There is little prior art that attempts to understand the fundamental limitations to time-of-arrival estimation from Poisson statistics. This work establishes both a simple mathematical description of the error behavior, and the associated physical processes that yield this behavior. Previous work on mean-square error characterization for ML estimators has predominantly focused on additive Gaussian noise. This work demonstrates that the discrete nature of the Poisson noise process leads to a distinctly different error behavior.

Erkmen, Baris I.↗

EMI from Spacecraft Docking Systems Spacecraft Charging - Plasma Contact Potentials

The plasma contact potential of a visiting vehicle (VV), such as the Orion Service Module (SM), is determined while docking at the Orion Crew Exploration Vehicle (CEV). Due to spacecraft charging effects on-orbit, the potential difference between the CEV and the VV can be large at docking, and an electrostatic discharge (ESD) could occur at capture, which could degrade, disrupt, damage, or destroy sensitive electronic equipment on the CEV and/or VV. Analytical and numerical models of the CEV are simulated to predict the worst-case potential difference between the CEV and the VV when the CEV is unbiased (solar panels unlit: eclipsed in the dark and inactive) or biased (solar panels sunlit: in the light and active).

Norgard, John D.↗

Stochastic Modelling of Wireless Energy Transfer

This study investigates the efficiency of a new method of powering remote sensors by the means of wireless energy transfer. The increased use of sensors for data collection comes with the inherent cost of supplying power from sources such as power cables or batteries. Wireless energy transfer technology eliminates the need for power cables or periodic battery replacement. The time and cost of setting up or expanding a sensor network will be reduced while allowing sensors to be placed in areas where running power cables or battery replacement is not feasible. This paper models wireless channels for power and data separately. Smart scheduling for the data channel is proposed to avoid transmitting data on a noisy channel where the probability of data loss is high to improve power efficiency. Analytical models have been developed and verified using simulations.

Veilleux, Shaun↗

Reduced Crew/Single Pilot Operations for Commercial Aircraft - Concept of Operations and Technology Needs

In 2012, NASA began exploring the feasibility of single pilot operations (SPO) and reduced crew operations (RCO) in the context of scheduled passenger air carrier operations (i.e., Parts 121 and 125). Technology and automation, especially aircraft automation, have significantly advanced in the 21st century and may be enabling to SPO. However, a move to SPO also has significant challenges. The purpose of a three-year NASA effort was to identify those challenges through workshops, analytic studies, and human-in-the-loop simulations assessing promising concepts and prototype solutions. This presentation will describe the progress that has been made in that three year effort.

reduced crew operations↗

Lightweight Thrust Chamber Assemblies using Multi-Alloy Additive Manufacturing and Composite Overwrap

Additive Manufacturing (AM) has brought significant design and fabrication opportunities for complex components with internal features such as liquid rocket engine thrust chambers not previously possible. This technology allows for significant cost savings and schedule reductions in addition to new performance optimization through weight reduction and increased margins. Specific to regeneratively-cooled combustion chambers and nozzles for liquid rocket engines, additive manufacturing offers the ability to form the complex internal coolant channels and the closeout of the channels to contain the high pressure liquid propellants with a single operation. Much of additive manufacturing development has focused on monolithic alloys using Laser Powder Bed Fusion (L-PBF), which do not allow for complete optimization of the structure. The National Aeronautics and Space Administration (NASA) completed feasibility of an AM bimetallic L-PBF GRCop-84 copper-alloy combustion chamber with an AM electron beam freeform Inconel 625 structural jacket under the Low Cost Upper Stage Propulsion (LCUSP) Project. A follow-on project called Rapid Analysis and Manufacturing Propulsion Technology (RAMPT) is under development to further expand large-scale multi-alloy thrust chambers while maturing composite overwrap technology for significant weight savings opportunities. The RAMPT project has three primary objectives: 1) Advancing blown powder Directed Energy Deposition (DED) to fabricate integral-channel large scale nozzles, 2) Develop composite overwrap technology to reduce weight and provide structural capability for thrust chamber assemblies, and 3) Develop bimetallic and multi-metallic additively manufactured radial and axial joints to optimize material performance. In addition to these primary manufacturing developments, analytical modeling efforts compliment the process development to simulate the AM processes to reduce build failures and distortions. The RAMPT project is also maturing the supply chain for various manufacturing processes described above in addition to L-PBF of GRCop-42. This paper will present an overview of the RAMPT project, the process development and hardware progress to date, material and hot-fire testing results, along with future developments.

Additive Manufacturing↗

The Next Generation X-Ray Galaxy Survey with eROSITA

eROSITA, launched on 13 July 2019, will be completing the first all-sky survey in the soft and medium X-ray band in nearly three decades. This 4-year survey, finishing in late 2023, will present a rich legacy for the entire astrophysics community and complement upcoming multi-wavelength surveys (with, e.g. the Large Synoptic Survey Telescope and the Dark Energy Survey). Besides the major scientific aim to study active galactic nuclei (AGN) and galaxy clusters, eROSITA will contribute significantly to X-ray studies of normal (i.e. not AGN) galaxies. Starting from multi-wavelength catalogues, we measure star formation rates and stellar masses for 60 212 galaxies constrained to distances of 50–200 Mpc. We chose this distance range to focus on the relatively unexplored volume outside the local Universe, where galaxies will be largely spatially unresolved and probe a range of X-ray luminosities that overlap with the low luminosity and/or highly obscured AGN population. We use the most recent X-ray scaling relations as well as the on-orbit eROSITA instrument performance to predict the X-ray emission from XRBs and diffuse hot gas and to perform both an analytic prediction and an end-to-end simulation using the mission simulation software, SIXTE. We consider potential contributions from hidden AGN and comment on the impact of normal galaxies on the measurement of the faint end of the AGN luminosity function. We predict that the eROSITA 4-year survey, will detect ≳15 000 galaxies (3σ significance) at 50–200 Mpc, which is ~ 100× more normal galaxies than detected in any X-ray survey to date.

Surveys↗

A simulation technique for predicting thickness of thermal sprayed coatings

The complexity of many of the components being coated today using the thermal spray process makes the trial and error approach traditionally followed in depositing a uniform coating inadequate, thereby necessitating a more analytical approach to developing robotic trajectories. A two dimensional finite difference simulation model has been developed to predict the thickness of coatings deposited using the thermal spray process. The model couples robotic and component trajectories and thermal spraying parameters to predict coating thickness. Simulations and experimental verification were performed on a rotating disk to evaluate the predictive capabilities of the approach.

Goedjen, John G.↗

An investigation of the information propagation and entropy transport aspects of Stirling machine numerical simulation

Aspects of the information propagation modeling behavior of integral machine computer simulation programs are investigated in terms of a transmission line. In particular, the effects of pressure-linking and temporal integration algorithms on the amplitude ratio and phase angle predictions are compared against experimental and closed-form analytic data. It is concluded that the discretized, first order conservation balances may not be adequate for modeling information propagation effects at characteristic numbers less than about 24. An entropy transport equation suitable for generalized use in Stirling machine simulation is developed. The equation is evaluated by including it in a simulation of an incompressible oscillating flow apparatus designed to demonstrate the effect of flow oscillations on the enhancement of thermal diffusion. Numerical false diffusion is found to be a major factor inhibiting validation of the simulation predictions with experimental and closed-form analytic data. A generalized false diffusion correction algorithm is developed which allows the numerical results to match their analytic counterparts. Under these conditions, the simulation yields entropy predictions which satisfy Clausius' inequality.

Goldberg, Louis F.↗

Analytical and experimental study of control effort associated with model reference adaptive control

Numerical simulation results presently obtained for the performance of model reference adaptive control (MRAC) are experimentally verified, with a view to accounting for differences between the plant and the reference model after the control function has been brought to bear. MRAC is both experimentally and analytically applied to a single-degree-of-freedom system, as well as analytically to a MIMO system having controlled differences between the reference model and the plant. The control effort is noted to be sensitive to differences between the plant and the reference model.

Messer, R. S.↗

Boundary value problems in magnetohydrodynamics (and fluid dynamics). I - Radiation boundary condition

In view of the difficulties that emerge in numerical simulations of time-dependent phenomena involving MHD waves due to waves' artificial reflection at numerical boundaries, where time-dependent magnetic fields outside those boundaries are unknown, a method previously developed for the determination of radiative boundary conditions is generalized for the treatment of isothermal MHD waves. The illustrative problem thus treated is that of the magnetic braking of a rigidly rotating disk which is embedded in a differentially-rotating fluid. The analytic results obtained are compared with those from numerical simulations with and without the present boundary treatment.

Vanajakshi, T. C.↗

Numerical investigations in three-dimensional internal flows

Renewed interest in simulation of high enthalpy flows was brought about by interest in the design and testing of hypersonic propulsion components. This renewed interest has led to refurbishment of high speed ground based wind tunnel facilities. Simulation of the high enthalpy flows occurring in the NASP propulsion system is of particular interest. Studying the combustor presentes some of the more challenging problems in the ground based simulation of the fluid flow. A difficulty in simulating the required conditions at the entrance to the combustor is related to achieveing a sufficiently high static pressure such that combustion would be possible at a high Mach number. An analysis of the flow within the original nozzle was not performed previously using numerical simulations, although it was operated experimentally. A Navier-Stokes code was chosen in the present study to examine analytically the flow within representative selected nozzles for use in combustion experiments. The numerical simulation of these fluid flows is popular and may be able to shed light on the direction that redesign of the shock tube nozzles should take. The numerical codes were modified to account for the axially symmetric flow field, for the inclusion of equilibrium real gas effects and for ideal gases whose thermodynamic properties are frozen at the stagnation conditions. Numerical simulation of flows within the original and two proposed nozzles, including the energy lost to the relatively cold walls, are investigated.

Rose, William C.↗

Double Power-law Formation by Sequential Particle Acceleration

Spectral double power laws are common in solar high-energy phenomena such as flares and interplanetary energetic-electron events. However, the physical mechanism that produces the changes in power-law index within a single spectrum is unclear. We developed a fully analytical method of forming single power-law spectra from sequential acceleration of particles orbiting inside and hopping between simulated large-scale magnetic islands formed by flare reconnection. Here, we extend the analytical method to the formation of double power-law spectra by assuming sequential acceleration in two successive regions with different acceleration and particle-transport rates. The resulting spectral distribution is continuous and smooth, with a flattening at low energies, two power-law regions at mid-energies, and a steep rollover at high energies. The model provides analytical expressions for the spectral indices, all energy breaks, and normalization constants as functions of just three physical parameters of each acceleration region: (1) the energy gain in each accelerator, (2) the percentage of particles transferred between accelerators, and (3) the number of accelerators visited. One of the most salient predictions of our work is that the spectral index at high (low) energies is determined by the parameters of the first "seed" (second) acceleration region. By constructing the spectral distribution through an iterative analytical process, the evolution toward a double power law is easily characterized and explained. Our analytical model provides tools to interpret space- and ground-based observations from RHESSI, FOXSI, NuSTAR, Solar Orbiter/STIX, EOVSA, and future high-energy missions.

Solar flares↗

Axisymmetric Plume Simulations with NASA's DSMC Analysis Code

A comparison of axisymmetric Direct Simulation Monte Carlo (DSMC) Analysis Code (DAC) results to analytic and Computational Fluid Dynamics (CFD) solutions in the near continuum regime and to 3D DAC solutions in the rarefied regime for expansion plumes into a vacuum is performed to investigate the validity of the newest DAC axisymmetric implementation. This new implementation, based on the standard DSMC axisymmetric approach where the representative molecules are allowed to move in all three dimensions but are rotated back to the plane of symmetry by the end of the move step, has been fully integrated into the 3D-based DAC code and therefore retains all of DAC s features, such as being able to compute flow over complex geometries and to model chemistry. Axisymmetric DAC results for a spherically symmetric isentropic expansion are in very good agreement with a source flow analytic solution in the continuum regime and show departure from equilibrium downstream of the estimated breakdown location. Axisymmetric density contours also compare favorably against CFD results for the R1E thruster while temperature contours depart from equilibrium very rapidly away from the estimated breakdown surface. Finally, axisymmetric and 3D DAC results are in very good agreement over the entire plume region and, as expected, this new axisymmetric implementation shows a significant reduction in computer resources required to achieve accurate simulations for this problem over the 3D simulations.

Stewart, B. D.↗