Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “simulation methods and programs”

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,117 records · Page 62

Perform a gyro test of general relativity in a satellite and develop associated control technology

The progress accomplished in the Stanford Gyro Relativity program during the period November 1974 to October 1975 was described. Gyro developments were continued in the main laboratory dewar, concentrating on the operation of a three axis gyro readout and on improvements to the methods of canceling trapped fields in the rotor; these efforts culminated in the first successful observation of the London moment in the spinning gyro rotor in March 1975. Following a review meeting at that time, a new goal was formulated for the next 12 to 18 months, namely to operate a gyroscope in the new ultra-low field facility with readout resolution approaching 1 arc-second. The following other tasks were also completed: (1) sputtering work, (2) magnetometry, (3) construction and installation of the North Star simulator, (4) analysis of torques on the gyro, especially in inclined orbits, (5) equivalence principle accelerometer, and (6) analysis of a twin-satellite test of relativity.

Fairbank, W. M.↗

Automated problem scheduling and reduction of synchronization delay effects

It is anticipated that in order to make effective use of many future high performance architectures, programs will have to exhibit at least a medium grained parallelism. A framework is presented for partitioning very sparse triangular systems of linear equations that is designed to produce favorable preformance results in a wide variety of parallel architectures. Efficient methods for solving these systems are of interest because: (1) they provide a useful model problem for use in exploring heuristics for the aggregation, mapping and scheduling of relatively fine grained computations whose data dependencies are specified by directed acrylic graphs, and (2) because such efficient methods can find direct application in the development of parallel algorithms for scientific computation. Simple expressions are derived that describe how to schedule computational work with varying degrees of granularity. The Encore Multimax was used as a hardware simulator to investigate the performance effects of using the partitioning techniques presented in shared memory architectures with varying relative synchronization costs.

Saltz, Joel H.↗

Narrow Scale Flow and a Weak Field by the Top of Earth's Core: Evidence from Orsted, Magsat and Secular Variation

As Earth's main magnetic field weakens, our magnetic shield against the onslaught of the solar wind thins. And the field strength needed to fend off battering by solar coronal mass ejections is decreasing, just when the delicate complexity of modem, vulnerable, electro-technological systems is increasing at an unprecedented rate. Recently, a working group of distinguished scientist from across the nation has asked NASA's Solid Earth and Natural Hazards program a key question: What are the dynamics of Earth s magnetic field and its interactions with the Earth system? Paleomagnetic studies of crustal rocks magnetized in the geologic past reveal that polarity reversals have occurred many times during Earth s history. Networked super-computer simulations of core field and flow, including effects of gravitational, pressure, rotational Coriolis, magnetic and viscous forces, suggest how this might happen in detail. And space-based measurements of the real, time-varying magnetic field help constrain estimates of the speed and direction of fluid iron flowing near the top of the core and enable tests of some hypotheses about such flow. Now scientists at NASA s Goddard Space Flight Center have developed and applied methods to test the hypotheses of narrow scale flow and of a dynamically weak magnetic field near the top of Earth s core. Using two completely different methods, C. V. Voorhies has shown these hypotheses lead to specific theoretical forms for the "spectrum" of Earth s main magnetic field and the spectrum of its rate of change. Much as solar physicists use a prism to separate sunlight into its spectrum, from long wavelength red to short wavelength blue light, geophysicists use a digital prism, spherical harmonic analysis, to separate the measured geomagnetic field into its spectrum, from long to short wavelength fields. They do this for the rate of change of the field as well.

Voorhies, Coerte V.↗

Squash-Box Feasibility Driven Differential Dynamic Programming

Recently, Differential Dynamic Programming (DDP) and other similar algorithms have become the solvers of choice when performing non-linear Model Predictive Control (nMPC) with modern robotic devices. The reason is that they have a lower computational cost per iteration when compared with off-the-shelf Non-Linear Programming (NLP) solvers, which enables its online operation. However, they cannot handle constraints, and are known to have poor convergence capabilities. In this paper, we propose a method to solve the optimal control problem with control bounds through a squashing function (i.e., a sigmoid, which is bounded by construction). It has been shown that a naive use of squashing functions damage the convergence rate. To tackle this, we first propose to add a quadratic barrier that avoids the difficulty of the plateau produced by the sigmoid. Second, we add an outer loop that adapts both the sigmoid and the barrier; it makes the optimal control problem with the squashing function converge to the original control-bounded problem. To validate our method, we present simulation results for different types of platforms including a multi-rotor, a biped, a quadruped and a humanoid robot.

