Search NASA⌕ Search

SEARCH · Search NASA

Results for “simulation codes”

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 145 records · Page 8

A Magnetohydrodynamic Boost for Relativistic Jets

We have performed relativistic magnetohydrodynamic simulations of the hydrodynamic boosting mechanism for relativistic jets explored by Aloy & Rezzolla (2006) using the RAISHIN code. Simulation results show that the presence of a magnetic field may change the properties of the shock interface between the tenuous, overpressured jet (V(sub j) (sup z)) flowing tangentially to a dense external medium. Magnetic fields can lead to more efficient acceleration of the jet, in comparison to the pure-hydrodynamic case. A poloidal magnetic field (B(sup z)), tangent to the interface and parallel to the jet flow, produces both a stronger outward moving shock and inward moving rarefaction wave. This leads to a large velocity component normal to the interface in addition to acceleration tangent to the interface, and the jet is thus accelerated to a larger Lorentz factors than those obtained in the pure-hydrodynamic case. In contrast, a strong toroidal magnetic field (B(sup y)), tangent to the interface but perpendicular to the jet flow, also leads to stronger acceleration tangent to the shock interface relative to the pure-hydrodynamic case, but to a lesser extent than found for the poloidal case due to the fact that the velocity component normal to the shock interface is now much smaller. Overall, the acceleration efficiency in the toroidal case is less than that of the poloidal case but both geometries still result in higher Lorentz factors than the pure-hydrodynamic case. Thus, the presence and relative orientation of a magnetic field in relativistic jets can have a significant influence on the hydrodynamic boost mechanism studied by Aloy & Rezzolla (2006).

Mizuno, Yosuke↗

A Magnetohydrodynamic Boost for Relativistic Jets

We performed relativistic magnetohydrodynamic simulations of the hydrodynamic boosting mechanism for relativistic jets explored by Aloy & Rezzolla (2006) using the RAISHIN code. Simulation results show that the presence of a magnetic field changes the properties of the shock interface between the tenuous, overpressured jet (V^z j) flowing tangentially to a dense external medium. We find that magnetic fields can lead to more efficient acceleration of the jet, in comparison to the pure-hydrodynamic case. A "poloidal" magnetic field (B^z), tangent to the interface and parallel to the jet flow, produces both a stronger outward moving shock and a stronger inward moving rarefaction wave. This leads to a large velocity component normal to the interface in addition to acceleration tangent to the interface, and the jet is thus accelerated to larger Lorentz factors than those obtained in the pure-hydrodynamic case. Likewise, a strong "toroidal" magnetic field (B^y), tangent to the interface but perpendicular to the jet flow, also leads to stronger acceleration tangent to the shock interface relative to the pure-hydrodynamic case. Thus. the presence and relative orientation of a magnetic field in relativistic jets can significant modify the hydrodynamic boost mechanism studied by Aloy & Rezzolla (2006).

Mizuno, Yosuke↗

Simulation of spacecraft attitude dynamics using TREETOPS and model-specific computer Codes

The simulation of spacecraft attitude dynamics and control using the generic, multi-body code called TREETOPS and other codes written especially to simulate particular systems is discussed. Differences in the methods used to derive equations of motion--Kane's method for TREETOPS and the Lagrangian and Newton-Euler methods, respectively, for the other two codes--are considered. Simulation results from the TREETOPS code are compared with those from the other two codes for two example systems. One system is a chain of rigid bodies; the other consists of two rigid bodies attached to a flexible base body. Since the computer codes were developed independently, consistent results serve as a verification of the correctness of all the programs. Differences in the results are discussed. Results for the two-rigid-body, one-flexible-body system are useful also as information on multi-body, flexible, pointing payload dynamics.

Cochran, John E.↗

Nonlinear Waves and Particle Heating in the Auroral Ionosphere

In the two year period, the researchers have developed innovative 2-D simulation codes for modeling the interaction between electron beams and plasma waves and particles in the auroral ionosphere. These new simulations include kinetic wave particles as well as wave-wave interactions near and above approximately 500 km. The principal numerical model consists of partial-differential equations which evolve large amplitude magnetized 2-D Langmuir waves self-consistently with the electron distribution function. Progress was also made in developing 2-D Particle in Cell (PIC) codes.

Goldman, Martin↗

