Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simplified Simulator”

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 397 records · Page 22

M51's spiral structure

The M51 system (NGC 5194/5195) provides an excellent problem both in spiral structure and in galaxy interactions. The authors present an analytic study of a computer experiment on the excitation mechanisms for M51's spiral arms and whether or not a halo is important for these mechanisms. This work extends previous numerical studies of the M51 system by including self-gravitation in a two component disk: gas and stars, and a dark halo. The analytic study provides two new observational constraints: the time (approx. 70 to 84 million years ago) and position angle of perigalacticon (300 degrees). By using these constraints and a simple conic approximation, the search for the companion's possible orbit is greatly simplified. This requires fewer N-body simulations than a fully self-gravitating orbit search.

Howard, S.↗

Analysis of fault-tolerant neurocontrol architectures

The fault-tolerance of analog parallel distributed implementations of a multivariable aircraft neurocontroller is analyzed by simulating weight and neuron failures in a simplified scheme of analog processing based on the functional architecture of the ETANN chip (Electrically Trainable Artificial Neural Network). The neural information processing is found to be only partially distributed throughout the set of weights of the neurocontroller synthesized with the backpropagation algorithm. Although the degree of distribution of the neural processing, and consequently the fault-tolerance of the neurocontroller, could be enhanced using Locally Distributed Weight and Neuron Approaches, a satisfactory level of fault-tolerance could only be obtained by retraining the degrated VLSI neurocontroller. The possibility of maintaining neurocontrol performance and stability in the presence of single weight of neuron failures was demonstrated through an automated retraining procedure of the neurocontroller based on a pre-programmed choice and sequence of the training parameters.

Troudet, T.↗

Laboratory and theoretical models of planetary-scale instabilities and waves

Research work is proceeding in theoretical, numerical, and experimental geophysical fluid dynamics leading up to a reflight of the GFFC (Geophysical Fluid Flow Cell Experiment) on USML-2. The work is intended not only to generate ideas for future space experiments, but to provide fundamental results concerned with nonlinear and chaotic properties of thermal convection and baroclinic waves in terrestrial and planetary atmospheres. The major efforts are focussed on thermal convection in a rapidly rotating annulus relevant to Jovian atmospheric dynamics, and on the chaotic behavior of baroclinic waves relevant to the Earth's atmosphere. The approach, in preparation for USML-2, is primarily theoretical and numerical. Mechanistic process models are solved numerically in order to identify physical mechanisms that may be observed in the GFFC, and which are important in real geophysical applications. The results from numerical simulations of geophysical fluid flow (subject to rotation and stratification) are compared with previous GFFC experiments on Spacelab-3 and with existing and proposed terrestrial laboratory experiments of various types. Pattern recognition algorithms have been employed to generate low-dimensional descriptions of the highly nonlinear and turbulent numerical simulations. Such empirically truncated descriptions provide for simplified but robust physical interpretations of the dynamics, as well as yielding highly efficient computations of these chaotic flows.

Hart, John E.↗

Design for inadvertent damage in composite laminates

Simplified predictive methods and models to computationally simulate durability and damage in polymer matrix composite materials/structures are described. The models include (1) progressive fracture, (2) progressively damaged structural behavior, (3) progressive fracture in aggressive environments, (4) stress concentrations, and (5) impact resistance. Several examples are included to illustrate applications of the models and to identify significant parameters and sensitivities. Comparisons with limited experimental data are made.

Singhal, Surendra N.↗

On the application of Chimera/unstructured hybrid grids for conjugate heat transfer

A hybrid grid system that combines the Chimera overset grid scheme and an unstructured grid method is developed to study fluid flow and heat transfer problems. With the proposed method, the solid structural region, in which only the heat conduction is considered, can be easily represented using an unstructured grid method. As for the fluid flow region external to the solid material, the Chimera overset grid scheme has been shown to be very flexible and efficient in resolving complex configurations. The numerical analyses require the flow field solution and material thermal response to be obtained simultaneously. A continuous transfer of temperature and heat flux is specified at the interface, which connects the solid structure and the fluid flow as an integral system. Numerical results are compared with analytical and experimental data for a flat plate and a C3X cooled turbine cascade. A simplified drum-disk system is also simulated to show the effectiveness of this hybrid grid system.

