Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Coupled Simulation”

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 649 records · Page 36

LADEE Multi-Domain Simulation

The Lunar Atmosphere Dust Environment Explorer (LADEE) was a small explorer class spacecraft that was launched on Sept 7, 2013 and that was de-orbited and successfully impacted the Moons surface on April 17, 2014 after completing all of the mission objectives. The low-cost rapidly prototyped hardware design used for the spacecraft was extend to the development of the software base. To achieve this goal, a Model Based Design approach was utilized to develop the onboard flight software, and out of this development a model based multipurpose simulator was created of the LADEE spacecraft and its mission environment. This simulator extended the traditional function of propagating the vehicle's kinematic and rotational states and included the electrical and thermal states propagation. Traditionally, these domains are handled by domain specific high fidelity simulations that use the states histories from other domains as input. By reducing the fidelity and abstracting the relevant features being monitored and controlled by the flight software, it was possible to model the coupling across these domains resulting in more accurate overall system behavior. A faster than real-time workstation (WSIM) version of the LADEE simulator was used to develop and test the software control algorithms in the Simulink environment. To maximize the performance of the simulation, modeling knobs were introduced to reduce the resolution of the some of the domains models when the effects of that domain were not significant for the scope of that simulation. The automatic code generation feature in Simulink was used to port the simulation to several real-time environments to support Processor-in-the-Loop (PIL) and Hardware-in-the-Loop (HIL) testing, verification and validation. The real-time environment required that the design of each of the domain models be deterministic as possible in the time required to perform all of the calculations to update its states. The simulation interface was designed to be compatible with the command interface employed by the LADEE mission operation team. The WSIM, PIL, and HIL simulators thus used a common interface and thus were used for flight software testing, for mission operations personnel training (nominal and off-nominal operations) prior to the mission and to perform command sequence verification during the mission. This presentation will look at the modeling strategies used to create a common interface to the simulator and to model and couple multiple domains within the simulation, the results of those strategies, and the lessons learned.

Multi-Domain↗

Enabling Parametric Optimal Ascent Trajectory Modeling During Early Phases of Design

During the early phases of engineering design, the costs committed are high, costs incurred are low, and the design freedom is high. It is well documented that decisions made in these early design phases drive the entire design's life cycle. In a traditional paradigm, key design decisions are made when little is known about the design. As the design matures, design changes become more difficult -- in both cost and schedule -- to enact. Indeed, the current capability-based paradigm that has emerged because of the constrained economic environment calls for the infusion of knowledge acquired during later design phases into earlier design phases, i.e. bring knowledge acquired during preliminary and detailed design into pre-conceptual and conceptual design. An area of critical importance to launch vehicle design is the optimization of its ascent trajectory, as the optimal trajectory will be able to take full advantage of the launch vehicle's capability to deliver a maximum amount of payload into orbit. Hence, the optimal ascent trajectory plays an important role in the vehicle's affordability posture as the need for more economically viable access to space solutions are needed in today's constrained economic environment. The problem of ascent trajectory optimization is not a new one. There are several programs that are widely used in industry that allows trajectory analysts to, based on detailed vehicle and insertion orbit parameters, determine the optimal ascent trajectory. Yet, little information is known about the launch vehicle early in the design phase - information that is required of many different disciplines in order to successfully optimize the ascent trajectory. Thus, the current paradigm of optimizing ascent trajectories involves generating point solutions for every change in a vehicle's design parameters. This is often a very tedious, manual, and time-consuming task for the analysts. Moreover, the trajectory design space is highly non-linear and multi-modal due to the interaction of various constraints. Additionally, when these obstacles are coupled with The Program to Optimize Simulated Trajectories [1] (POST), an industry standard program to optimize ascent trajectories that is difficult to use, it requires expert trajectory analysts to effectively optimize a vehicle's ascent trajectory. As it has been pointed out, the paradigm of trajectory optimization is still a very manual one because using modern computational resources on POST is still a challenging problem. The nuances and difficulties involved in correctly utilizing, and therefore automating, the program presents a large problem. In order to address these issues, the authors will discuss a methodology that has been developed. The methodology is two-fold: first, a set of heuristics will be introduced and discussed that were captured while working with expert analysts to replicate the current state-of-the-art; secondly, leveraging the power of modern computing to evaluate multiple trajectories simultaneously, and therefore, enable the exploration of the trajectory's design space early during the pre-conceptual and conceptual phases of design. When this methodology is coupled with design of experiments in order to train surrogate models, the authors were able to visualize the trajectory design space, enabling parametric optimal ascent trajectory information to be introduced with other pre-conceptual and conceptual design tools. The potential impact of this methodology's success would be a fully automated POST evaluation suite for the purpose of conceptual and preliminary design trade studies. This will enable engineers to characterize the ascent trajectory's sensitivity to design changes in an arbitrary number of dimensions and for finding settings for trajectory specific variables, which result in optimal performance for a "dialed-in" launch vehicle design. The effort described in this paper was developed for the Advanced Concepts Office [2] at NASA Marshall Space Flight Center