Three-dimensional numerical simulation of interplanetary magnetic field changes at 1 AU as a consequence of simulated solar flares

In order to study the relationship between different forms of activity and transient variations of the north-south component, B(sub z), of the Interplanetary Magnetic Field (IMF), at 1 AU, a three dimensional numerical simulation code was employed to study several aspects of this problem. One form of solar activity, the flare, is simulated by using a pressure pulse at different locations near the solar surface and observing the simulated IMF evolution of B(sub theta)(-B(sub z)) at 1 AU. It was found, for a given pressure pulse, that the orientation of the corresponding transient variation of B(sub z) has a strong relationship with the location of the pressure pulse and the initial condition of IMF.

Wu, S. T.↗

ISS Destiny Laboratory Smoke Detection Model

Smoke transport and detection were modeled numerically in the ISS Destiny module using the NIST, Fire Dynamics Simulator code. The airflows in Destiny were modeled using the existing flow conditions and the module geometry included obstructions that simulate the currently installed hardware on orbit. The smoke source was modeled as a 0.152 by 0.152 m region that emitted smoke particulate ranging from 1.46 to 8.47 mg/s. In the module domain, the smoke source was placed in the center of each Destiny rack location and the model was run to determine the time required for the two smoke detectors to alarm. Overall the detection times were dominated by the circumferential flow, the axial flow from the intermodule ventilation and the smoke source strength.

Brooker, John E.↗

Benchmark Active Controls Technology (BACT) Wing CFD Results

The Benchmark Active Controls Technology (BACT) wing test (see chapter 8E) provides data for the validation of aerodynamic, aeroelastic, and active aeroelastic control simulation codes. These data provide a rich database for development and validation of computational aeroelastic and aeroservoelastic methods. In this vein, high-level viscous CFD analyses of the BACT wing have been performed for a subset of the test conditions available in the dataset. The computations presented in this section investigate the aerodynamic characteristics of the rigid clean wing configuration as well as simulations of the wing with a static and oscillating aileron and spoiler deflection. Two computational aeroelasticity codes extensively used at NASA Langley Research Center are implemented in this simulation. They are the ENS3DAE and CFL3DAE computational aeroelasticity programs. Both of these methods solve the three-dimensional compressible Navier-Stokes equations for both rigid and flexible vehicles, but they use significantly different approaches to the solution 6f the aerodynamic equations of motion. Detailed descriptions of both methods are presented in the following section.

Schuster, David M.↗

Trajectories for High Specific Impulse High Specific Power Deep Space Exploration

Flight times and deliverable masses for electric and fusion propulsion systems are difficult to approximate. Numerical integration is required for these continuous thrust systems. Many scientists are not equipped with the tools and expertise to conduct interplanetary and interstellar trajectory analysis for their concepts. Several charts plotting the results of well-known trajectory simulation codes were developed and are contained in this paper. These charts illustrate the dependence of time of flight and payload ratio on jet power, initial mass, specific impulse and specific power. These charts are intended to be a tool by which people in the propulsion community can explore the possibilities of their propulsion system concepts. Trajectories were simulated using the tools VARITOP and IPOST. VARITOP is a well known trajectory optimization code that involves numerical integration based on calculus of variations. IPOST has several methods of trajectory simulation; the one used in this paper is Cowell's method for full integration of the equations of motion. An analytical method derived in the companion paper was also evaluated. The accuracy of this method is discussed in the paper.

Polsgrove, Tara↗

Simulation of Trajectories for High Specific Impulse Deep Space Exploration

Difficulties in approximating flight times and deliverable masses for continuous thrust propulsion systems have complicated comparison and evaluation of proposed propulsion concepts. These continuous thrust propulsion systems are of interest to many groups, not the least of which are the electric propulsion and fusion communities. Several charts plotting the results of well-known trajectory simulation codes were developed and are contained in this paper. These charts illustrate the dependence of time of flight and payload ratio on jet power, initial mass, specific impulse and specific power. These charts are intended to be a tool by which people in the propulsion community can explore the possibilities of their propulsion system concepts. Trajectories were simulated using the tools VARITOP and IPOST. VARITOP is a well known trajectory optimization code that involves numerical integration based on calculus of variations. IPOST has several methods of trajectory simulation; the one used in this paper is Cowell's method for full integration of the equations of motion. The analytical method derived in the companion paper was also used to simulate the trajectory. The accuracy of this method is discussed in the paper.