Kao, Kai-Hsiung↗

Numerical Modeling of Solidification in Space With MEPHISTO-4

A pre-flight analysis of the directional solidification of BiSn with MEPHISTO-4 is presented. Simplified Bridgman growth under microgravity conditions is simulated using a two dimensional finite element model. This numerical model is a single domain, pseudo-steady state model, and includes the effects of both thermal and solutal convection. The results show that for all orientations of the applied steady state gravity vector, of magnitude 1 micro-g, the directional solidification process remains diffusion controlled. The maximum convective velocity was found to be 4.424 x 10(exp -5) cm/s for the horizontal Bridgman growth configuration. This value is an order of magnitude lower than the growth velocity. The maximum and minimum values or solute concentration in the liquid at the crystal-melt interface were 13.867 at.% and 13.722 at.%, respectively. This gives a radial segregation value of xi = 1.046% at the interface. A secondary objective of this work was to compare the results obtained to those that consider thermal convection only (no solutal convection). It was found that the convective flow patterns in simulations which included solutal convection were significantly different from those which ignored solutal convection. The level of radial segregation predicted by the current simulations is an order of magnitude lower than that found in simulations which ignore solutal convection. The final aim was to investigate the effect of g-jitter on the crystal growth process. A simulation was performed to calculate the system response to a 1 second, 100 micro-g gravity impulse acting normal to the direction of growth. This pulse is consistent with that induced by Orbiter thruster firings. The results obtained indicate that such a gravity pulse causes an increase in the level of radial solute segregation at the interface from the steady state values. The maximum value of solute concentration in the liquid was found to be 13.888 at.%, the minimum value calculated was 13.706 at.%, yielding a radial segregation value of xi = 1.31% at the interface. These values occurred 126 seconds after the pulse terminated. Thus it is anticipated that the process will remain diffusion controlled even when subjected to such g-jitter.

Simpson, James E.↗

Flame-Vortex Studies to Quantify Markstein Numbers Needed to Model Flame Extinction Limits

This has quantified a database of Markstein numbers for unsteady flames; future work will quantify a database of flame extinction limits for unsteady conditions. Unsteady extinction limits have not been documented previously; both a stretch rate and a residence time must be measured, since extinction requires that the stretch rate be sufficiently large for a sufficiently long residence time. Ma was measured for an inwardly-propagating flame (IPF) that is negatively-stretched under microgravity conditions. Computations also were performed using RUN-1DL to explain the measurements. The Markstein number of an inwardly-propagating flame, for both the microgravity experiment and the computations, is significantly larger than that of an outwardy-propagating flame. The computed profiles of the various species within the flame suggest reasons. Computed hydrogen concentrations build up ahead of the IPF but not the OPF. Understanding was gained by running the computations for both simplified and full-chemistry conditions. Numerical Simulations. To explain the experimental findings, numerical simulations of both inwardly and outwardly propagating spherical flames (with complex chemistry) were generated using the RUN-1DL code, which includes 16 species and 46 reactions.

Driscoll, James F.↗

A GCM Recent History of the Northern Martian Polar Layered Deposits

The polar layered deposits are thought to contain alternate layers of water and dust in different proportions resulting from the astronomical forcing of the martian climate. In particular, longterm variations in the orbital and axial elements of Mars are presumed to generate variations of the latitudes of surface water ice stability and of the amount of water exchanged in the polar areas. At high obliquity, simplified climate models and independent general circulation simulations suggest a transfer of water ice from the north polar region to tropical areas, whereas at lower and present obliquities, water ice is expected to be stable only at the poles. If so, over obliquity cycles, water ice may be redistributed between the surface water reservoirs leading to their incremental building or disintegration depending on the rates of water transfer. If only a relative limited amount of the available water is exchanged on orbital timescales, this may provide an efficient mechanism for the formation of the observed polar deposits. Within this context, GCM simulations of the martian water cycle have been performed for various obliquities ranging from 15 degrees to 45 degrees and for a large set of initial water ice locations to determine the rate of water exchange between the surface water reservoirs as a function of the obliquity. Propagating these rates over the last 10 Ma orbital history gives a possible recent evolution of these reservoirs.