Holt, James B.↗

Numerical simulation of rarefied gas flow through a slit

Two different approaches, the finite-difference method coupled with the discrete-ordinate method (FDDO), and the direct-simulation Monte Carlo (DSMC) method, are used in the analysis of the flow of a rarefied gas from one reservoir to another through a two-dimensional slit. The cases considered are for hard vacuum downstream pressure, finite pressure ratios, and isobaric pressure with thermal diffusion, which are not well established in spite of the simplicity of the flow field. In the FDDO analysis, by employing the discrete-ordinate method, the Boltzmann equation simplified by a model collision integral is transformed to a set of partial differential equations which are continuous in physical space but are point functions in molecular velocity space. The set of partial differential equations are solved by means of a finite-difference approximation. In the DSMC analysis, three kinds of collision sampling techniques, the time counter (TC) method, the null collision (NC) method, and the no time counter (NTC) method, are used.

Keith, Theo G., Jr.↗

Computational Analysis of a Quiet Single-Main Rotor Helicopter for Air Taxi Operations

A computational study has been done on NASA’s quiet single-main rotor helicopter concept for urban air mobility using high-fidelity computational fluid dynamics, rotorcraft comprehensive analysis tools, and computational aeroacoustics. High-order accurate schemes, dual-time stepping, and the delayed detached-eddy simulation model have been employed. A loose-coupling approach between the flow solver and the rotorcraft comprehensive code is implemented to include vehicle trim and blade motions. Acoustic simulations based on the Ffowcs-Williams and Hawkings equations have been performed to compute the rotor noise. Different blade geometries and tip speeds are analyzed, and performance and acoustics results are compared. The vehicle has been simulated in hover and cruise, with flight conditions representative of an air taxi mission. The quiet single-main rotor helicopter is one of the conceptual designs intended to focus and guide NASA’s research activities in support of aircraft development for vertical take-off and landing air taxi operations.

ARMD↗

(abstract) Role of Long Equatorial Wave Reflection in the Low-Frequency Variability Observed During the TOPEX/POSEIDON Period: Data Analysis and Coupled Model Study

The analysis of TOPEX/POSEIDON sea level and ERS-1 zonal wind stress data gives evidence that first mode meridonal long Rossby waves fully reflected into Kelvin waves at the equatorial western Pacific boundary at different periods including the January-June 1994 period. The evolution of the conditions (zonal wind stress, sea surface temperature and sea level anomalies) in the central Pacific suggests that this reflection may have played a role in the reversal of cold to warm anomalies observed in the central Pacific in June-July 1994. To investigate the actual role of reflected Kelvin waves observed by TOPEX/POSEIDON, a simple ocean-atmosphere coupled model is designed. Starting from these initial conditions, several simulations are run for testing the role western bundary reflection evidenced in TOPEX/POSEIDON data. Coupled mechanisms are suggested to explain the low-frequency variability observed during this period.

Rossby waves Kelvin waves zonal wind stress TOPEX/↗

Modeling of Convective-Stratiform Precipitation Processes: Sensitivity to Partitioning Methods and Numerical Advection Schemes

