Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Simulation Results”

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 127 records · Page 7

Towards Full Aircraft Airframe Noise Prediction: Detached Eddy Simulations

Results from a computational study on the aeroacoustic characteristics of an 18%-scale, semi-span Gulf-stream aircraft model are presented in this paper. NASA's FUN3D unstructured compressible Navier-Stokes solver was used to perform steady and unsteady simulations of the flow field associated with this high-fidelity aircraft model. Solutions were obtained for free-air at a Mach number of 0.2 with the flap deflected at 39 deg, with the main gear off and on (the two baseline configurations). Initially, the study focused on accurately predicting the prominent noise sources at both flap tips for the baseline configuration with deployed flap only. Building upon the experience gained from this initial effort, subsequent work involved the full landing configuration with both flap and main landing gear deployed. For the unsteady computations, we capitalized on the Detached Eddy Simulation capability of FUN3D to capture the complex time-dependent flow features associated with the flap and main gear. To resolve the noise sources over a broad frequency range, the tailored grid was very dense near the flap inboard and outboard tips and the region surrounding the gear. Extensive comparison of the computed steady and unsteady surface pressures with wind tunnel measurements showed good agreement for the global aerodynamic characteristics and the local flow field at the flap inboard tip. However, the computed pressure coefficients indicated that a zone of separated flow that forms in the vicinity of the outboard tip is larger in extent along the flap span and chord than measurements suggest. Computed farfield acoustic characteristics from a FW-H integral approach that used the simulated pressures on the model solid surface were in excellent agreement with corresponding measurements.

Khorrami, Mehdi R.↗

Direct Simulation of Evolution and Control of Three-Dimensional Instabilities in Attachment-Line Boundary Layers

The spatial evolution of three-dimensional disturbances in an attachment-line boundary layer is computed by direct numerical simulation of the unsteady, incompressible Navier-Stokes equations. Disturbances are introduced into the boundary layer by harmonic sources that involve unsteady suction and blowing through the wall. Various harmonic- source generators are implemented on or near the attachment line, and the disturbance evolutions are compared. Previous two-dimensional simulation results and nonparallel theory are compared with the present results. The three-dimensional simulation results for disturbances with quasi-two-dimensional features indicate growth rates of only a few percent larger than pure two-dimensional results; however, the results are close enough to enable the use of the more computationally efficient, two-dimensional approach. However, true three-dimensional disturbances are more likely in practice and are more stable than two-dimensional disturbances. Disturbances generated off (but near) the attachment line spread both away from and toward the attachment line as they evolve. The evolution pattern is comparable to wave packets in at-plate boundary-layer flows. Suction stabilizes the quasi-two-dimensional attachment-line instabilities, and blowing destabilizes these instabilities; these results qualitatively agree with the theory. Furthermore, suction stabilizes the disturbances that develop off the attachment line. Clearly, disturbances that are generated near the attachment line can supply energy to attachment-line instabilities, but suction can be used to stabilize these instabilities.

Joslin, Ronald D.↗

Predictive Simulation of Gas Adsorption in Fixed-Beds and Limitations due to the Ill-Posed Danckwerts Boundary Condition

The 1-D axially dispersed plug flow model is a mathematical model widely used for the simulation of adsorption processes. Lumped mass transfer coefficients such as the Glueckauf linear driving force (LDF) term and the axial dispersion coefficient are generally obtained by fitting simulation results to the experimental breakthrough test data. An approach is introduced where these parameters, along with the only free parameter in the energy balance equations, are individually fit to specific test data that isolates the appropriate physics. It is shown that with this approach this model provides excellent simulation results for the C02 on zeolite SA sorbent/sorbate system; however, for the H20 on zeolite SA system, non-physical deviations from constant pattern behavior occur when fitting dispersive experimental results with a large axial dispersion coefficient. A method has also been developed that determines a priori what values of the LDF and axial dispersion terms will result in non-physical simulation results for a specific sorbent/sorbate system when using the one-dimensional axially dispersed plug flow model. A relationship between the steepness of the adsorption equilibrium isotherm as indicated by the distribution factor, the magnitude of the axial dispersion and mass transfer coefficient, and the resulting non-physical behavior is derived. This relationship is intended to provide a guide for avoiding non-physical behavior by limiting the magnitude of the axial dispersion term on the basis of the mass transfer coefficient and distribution factor.