Polsgrove, Tara↗

Structures of Alfven shocks - S-shaped magnetic hodogram observed at the magnetopause

The Alfven shock structure is obtained by a two-fluid simulation code including dissipations due to resistivity, viscosity and heat conductivity. The hodogram of the tangential magnetic field component of the simulated Alfven shock is S-shaped with a 180-deg angular rotation. This is the first simulation result to show the S-shaped hodogram associated with Alfven shocks. The S-shaped hodograms have been observed at the magnetopause (Berchem and Russell, 1982). On the basis of this result, it is proposed that the magnetopause is an Alfven shock where and when the hodogram is S-shaped and rotates 180 deg.

Lyu, L. H.↗

Monte Carlo Calculations of Neutron Number Spectra and Buildup Factors in Infinite Conical Configurations

A Monte Carlo code simulating neutron transport in infinite cones of water and water-equivalent hydrogen was prepared for an IBM 704 computer. The code was essentially a modification of the point-source, infinite-medium code used in NASA TN D-850. Studies were made of differential neutron number spectra and associated buildup factors for infinite cones having apex half-angles of 15 degrees, 30 degrees, 45 degrees, and 60 degrees. The buildup factors obtained were compared with those for the appropriate infinite medium, which allowed an examination of the effect of solid angle subtended by material on the transport of 6-Mev source neutrons emanating from the cone apex. The variation of number buildup factor with distance for the various cones shows that neutron scattering out of the cones is predominant in the first 30 t o 40 centimeters of material, and that transport beyond this distance is of a similar nature in all the cones.

CONE↗

Some aspects of robotics calibration, design and control

The main objective is to introduce techniques in the areas of testing and calibration, design, and control of robotic systems. A statistical technique is described that analyzes a robot's performance and provides quantitative three-dimensional evaluation of its repeatability, accuracy, and linearity. Based on this analysis, a corrective action should be taken to compensate for any existing errors and enhance the robot's overall accuracy and performance. A comparison between robotics simulation software packages that were commercially available (SILMA, IGRIP) and that of Kennedy Space Center (ROBSIM) is also included. These computer codes simulate the kinematics and dynamics patterns of various robot arm geometries to help the design engineer in sizing and building the robot manipulator and control system. A brief discussion on an adaptive control algorithm is provided.

Tawfik, Hazem↗

LOOPREF: A Fluid Code for the Simulation of Coronal Loops

This report documents the code LOOPREF. LOOPREF is a semi-one dimensional finite element code that is especially well suited to simulate coronal-loop phenomena. It has a full implementation of adaptive mesh refinement (AMR), which is crucial for this type of simulation. The AMR routines are an improved version of AMR1D. LOOPREF's versatility makes is suitable to simulate a wide variety of problems. In addition to efficiently providing very high resolution in rapidly changing regions of the domain, it is equipped to treat loops of variable cross section, any non-linear form of heat conduction, shocks, gravitational effects, and radiative loss.

deFainchtein, Rosalinda↗

Advanced Combustion Modeling for Complex Turbulent Flows

The next generation of aircraft engines will need to pass stricter efficiency and emission tests. NASA's Ultra-Efficient Engine Technology (UEET) program has set an ambitious goal of 70% reduction of NO(x) emissions and a 15% increase in fuel efficiency of aircraft engines. We will demonstrate the state-of-the-art combustion tools developed a t Stanford's Center for Turbulence Research (CTR) as part of this program. In the last decade, CTR has spear-headed a multi-physics-based combustion modeling program. Key technologies have been transferred to the aerospace industry and are currently being used for engine simulations. In this demo, we will showcase the next-generation combustion modeling tools that integrate a very high level of detailed physics into advanced flow simulation codes. Combustor flows involve multi-phase physics with liquid fuel jet breakup, evaporation, and eventual combustion. Individual components of the simulation are verified against complex test cases and show excellent agreement with experimental data.

Ham, Frank Stanford↗

InSight Entry, Descent and Landing Pre-Flight Performance Predictions