Navarro, Angel Santamaria↗

Hybrid data-driven and model-informed online tool wear detection in milling machines

Precision machining tool wear is responsible for low product throughput and quality. Monitoring the tool wear online is vital to prevent degradation in machining quality. However, direct real-time tool wear measurement is not practical. This paper presents residual-based anomaly detection models, combining a hybrid model comprised of a physics-based model and a data-driven model (a decision tree or a neural network) to predict signals of interest (e.g., power or forces) under nominal conditions, followed by Page’s cumulative sum test for detecting tool wear on-line using the computer numerical control machine measurements. The most informative features are ranked using dynamic programming and its approximation variants from real-time measurements and machine settings, such as the width of cut, depth of cut, feed rate and spindle speed, that serve as inputs to the predictive models. The baseline nominal model is incrementally updated with experimental data via a gradient boosted adaptation model to generate the residuals that account for discrepancies between the actual machine data under normal conditions and the baseline nominal model predictions. The hybrid model is validated against 20 Mazak milling machine experimental tests and one Haas run-to-failure experiment. The proposed anomaly detector is applied to synthetic data from simulations of the physics-based model at different operating conditions, measurement noise levels, and tool wear levels, and the methods were able to achieve an overall 92% accuracy in data with 1% noise. The anomaly detection methods based on hybrid model reduced the false alarms of either the data-driven or physical-based models alone, and are found to be capable of good online detection of tool wear.

Online anomaly detection↗

Two-Dimensional Network Flow Modeling of No-Vent Tank Filling of A Cryogenic Tank With Thermo-Dynamic Vent System Assisted Injector

A finite volume-based network flow modeling tool, Generalized Fluid System Simulation Program (GFSSP), was used to simulate the no-vent filling of a cryogenic tank with a thermodynamic vent system assisted injector. The tank was discretized into nodes in the radial and axial directions to calculate two-dimensional axisymmetric flow with heat transfer between the solid wall and fluid. The nodes are connected by branches where momentum equations are solved to calculate the flowrates. The mass and energy conservation equations and the equation of state are solved to calculate pressure, temperature, and resident mass at the nodes. The energy conservation equations at the solid nodes were solved to calculate the temperature of the solid wall. The system of equations was solved by a combination of successive substitution and simultaneous Newton-Raphson method. The numerical model accounts for a) different regimes of pool boiling heat transfer, b) condensation of vapor around the spray droplets, c) condensation of vapor at the interface of the cooled injector and ullage, and d) condensation at the liquid-vapor interphase. The numerical predictions of tank pressure, filling rate, and wall temperatures were compared with the test data. The paper also presents the effect of nodal discretization by comparing the results of the 1D and 2D models and examines the effect of the droplet diameter on tank filling.

Nodal Model↗

Two-Dimensional Network Flow Modeling of No-Vent Tank Filling of A Cryogenic Tank With Thermo-Dynamic Vent System Assisted Injector

A finite volume-based network flow modeling tool, Generalized Fluid System Simulation Program (GFSSP), was used to simulate the no-vent filling of a cryogenic tank with a thermodynamic vent system assisted injector. The tank was discretized into nodes in the radial and axial directions to calculate two-dimensional axisymmetric flow with heat transfer between the solid wall and fluid. The nodes are connected by branches where momentum equations are solved to calculate the flowrates. The mass and energy conservation equations and the equation of state are solved to calculate pressure, temperature, and resident mass at the nodes. The energy conservation equations at the solid nodes were solved to calculate the temperature of the solid wall. The system of equations was solved by a combination of successive substitution and simultaneous Newton-Raphson method. The numerical model accounts for a) different regimes of pool boiling heat transfer, b) condensation of vapor around the spray droplets, c) condensation of vapor at the interface of the cooled injector and ullage, and d) condensation at the liquid-vapor interphase. The numerical predictions of tank pressure, filling rate, and wall temperatures were compared with the test data. The paper also presents the effect of nodal discretization by comparing the results of the 1D and 2D models and examines the effect of the droplet diameter on tank filling.