Knox, James Clinton↗

Airborne measurements of solar and planetary near ultraviolet radiation during the NASA/ESA CV-900 spacelab simulation

Results from a comparative study of the feasibility of employing experiment operators on the space shuttle to acquire scientifically worthwhile data are presented. The experiments performed during these tests included spectral observations of the Sun and Venus in the near ultraviolet region. The solar measurements were analyzed to determine ozone abundance in the terrestrial atmosphere. Using a detailed spectral matching technique to compare airborne solar UV measurements with synthetic spectral profiles of sunlight, it is deduced that in winter the total atmospheric ozone abundance is about 0.33 atm/cm at midlatitudes in the northern hemisphere.

Sivjee, G. G.↗

Automated Unstructured Grid Adaptation on a Strut Fuel Injector at Hypervelocity Flow Conditions

Computational fluid dynamics (CFD) analysis is presented with the use of an automated unstructured grid adaptation tool on a strut fuel injector at hypervelocity flow conditions. The analysis was carried out with the VULCAN-CFD solver using Reynolds-averaged simulations (RAS). The hypervelocity flow conditions match the high Mach number flow of the experiments conducted as part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). The current work uses an automated grid adaptation tool recently implemented in VULCAN-CFD, and explores this tool’s ability to solve highspeed mixing problems. Simulation results obtained using the unstructured adaptive grid approach are compared to those on a user-generated structured grid. These results are evaluated by analyzing how efficiently comparable fidelity results are obtained from both adapted and structured simulations. In addition, two adaptation strategies were used to explore the impact on the final solution. In the current work, the unstructured grid adaptation tool automatically generates unstructured grids and performs adaptation of the grid based on a Hessian error estimate of a specified flowfield parameter. Multiple adaptations were executed using each run strategy with the one-dimensional values of the mixing efficiency used to determine grid convergence and for comparison with the structured grid simulation results. It was found that the unstructured adaptive grid simulations were able to produce results that matched closely with those on structured grids using far fewer grid cells, and thus, requiring far less computational time to reach the solution. It was also discovered that the adaptation run strategy influenced the total number of grid cells and the efficiency with which a final grid-adapted solution was reached. Furthermore, motivated by the grid convergence index (GCI) used for structured grid simulations, a grid convergence estimate (GCE) was developed and demonstrated for the grid adaptation. Overall, the investigation demonstrated that the automated unstructured grid adaptation tool implemented in VULCAN-CFD is capable of accurately and efficiently solving complex high-speed mixing problems with only a fraction of the grid cells required to obtain comparable results on a user-generated structured grid.

hypersonics↗

Automated Unstructured Grid Adaptation on a Strut Fuel Injector at Hypervelocity Flow Conditions

Computational fluid dynamics (CFD) analysis is presented with the use of an automated unstructured grid adaptation tool on a strut fuel injector at hypervelocity flow conditions. The analysis was carried out with the VULCAN-CFD solver using Reynolds-averaged simulations (RAS). The hypervelocity flow conditions match the high Mach number flow of the experiments conducted as part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). The current work uses an automated grid adaptation tool recently implemented in VULCAN-CFD, and explores this tool’s ability to solve highspeed mixing problems. Simulation results obtained using the unstructured adaptive grid approach are compared to those on a user-generated structured grid. These results are evaluated by analyzing how efficiently comparable fidelity results are obtained from both adapted and structured simulations. In addition, two adaptation strategies were used to explore the impact on the final solution. In the current work, the unstructured grid adaptation tool automatically generates unstructured grids and performs adaptation of the grid based on a Hessian error estimate of a specified flowfield parameter. Multiple adaptations were executed using each run strategy with the one-dimensional values of the mixing efficiency used to determine grid convergence and for comparison with the structured grid simulation results. It was found that the unstructured adaptive grid simulations were able to produce results that matched closely with those on structured grids using far fewer grid cells, and thus, requiring far less computational time to reach the solution. It was also discovered that the adaptation run strategy influenced the total number of grid cells and the efficiency with which a final grid-adapted solution was reached. Furthermore, motivated by the grid convergence index (GCI) used for structured grid simulations, a grid convergence estimate (GCE) was developed and demonstrated for the grid adaptation. Overall, the investigation demonstrated that the automated unstructured grid adaptation tool implemented in VULCAN-CFD is capable of accurately and efficiently solving complex high-speed mixing problems with only a fraction of the grid cells required to obtain comparable results on a user-generated structured grid.