Levrard, B.↗

Orbit Stability of OSIRIS-REx in the Vicinity of Bennu Using a High-Fidelity Solar Radiation Model

Solar radiation pressure is one of the largest perturbing forces on the OSIRISRex trajectory as it orbits the asteroid Bennu. In this work, we investigate how forces due to solar radiation perturb the OSIRIS-REx trajectory in a high-fidelity model. The model accounts for Bennu's non-spherical gravity field, third-body gravity forces from the Sun and Jupiter, as well as solar radiation forces acting on a simplified spacecraft model. Such high-fidelity simulations indicate significant solar radiation pressure perturbations from the nominal orbit. Modifications to the initial design of the nominal orbit are found using a variation of parameters approach that reduce the perturbation in eccentricity by a factor of one-half.

Williams, Trevor W.↗

Computational Assessment of the Impact of Wave Count on Rotating Detonation Engine Performance

The impact of the number of azimuthally propagating waves on the performance of a Rotating Detonation Engine (RDE) is investigated using a simplified two-dimensional computational fluid dynamic simulation. The basic RDE configuration examined has no exit throat. The inlet is assumed lossless and does not allow backflow. The adiabatic, inviscid, and premixed simulation utilizes a particularly simple finite rate reaction mechanism that allows user control over the relative amounts of deflagration and detonation that occur, the presence and extent of a reaction delay associated with fuel and oxidizer mixing, and the number of waves present in the domain. The simplifications and idealizations decouple the effect of wave count from other potential loss mechanisms. Performance is measured using exhaust flow total pressure gain relative to the inlet total pressure. One, two, and three wave solutions are computed under identical boundary conditions and grid resolution. For simulations with no mixing delay and minimal deflagration, the number of waves present has negligible impact on performance. With a mixing delay, performance decreases with increasing wave number. With increased deflagration, performance increases with increasing wave number. With both effects simulated it is found that the two wave solution performs better than either the one or three wave solutions. The causes of these trends are explored. They imply that for practical, RDE’s, the number of waves that are present can impact performance.

Detonation↗

Computational Assessment of the Impact of Wave Count on Rotating Detonation Engine Performance

The impact of the number of azimuthally propagating waves on the performance of a Rotating Detonation Engine (RDE) is investigated using a simplified two-dimensional computational fluid dynamic simulation. The basic RDE configuration examined has no exit throat. The inlet is assumed lossless and does not allow backflow. The adiabatic, inviscid, and premixed simulation utilizes a particularly simple finite rate reaction mechanism that allows user control over the relative amounts of deflagration and detonation that occur, the presence and extent of a reaction delay associated with fuel and oxidizer mixing, and the number of waves present in the domain. The simplifications and idealizations decouple the effect of wave count from other potential loss mechanisms. Performance is measured using exhaust flow total pressure gain relative to the inlet total pressure. One, two, and three wave solutions are computed under identical boundary conditions and grid resolution. For simulations with no mixing delay and minimal deflagration, the number of waves present has negligible impact on performance. With a mixing delay, performance decreases with increasing wave number. With increased deflagration, performance increases with increasing wave number. With both effects simulated it is found that the two wave solution performs better than either the one or three wave solutions. The causes of these trends are explored. They imply that for practical, RDE’s, the number of waves that are present can impact performance.

Detonation↗

Simplified Applied Loads of a Transonic Truss-Braced Wing Derived from Computational Fluid Dynamics