Nodal Model↗

Perturbation Modeling and Navigation Tuning for a Crewed Station in Near Rectilinear Halo Orbit

NASA's Gateway program will build a crew-tended station in an Earth-Moon Near Rectilinear Halo Orbit (NRHO) to support deep space missions to the lunar surface and heliocentric space. The station in the NRHO will be visited by crewed vehicles including lunar landers that are high mass relative to Gateway. Docking events with large vehicles induce significant thruster plume and docking contact velocity perturbations. Crewed vehicles induce significant venting perturbations that can stress the navigation filter's ability to remain converged about an estimated state. The torques that arise from venting and lunar gravity gradient require reaction control system (RCS) thruster intervention, which themselves induce a velocity perturbation. The best method to model RCS thruster perturbations is currently in development. This paper simulates different RCS perturbation models, schedules, and an increasing magnitude to investigate the navigation sensitivity to different RCS perturbation models. From the RCS perturbation models simulated, a worst-case scenario is chosen for further scrutiny. The navigation filter is further tuned to handle the most disruptive perturbation modeling and the resulting performance is analyzed.

Navigation↗

A multi-purpose method for analysis of spur gear tooth loading

A large digitized approach was developed for the static and dynamic load analysis of spur gearing. An iterative procedure was used to calculate directly the "variable-variable" gear mesh stiffness as a function of transmitted load, gear tooth profile errors, gear tooth deflections and gear hub torsional deformation, and position of contacting profile points. The developed approach can be used to analyze the loads, Hertz stresses, and PV for the normal and high contrast ratio gearing, presently the modeling is limited to the condition that for a given gear all teeth have identical spacing and profiles (with or without surface imperfections). Certain types of simulated sinusoidal profile errors and pitting can cause interruptions of the gear mesh stiffness function and, thus, increase the dynamic loads in spur gearing. In addition, a finite element stress and mesh subprogram was developed for future introduction into the main program for calculating the gear tooth bending stresses under dynamic loads.

Kasuba, R.↗

NASA/FAA general aviation crash dynamics program

The program involves controlled full scale crash testing, nonlinear structural analyses to predict large deflection elastoplastic response, and load attenuating concepts for use in improved seat and subfloor structure. Both analytical and experimental methods are used to develop expertise in these areas. Analyses include simplified procedures for estimating energy dissipating capabilities and comprehensive computerized procedures for predicting airframe response. These analyses are developed to provide designers with methods for predicting accelerations, loads, and displacements on collapsing structure. Tests on typical full scale aircraft and on full and subscale structural components are performed to verify the analyses and to demonstrate load attenuating concepts. A special apparatus was built to test emergency locator transmitters when attached to representative aircraft structure. The apparatus is shown to provide a good simulation of the longitudinal crash pulse observed in full scale aircraft crash tests.

Thomson, R. G.↗

Geological terrain models

The initial phase of a program to determine the best interpretation strategy and sensor configuration for a radar remote sensing system for geologic applications is discussed. In this phase, terrain modeling and radar image simulation were used to perform parametric sensitivity studies. A relatively simple computer-generated terrain model is presented, and the data base, backscatter file, and transfer function for digital image simulation are described. Sets of images are presented that simulate the results obtained with an X-band radar from an altitude of 800 km and at three different terrain-illumination angles. The simulations include power maps, slant-range images, ground-range images, and ground-range images with statistical noise incorporated. It is concluded that digital image simulation and computer modeling provide cost-effective methods for evaluating terrain variations and sensor parameter changes, for predicting results, and for defining optimum sensor parameters.

Kaupp, V. H.↗

Heat Transfer and Thermal Stability Research for Advanced Hydrocarbon Fuel Technologies