On November 26, 2018, the Interior Exploration using Seismic Investigations, Geodesy and Heat Transport (InSight) lander successfully touched down on the surface of Mars. Over its seven-plus year development, NASA Langley Research Center’s (LaRC) Program to Optimize Simulated Trajectories II (POST2) was used to assess the mission’s Entry, Descent and Landing (EDL) vehicle system performance against related requirements across the full range of possible environmental and spacecraft conditions. Much of the simulation code was derived from the Phoenix mission, for which this vehicle is very similar. The InSight six degree-of-freedom simulation included models for Mars atmosphere, gravity and digital elevation maps of the landing location. Additionally, vehicle specific aerodynamic, parachute, engine, navigation sensor, flight software and landing radar models were also included. A set of dispersions for each model, as well as for additional simulation input parameters, were also included in order to provide a statistical, Monte Carlo prediction of the EDL system performance. An overview of the pre-flight performance assessments completed, including the various simulation campaigns used, will be provided. Ultimately, this work was critical in the assessment of readiness for InSight launch. A brief description of the use of this simulation in support of flight operations is also discussed.

Robert W Maddock↗

NASCAP user's manual

The NASCAP (NASA Charging Analyzer Program) code simulates the charging process for a complex object in either tenuous plasma or ground test environment. Detailed specifications needed to run the code are presented. The object definition section, OBJDEF, allows the test object to be easily defined in the cubic mesh. The test object is composed of conducting sections which may be wholly or partially covered with thin dielectric coatings. The potential section, POTENT, obtains the electrostatic potential in the space surrounding the object. It uses the conjugate gradient method to solve the finite element formulation of Poisson's equation. The CHARGE section of NASCAP treats charge redistribution among the surface cells of the object as well as charging through radiation bombardment. NASCAP has facilities for extensive graphical output, including several types of object display plots, potential contour plots, space charge density contour plots, current density plots, and particle trajectory plots.

Mandell, M. J.↗

Computational Icing Analysis on NASA’s SIDRM Geometry to Investigate Collection Efficiency

Computational icing analysis results were compared to experimental icing tunnel data including aerothermal (e.g. dry air) and supercooled water droplet rime-ice conditions from tests conducted in early 2022 at the NASA Icing Research Tunnel (IRT). The Simulated Intercompressor Duct Research Model (SIDRM) test article was used in this study, and its geometry represents the inter-compressor duct region of a turbofan engine. The test article’s purpose is to study the physics of supercooled water icing and ice crystal icing. This study compared three different icing codes: FENSAP-ICE (Eulerian approach), LEWICE3D (Lagrangian approach), and GlennICE (Lagrangian approach). All three icing codes were conducted on SIDRM’s complex body flow-field and compared to different experimental supercooled water rime runs. The test article instrumentation (pressure taps, thermocouples, etc.) and 3D laser scans of final ice shapes were used to compare against the different icing code simulations. The overall objectives are to understand how the icing codes handle capturing collection efficiency on the complex test article’s unheated surfaces. In the aerothermal cases, pressure tap readings matched the CFD results, but dry air CFD underpredicted thermocouple readings. Collection efficiency results from all three icing codes matched well together on the main body leading edge, main body slope, and the strut leading edges. All three icing codes underpredicted the final ice shape using a single-shot constant ice density approach, with more difficulty coming from the strut leading edge ice shape due to the swept wing like flow field. NASA’s overall goal for this effort is to develop computational icing tools to assist in the design and certification of engines for flight in icing conditions.

Supercooled icing, SIDRM, aerodynamics, computatio↗

Numerical simulations of the electrodynamic interactions between the Tethered-Satellite-System and space plasma

The first Tethered-Satellite-System (TSS-1), scheduled for a flight in late 1992, is expected to provide relevant information related to the concept of generating an emf in a 20-km-long (or longer) conducting wire. This paper presents numerical simulations of the electrodynamic interactions between the TSS system and space plasma, using a 2D and 3D models of the system. The 2D case code simulates the motion of a long cylinder past a plasma, which is composed of electrons and H(+) ions. The system is solved by allowing the plasma to flow past the cylinder with an imposed magnetic field. The more complex 3D case is considered to study the dynamics in great detail. Results of 2D simulation show that the interaction of a satellite with plasma flowing perpendicularly to the magnetic field results in an enhancement in the current collection.

Vashi, Bharat I.↗