The transonic truss-braced wing (TTBW) concept is under consideration by the NASA Advanced Air Transport Technology (AATT) project due to its potential to reduce fuel consumption of transport aircraft. To properly evaluate the TTBW concept, detailed structural modeling and subsequent test planning is underway, which requires a simplified set of discrete applied loads that simulate the distributed aerodynamic pressures. Given the unique architecture of the TTBW and its intended operation in the transonic regime, computational fluid dynamics (CFD) is the ideal approach to employ to ensure the aerodynamics of this configuration are sufficiently captured. To convert aerodynamic results into input for structural models, application of a data transfer tool, Discrete Data Transfer Between Dissimilar Meshes (DDTBDM), is presented to resolve aerodynamic pressures from CFD into single load vectors at hardpoints along both the wing and the strut. Evaluation of mapping the aerodynamic forces to hardpoints at rib and spar intersections in a structural model is performed, and the results are also compared with a linear static aeroelastic analysis. The structural responses indicate that distributing forces along only the wing and strut forward spars is insufficient for simulating the distributed aerodynamic loads. Adding moments along the forward spars improved the global response of the wing, but the point moments created unrealistic local rotations within the structure. When distributing loads along both the forward and aft spars of the wing and strut, the global response was less sensitive to the inclusion of point moments, which indicated the arrangement of point forces alone was sufficient input for structural design work and testing and provided a more continuous loading profile.

Finite Element↗

Simplified Applied Loads of a Transonic Truss-Braced Wing Derived from Computational Fluid Dynamic

The transonic truss-braced wing (TTBW) concept is under consideration by the NASA Advanced Air Transport Technology (AATT) project due to its potential to reduce fuel consumption of transport aircraft. To properly evaluate the TTBW concept, detailed structural modeling and subsequent test planning is underway, which requires a simplified set of discrete applied loads that simulate the distributed aerodynamic pressures. Given the unique architecture of the TTBW and its intended operation in the transonic regime, computational fluid dynamics (CFD) is the ideal approach to employ to ensure the aerodynamics of this configuration are sufficiently captured. To convert aerodynamic results into input for structural models, application of a data transfer tool, Discrete Data Transfer Between Dissimilar Meshes (DDTBDM), is presented to resolve aerodynamic pressures from CFD into single load vectors at hardpoints along both the wing and the strut. Evaluation of mapping the aerodynamic forces to hardpoints at rib and spar intersections in a structural model is performed, and the results are also compared with a linear static aeroelastic analysis. The structural responses indicate that distributing forces along only the wing and strut forward spars is insufficient for simulating the distributed aerodynamic loads. Adding moments along the forward spars improved the global response of the wing, but the point moments created unrealistic local rotations within the structure. When distributing loads along both the forward and aft spars of the wing and strut, the global response was less sensitive to the inclusion of point moments, which indicated the arrangement of point forces alone was sufficient input for structural design work and testing and provided a more continuous loading profile.

Finite Element↗

Measuring and Utilizing Gravity-Gradient Induced Torques on Future Gravity Recovery Missions

This research is a novel investigation into the use of newly-available relative angular acceleration measurements between a spacecraft utilized for gravity recovery missions and an internally located test mass. The gravity-gradient torque equation for zonal spherical harmonic order n is formulated for a known gravitational potential field, and through simulations it is proven that the Simplified-Gravitational Reference Sensor will be sensitive to the gravity-gradient induced torques acting on its test mass. The full conference paper will demonstrate how the presented gravity-gradient torque equation and the measured relative angular acceleration between the spacecraft and test mass will improve the accuracy of the gravity field models acquired by future gravity recovery missions by directly measuring the drag acting on the spacecraft with a single accelerometer.

gravity recovery↗

Numerical methods for the simulation of complex multi-body flows with applications for the integrated Space Shuttle vehicle