In recent years there has been increased interest in the development of a new generation of high performance boost rocket engines. These efforts, which will represent a substantial advancement in boost engine technology over that developed for the Space Shuttle Main Engines in the early 1970s, are being pursued both at NASA and the United States Air Force. NASA, under its Space Launch Initiative s Next Generation Launch Technology Program, is investigating the feasibility of developing a highly reliable, long-life, liquid oxygen/kerosene (RP-1) rocket engine for launch vehicles. One of the top technical risks to any engine program employing hydrocarbon fuels is the potential for fuel thermal stability and material compatibility problems to occur under the high-pressure, high-temperature conditions required for regenerative fuel cooling of the engine combustion chamber and nozzle. Decreased heat transfer due to carbon deposits forming on wetted fuel components, corrosion of materials common in engine construction (copper based alloys), and corrosion induced pressure drop increases have all been observed in laboratory tests simulating rocket engine cooling channels. To mitigate these risks, the knowledge of how these fuels behave in high temperature environments must be obtained. Currently, due to the complexity of the physical and chemical process occurring, the only way to accomplish this is empirically. Heated tube testing is a well-established method of experimentally determining the thermal stability and heat transfer characteristics of hydrocarbon fuels. The popularity of this method stems from the low cost incurred in testing when compared to hot fire engine tests, the ability to have greater control over experimental conditions, and the accessibility of the test section, facilitating easy instrumentation. These benefits make heated tube testing the best alternative to hot fire engine testing for thermal stability and heat transfer research. This investigation used the Heated Tube Facility at the NASA Glenn Research Center to perform a thermal stability and heat transfer characterization of RP-1 in an environment simulating that of a high chamber pressure, regenerative cooled rocket engine. The first step in the research was to investigate the carbon deposition process of previous heated tube experiments by performing scanning electron microscopic analysis in conjunction with energy dispersive spectroscopy on the tube sections. This analysis gave insight into the carbon deposition process and the effect that test conditions played in the formation of deleterious coke. Furthermore, several different formations were observed and noted. One other crucial finding of this investigation was that in sulfur containing hydrocarbon fuels, the interaction of the sulfur components with copper based wall materials presented a significant corrosion problem. This problem in many cases was more life limiting than those posed by the carbon deposition process. The results of this microscopic analysis was detailed and presented at the December 2003 JANNAF Air-Breathing Propulsion Meeting as a Materials Compatibility and Thermal Stability Analysis of common Hydrocarbon Fuels (reference 1).

DeWitt, Kenneth↗

A Monte Carlo/response surface strategy for sensitivity analysis: application to a dynamic model of vegetative plant growth

We describe the application of a strategy for conducting a sensitivity analysis for a complex dynamic model. The procedure involves preliminary screening of parameter sensitivities by numerical estimation of linear sensitivity coefficients, followed by generation of a response surface based on Monte Carlo simulation. Application is to a physiological model of the vegetative growth of soybean plants. The analysis provides insights as to the relative importance of certain physiological processes in controlling plant growth. Advantages and disadvantages of the strategy are discussed.

Non-NASA Center↗

Finding curvilinear features in speckled images

A method for finding curves in digital images with speckle noise is described. The solution method differs from standard linear convolutions followed by thresholds in that it explicitly allows curvature in the features. Maximum a posteriori (MAP) estimation is used, together with statistical models for the speckle noise and for the curve-generation process, to find the most probable estimate of the feature, given the image data. The estimation process is first described in general terms. Then, incorporation of the specific neighborhood system and a multiplicative noise model for speckle allows derivation of the solution, using dynamic programming, of the estimation problem. The detection of curvilinear features is considered separately. The detection results allow the determination of the minimal size of detectable feature. Finally, the estimation of linear features, followed by a detection step, is shown for computer-simulated images and for a SAR image of sea ice.

Samadani, Ramin↗

Plan recognition and generalization in command languages with application to telerobotics

A method for pragmatic inference as a necessary accompaniment to command languages is proposed. The approach taken focuses on the modeling and recognition of the human operator's intent, which relates sequences of domain actions ('plans') to changes in some model of the task environment. The salient feature of this module is that it captures some of the physical and linguistic contextual aspects of an instruction. This provides a basis for generalization and reinterpretation of the instruction in different task environments. The theoretical development is founded on previous work in computational linguistics and some recent models in the theory of action and intention. To illustrate these ideas, an experimental command language to a telerobot is implemented. The program consists of three different components: a robot graphic simulation, the command language itself, and the domain-independent pragmatic inference module. Examples of task instruction processes are provided to demonstrate the benefits of this approach.