hypersonics↗

Coupling of Charged Particles Via Coulombic Interactions: Numerical Simulations and Resultant Kappa-Like Velocity Space Distribution Functions

A parametric study is performed using the electrostatic simulations of Randol and Christian (2014) in which the number density, n, and initial thermal speed, theta, are varied. The range of parameters covers an extremely broad plasma regime, all the way from the very weak coupling of space plasmas to the very strong coupling of solid plasmas. The first result is that simulations at the same Lambda(sub D), where Lambda(sub D) is the plasma coupling parameter, but at different combinations of n and theta, behave exactly the same. As a function of Lambda(sub D), the form of p(v), the velocity distribution function of v, the magnitude of v, the velocity vector, is studied. For intermediate to high D, heating is observed in p(v) that obeys conservation of energy, and a suprathermal tail is formed, with a spectral index that depends on Lambda(sub D). For strong coupling (Lambda(sub D) much > 1), the form of the tail is v5, consistent with the findings of Randol and Christian (2014). For weak coupling (Lambda(sub D much <1), no acceleration or heating occurs, as there is no free energy. The dependence on N, the number of particles in the simulation, is also explored. There is a subtle dependence in the index of the tail, such that v5 appears to be the N approaches infinity limit.

Randol, Brent M.↗

Automated Unstructured Grid Adaptation on a Strut Fuel Injector at Hypervelocity Flow Conditions

Computational fluid dynamics (CFD) analysis is presented with the use of an automated unstructured grid adaptation tool on a strut fuel injector at hypervelocity flow conditions. The analysis was carried out with the VULCAN-CFD solver using Reynolds-averaged simulations (RAS). The hypervelocity flow conditions match the high Mach number flow of the experiments conducted as part of the Enhanced Injection and Mixing Project (EIMP) at the NASA Langley Research Center (LaRC). The current work utilizes an automated grid adaptation tool recently implemented into VULCAN-CFD, and explores this tool’s ability to solve high-speed mixing problems. Simulation results obtained using the unstructured adaptive grid approach are compared to those on a user generated structured grid. These results are evaluated by analyzing how efficiently comparable fidelity results are obtained from both adapted and structured simulations. In addition, two adaptation strategies were used to explore the impact on the final solution. In the current work, the unstructured grid adaptation tool automatically generates unstructured grids and performs adaptation of the grid based on a Hessian error estimate of a specified flow field parameter. Multiple adaptations were executed using each run strategy with the one-dimensional values of the mixing efficiency used to determine grid convergence and for comparison with the structured grid simulation results. It was found that the unstructured adaptive grid simulations were able to produce results that matched closely with those on structured grids using far fewer grid cells, and thus, requiring far less computational time to reach the solution. It was also discovered that the adaptation run strategy influenced the total number of grid cells and the efficiency with which a final grid-adapted solution was reached. Overall, the investigation demonstrated that the automated unstructured grid adaptation tool implemented in VULCAN-CFD is capable of accurately and efficiently solving complex highspeed mixing problems using only a fraction of the grid cells required to obtain comparable results using a user-generated structured grid.

grid adaptation↗

GRMHD and GRPIC Simulations

We have developed a new three-dimensional general relativistic magnetohydrodynamic (GRMHD) code by using a conservative, high-resolution shock-capturing scheme. The numerical fluxes are calculated using the HLL approximate Riemann solver scheme. The flux-interpolated constrained transport scheme is used to maintain a divergence-free magnetic field. We have performed various 1-dimensional test problems in both special and general relativity by using several reconstruction methods and found that the new 3D GRMHD code shows substantial improvements over our previous code. The simulation results show the jet formations from a geometrically thin accretion disk near a nonrotating and a rotating black hole. We will discuss the jet properties depended on the rotation of a black hole and the magnetic field configuration including issues for future research. A General Relativistic Particle-in-Cell Code (GRPIC) has been developed using the Kerr-Schild metric. The code includes kinetic effects, and is in accordance with GRMHD code. Since the gravitational force acting on particles is extreme near black holes, there are some difficulties in numerically describing these processes. The preliminary code consists of an accretion disk and free-falling corona. Results indicate that particles are ejected from the black hole. These results are consistent with other GRMHD simulations. The GRPIC simulation results will be presented, along with some remarks and future improvements. The emission is calculated from relativistic flows in black hole systems using a fully general relativistic radiative transfer formulation, with flow structures obtained by GRMHD simulations considering thermal free-free emission and thermal synchrotron emission. Bright filament-like features protrude (visually) from the accretion disk surface, which are enhancements of synchrotron emission where the magnetic field roughly aligns with the line-of-sight in the co-moving frame. The features move back and forth as the accretion flow evolves, but their visibility and morphology are robust. We would like to extend this research using GRPIC simulations and examine a possible new mechanism for certain X-ray quasi-periodic oscillations (QPOs) observed in blackhole X-ray binaries.

Nishikawa, K.-I.↗

Hierarchical Control of Distributed Battery Energy Storage Systems in a DC Microgrid

Contents - NASA Future Power Systems - DC Microgrid with Distributed Battery Energy Storage - A simplified gateway power system - Droop Control (Primary Control) - Distributed battery units, droop control, and load sharing characteristics - Simulation results - Three-Level Hierarchical Control of Distributed Battery Units During Eclipse - Secondary Control - Objective and definition of unit control error - Implementation of secondary control - Simulation results - Tertiary Control - Schedule BU current sharing with the objective of equal SoC - Simulation results•Conclusions

Jing Zhang↗

Boundedness regions of discrete-time dynamic systems

Techniques for obtaining quantitative information about boundedness properties are developed and applied to the sampled-data control of satellite attitude with quantization. Relevant stability concepts are introduced as a series of definitions, and interrelationships between various definitions are discussed. The boundedness regions are estimated by means of quadratic Liapunov functions, and a sufficient condition for the existence of a boundedness region is given for a certain class of systems. A quadratic Liapunov function is applied to the Lur'e-Postinkov class of systems, where the linear part of the system is not asymptotically stable and the quantizer represents the nonlinear characteristic. A numerical calculation of the region of boundedness estimates is performed for satellite attitude control and is compared with simulation results. It is tentatively concluded that the Liapunov results may be good and that simulation results may be difficult to interpret and time-consuming to generate. The Lur'e-based technique yields estimates of regions of absolute boundedness, but at the cost of greater analytical complexity.

Siljak, D.↗

MIA computer simulation test results report

Results of the first noise susceptibility computer simulation tests of the complete MIA receiver analytical model are presented. Computer simulation tests were conducted with both Gaussian and pulse noise inputs. The results of the Gaussian noise tests were compared to results predicted previously and were found to be in substantial agreement. The results of the pulse noise tests will be compared to the results of planned analogous tests in the Data Bus Evaluation Laboratory at a later time. The MIA computer model is considered to be fully operational at this time.

Unger, G. E.↗

Validation of a High-Fidelity Supersonic Parachute Inflation Dynamics Model and Best Practice

The parachute inflation dynamics (PID) of a Disk-Gap-Band (DGB) parachute system travelling at supersonic velocities in the high-altitude Earth atmosphere is simulated using a high-fidelity, multiphysics, massively-parallel computational framework. The sensitivity of the simulation results with respect to critical modelling and simulation parameters, including the fluid mechanical governing equations (Euler vs Navier-Stokes with RANS turbulence closure vs Navier-Stokes with LES turbulence closer), the interaction between the flow and the suspension lines (accounted for vs not accounted for), the dynamics of the payload (restrained vs free), and the material model of the heterogeneous fabric canopy (Saint-Venant Kirchoff vs mechanics-informed multiscale neural network), is extensively investigated. The parameter investigation is enabled by advancements to the state-of-the-art in the modelling and simulation of the supersonic parachute inflation process, which are presented as part of the computational framework. The predictive quality of the computational framework in the context of supersonic PID is validated by comparison of the simulation results with the experimental results from NASA’s Advanced Supersonic Parachute Research (ASPIRE) project. Through the validation, best practices for the modelling and simulation of supersonic PID are developed and proposed.

Rabinovitch, Jason↗

Functional models of power electronic components for system studies

A novel approach to model power electronic circuits has been developed to facilitate simulation studies of system-level issues. The underlying concept for this approach is to develop an equivalent circuit, the functional model, that performs the same functions as the actual circuit but whose operation can be simulated by using larger time step size and the reduction in model complexity, the computation time required by a functional model is significantly shorter than that required by alternative approaches. The authors present this novel modeling approach and discuss the functional models of two major power electronic components, the DC/DC converter unit and the load converter, that are being considered by NASA for use in the Space Station Freedom electric power system. The validity of these models is established by comparing the simulation results with available experimental data and other simulation results obtained by using a more established modeling approach. The usefulness of this approach is demonstrated by incorporating these models into a power system model and simulating the system responses and interactions between components under various conditions.

Tam, Kwa-Sur↗

Electron injection and acceleration at nonlinear shocks: Results of numerical simulations

We present results of numerical simulations of electron injection and acceleration at nonlinear high Mach number shocks. The electrons are assumed to be heated at the thermal subshock to an energy E(sub inj), which is treated as a free parameter, above which they are injected by self-generated whistlers to momentum m(sub p) x V(sub A). This injection mechanism requires Mach numbers greater than (43/(beta(sub -)))(((k T(sub e))/(E(sub inj)))(exp 1/2)), where T(sub e) and beta (sub -) are the upstream electron temperature and plasma beta parameter. Above m(sub p) x V(sub A) electrons are trapped in the shock by Alfven waves. In the proton precursor region the Alfven waves are assumed to be generated by protons accelerated at the shock, and have nonlinear intensities. Below GeV, however, electrons of a given rigidity propagate faster than protons with a similar rigidity and therefore diffuse to regions ahead of the proton precursor. In those regions the Alfven waves are generated by the electrons themselves. The diffusion coefficient appears to increase with decreasing acceleration efficiency. As a result, the number of electrons accelerated to energies GeV and above and, hence, the electron to proton ratio, depend only weakly on the extent of electron heating at the subshock. The negative feedback also renders the electron spectra insensitive to shock compression ratio and smoothing length scale. The estimated e/p ratio at GeV is between approximately 1%-10%.

Levinson, Amir↗

Simulation test results for lift/cruise fan research and technology aircraft

A flight simulation program was conducted on the flight simulator for advanced aircraft (FSAA). The flight simulation was a part of a contracted effort to provide a lift/cruise fan V/STOL aircraft mathematical model for flight simulation. The simulated aircraft is a configuration of the Lift/Cruise Fan V/STOL research technology aircraft (RTA). The aircraft was powered by three gas generators driving three fans. One lift fan was installed in the nose of the aircraft, and two lift/cruise fans at the wing root. The thrust of these fans was modulated to provide pitch and roll control, and vectored to provide yaw, side force control, and longitudinal translation. Two versions of the RTA were defined. One was powered by the GE J97/LF460 propulsion system which was gas-coupled for power transfer between fans for control. The other version was powered by DDA XT701 gas generators driving 62 inch variable pitch fans. The flight control system in both versions of the RTA was the same.

Bland, M. P.↗

An analysis of the 70-meter antenna hydrostatic bearing by means of computer simulation

Recently, the computer program 'A Computer Solution for Hydrostatic Bearings with Variable Film Thickness,' used to design the hydrostatic bearing of the 70-meter antennas, was modified to improve the accuracy with which the program predicts the film height profile and oil pressure distribution between the hydrostatic bearing pad and the runner. This article presents a description of the modified computer program, the theory upon which the computer program computations are based, computer simulation results, and a discussion of the computer simulation results.

Bartos, R. D.↗

Crossflow-Vortex Breakdown on Swept Wings: Correlation of Nonlinear Physics

The spatial evolution of cross flow-vortex packets in a laminar boundary layer on a swept wing are computed by the direct numerical simulation of the incompressible Navier- Stokes equations. A wall-normal velocity distribution of steady suction and blowing at the wing surface is used to generate a strip of equally spaced and periodic disturbances along the span. Three simulations are conducted to study the effect of initial amplitude on the disturbance evolution, to determine the role of traveling cross ow modes in transition, and to devise a correlation function to guide theories of transition prediction. In each simulation, the vortex packets first enter a chordwise region of linear independent growth, then, the individual packets coalesce downstream and interact with adjacent packets, and, finally, the vortex packets nonlinearly interact to generate inflectional velocity profiles. As the initial amplitude of the disturbance is increased, the length of the evolution to breakdown decreases. For this pressure gradient, stationary modes dominate the disturbance evolution. A two-coeffcient function was devised to correlate the simulation results. The coefficients, combined with a single simulation result, provide sufficient information to generate the evolution pattern for disturbances of any initial amplitude.

Joslin, R. D.↗