Six different convective-stratiform separation techniques, including a new technique that utilizes the ratio of vertical and terminal velocities, are compared and evaluated using two-dimensional numerical simulations of a tropical [Tropical Ocean Global Atmosphere Coupled Ocean-Atmosphere Response Experiment (TOGA COARE)] and midlatitude continental [Preliminary Regional Experiment for STORM-Central (PRESTORM)] squall line. The simulations are made using two different numerical advection schemes: 4th order and positive definite advection. Comparisons are made in terms of rainfall, cloud coverage, mass fluxes, apparent heating and moistening, mean hydrometeor profiles, CFADs (Contoured Frequency with Altitude Diagrams), microphysics, and latent heating retrieval. Overall, it was found that the different separation techniques produced results that qualitatively agreed. However, the quantitative differences were significant. Observational comparisons were unable to conclusively evaluate the performance of the techniques. Latent heating retrieval was shown to be sensitive to the use of separation technique mainly due to the stratiform region for methods that found very little stratiform rain. The midlatitude PRESTORM simulation was found to be nearly invariant with respect to advection type for most quantities while for TOGA COARE fourth order advection produced numerous shallow convective cores and positive definite advection fewer cells that were both broader and deeper penetrating above the freezing level.

Lang, Steve↗

Real-Time Modeling of Thermal Interactions in Cryogenic Ball Bearings

Based on the integration of classical differential equations of motion of bearing elements, a real-time dynamic simulation of ball bearing performance is coupled with simultaneous modeling of thermal interactions between the bearing elements. The transient heat generations are time averaged over a significantly larger time step before they are used to compute the changes in bearing temperature fields. Such an averaging algorithm eliminates the numerical difficulties associated with simultaneous integration of mechanical thermal differential equations with vastly different time scales. Under stable operation, the step change in temperatures at each thermal time step converges to a steady-state solution. Following is a list of accomplishments in the present investigation: (1) A real-time dynamic performance of both 440C and hybrid ball bearings in a LOX environment is modeled in both transient and steady-state time domains, as the inner race accelerates to a prescribed operating speed. (2) The frictional interactions at the ball to race contacts are modeled by experimentally measured traction coefficients as a function of slide-to-roll ratios. (3) Bearing heat generation, as defined by the heat transferred to circulating LOX, is validated against experimental values determined by measured LOX temperatures under prescribed pressures and flow rates. The predicted heat generation is in good agreement with the experimental values. (4) The steady-state solutions are shown to be independent of initial conditions. This establishes convergence of the numerical integration of a differential equation of motion of the bearing elements. (5) Churning and drag losses in the circulating LOX constitutes the majority of bearing heat generation under the experimental test conditions. (6) The heat dissipated in frictional interactions in hybrid bearings is significantly lower in hybrid bearings in comparison to that simulated in all-steel bearings, while the churning and drag losses are relatively unchanged.

Gupta, P. K.↗

A helicopter flight investigation of roll-control sensitivity, damping and cross coupling in a low altitude lateral maneuvering task

A helicopter in-flight simulation was conducted to determine the effects of variations in roll damping, roll sensitivity, and pitch and roll rate cross-coupling on helicopter flying qualities in a low altitude maneuver. The experiment utilized the UH-1H helicopter in-flight simulator, which is equipped with the V/STOLAND avionics system. The response envelope of this vehicle allowed simulation of configurations with low to moderate damping and sensitivity. A visual, low level slalom course was set up, consisting of constant speed and constant altitude S-turns around the 1000 ft makers of an 8000 ft runway. Results are shown in terms of Cooper-Harper pilot ratings, pilot commentary, and statistical and frequency analyses of the lateral characteristics. These results show good consistency with previous ground simulator results and are compared with existing flying qualities criteria.

Corliss, L. D.↗

An Ocean-Atmosphere Simulation for Studying Air-Sea Interactions

During the past few years the Goddard Earth Observing System (GEOS) and Massachusetts Institute of Technology (MIT) modeling groups have produced, respectively, global atmosphere-only and ocean-only simulations with km-scale grid spacing. These simulations have proved invaluable for process studies and for the development of satellite and in-situ sampling strategies. Nevertheless, a key limitation of these "nature" simulations is the lack of interactivity between the ocean and the atmosphere, which limits their usefulness for studying air-sea interactions and for designing observing missions to study these interactions. To remove this limitation, we aim to couple the km-scale GEOS atmosphere simulation to the km-scale MIT ocean simulation.

Strobach, Ehud↗

Interpreting the Ultraviolet Aerosol Index Observed with the OMI Satellite Instrument to Understand Absorption by Organic Aerosols: Implications for Atmospheric Oxidation and Direct Radiative Effects