Yared, Wael I.↗

Generalized Fluid System Simulation Program (GFSSP) - Version 6

The Generalized Fluid System Simulation Program (GFSSP) is a finite-volume based general-purpose computer program for analyzing steady state and time-dependent flow rates, pressures, temperatures, and concentrations in a complex flow network. The program is capable of modeling real fluids with phase changes, compressibility, mixture thermodynamics, conjugate heat transfer between solid and fluid, fluid transients, pumps, compressors, flow control valves and external body forces such as gravity and centrifugal. The thermo-fluid system to be analyzed is discretized into nodes, branches, and conductors. The scalar properties such as pressure, temperature, and concentrations are calculated at nodes. Mass flow rates and heat transfer rates are computed in branches and conductors. The graphical user interface allows users to build their models using the 'point, drag, and click' method; the users can also run their models and post-process the results in the same environment. The integrated fluid library supplies thermodynamic and thermo-physical properties of 36 fluids, and 24 different resistance/source options are provided for modeling momentum sources or sinks in the branches. Users can introduce new physics, non-linear and time-dependent boundary conditions through user-subroutine.

Majumdar, Alok↗

Experimental and analytical investigations of wave enhanced supersonic combustors

Supersonic combustion ramjet (scramjet) engines rely on rapid mixing and combustion of fuel for good performance. However, both of these processes are relatively slow compared to the residence times in supersonic combusters. Methods of improving the mixing process include subjecting the fuel streams to shock waves and the generation of vortices by spiriling, bifurcating struts. The combustion process can also be enhanced by shock or detonation waves. An oblique shock wave can act as a flameholder by increasing the pressure and temperature of the air-fuel mixture and thereby decreasing the ignition delay. If the oblique shock is sufficiently strong, then the combustion front and the shock wave can couple into a detonation wave. In this case, combustion occurs almost instantaneously in a thin zone behind the wave front. While the existence of standing oblique detonation waves has been proven in computer simulations, the experimental validation is still in progress. Currently, there is an experimental and analytical program at NASA-Ames Research Center to study detonations and other means of wave enhanced supersonic mixing and combustion.

Adelman, Henry G.↗

Improving Search Properties in Genetic Programming

With the advancing computer processing capabilities, practical computer applications are mostly limited by the amount of human programming required to accomplish a specific task. This necessary human participation creates many problems, such as dramatically increased cost. To alleviate the problem, computers must become more autonomous. In other words, computers must be capable to program/reprogram themselves to adapt to changing environments/tasks/demands/domains. Evolutionary computation offers potential means, but it must be advanced beyond its current practical limitations. Evolutionary algorithms model nature. They maintain a population of structures representing potential solutions to the problem at hand. These structures undergo a simulated evolution by means of mutation, crossover, and a Darwinian selective pressure. Genetic programming (GP) is the most promising example of an evolutionary algorithm. In GP, the structures that evolve are trees, which is a dramatic departure from previously used representations such as strings in genetic algorithms. The space of potential trees is defined by means of their elements: functions, which label internal nodes, and terminals, which label leaves. By attaching semantic interpretation to those elements, trees can be interpreted as computer programs (given an interpreter), evolved architectures, etc. JSC has begun exploring GP as a potential tool for its long-term project on evolving dextrous robotic capabilities. Last year we identified representation redundancies as the primary source of inefficiency in GP. Subsequently, we proposed a method to use problem constraints to reduce those redundancies, effectively reducing GP complexity. This method was implemented afterwards at the University of Missouri. This summer, we have evaluated the payoff from using problem constraints to reduce search complexity on two classes of problems: learning boolean functions and solving the forward kinematics problem. We have also developed and implemented methods to use additional problem heuristics to fine-tune the searchable space, and to use typing information to further reduce the search space. Additional improvements have been proposed, but they are yet to be explored and implemented.

Janikow, Cezary Z.↗