This project forms part of the long term computational effort to simulate the time dependent flow over the integrated Space Shuttle vehicle (orbiter, solid rocket boosters (SRB's), external tank (ET), and attach hardware) during its ascent mode for various nominal and abort flight conditions. Due to the limitations of experimental data such as wind tunnel wall effects and the difficulty of safely obtaining valid flight data, numerical simulations are undertaken to supplement the existing data base. This data can then be used to predict the aerodynamic behavior over a wide range of flight conditions. Existing computational results show relatively good overall comparison with experiments but further refinement is required to reduce numerical errors and to obtain finer agreements over a larger parameter space. One of the important goals of this project is to obtain better comparisons between numerical simulations and experiments. In the simulations performed so far, the geometry has been simplified in various ways to reduce the complexity so that useful results can be obtained in a reasonable time frame due to limitations in computer resources. In this project, the finer details of the major components of the Space Shuttle are modeled better by including more complexity in the geometry definition. Smaller components not included in early Space Shuttle simulations will now be modeled and gridded.

Chan, William M.↗

Random Dopant Induced Threshold Voltage Lowering and Fluctuations in Sub 50 nm MOSFETs: a Statistical 3D 'Atomistic' Simulation Study

A 3D 'atomistic' simulation study of random dopant induced threshold voltage fluctuations and lowering in sub 50 nm MOSFETs is presented. The attention is focused mainly on devices with 30 nm effective channel length which represent the expected level of scaling at the end of the Silicon Roadmap. An efficient algorithm, based on a single 3D ap solution of the Poisson equation and a simplified current continuity equation, is used in the simulations. Large samples of microscopically different devices (typically 200) arc used in order to obtain statistically reliable results. The influence of different aspects of the conventional MOSFET design on the threshold voltage fluctuations and lowering are investigated. Results for fluctuation resistant device architectures based on low-doped epitaxial channel MOSFETs are also presented.

Asenov, Asen↗

High-Fidelity Computational Aerodynamics of Multi-Rotor Unmanned Aerial Vehicles

High-fidelity Computational Fluid Dynamics (CFD) simulations have been carried out for several multi-rotor Unmanned Aerial Vehicles (UAVs). Three vehicles have been studied: the classic quadcopter DJI Phantom 3, an unconventional quadcopter specialized for forward flight, the SUI Endurance, and an innovative concept for Urban Air Mobility (UAM), the Elytron 4S UAV. The three-dimensional unsteady Navier-Stokes equations are solved on overset grids using high-order accurate schemes, dual-time stepping, and a hybrid turbulence model. The DJI Phantom 3 is simulated with different rotors and with both a simplified airframe and the real airframe including landing gear and a camera. The effects of weather are studied for the DJI Phantom 3 quadcopter in hover. The SUI En- durance original design is compared in forward flight to a new configuration conceived by the authors, the hybrid configuration, which gives a large improvement in forward thrust. The Elytron 4S UAV is simulated in helicopter mode and in airplane mode. Understanding the complex flows in multi-rotor vehicles will help design quieter, safer, and more efficient future drones and UAM vehicles.

High-Fidelity↗

Simplified Modeling of Oxidation of Hydrocarbons

A method of simplified computational modeling of oxidation of hydrocarbons is undergoing development. This is one of several developments needed to enable accurate computational simulation of turbulent, chemically reacting flows. At present, accurate computational simulation of such flows is difficult or impossible in most cases because (1) the numbers of grid points needed for adequate spatial resolution of turbulent flows in realistically complex geometries are beyond the capabilities of typical supercomputers now in use and (2) the combustion of typical hydrocarbons proceeds through decomposition into hundreds of molecular species interacting through thousands of reactions. Hence, the combination of detailed reaction- rate models with the fundamental flow equations yields flow models that are computationally prohibitive. Hence, further, a reduction of at least an order of magnitude in the dimension of reaction kinetics is one of the prerequisites for feasibility of computational simulation of turbulent, chemically reacting flows. In the present method of simplified modeling, all molecular species involved in the oxidation of hydrocarbons are classified as either light or heavy; heavy molecules are those having 3 or more carbon atoms. The light molecules are not subject to meaningful decomposition, and the heavy molecules are considered to decompose into only 13 specified constituent radicals, a few of which are listed in the table. One constructs a reduced-order model, suitable for use in estimating the release of heat and the evolution of temperature in combustion, from a base comprising the 13 constituent radicals plus a total of 26 other species that include the light molecules and related light free radicals. Then rather than following all possible species through their reaction coordinates, one follows only the reduced set of reaction coordinates of the base. The behavior of the base was examined in test computational simulations of the combustion of heptane in a stirred reactor at various initial pressures ranging from 0.1 to 6 MPa. Most of the simulations were performed for stoichiometric mixtures; some were performed for fuel/oxygen mole ratios of 1/2 and 2.

Bellan, Josette↗