Satellite observations of the ultraviolet aerosol index (UVAI) are sensitive to absorption of solar radiation by aerosols; this absorption affects photolysis frequencies and radiative forcing. We develop a global simulation of the UVAI using the 3-D chemical transport model GEOSChem coupled with the Vector Linearized Discrete Ordinate Radiative Transfer model (VLIDORT). The simulation is applied to interpret UVAI observations from the Ozone Monitoring Instrument (OMI) for the year 2007. Simulated and observed values are highly consistent in regions where mineral dust dominates the UVAI, but a large negative bias (-0.32 to -0.97) exists between simulated and observed values in biomass burning regions. We determine effective optical properties for absorbing organic aerosol, known as brown carbon (BrC), and implement them into GEOS-Chem to better represent observed UVAI values over biomass burning regions. The inclusion of absorbing BrC decreases the mean bias between simulated and OMI UVAI values from -0.57 to -0.09 over West Africa in January, from -0.32 to +0.0002 over South Asia in April, from -0.97 to -0.22 over southern Africa in July, and from -0.50 to +0.33 over South America in September. The spectral dependence of absorption after including BrC in the model is broadly consistent with reported observations for biomass burning aerosol, with absorbing Angstrom exponent (AAE) values ranging from 2.9 in the ultraviolet (UV) to 1.3 across the UV-Near IR spectrum. We assess the effect of the additional UV absorption by BrC on atmospheric photochemistry by examining tropospheric hydroxyl radical (OH) concentrations in GEOS-Chem. The inclusion of BrC decreases OH by up to 30% over South America in September, up to 20% over southern Africa in July, and up to 15% over other biomass burning regions. Global annual mean OH concentrations in GEOS-Chem decrease due to the presence of absorbing BrC, increasing the methyl chloroform lifetime from 5.62 to 5.68 years, thus reducing the bias against observed values. We calculate the direct radiative effect (DRE) of BrC using GEOS-Chem coupled with the radiative transfer model RRTMG (GC-RT). Treating organic aerosol as containing more strongly absorbing BrC changes the global annual mean all-sky top of atmosphere (TOA) DRE by +0.03Wm(exp -2) and all-sky surface DRE by -0.08Wm(exp -2). Regional changes of up to +0.3Wm(exp -2) at TOA and down to -1.5Wm(exp -2) at the surface are found over major biomass burning regions.

Air Pollution↗

Computational Non-Destructive Evaluation Improving Ultrasonic Interrogation of Complex Geometry Composite Parts

Finite element simulation was employed in modeling the ultrasound (UT) pressure pulse propagation through a coupled liquid-composite medium to reproduce experimental data. From the simulation point of view, the proposed approach is challenging when there is a large simulation domain. For example, it is shown that a sub-micron wavelength of an ultrasonic wave requires a mesh size of several microns and this in turn requires significant computational resources, as well as special care in modeling. Some of the simulation results are presented considering that such modeling should reproduce experimental data for a healthy and faulty composite structure with complex geometry. Many possible experimental setups are simulated to demonstrate the non-destructive testing technique. This setup includes the generation of pressure pulse propagating through the tested composite plate and possible scattering by discontinuities (area of different impedance) that may be present in the panel. This scattered pulse together with the baseline pressure pulse generates a signature on the probe element which can be used to locate the position of defects in the structures.

modeling↗

Analysis of internal flows relative to the space shuttle main engine

Cooperative efforts between the Lockheed-Huntsville Computational Mechanics Group and the NASA-MSFC Computational Fluid Dynamics staff has resulted in improved capabilities for numerically simulating incompressible flows generic to the Space Shuttle Main Engine (SSME). A well established and documented CFD code was obtained, modified, and applied to laminar and turbulent flows of the type occurring in the SSME Hot Gas Manifold. The INS3D code was installed on the NASA-MSFC CRAY-XMP computer system and is currently being used by NASA engineers. Studies to perform a transient analysis of the FPB were conducted. The COBRA/TRAC code is recommended for simulating the transient flow of oxygen into the LOX manifold. Property data for modifying the code to represent LOX/GOX flow was collected. The ALFA code was developed and recommended for representing the transient combustion in the preburner. These two codes will couple through the transient boundary conditions to simulate the startup and/or shutdown of the fuel preburner. A study, NAS8-37461, is currently being conducted to implement this modeling effort.

Source record↗

The Nike canted-fin Apache

Nike-Apache sounding rocket vehicle with canted fins, examining pitch-roll coupling characteristics by equilibrium and dynamic simulation methods

Mayo, E. E.↗

The motion of interconnected flexible bodies

The equations of motion for an arbitrarily interconnected collection of substructures are derived. The substructures are elastic bodies which may be idealized as finite element assemblies and are subject to small deformations relative to a nominal state. Interconnections between the elastic substructures permit large relative translations and rotations between substructures, governed by Pfaffian constraints describing the connections. Screw connections (permitting rotation about and translation along a single axis) eliminate constraint forces and incorporate modal coupling. The problem of flexible spacecraft simulation is discussed. Hurty's component mode approach is extended by permitting interconnected elastic substructures large motions relative to each other and relative to inertial space. The hybrid coordinate methods are generalized by permitting all substructures to be flexible (rather than only the terminal members of a topological tree of substructures). The basic relationships of continuum mechanics are developed.

Hopkins, A. S.↗

An AI approach for scheduling space-station payloads at Kennedy Space Center

The Payload Processing for Space-Station Operations (PHITS) is a prototype modeling tool capable of addressing many Space Station related concerns. The system's object oriented design approach coupled with a powerful user interface provide the user with capabilities to easily define and model many applications. PHITS differs from many artificial intelligence based systems in that it couples scheduling and goal-directed simulation to ensure that on-orbit requirement dates are satisfied.

Castillo, D.↗

Comprehensive rotorcraft analysis methods

The development and application of comprehensive rotorcraft analysis methods in the field of rotorcraft technology are described. These large scale analyses and the resulting computer programs are intended to treat the complex aeromechanical phenomena that describe the behavior of rotorcraft. They may be used to predict rotor aerodynamics, acoustic, performance, stability and control, handling qualities, loads and vibrations, structures, dynamics, and aeroelastic stability characteristics for a variety of applications including research, preliminary and detail design, and evaluation and treatment of field problems. The principal comprehensive methods developed or under development in recent years and generally available to the rotorcraft community because of US Army Aviation Research and Technology Activity (ARTA) sponsorship of all or part of the software systems are the Rotorcraft Flight Simulation (C81), Dynamic System Coupler (DYSCO), Coupled Rotor/Airframe Vibration Analysis Program (SIMVIB), Comprehensive Analytical Model of Rotorcraft Aerodynamics and Dynamics (CAMRAD), General Rotorcraft Aeromechanical Stability Program (GRASP), and Second Generation Comprehensive Helicopter Analysis System (2GCHAS).

Stephens, Wendell B.↗

Cryogenic Limb Array Etalon Spectrometer (CLAES) S/N001 cryostat test results

The results of a test program designed to verify the operational and thermal performance of the Cryogenic Limb Array Etalon Spectrometer (CLAES) S/N001 Ne/CO2 cryostat, as monitored and serviced by the Data Acquisition System and Ground Support Equipment, are presented. The S/N001 model test assembly consists of the S/N001 configuration coupled with an instrument thermal mass simulator (ITMS) for simulating the thermal behavior of the instrument package. The test data for the cryostat include the rate of temperature decrease of the various parts of the cryostat/ITMS and the final equilibrium temperatures of the shields, thermal links, and tank. The test results showed that system heat rates, operating temperatures, and expected end-of-life temperature predictions all agree with predictions within estimated margins.

Horsley, W. J.↗

Free wake analysis of hover performance using a new influence coefficient method

A new approach to the prediction of helicopter rotor performance using a free wake analysis was developed. This new method uses a relaxation process that does not suffer from the convergence problems associated with previous time marching simulations. This wake relaxation procedure was coupled to a vortex-lattice, lifting surface loads analysis to produce a novel, self contained performance prediction code: EHPIC (Evaluation of Helicopter Performance using Influence Coefficients). The major technical features of the EHPIC code are described and a substantial amount of background information on the capabilities and proper operation of the code is supplied. Sample problems were undertaken to demonstrate the robustness and flexibility of the basic approach. Also, a performance correlation study was carried out to establish the breadth of applicability of the code, with very favorable results.

Quackenbush, Todd R.↗