User's Manual for the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA)
NASA Technical Reports Server (NTRS)
Gnoffo, Peter A.; Cheatwood, F. McNeil
1996-01-01
This user's manual provides detailed instructions for the installation and the application of version 4.1 of the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA). Also provides simulation of flow field in thermochemical nonequilibrium around vehicles traveling at hypersonic velocities through the atmosphere. Earlier versions of LAURA were predominantly research codes, and they had minimal (or no) documentation. This manual describes UNIX-based utilities for customizing the code for special applications that also minimize system resource requirements. The algorithm is reviewed, and the various program options are related to specific equations and variables in the theoretical development.
Distributed-Memory Computing With the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA)
NASA Technical Reports Server (NTRS)
Riley, Christopher J.; Cheatwood, F. McNeil
1997-01-01
The Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA), a Navier-Stokes solver, has been modified for use in a parallel, distributed-memory environment using the Message-Passing Interface (MPI) standard. A standard domain decomposition strategy is used in which the computational domain is divided into subdomains with each subdomain assigned to a processor. Performance is examined on dedicated parallel machines and a network of desktop workstations. The effect of domain decomposition and frequency of boundary updates on performance and convergence is also examined for several realistic configurations and conditions typical of large-scale computational fluid dynamic analysis.
NASA Technical Reports Server (NTRS)
Gnoffo, Peter A.
1989-01-01
An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for three-dimensional, viscous, hypersonic flows in chemical and thermal nonequilibrium is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. The relaxation strategy is well suited for computers employing either vector or parallel architectures, and the relation between computer architecture and algorithm is emphasized. It is also well suited to the numerical solution of the governing equations on unstructured grids. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large. block tri-diagonal systems. A single relaxation step depends only on information from nearest neighbors. Predictions for pressure distributions, surface heating, and aerodynamic coefficients compare well with experimental data for Mach 10 flow over a blunt body. Predictions for the hypersonic flow of air in chemical and thermal nonequilibrium (velocity = 8917 m/s, altitude = 78 km.) over the Aeroassist Flight Experiment (AFE) configuration obtained on a multi-domain grid are discussed.
Numerical Study of Three-Dimensional Flows Using Unfactored Upwind-Relaxation Sweeping Algorithm
NASA Astrophysics Data System (ADS)
Zha, G.-C.; Bilgen, E.
1996-05-01
The linear stability analysis of the unfactored upwind relaxation-sweeping (URS) algorithm for 3D flow field calculations has been carried out and it is shown that the URS algorithm is unconditionally stable. The algorithm is independent of the global sweeping direction selection. However, choosing the direction with relatively low variable gradient as the global sweeping direction results in a higher degree of stability. Three-dimensional compressible Euler equations are solved by using the implicit URS algorithm to study internal flows of a non-axisymmetric nozzle with a circular-to-rectangular transition duct and complex shock wave structures for a 3D channel flow. The efficiency and robustness of the URS algorithm has been demonstrated.
An upwind-biased, point-implicit relaxation algorithm for viscous, compressible perfect-gas flows
NASA Technical Reports Server (NTRS)
Gnoffo, Peter A.
1990-01-01
An upwind-biased, point-implicit relaxation algorithm for obtaining the numerical solution to the governing equations for three-dimensional, viscous, compressible, perfect-gas flows is described. The algorithm is derived using a finite-volume formulation in which the inviscid components of flux across cell walls are described with Roe's averaging and Harten's entropy fix with second-order corrections based on Yee's Symmetric Total Variation Diminishing scheme. Viscous terms are discretized using central differences. The relaxation strategy is well suited for computers employing either vector or parallel architectures. It is also well suited to the numerical solution of the governing equations on unstructured grids. Because of the point-implicit relaxation strategy, the algorithm remains stable at large Courant numbers without the necessity of solving large, block tri-diagonal systems. Convergence rates and grid refinement studies are conducted for Mach 5 flow through an inlet with a 10 deg compression ramp and Mach 14 flow over a 15 deg ramp. Predictions for pressure distributions, surface heating, and aerodynamics coefficients compare well with experiment data for Mach 10 flow over a blunt body.
NASA Technical Reports Server (NTRS)
Newsome, Richard W.; Walters, Robert W.; Thomas, James L.
1987-01-01
A previously developed upwind/relaxation algorithm for solving the unsteady, compressible, thin-layer Navier-Stokes equations is presently modified so that the downstream influence of the subsonic part of the boundary layer in an otherwise supersonic flow is suppressed by restricting the streamwise pressure gradient. A 'parabolized' solution is then efficiently obtained by marching downstream and iterating locally in each crossflow plane until achieving convergence. This parabolized solution is an excellent final one for problems without large adverse streamwise pressure gradients.
Computational aerothermodynamics
NASA Technical Reports Server (NTRS)
Deiwert, George S.; Green, Michael J.
1987-01-01
Computational aerothermodynamics (CAT) has in the past contributed to the understanding of real-gas flows encountered by hypervelocity reentry vehicles. With advances in computational fluid dynamics, in the modeling of high temperature phenomena, and in computer capability, CAT is an enabling technology for the design of many future space vehicles. An overview of the current capabilities of CAT is provided by describing available methods and their applications. Technical challenges that need to be met are discussed.
Aerothermodynamic methods for a Mars environmental survey Mars entry
NASA Technical Reports Server (NTRS)
Mitcheltree, R. A.
1994-01-01
Computational fluid dynamics models for the thermodynamics and transport properties used in an equilibrium version of the Langley aerothermodynamics upwind relaxation algorithm (LAURA) for Mars atmospheric entries are described. In addition, the physical models used in a nonequilibrium version of LAURA for Mars-entry flows are described. Uncertainties exist in defining constants used in the transport properties for the equilibrium model and in many of the physical models for the nonequilibrium version. Solutions from the two codes using the best available constants are examined at the Mars-entry conditions characteristics of the Mars environmental survey mission. While the flowfields are near thermal equilibrium, chemical nonequilibrium effects are present in the entry cases examined. Convective heating at the stagnation point for these flows (assuming fully catalytic wall boundary conditions) is approximately 100 W/cm(exp 2). Radiative heating is negligible.
Aerothermodynamic Analysis of Commercial Experiment Transporter (COMET) Reentry Capsule
NASA Technical Reports Server (NTRS)
Wood, William A.; Gnoffo, Peter A.; Rault, Didier F. G.
1996-01-01
An aerothermodynamic analysis of the Commercial Experiment Transporter (COMET) reentry capsule has been performed using the laminar thin-layer Navier-Stokes solver Langley Aerothermodynamic Upwind Relaxation Algorithm. Flowfield solutions were obtained at Mach numbers 1.5, 2, 5, 10, 15, 20, 25, and 27.5. Axisymmetric and 5, 10, and 20 degree angles of attack were considered across the Mach-number range, with the Mach 25 conditions taken to 90 degrees angle of attack and the Mach 27.5 cases taken to 60 degrees angle of attack. Detailed surface heat-transfer rates were computed at Mach 20 and 25, revealing that heating rates on the heat-shield shoulder ,can exceed the stagnation-point heating by 230 percent. Finite-rate chemistry solutions were performed above Mach 10, otherwise perfect gas computations were made. Drag, lift, and pitching moment coefficients are computed and details of a wake flow are presented. The effect of including the wake in the solution domain was investigated and base pressure corrections to forebody drag coefficients were numerically determined for the lower Mach numbers. Pitching moment comparisons are made with direct simulation Monte Carlo results in the more rarefied flow at the highest Mach numbers, showing agreement within two-percent. Thin-layer Navier-Stokes computations of the axial force are found to be 15 percent higher across the speed range than the empirical/Newtonian based results used during the initial trajectory analyses.
An upwind-biased space marching algorithm for supersonic viscous flow
NASA Technical Reports Server (NTRS)
Greene, Francis A.
1991-01-01
The modifications are documented which were made to the Langley Aerothermodynamic Upwind Relaxation Algorithm which allow it to compute solutions in a space marching manner. The space marching flux is formulated to be either first- or second-order accurate, using Roe's upwind differencing or Symmetric Total Variation Diminishing differencing, respectively. The algorithm solves the thin layer Navier-Stokes equations, and is subject to the same flow restrictions as a parabolized Navier-Stokes solver. Each cross flow plane is locally iterated in time until converged, then marched in space to the next station. The algorithm is tested on a sphere-cone geometry and a geometry which models the windward side of the Space Shuttle Orbiter. Computational results for surface heating are compared to ground based experimental data. In addition, space marching predictions for surface pressure are compared against values from the original algorithm.
Multi-Component Diffusion with Application To Computational Aerothermodynamics
NASA Technical Reports Server (NTRS)
Sutton, Kenneth; Gnoffo, Peter A.
1998-01-01
The accuracy and complexity of solving multicomponent gaseous diffusion using the detailed multicomponent equations, the Stefan-Maxwell equations, and two commonly used approximate equations have been examined in a two part study. Part I examined the equations in a basic study with specified inputs in which the results are applicable for many applications. Part II addressed the application of the equations in the Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) computational code for high-speed entries in Earth's atmosphere. The results showed that the presented iterative scheme for solving the Stefan-Maxwell equations is an accurate and effective method as compared with solutions of the detailed equations. In general, good accuracy with the approximate equations cannot be guaranteed for a species or all species in a multi-component mixture. 'Corrected' forms of the approximate equations that ensured the diffusion mass fluxes sum to zero, as required, were more accurate than the uncorrected forms. Good accuracy, as compared with the Stefan- Maxwell results, were obtained with the 'corrected' approximate equations in defining the heating rates for the three Earth entries considered in Part II.
NASA Technical Reports Server (NTRS)
1984-01-01
Space shuttle aerothermodynamic data, collected from a continuing series of wind tunnel tests, are permanently stored with the Data Management Services (DMS) system. Information pertaining to current baseline configuration definition is also stored. A list of documentation of DMS processed data arranged sequentially and by space shuttle configuration is presented. The listing provides an up to date record of all applicable aerothermodynamic data collected, processed, or summarized during the space shuttle program. Tables are designed to provide survey information to the various space shuttle managerial and technical levels.
Aerothermodynamic Measurement and Prediction for Modified Orbiter at Mach 6 and 10
NASA Technical Reports Server (NTRS)
Micol, John R.
1995-01-01
Detailed heat-transfer rate distributions measured laterally over the windward surface of an orbiter-like configuration using thin-film resistance heat-transfer gauges and globally using the newly developed relative intensity, two-color thermographic phosphor technique are presented for Mach 6 and 10 in air. The angle of attack was varied from 0 to 40 deg, and the freestream Reynolds number based on the model length was varied from 4 x 10(exp 5) to 6 x 10(exp 6) at Mach 6, corresponding to laminar, transitional, and turbulent boundary layers; the Reynolds number at Mach 10 was 4 x 10(exp 5), corresponding to laminar flow. The primary objective of the present study was to provide detailed benchmark heat-transfer data for the calibration of computational fluid-dynamics codes. Predictions from a Navier-Stokes solver referred to as the Langley aerothermodynamic upwind relaxation algorithm and an approximate boundary-layer solving method known as the axisymmetric analog three-dimensional boundary layer code are compared with measurement. In general, predicted laminar heat-transfer rates are in good agreement with measurements.
HEART Aerothermodynamic Analysis
NASA Technical Reports Server (NTRS)
Mazaheri, Alireza
2012-01-01
This paper presents an assessment of the aerothermodynamic environment around an 8.3 meter High Energy Atmospheric Reentry Test (HEART) vehicle. This study generated twelve nose shape configurations and compared their responses at the peak heating trajectory point against the baseline nose shape. The heat flux sensitivity to the angle of attack variations are also discussed. The possibility of a two-piece Thermal Protection System (TPS) design at the nose is also considered, as are the surface catalytic affects of the aeroheating environment of such configuration. Based on these analyses, an optimum nose shape is proposed to minimize the surface heating. A recommendation is also made for a two-piece TPS design, for which the surface catalytic uncertainty associated with the jump in heating at the nose-IAD juncture is reduced by a minimum of 93%. In this paper, the aeroshell is assumed to be rigid and the inflatable fluid interaction effect is left for future investigations.
X-38 Experimental Aerothermodynamics
NASA Technical Reports Server (NTRS)
Horvath, Thomas J.; Berry, Scott A.; Merski, N. Ronald; Fitzgerald, Steve M.
2000-01-01
The X-38 program seeks to demonstrate an autonomously returned orbital test flight vehicle to support the development of an operational Crew Return Vehicle for the International Space Station. The test flight, anticipated in 2002, is intended to demonstrate the entire mission profile of returning Space Station crew members safely back to earth in the event of medical or mechanical emergency. Integral to the formulation of the X-38 flight data book and the design of the thermal protection system, the aerothermodynamic environment is being defined through a synergistic combination of ground based testing and computational fluid dynamics. This report provides an overview of the hypersonic aerothermodynamic wind tunnel program conducted at the NASA Langley Research Center in support of the X-38 development. Global and discrete surface heat transfer force and moment, surface streamline patterns, and shock shapes were measured on scaled models of the proposed X-38 configuration in different test gases at Mach 6, 10 and 20. The test parametrics include angle of attack from 0 to 50 degs, unit Reynolds numbers from 0.3 x 10 (exp 6) to 16 x 10 (exp 6)/ ft, rudder deflections of 0, 2, and 5 deg. and body flap deflections from 0 to 30 deg. Results from hypersonic aerodynamic screening studies that were conducted as the configuration evolved to the present shape at, presented. Heavy gas simulation tests have indicated that the primary real gas effects on X-38 aerodynamics at trim conditions are expected to favorably influence flap effectiveness. Comparisons of the experimental heating and force and moment data to prediction and the current aerodynamic data book are highlighted. The effects of discrete roughness elements on boundary layer transition were investigated at Mach 6 and the development of a transition correlation for the X-38 vehicle is described. Extrapolation of ground based heating measurements to flight radiation equilibrium wall temperatures at Mach 6 and 10 were
Aerothermodynamic Flight Simulation Capabilities for Aerospace Vehicles
NASA Technical Reports Server (NTRS)
Miller, Charles G.
1998-01-01
Aerothermodynamics, encompassing aerodynamics, aeroheating, and fluid dynamics and physical processes, is the genesis for the design and development of advanced space transportation vehicles and provides crucial information to other disciplines such as structures, materials, propulsion, avionics, and guidance, navigation and control. Sources of aerothermodynamic information are ground-based facilities, Computational Fluid Dynamic (CFD) and engineering computer codes, and flight experiments. Utilization of this aerothermodynamic triad provides the optimum aerothermodynamic design to safely satisfy mission requirements while reducing design conservatism, risk and cost. The iterative aerothermodynamic process for initial screening/assessment of aerospace vehicle concepts, optimization of aerolines to achieve/exceed mission requirements, and benchmark studies for final design and establishment of the flight data book are reviewed. Aerothermodynamic methodology centered on synergism between ground-based testing and CFD predictions is discussed for various flow regimes encountered by a vehicle entering the Earth s atmosphere from low Earth orbit. An overview of the resources/infrastructure required to provide accurate/creditable aerothermodynamic information in a timely manner is presented. Impacts on Langley s aerothermodynamic capabilities due to recent programmatic changes such as Center reorganization, downsizing, outsourcing, industry (as opposed to NASA) led programs, and so forth are discussed. Sample applications of these capabilities to high Agency priority, fast-paced programs such as Reusable Launch Vehicle (RLV)/X-33 Phases I and 11, X-34, Hyper-X and X-38 are presented and lessons learned discussed. Lastly, enhancements in ground-based testing/CFD capabilities necessary to partially/fully satisfy future requirements are addressed.
Aerothermodynamic systems engineering and design
NASA Astrophysics Data System (ADS)
A reference source for various aspects of aerothermodynamic systems engineering and design is presented. Air conditioning load analysis is addressed, including physiological requirements, heat and cooling load equations, skin temperature computational methods, cooling loads due to radiation through transparent areas, heating and cooling loads due to internal sources, and practical considerations in the determination of overall heating and cooling loads. Refrigeration system design is considered, including air cycle systems, vapor cycle systems, combined vapor cycle and air cycle systems, and thermoelectric cooling. Heating methods is heating system design and low pressure and high pressure systems in air distribution system design are addressed. Procedures and equations commonly used for aerospace applications of these technologies are included.
Aerothermodynamics of manned Mars missions
NASA Technical Reports Server (NTRS)
Park, Chul; Davies, Carol B.
1989-01-01
The aerothermodynamic problems associated with the aerobraking of the spacecraft proposed for the manned Mars mission are studied. The propulsive Delta V necessary at departure from earth and Mars and the velocities of the atmospheric entries into the two planets are deduced. It is shown that the propulsive Delta V can be reduced by increasing the entry velocities and that entry velocities up to about 15 km/sec are appropriate at both earth and Mars. L/D values of 0.8 and 2.0 are found to be necessary at earth and Mars, respectively. Density, pressure, and stagnation-point convective-heat-transfer rates are calculated for the typical aerobraking flights. Assuming the shock layer flow to be in equilibrium, the stagnation-point radiative-heat-transfer rates are calculated to be larger than the convective-heat-transfer rates. The possible impact of ablation, turbulence, and nonequilibrium are discussed.
Opportunities for research in aerothermodynamics
NASA Technical Reports Server (NTRS)
Graham, R. W.
1983-01-01
"Aerothermodynamics' involves the disciplines of chemistry, thermodynamics, fluid mechanics and heat transfer which have collaborative importance in propulsion systems. There are growing opportunities for the further application of these disciplines to improve the methodology for the design of advanced gas turbines; particularly, the combustor and turbine. Design procedures follow empirical or cut and try guidelines. The tremendous advances in computational analysis and in instrumentation techniques hold promise for research answers to complex physical processes that are currently not well understood. The transfer of basic research understanding to engineering design should result in shorter, less expensive development commitments for engines. The status and anticipated opportunities in research topics relevant to combustors and turbines is reviewed.
Aerothermodynamics Overview and Prediction Assessment
NASA Technical Reports Server (NTRS)
Heidmann, James D.
2007-01-01
An overview of the Aerothermodynamics Discipline within NASA s Subsonic Fixed Wing Project is given. The primary focus of the presentation is on the research efforts conducted in fiscal year 2007. This year (2007), the work primarily consisted of efforts under level 1 (foundational research) and level 2 (tools and technology development). Examples of work under level 1 are large eddy simulation development, advanced turbine cooling concept development, and turbomachinery flow control development. Examples of level 2 research are the development of highly-loaded compressor and turbine test programs and advanced turbomachinery simulation development, including coupled inlet-fan simulations. An overview of the NRA research activity is also provided. This NRA focused on plasma and aspiration flow control for low pressure turbine application. Finally, a status report on the turbomachinery CFD code assessment activity is provided. This activity focuses on the use of several NASA in-house codes for the NASA rotor 37 and stage 35 test cases.
A Perspective on Computational Aerothermodynamics at NASA
NASA Technical Reports Server (NTRS)
Gnoffo, Peter A.
2007-01-01
The evolving role of computational aerothermodynamics (CA) within NASA over the past 20 years is reviewed. The presentation highlights contributions to understanding the Space Shuttle pitching moment anomaly observed in the first shuttle flight, prediction of a static instability for Mars Pathfinder, and the use of CA for damage assessment in post-Columbia mission support. In the view forward, several current challenges in computational fluid dynamics and aerothermodynamics for hypersonic vehicle applications are discussed. Example simulations are presented to illustrate capabilities and limitations. Opportunities to advance the state-of-art in algorithms, grid generation and adaptation, and code validation are identified.
Overview of Aerothermodynamic Loads Definition Study
NASA Technical Reports Server (NTRS)
Povinelli, L. A.
1985-01-01
The Aerothermodynamic Loads Definition were studied to develop methods to more accurately predict the operating environment in the space shuttle main engine (SSME) components. Development of steady and time-dependent, three-dimensional viscous computer codes and experimental verification and engine diagnostic testing are considered. The steady, nonsteady, and transient operating loads are defined to accurately predict powerhead life. Improvements in the structural durability of the SSME turbine drive systems depends on the knowledge of the aerothermodynamic behavior of the flow through the preburner, turbine, turnaround duct, gas manifold, and injector post regions.
Space Shuttle aerothermodynamic data report, phase C
NASA Technical Reports Server (NTRS)
1985-01-01
Space shuttle aerothermodynamic data, collected from a continuing series of wind tunnel tests, are permanently stored with the Data Management Services (DMS) system. Information pertaining to current baseline configuration definition is also stored. Documentation of DMS processed data arranged sequentially and by space shuttle configuration are included. An up-to-date record of all applicable aerothermodynamic data collected, processed, or summarized during the space shuttle program is provided. Tables are designed to provide suvery information to the various space shuttle managerial and technical levels.
Aerothermodynamic Analyses of Towed Ballutes
NASA Technical Reports Server (NTRS)
Gnoffo, Peter A.; Buck, Greg; Moss, James N.; Nielsen, Eric; Berger, Karen; Jones, William T.; Rudavsky, Rena
2006-01-01
A ballute (balloon-parachute) is an inflatable, aerodynamic drag device for application to planetary entry vehicles. Two challenging aspects of aerothermal simulation of towed ballutes are considered. The first challenge, simulation of a complete system including inflatable tethers and a trailing toroidal ballute, is addressed using the unstructured-grid, Navier-Stokes solver FUN3D. Auxiliary simulations of a semi-infinite cylinder using the rarefied flow, Direct Simulation Monte Carlo solver, DSV2, provide additional insight into limiting behavior of the aerothermal environment around tethers directly exposed to the free stream. Simulations reveal pressures higher than stagnation and corresponding large heating rates on the tether as it emerges from the spacecraft base flow and passes through the spacecraft bow shock. The footprint of the tether shock on the toroidal ballute is also subject to heating amplification. Design options to accommodate or reduce these environments are discussed. The second challenge addresses time-accurate simulation to detect the onset of unsteady flow interactions as a function of geometry and Reynolds number. Video of unsteady interactions measured in the Langley Aerothermodynamic Laboratory 20-Inch Mach 6 Air Tunnel and CFD simulations using the structured grid, Navier-Stokes solver LAURA are compared for flow over a rigid spacecraft-sting-toroid system. The experimental data provides qualitative information on the amplitude and onset of unsteady motion which is captured in the numerical simulations. The presence of severe unsteady fluid - structure interactions is undesirable and numerical simulation must be able to predict the onset of such motion.
Aerothermodynamics at NASA-Langley Research Center
NASA Technical Reports Server (NTRS)
Weilmuenster, K. James
2001-01-01
The Aerothermodynamics Branch at NASA - Langley Research Center is tasked with developing, assessing and applying aerothermodynamic technologies to enable the development of hypersonic aircraft, launch vehicles, and planetary/earth entry systems. To accomplish this mission, the Branch capitalizes on the synergism between the experimental and computational facilities/tools which reside in the branch and a staff that can draw on five decades of experience in aerothermodynamics. The Aerothermodynamics Branch is staffed by 30 scientists/engineers. The staff, of which two-thirds are less than 40 years old, is split evenly between experimentalists and computationalists. Approximately 90 percent of the staff work on space transportation systems while the remainder work on planetary missions. The Branch manages 5 hypersonic wind tunnels which are staffed by 14 technicians, numerous high end work stations and a SGI Origin 2000 system. The Branch also utilizes other test facilities located at Langley as well as other national and international test sites. Large scale computational requirements are met by access to Agency resources.
Preliminary aerothermodynamic design method for hypersonic vehicles
NASA Technical Reports Server (NTRS)
Harloff, G. J.; Petrie, S. L.
1987-01-01
Preliminary design methods are presented for vehicle aerothermodynamics. Predictions are made for Shuttle orbiter, a Mach 6 transport vehicle and a high-speed missile configuration. Rapid and accurate methods are discussed for obtaining aerodynamic coefficients and heat transfer rates for laminar and turbulent flows for vehicles at high angles of attack and hypersonic Mach numbers.
Computational Aerothermodynamic Design Issues for Hypersonic Vehicles
NASA Technical Reports Server (NTRS)
Gnoffo, Peter A.; Weilmuenster, K. James; Hamilton, H. Harris, II; Olynick, David R.; Venkatapathy, Ethiraj
2005-01-01
A brief review of the evolutionary progress in computational aerothermodynamics is presented. The current status of computational aerothermodynamics is then discussed, with emphasis on its capabilities and limitations for contributions to the design process of hypersonic vehicles. Some topics to be highlighted include: (1) aerodynamic coefficient predictions with emphasis on high temperature gas effects; (2) surface heating and temperature predictions for thermal protection system (TPS) design in a high temperature, thermochemical nonequilibrium environment; (3) methods for extracting and extending computational fluid dynamic (CFD) solutions for efficient utilization by all members of a multidisciplinary design team; (4) physical models; (5) validation process and error estimation; and (6) gridding and solution generation strategies. Recent experiences in the design of X-33 will be featured. Computational aerothermodynamic contributions to Mars Path finder, METEOR, and Stardust (Comet Sample return) will also provide context for this discussion. Some of the barriers that currently limit computational aerothermodynamics to a predominantly reactive mode in the design process will also be discussed, with the goal of providing focus for future research.
Computational Aerothermodynamic Design Issues for Hypersonic Vehicles
NASA Technical Reports Server (NTRS)
Gnoffo, Peter A.; Weilmuenster, K. James; Hamilton, H. Harris, II; Olynick, David R.; Venkatapathy, Ethiraj
1997-01-01
A brief review of the evolutionary progress in computational aerothermodynamics is presented. The current status of computational aerothermodynamics is then discussed, with emphasis on its capabilities and limitations for contributions to the design process of hypersonic vehicles. Some topics to be highlighted include: (1) aerodynamic coefficient predictions with emphasis on high temperature gas effects; (2) surface heating and temperature predictions for thermal protection system (TPS) design in a high temperature, thermochemical nonequilibrium environment; (3) methods for extracting and extending computational fluid dynamic (CFD) solutions for efficient utilization by all members of a multidisciplinary design team; (4) physical models; (5) validation process and error estimation; and (6) gridding and solution generation strategies. Recent experiences in the design of X-33 will be featured. Computational aerothermodynamic contributions to Mars Pathfinder, METEOR, and Stardust (Comet Sample return) will also provide context for this discussion. Some of the barriers that currently limit computational aerothermodynamics to a predominantly reactive mode in the design process will also be discussed, with the goal of providing focus for future research.
Computational Aerothermodynamic Design Issues for Hypersonic Vehicles
NASA Technical Reports Server (NTRS)
Olynick, David R.; Venkatapathy, Ethiraj
2004-01-01
A brief review of the evolutionary progress in computational aerothermodynamics is presented. The current status of computational aerothermodynamics is then discussed, with emphasis on its capabilities and limitations for contributions to the design process of hypersonic vehicles. Some topics to be highlighted include: (1) aerodynamic coefficient predictions with emphasis on high temperature gas effects; (2) surface heating and temperature predictions for thermal protection system (TPS) design in a high temperature, thermochemical nonequilibrium environment; (3) methods for extracting and extending computational fluid dynamic (CFD) solutions for efficient utilization by all members of a multidisciplinary design team; (4) physical models; (5) validation process and error estimation; and (6) gridding and solution generation strategies. Recent experiences in the design of X-33 will be featured. Computational aerothermodynamic contributions to Mars Pathfinder, METEOR, and Stardust (Comet Sample return) will also provide context for this discussion. Some of the barriers that currently limit computational aerothermodynamics to a predominantly reactive mode in the design process will also be discussed, with the goal of providing focus for future research.
Development of upwind schemes for the Euler equations
NASA Technical Reports Server (NTRS)
Chakravarthy, Sukumar R.
1987-01-01
Described are many algorithmic and computational aspects of upwind schemes and their second-order accurate formulations based on Total-Variation-Diminishing (TVD) approaches. An operational unification of the underlying first-order scheme is first presented encompassing Godunov's, Roe's, Osher's, and Split-Flux methods. For higher order versions, the preprocessing and postprocessing approaches to constructing TVD discretizations are considered. TVD formulations can be used to construct relaxation methods for unfactored implicit upwind schemes, which in turn can be exploited to construct space-marching procedures for even the unsteady Euler equations. A major part of the report describes time- and space-marching procedures for solving the Euler equations in 2-D, 3-D, Cartesian, and curvilinear coordinates. Along with many illustrative examples, several results of efficient computations on 3-D supersonic flows with subsonic pockets are presented.
Overview of aerothermodynamic loads definition study
NASA Technical Reports Server (NTRS)
Gaugler, Raymond E.
1989-01-01
Over the years, NASA has been conducting the Advanced Earth-to-Orbit (AETO) Propulsion Technology Program to provide the knowledge, understanding, and design methodology that will allow the development of advanced Earth-to-orbit propulsion systems with high performance, extended service life, automated operations, and diagnostics for in-flight health monitoring. The objective of the Aerothermodynamic Loads Definition Study is to develop methods to more accurately predict the operating environment in AETO propulsion systems, such as the Space Shuttle Main Engine (SSME) powerhead. The approach taken consists of 2 parts: to modify, apply, and disseminate existing computational fluid dynamics tools in response to current needs and to develop new technology that will enable more accurate computation of the time averaged and unsteady aerothermodynamic loads in the SSME powerhead. The software tools are detailed. Significant progress was made in the area of turbomachinery, where there is an overlap between the AETO efforts and research in the aeronautical gas turbine field.
Aerothermodynamic Insight From The HIFIRE Program
NASA Astrophysics Data System (ADS)
Kimmel, Roger L.; Adamczak, David; Dolvin, Douglas; Borg, Matthew; Stanfield, Scott
2011-05-01
The HIFiRE (Hypersonic International Flight Research and Experimentation) program is a joint venture of the United States Air Force Research Laboratory and Australian Defence Science and Technology Organisation to utilize economical flight research opportunities in the exploration of flight science issues for space access systems. Flights 1 and 5 focus on collecting high-resolution experimental data on critical aerothermodynamic phenomena, including laminar-turbulent transition and shock/boundary layer interactions. Flight 1, successfully flown in March 2010, employed a test article composed of a 7-deg right angle cone, followed by a cylinder and flare. The test article remained attached to the second-stage booster throughout the ballistic trajectory. Flight 5, to be launched in a similar fashion, will feature a 2:1 elliptic cross-section cone as the test article. For both flights significant resources have been invested in pre-flight aerothermodynamic analysis and testing. This manuscript will summarize the overall strategy of the HIFiRE program, review the pre-flight aerothermodynamic analysis for Flights 1 and 5, and present a brief look at preliminary results from the post-flight analysis of Flight 1.
Computational Aerothermodynamics in Aeroassist Applications
NASA Technical Reports Server (NTRS)
Gnoffo, Peter A.
2001-01-01
Aeroassisted planetary entry uses atmospheric drag to decelerate spacecraft from super-orbital to orbital or suborbital velocities. Numerical simulation of flow fields surrounding these spacecraft during hypersonic atmospheric entry is required to define aerothermal loads. The severe compression in the shock layer in front of the vehicle and subsequent, rapid expansion into the wake are characterized by high temperature, thermo-chemical nonequilibrium processes. Implicit algorithms required for efficient, stable computation of the governing equations involving disparate time scales of convection, diffusion, chemical reactions, and thermal relaxation are discussed. Robust point-implicit strategies are utilized in the initialization phase; less robust but more efficient line-implicit strategies are applied in the endgame. Applications to ballutes (balloon-like decelerators) in the atmospheres of Venus, Mars, Titan, Saturn, and Neptune and a Mars Sample Return Orbiter (MSRO) are featured. Examples are discussed where time-accurate simulation is required to achieve a steady-state solution.
Upwind dynamic soaring of albatrosses and UAVs
NASA Astrophysics Data System (ADS)
Richardson, Philip L.
2015-01-01
Albatrosses have been observed to soar in an upwind direction using what is called here an upwind mode of dynamic soaring. The upwind mode was modeled using the dynamics of a two-layer Rayleigh cycle in which the lower layer has zero velocity and the upper layer has a uniform wind speed of W. The upwind mode consists of a climb across the wind-shear layer headed upwind, a 90° turn and descent across the wind-shear layer perpendicular to the wind, followed by a 90° turn into the wind. The increase of airspeed gained from crossing the wind-shear layer headed upwind was balanced by the decrease of airspeed caused by drag. Results show that a wandering albatross can soar over the ocean in an upwind direction at a mean speed of 8.4 m/s in a 3.6 m/s wind, which is the minimum wind speed necessary for sustained dynamic soaring. A main result is that albatrosses can soar upwind much faster than the wind speed. Furthermore, albatrosses were found to be able to increase upwind speeds in winds greater than 3.6 m/s, reaching an upwind speed of 12.1 m/s in a wind speed of 7 m/s (for example). The upwind dynamic soaring mode of a possible robotic albatross UAV (Unmanned Aerial Vehicle) was modeled using a Rayleigh cycle and characteristics of a high-performance glider. Maximum possible airspeeds are equal to approximately 9.5 times the wind speed of the upper layer. In a wind of 10 m/s, the maximum possible upwind (56 m/s) and across-wind (61 m/s) components of UAV velocity over the ocean result in a diagonal upwind velocity of 83 m/s. In sufficient wind, a UAV could, in principle, use fast diagonal speeds to rapidly survey large areas of the ocean surface and the marine boundary layer. In practice, the maximum speeds of a UAV soaring over the ocean could be significantly less than these predictions. Some limitations to achieving fast travel velocities over the ocean are discussed and suggestions are made for further studies to test the concept of a robotic albatross.
Atmospheric Entry Aerothermodynamics Flight Test on CubeSat Platform
NASA Astrophysics Data System (ADS)
Sakraker, I.; Umit, E.; van der Haegen, V.; Chazot, O.
2014-06-01
The challenging aerothermochemistry of atmospheric entry is aimed to be experimented on a triple CubeSat platform having ablative TPS in the front unit and ceramic TPS on the side panels. Five aerothermodynamics payloads are presented in this paper.
Aerothermodynamic data base. Data file contents report, phase C
NASA Technical Reports Server (NTRS)
Lutz, G. R.
1983-01-01
Space shuttle aerothermodynamic data, collected from a continuing series of wind tunnel tests, are permanently stored with the Data Management Services (DMS) system. Information pertaining to current baseline configuration definition is also stored. Documentation of DMS processed data arranged sequentially and by space shuttle configuration is listed to provide an up-to-date record of all applicable aerothermodynamic data collected, processed, or summarized during the space shuttle program. Tables provide survey information to the various space shuttle managerial and technical levels.
Shuttle Tethered Aerothermodynamics Research Facilty (STARFAC) instrumentation requirements
NASA Technical Reports Server (NTRS)
Wood, G. M.; Siemers, P. M.; Carlomagno, G. M.; Hoffman, J.
1986-01-01
The instrumentation requirements for the Shuttle Tethered Aerothermodynamic Research Facility (STARFAC) are presented. The typical physical properties of the terrestrial atmosphere are given along with representative atmospheric daytime ion concentrations and the equilibrium and nonequilibrium gas property comparison from a point away from a wall. STARFAC science and engineering measurements are given as are the TSS free stream gas analysis. The potential nonintrusive measurement techniques for hypersonic boundary layer research are outlined along with the quantitative physical measurement methods for aerothermodynamic studies.
Aerothermodynamics of the Mars Global Surveyor Spacecraft
NASA Technical Reports Server (NTRS)
Shane, Russell W.; Tolson, Robert H.
1998-01-01
The aerothermodynamics characteristics of the Mars Global Surveyor spacecraft are investigated and reported. These results have been used by the Mars Global Surveyor mission planners to design the aerobraking phase of the mission. Analytical and Direct Simulation Monte Carlo computer codes were used with a detailed, three dimensional model of the spacecraft to evaluate spacecraft aerobraking characteristics for flight in free molecular and transitional flow regimes. The spacecraft is found to be aerodynamically stable in aerobraking and planned contingency configurations. Aerodynamic forces, moments, and heating are found to be highly dependent on atmospheric density. Accommodation coefficient. is seen to strongly influence drag coefficient. Transitional flow effects are found to reduce overall solar panel heating. Attitude control thruster plumes are shown to interact with the freestream, diminishing the effectiveness of the attitude control system and even leading to thrust reversal. These plume-freestream interaction effects are found to be highly dependent on freestream density.
Shuttle entry technology payloads. [for aerothermodynamic research
NASA Technical Reports Server (NTRS)
Siemers, P. M., III
1975-01-01
The flight frequency of the Space Transportation System (STS) coupled with its large payload-carrying capability will provide an unprecedented opportunity for conducting aerothermodynamic/entry technology research. This STS research opportunity can be characterized into two distinct categories: (1) that research which will utilize the STS orbiter as the test vehicle, and (2) that research which will utilize a vehicle launched from the orbiter for entry. To date, on-going studies have defined experiments as well as the support systems required for the shuttle launched research program. The proposed Entry Technology Program will provide a flight data base from which accurate correlations can be performed relative to ground test and analysis data. These correlations will result in optimized designs for future flight systems.
Overview of aerothermodynamic loads definition study
NASA Technical Reports Server (NTRS)
Gaugler, Raymond E.
1991-01-01
The objective of the Aerothermodynamic Loads Definition Study is to develop methods of accurately predicting the operating environment in advanced Earth-to-Orbit (ETO) propulsion systems, such as the Space Shuttle Main Engine (SSME) powerhead. Development of time averaged and time dependent three dimensional viscous computer codes as well as experimental verification and engine diagnostic testing are considered to be essential in achieving that objective. Time-averaged, nonsteady, and transient operating loads must all be well defined in order to accurately predict powerhead life. Described here is work in unsteady heat flow analysis, improved modeling of preburner flow, turbulence modeling for turbomachinery, computation of three dimensional flow with heat transfer, and unsteady viscous multi-blade row turbine analysis.
On symmetric and upwind TVD schemes
NASA Technical Reports Server (NTRS)
Yee, H. C.
1985-01-01
A class of explicit and implicit total variation diminishing (TVD) schemes for the compressible Euler and Navier-Stokes equations was developed. They do not generate spurious oscillations across shocks and contact discontinuities. In general, shocks can be captured within 1 to 2 grid points. For the inviscid case, these schemes are divided into upwind TVD schemes and symmetric (nonupwind) TVD schemes. The upwind TVD scheme is based on the second-order TVD scheme. The symmetric TVD scheme is a generalization of Roe's and Davis' TVD Lax-Wendroff scheme. The performance of these schemes on some viscous and inviscid airfoil steady-state calculations is investigated. The symmetric and upwind TVD schemes are compared.
Experimental aerothermodynamic research of hypersonic aircraft
NASA Technical Reports Server (NTRS)
Cleary, Joseph W.
1990-01-01
Wind tunnel tests were conducted to establish a benchmark experimental data base for a genetic hypersonic aircraft vehicle. Comprehensive measurements were made at Mach 7 to give flow visualization, surface pressure, surface convective heat transfer, and flow field Pitot pressure for a delta platform all-body vehicle. The tests were conducted in the NASA/Ames 3.5-Foot Hypersonic Wind Tunnel at Reynolds numbers sufficient to give turbulent flow. Comparisons are made of the experimental results with computational solutions of the flow by an upwind parabolized Navier-Stokes code developed at Ames. Good agreement of experiment with solutions by the code is demonstrated.
Uncertainty Assessment of Hypersonic Aerothermodynamics Prediction Capability
NASA Technical Reports Server (NTRS)
Bose, Deepak; Brown, James L.; Prabhu, Dinesh K.; Gnoffo, Peter; Johnston, Christopher O.; Hollis, Brian
2011-01-01
The present paper provides the background of a focused effort to assess uncertainties in predictions of heat flux and pressure in hypersonic flight (airbreathing or atmospheric entry) using state-of-the-art aerothermodynamics codes. The assessment is performed for four mission relevant problems: (1) shock turbulent boundary layer interaction on a compression corner, (2) shock turbulent boundary layer interaction due a impinging shock, (3) high-mass Mars entry and aerocapture, and (4) high speed return to Earth. A validation based uncertainty assessment approach with reliance on subject matter expertise is used. A code verification exercise with code-to-code comparisons and comparisons against well established correlations is also included in this effort. A thorough review of the literature in search of validation experiments is performed, which identified a scarcity of ground based validation experiments at hypersonic conditions. In particular, a shortage of useable experimental data at flight like enthalpies and Reynolds numbers is found. The uncertainty was quantified using metrics that measured discrepancy between model predictions and experimental data. The discrepancy data is statistically analyzed and investigated for physics based trends in order to define a meaningful quantified uncertainty. The detailed uncertainty assessment of each mission relevant problem is found in the four companion papers.
Aero-Thermo-Dynamic Mass Analysis
Shiba, Kota; Yoshikawa, Genki
2016-01-01
Each gas molecule has its own molecular weight, while such a microscopic characteristic is generally inaccessible, and thus, it is measured indirectly through e.g. ionization in conventional mass analysis. Here, we present a novel approach to the direct measurement of molecular weight through a nanoarchitectonic combination of aerodynamics, thermodynamics, and mechanics, transducing microscopic events into macroscopic phenomena. It is confirmed that this approach can provide molecular weight of virtually any gas or vaporizable liquid sample in real-time without ionization. Demonstrations through analytical calculations, numerical simulations, and experiments verify the validity and versatility of the novel mass analysis realized by a simple setup with a flexible object (e.g. with a bare cantilever and even with a business card) placed in a laminar jet. Owing to its unique and simple working principle, this aero-thermo-dynamic mass analysis (AMA) can be integrated into various analytical devices, production lines, and consumer mobile platforms, opening new chapters in aerodynamics, thermodynamics, mechanics, and mass analysis. PMID:27412335
Aeroassist flight experiment aerodynamics and aerothermodynamics
NASA Technical Reports Server (NTRS)
Brewer, Edwin B.
1989-01-01
The problem is to determine the transitional flow aerodynamics and aerothermodynamics, including the base flow characteristics, of the Aeroassist Flight Experiment (AFE). The justification for the computational fluid dynamic (CFD) Application stems from MSFC's system integration responsibility for the AFE. To insure that the AFE objectives are met, MSFC must understand the limitations and uncertainties of the design data. Perhaps the only method capable of handling the complex physics of the rarefied high energy AFE trajectory is Bird's Direct Simulation Monte Carlo (DSMC) technique. The 3-D code used in this analysis is applicable only to the AFE geometry. It uses the Variable Hard Sphere (VHS) collision model and five specie chemistry model available from Langley Research Center. The code is benchmarked against the AFE flight data and used as an Aeroassisted Space Transfer Vehicle (ASTV) design tool. The code is being used to understand the AFE flow field and verify or modify existing design data. Continued application to lower altitudes is testing the capability of the Numerical Aerodynamic Simulation Facility (NASF) to handle 3-D DSMC and its practicality as an ASTV/AFE design tool.
Aero-Thermo-Dynamic Mass Analysis
NASA Astrophysics Data System (ADS)
Shiba, Kota; Yoshikawa, Genki
2016-07-01
Each gas molecule has its own molecular weight, while such a microscopic characteristic is generally inaccessible, and thus, it is measured indirectly through e.g. ionization in conventional mass analysis. Here, we present a novel approach to the direct measurement of molecular weight through a nanoarchitectonic combination of aerodynamics, thermodynamics, and mechanics, transducing microscopic events into macroscopic phenomena. It is confirmed that this approach can provide molecular weight of virtually any gas or vaporizable liquid sample in real-time without ionization. Demonstrations through analytical calculations, numerical simulations, and experiments verify the validity and versatility of the novel mass analysis realized by a simple setup with a flexible object (e.g. with a bare cantilever and even with a business card) placed in a laminar jet. Owing to its unique and simple working principle, this aero-thermo-dynamic mass analysis (AMA) can be integrated into various analytical devices, production lines, and consumer mobile platforms, opening new chapters in aerodynamics, thermodynamics, mechanics, and mass analysis.
Aero-Thermo-Dynamic Mass Analysis.
Shiba, Kota; Yoshikawa, Genki
2016-01-01
Each gas molecule has its own molecular weight, while such a microscopic characteristic is generally inaccessible, and thus, it is measured indirectly through e.g. ionization in conventional mass analysis. Here, we present a novel approach to the direct measurement of molecular weight through a nanoarchitectonic combination of aerodynamics, thermodynamics, and mechanics, transducing microscopic events into macroscopic phenomena. It is confirmed that this approach can provide molecular weight of virtually any gas or vaporizable liquid sample in real-time without ionization. Demonstrations through analytical calculations, numerical simulations, and experiments verify the validity and versatility of the novel mass analysis realized by a simple setup with a flexible object (e.g. with a bare cantilever and even with a business card) placed in a laminar jet. Owing to its unique and simple working principle, this aero-thermo-dynamic mass analysis (AMA) can be integrated into various analytical devices, production lines, and consumer mobile platforms, opening new chapters in aerodynamics, thermodynamics, mechanics, and mass analysis. PMID:27412335
An Overview of the Space Shuttle Aerothermodynamic Design
NASA Technical Reports Server (NTRS)
Martin, Fred
2011-01-01
The Space Shuttle Thermal Protection System was one of the three areas that required the development of new technology. The talk discusses the pre-flight development of the aerothermodynamic environment which was based on Mach 8 wind tunnel data. A high level overview of the pre-flight heating rate predictions and comparison to the Orbiter Flight Test (OFT) data is presented, along with a discussion of the dramatic improvement in the state-of-the-art in aerothermodynamic capability that has been used to support the Shuttle Program. A high level review of the Orbiter aerothermodynamic design is discussed, along with improvements in Computational Fluid Dynamics and wind tunnel testing that was required for flight support during the last 30 years. The units have been removed from the plots, and the discussion is kept at a high level.
NASA's hypersonic fluid and thermal physics program (Aerothermodynamics)
NASA Technical Reports Server (NTRS)
Graves, R. A.; Hunt, J. L.
1985-01-01
This survey paper gives an overview of NASA's hypersonic fluid and thermal physics program (recently renamed aerothermodynamics). The purpose is to present the elements of, example results from, and rationale and projection for this program. The program is based on improving the fundamental understanding of aerodynamic and aerothermodynamic flow phenomena over hypersonic vehicles in the continuum, transitional, and rarefied flow regimes. Vehicle design capabilities, computational fluid dynamics, computational chemistry, turbulence modeling, aerothermal loads, orbiter flight data analysis, orbiter experiments, laser photodiagnostics, and facilities are discussed.
Prediction of the Aerothermodynamic Environment of the Huygens Probe
NASA Technical Reports Server (NTRS)
Hollis, Brian R.; Striepe, Scott A.; Wright, Michael J.; Bose, Deepak; Sutton, Kenneth; Takashima, Naruhisa
2005-01-01
An investigation of the aerothermodynamic environment of the Huygens entry probe has been conducted. A Monte Carlo simulation of the trajectory of the probe during entry into Titan's atmosphere was performed to identify a worst-case heating rate trajectory. Flowfield and radiation transport computations were performed at points along this trajectory to obtain convective and radiative heat-transfer distributions on the probe's heat shield. This investigation identified important physical and numerical factors, including atmospheric CH4 concentration, transition to turbulence, numerical diffusion modeling, and radiation modeling, which strongly influenced the aerothermodynamic environment.
Towards an "All Speed" Unstructured Upwind Scheme
NASA Technical Reports Server (NTRS)
Loh, Ching Y.; Jorgenson, Philip C.E.
2009-01-01
In the authors previous studies [1], a time-accurate, upwind finite volume method (ETAU scheme) for computing compressible flows on unstructured grids was proposed. The scheme is second order accurate in space and time and yields high resolution in the presence of discontinuities. The scheme features a multidimensional limiter and multidimensional numerical dissipation. These help to stabilize the numerical process and to overcome the annoying pathological behaviors of upwind schemes. In the present paper, it will be further shown that such multidimensional treatments also lead to a nearly all-speed or Mach number insensitive upwind scheme. For flows at very high Mach number, e.g., 10, local numerical instabilities or the pathological behaviors are suppressed, while for flows at very low Mach number, e.g., 0.02, computation can be directly carried out without invoking preconditioning. For flows in different Mach number regimes, i.e., low, medium, and high Mach numbers, one only needs to adjust one or two parameters in the scheme. Several examples with low and high Mach numbers are demonstrated in this paper. Thus, the ETAU scheme is applicable to a broad spectrum of flow regimes ranging from high supersonic to low subsonic, appropriate for both CFD (computational fluid dynamics) and CAA (computational aeroacoustics).
Accurate upwind methods for the Euler equations
NASA Technical Reports Server (NTRS)
Huynh, Hung T.
1993-01-01
A new class of piecewise linear methods for the numerical solution of the one-dimensional Euler equations of gas dynamics is presented. These methods are uniformly second-order accurate, and can be considered as extensions of Godunov's scheme. With an appropriate definition of monotonicity preservation for the case of linear convection, it can be shown that they preserve monotonicity. Similar to Van Leer's MUSCL scheme, they consist of two key steps: a reconstruction step followed by an upwind step. For the reconstruction step, a monotonicity constraint that preserves uniform second-order accuracy is introduced. Computational efficiency is enhanced by devising a criterion that detects the 'smooth' part of the data where the constraint is redundant. The concept and coding of the constraint are simplified by the use of the median function. A slope steepening technique, which has no effect at smooth regions and can resolve a contact discontinuity in four cells, is described. As for the upwind step, existing and new methods are applied in a manner slightly different from those in the literature. These methods are derived by approximating the Euler equations via linearization and diagonalization. At a 'smooth' interface, Harten, Lax, and Van Leer's one intermediate state model is employed. A modification for this model that can resolve contact discontinuities is presented. Near a discontinuity, either this modified model or a more accurate one, namely, Roe's flux-difference splitting. is used. The current presentation of Roe's method, via the conceptually simple flux-vector splitting, not only establishes a connection between the two splittings, but also leads to an admissibility correction with no conditional statement, and an efficient approximation to Osher's approximate Riemann solver. These reconstruction and upwind steps result in schemes that are uniformly second-order accurate and economical at smooth regions, and yield high resolution at discontinuities.
The definition of the Shuttle Tethered Aerothermodynamic Research Facility
NASA Technical Reports Server (NTRS)
Siemers, P. M., III; Wood, G. M., Jr.; Wolf, H.; Flanagan, P. F.; Henry, M. W.
1985-01-01
Studies have been conducted to define the feasibility and practical limitations of the Shuttle Orbiter Tethered 'wind-tunnel' concept. This concept, referred to as the Shuttle Tethered Aerothermodynamic Research Facility (STARFAC), is proposed to provide researchers access to altitudes above 90 km to accomplish aerothermodynamic research in the rarefied upper atmosphere. Determining the feasibility and limitations of the concept has required the enhancement and/or development of mission simulation analytical techniques and control laws; the accomplishment of candidate mission simulations; the definition of instrumentation requirements, both for science and engineering; and the establishment of tether and satellite design requirements to meet STARFAC objectives. The results of the study, to date, indicate that such a concept is both feasible and practical. Representative results are presented, as are recommendations for continued studies which would result in program implementation.
The tethered satellite system for low density aerothermodynamics studies
NASA Technical Reports Server (NTRS)
Carlomagno, Giovanni M.; De Luca, Luigi; Siemers, P. M., III; Wood, George M., Jr.
1986-01-01
The feasibility of the operation of the Tethered Satellite System (TSS) as a continuous open wind tunnel for low-density aerothermodynamic studies (applicable to the design of hypersonic space vehicles including STARFAC, AOTV, and ERV) is considered. The Shuttle Continuous Open Wind Tunnel (SCOWT) program, for the study of the energy and momentum transfer between the tethered satellite and its environmental medium during the TSS/2 mission, is described. Instrumentation and TSS design requirements to meet SCOWT objectives are also considered. SCOWT will provide information on the gasdynamic processes occurring downstream of the bow wave standing in front of the TS, the chemistry and physics of the upper atmosphere related to satellite aerothermodynamics, and TSS's overall experimental envelope of operation.
Iterated upwind schemes for gas dynamics
Smolarkiewicz, Piotr K. Szmelter, Joanna
2009-01-10
A class of high-resolution schemes established in integration of anelastic equations is extended to fully compressible flows, and documented for unsteady (and steady) problems through a span of Mach numbers from zero to supersonic. The schemes stem from iterated upwind technology of the multidimensional positive definite advection transport algorithm (MPDATA). The derived algorithms employ standard and modified forms of the equations of gas dynamics for conservation of mass, momentum and either total or internal energy as well as potential temperature. Numerical examples from elementary wave propagation, through computational aerodynamics benchmarks, to atmospheric small- and large-amplitude acoustics with intricate wave-flow interactions verify the approach for both structured and unstructured meshes, and demonstrate its flexibility and robustness.
Downward-deployed tethered platforms for high enthalpy aerothermodynamic research
NASA Technical Reports Server (NTRS)
Wood, George M.; Siemers, Paul M.; Squires, R. Kenneth; Wolf, Henry; Carlomagno, Giovanni M.
1988-01-01
The data on aerothermodynamic and aerodynamic interactions at altitudes above 50 km is extremely limited because of the relative inaccessibility of the region to research vehicles of any sort. This paper addresses the practicability of using downward deployed satellites tethered to an orbiting host vehicle in order to obtain steady-state data in the upper reaches of the region above 80 or 90 km.
Shuttle entry aerothermodynamic flight research - The Orbiter Experiments (OEX) Program
NASA Technical Reports Server (NTRS)
Throckmorton, David A.
1992-01-01
Results of the OEX program are summarized with emphasis on the information on entry aerothermodynamic phenomena derived from Space Shuttle operations. The discussion focuses on OEX experiment complement and operational history, freestream environment and vehicle attitude data, aerodynamic force and moment data, aerodynamic surface data, and vehicle configuration data. Attention is also given to orbiter aerodynamic performance, stability and control, high-altitude atmospheric density variability, direct simulation Monte Carlo validation, orbital drag variation, and computational fluid dynamic technique validation.
Numerical analysis and design of upwind sails
NASA Astrophysics Data System (ADS)
Shankaran, Sriram
The use of computational techniques that solve the Euler or the Navier-Stokes equations are increasingly being used by competing syndicates in races like the Americas Cup. For sail configurations, this desire stems from a need to understand the influence of the mast on the boundary layer and pressure distribution on the main sail, the effect of camber and planform variations of the sails on the driving and heeling force produced by them and the interaction of the boundary layer profile of the air over the surface of the water and the gap between the boom and the deck on the performance of the sail. Traditionally, experimental methods along with potential flow solvers have been widely used to quantify these effects. While these approaches are invaluable either for validation purposes or during the early stages of design, the potential advantages of high fidelity computational methods makes them attractive candidates during the later stages of the design process. The aim of this study is to develop and validate numerical methods that solve the inviscid field equations (Euler) to simulate and design upwind sails. The three dimensional compressible Euler equations are modified using the idea of artificial compressibility and discretized on unstructured tetrahedral grids to provide estimates of lift and drag for upwind sail configurations. Convergence acceleration techniques like multigrid and residual averaging are used along with parallel computing platforms to enable these simulations to be performed in a few minutes. To account for the elastic nature of the sail cloth, this flow solver was coupled to NASTRAN to provide estimates of the deflections caused by the pressure loading. The results of this aeroclastic simulation, showed that the major effect of the sail elasticity; was in altering the pressure distribution around the leading edge of the head and the main sail. Adjoint based design methods were developed next and were used to induce changes to the camber
Study of second order upwind differencing in a recirculating flow
NASA Technical Reports Server (NTRS)
Vanka, S. P.
1985-01-01
The accuracy and stability of the second order upwind differencing scheme was investigated. The solution algorithm employed is based on a coupled solution of the nonlinear finite difference equations by the multigrid technique. Calculations have been made of the driven cavity flow for several Reynolds numbers and finite difference grids. In comparison with the hybrid differencing, the second order upwind differencing is somewhat more accurate but it is not monotonically accurate with mesh refinement. Also, the convergence of the solution algorithm deteriorates with the use of the second order upwind differencing.
Kinetic theory based new upwind methods for inviscid compressible flows
NASA Technical Reports Server (NTRS)
Deshpande, S. M.
1986-01-01
Two new upwind methods called the Kinetic Numerical Method (KNM) and the Kinetic Flux Vector Splitting (KFVS) method for the solution of the Euler equations have been presented. Both of these methods can be regarded as some suitable moments of an upwind scheme for the solution of the Boltzmann equation provided the distribution function is Maxwellian. This moment-method strategy leads to a unification of the Riemann approach and the pseudo-particle approach used earlier in the development of upwind methods for the Euler equations. A very important aspect of the moment-method strategy is that the new upwind methods satisfy the entropy condition because of the Boltzmann H-Theorem and suggest a possible way of extending the Total Variation Diminishing (TVD) principle within the framework of the H-Theorem. The ability of these methods in obtaining accurate wiggle-free solution is demonstrated by applying them to two test problems.
Towards a genuinely multi-dimensional upwind scheme
NASA Technical Reports Server (NTRS)
Powell, Kenneth G.; Vanleer, Bram; Roe, Philip L.
1990-01-01
Methods of incorporating multi-dimensional ideas into algorithms for the solution of Euler equations are presented. Three schemes are developed and tested: a scheme based on a downwind distribution, a scheme based on a rotated Riemann solver and a scheme based on a generalized Riemann solver. The schemes show an improvement over first-order, grid-aligned upwind schemes, but the higher-order performance is less impressive. An outlook for the future of multi-dimensional upwind schemes is given.
SUPG Finite Element Simulations of Compressible Flows for Aerothermodynamic Applications
NASA Technical Reports Server (NTRS)
Kirk, Benjamin S.
2007-01-01
This viewgraph presentation reviews the Streamline-Upwind Petrov-Galerkin (SUPG) Finite Element Simulation. It covers the background, governing equations, weak formulation, shock capturing, inviscid flux discretization, time discretization, linearization, and implicit solution strategies. It also reviews some applications such as Type IV Shock Interaction, Forward-Facing Cavity and AEDC Sharp Double Cone.
Experimental Stage Separation Tool Development in NASA Langley's Aerothermodynamics Laboratory
NASA Technical Reports Server (NTRS)
Murphy, Kelly J.; Scallion, William I.
2005-01-01
As part of the research effort at NASA in support of the stage separation and ascent aerothermodynamics research program, proximity testing of a generic bimese wing-body configuration was conducted in NASA Langley's Aerothermodynamics Laboratory in the 20-Inch Mach 6 Air Tunnel. The objective of this work is the development of experimental tools and testing methodologies to apply to hypersonic stage separation problems for future multi-stage launch vehicle systems. Aerodynamic force and moment proximity data were generated at a nominal Mach number of 6 over a small range of angles of attack. The generic bimese configuration was tested in a belly-to-belly and back-to-belly orientation at 86 relative proximity locations. Over 800 aerodynamic proximity data points were taken to serve as a database for code validation. Longitudinal aerodynamic data generated in this test program show very good agreement with viscous computational predictions. Thus a framework has been established to study separation problems in the hypersonic regime using coordinated experimental and computational tools.
Aerothermodynamic Environments Definition for the Mars Science Laboratory Entry Capsule
NASA Technical Reports Server (NTRS)
Edquist, Karl T.; Dyakonov, Artem A.; Wright, Michael J.; Tang, Chun Y.
2007-01-01
An overview of the aerothermodynamic environments definition status is presented for the Mars Science Laboratory entry vehicle. The environments are based on Navier-Stokes flowfield simulations on a candidate aeroshell geometry and worst-case entry heating trajectories. Uncertainties for the flowfield predictions are based primarily on available ground data since Mars flight data are scarce. The forebody aerothermodynamics analysis focuses on boundary layer transition and turbulent heating augmentation. Turbulent transition is expected prior to peak heating, a first for Mars entry, resulting in augmented heat flux and shear stress at the same heatshield location. Afterbody computations are also shown with and without interference effects of reaction control system thruster plumes. Including uncertainties, analysis predicts that the heatshield may experience peaks of 225 W/sq cm for turbulent heat flux, 0.32 atm for stagnation pressure, and 400 Pa for turbulent shear stress. The afterbody heat flux without thruster plume interference is predicted to be 7 W/sq cm on the backshell and 10 W/sq cm on the parachute cover. If the reaction control jets are fired near peak dynamic pressure, the heat flux at localized areas could reach as high as 76 W/sq cm on the backshell and 38 W/sq cm on the parachute cover, including uncertainties. The final flight environments used for hardware design will be updated for any changes in the aeroshell configuration, heating design trajectories, or uncertainties.
NASA Astrophysics Data System (ADS)
Rufolo, Giuseppe C.; Di Benedetto, Sara; Walpot, Louis; Roncioni, Pietro; Marini, Marco
2011-05-01
In the frame of the Intermediate eXperimental Vehicle (IXV) project, the European Space Agency (ESA) is coordinating a series of technical assistance activities aimed at verifying and supporting the IXV industrial design and development process. The technical assistance is operated with the support of the Italian Space Agency (ASI), by means of the Italian Aerospace Research Center (CIRA), and the European Space Research and Technology Centre (ESTEC) under the super visioning and coordination of ESA IXV team. One of the purposes of the activity is to develop an independent capability for the assessment and verification of the industrial results with respect to the aerothermodynamic characterization of the IXV vehicle. To this aim CIRA is developing and independent AeroThermodynamics DataBase (ATDB), intended as a tool generating in output the time histories of local quantities (heat flux, pressure, skin friction) for each point of the IXV vehicle and for each trajectory (in a pre-defined envelope), together with an uncertainties model. The reference Computational Fluid Dynamics (CFD) solutions needed for the development of the tool have been provided by ESA-ESTEC (with the CFD code LORE) and CIRA (with the CFD code H3NS).
Development and application of computational aerothermodynamics flowfield computer codes
NASA Technical Reports Server (NTRS)
Venkatapathy, Ethiraj
1994-01-01
Research was performed in the area of computational modeling and application of hypersonic, high-enthalpy, thermo-chemical nonequilibrium flow (Aerothermodynamics) problems. A number of computational fluid dynamic (CFD) codes were developed and applied to simulate high altitude rocket-plume, the Aeroassist Flight Experiment (AFE), hypersonic base flow for planetary probes, the single expansion ramp model (SERN) connected with the National Aerospace Plane, hypersonic drag devices, hypersonic ramp flows, ballistic range models, shock tunnel facility nozzles, transient and steady flows in the shock tunnel facility, arc-jet flows, thermochemical nonequilibrium flows around simple and complex bodies, axisymmetric ionized flows of interest to re-entry, unsteady shock induced combustion phenomena, high enthalpy pulsed facility simulations, and unsteady shock boundary layer interactions in shock tunnels. Computational modeling involved developing appropriate numerical schemes for the flows on interest and developing, applying, and validating appropriate thermochemical processes. As part of improving the accuracy of the numerical predictions, adaptive grid algorithms were explored, and a user-friendly, self-adaptive code (SAGE) was developed. Aerothermodynamic flows of interest included energy transfer due to strong radiation, and a significant level of effort was spent in developing computational codes for calculating radiation and radiation modeling. In addition, computational tools were developed and applied to predict the radiative heat flux and spectra that reach the model surface.
Aerothermodynamic Design of the Mars Science Laboratory Heatshield
NASA Technical Reports Server (NTRS)
Edquist, Karl T.; Dyakonov, Artem A.; Wright, Michael J.; Tang, Chun Y.
2009-01-01
Aerothermodynamic design environments are presented for the Mars Science Laboratory entry capsule heatshield. The design conditions are based on Navier-Stokes flowfield simulations on shallow (maximum total heat load) and steep (maximum heat flux, shear stress, and pressure) entry trajectories from a 2009 launch. Boundary layer transition is expected prior to peak heat flux, a first for Mars entry, and the heatshield environments were defined for a fully-turbulent heat pulse. The effects of distributed surface roughness on turbulent heat flux and shear stress peaks are included using empirical correlations. Additional biases and uncertainties are based on computational model comparisons with experimental data and sensitivity studies. The peak design conditions are 197 W/sq cm for heat flux, 471 Pa for shear stress, 0.371 Earth atm for pressure, and 5477 J/sq cm for total heat load. Time-varying conditions at fixed heatshield locations were generated for thermal protection system analysis and flight instrumentation development. Finally, the aerothermodynamic effects of delaying launch until 2011 are previewed.
Development of X-33/X-34 Aerothermodynamic Data Bases: Lessons Learned and Future Enhancements
NASA Technical Reports Server (NTRS)
Miller, C. G.
1999-01-01
A synoptic of programmatic and technical lessons learned in the development of aerothermodynamic data bases for the X-33 and X-34 programs is presented in general terms and from the perspective of the NASA Langley Research Center Aerothermodynamics Branch. The format used is that of the aerothermodynamic chain, the links of which are personnel, facilities, models/test articles, instrumentation, test techniques, and computational fluid dynamics (CFD). Because the aerodynamic data bases upon which the X-33 and X-34 vehicles will fly are almost exclusively from wind tunnel testing, as opposed to CFD, the primary focus of the lessons learned is on ground-based testing.
Upwind impacts of ammonia from an intensive poultry unit.
Jones, L; Nizam, M S; Reynolds, B; Bareham, S; Oxley, E R B
2013-09-01
This study investigated potential ammonia impacts on a sand dune nature reserve 600 m upwind of an intensive poultry unit. Ammonia concentrations and total nitrogen deposition were measured over a calendar year. A series of ammonia and nitrogen exposure experiments using dune grassland species were conducted in controlled manipulations and in the field. Ammonia emissions from the intensive poultry unit were detected up to 2.8 km upwind, contributing to exceedance of critical levels of ammonia 800 m upwind and exceedance of critical loads of nitrogen 2.8 km upwind. Emissions contributed 30% of the total N load in parts of the upwind conservation site. In the nitrogen exposure experiments, plants showed elevated tissue nitrogen contents, and responded to ammonia concentrations and nitrogen deposition loads observed in the conservation site by increasing biomass. Estimated long-term impacts suggest an increase in the soil carbon pool of 9% over a 50-year timescale. PMID:23792381
Solving Upwind-Biased Discretizations. 2; Multigrid Solver Using Semicoarsening
NASA Technical Reports Server (NTRS)
Diskin, Boris
1999-01-01
This paper studies a novel multigrid approach to the solution for a second order upwind biased discretization of the convection equation in two dimensions. This approach is based on semi-coarsening and well balanced explicit correction terms added to coarse-grid operators to maintain on coarse-grid the same cross-characteristic interaction as on the target (fine) grid. Colored relaxation schemes are used on all the levels allowing a very efficient parallel implementation. The results of the numerical tests can be summarized as follows: 1) The residual asymptotic convergence rate of the proposed V(0, 2) multigrid cycle is about 3 per cycle. This convergence rate far surpasses the theoretical limit (4/3) predicted for standard multigrid algorithms using full coarsening. The reported efficiency does not deteriorate with increasing the cycle, depth (number of levels) and/or refining the target-grid mesh spacing. 2) The full multi-grid algorithm (FMG) with two V(0, 2) cycles on the target grid and just one V(0, 2) cycle on all the coarse grids always provides an approximate solution with the algebraic error less than the discretization error. Estimates of the total work in the FMG algorithm are ranged between 18 and 30 minimal work units (depending on the target (discretizatioin). Thus, the overall efficiency of the FMG solver closely approaches (if does not achieve) the goal of the textbook multigrid efficiency. 3) A novel approach to deriving a discrete solution approximating the true continuous solution with a relative accuracy given in advance is developed. An adaptive multigrid algorithm (AMA) using comparison of the solutions on two successive target grids to estimate the accuracy of the current target-grid solution is defined. A desired relative accuracy is accepted as an input parameter. The final target grid on which this accuracy can be achieved is chosen automatically in the solution process. the actual relative accuracy of the discrete solution approximation
Computational Aerodynamic Analysis of Offshore Upwind and Downwind Turbines
Zhao, Qiuying; Sheng, Chunhua; Afjeh, Abdollah
2014-01-01
Aerodynamic interactions of the model NREL 5 MW offshore horizontal axis wind turbines (HAWT) are investigated using a high-fidelity computational fluid dynamics (CFD) analysis. Four wind turbine configurations are considered; three-bladed upwind and downwind and two-bladed upwind and downwind configurations, which operate at two different rotor speeds of 12.1 and 16 RPM. In the present study, both steady and unsteady aerodynamic loads, such as the rotor torque, blade hub bending moment, and base the tower bending moment of the tower, are evaluated in detail to provide overall assessment of different wind turbine configurations. Aerodynamic interactions between the rotor and tower are analyzed,more » including the rotor wake development downstream. The computational analysis provides insight into aerodynamic performance of the upwind and downwind, two- and three-bladed horizontal axis wind turbines.« less
Aerothermodynamic design feasibility of a Mars aerocapture/aeromaneuver vehicle
NASA Technical Reports Server (NTRS)
Florence, D. E.
1981-01-01
Lifting aerodynamic configurations have been screened and selected for the Mars aerocapture mission that (1) meet the geometric packaging requirements of the various payloads and the Space Shuttle cargo bay and (2) provide the aerodynamic performance characteristics required to obtain the atmospheric exit steering accuracy and the parachute deployment conditions desired. Hypersonic heat transfer and aerodynamic loads to the vehicle in the CO2 atmosphere are evaluated. Contemporary low density ablative thermal protection materials were selected that meet all the atmospheric entry requirements and provide a minimum mass solution. Results are presented of the aerodynamic configuration and thermal protection materials screening and selection. It is concluded that the aerothermodynamic design of this concept is feasible using state-of-the-art technology.
High-Energy Atmospheric Reentry Test Aerothermodynamic Analysis
NASA Technical Reports Server (NTRS)
Mazaheri, Alireza
2013-01-01
This paper presents an assessment of the aerothermodynamic environment around an 8.3 meter High Energy Atmospheric Reentry Test (HEART) vehicle. This study generated twelve nose shape configurations and compared their responses at the peak heating trajectory point against the baseline nose shape. The heat flux sensitivity to the angle of attack variations are also discussed. The possibility of a two-piece Thermal Protection System (TPS) design at the nose is also considered, as are the surface catalytic affects of the aeroheating environment of such configuration. Based on these analyses, an optimum nose shape is proposed to minimize the surface heating. A recommendation is also made for a two-piece TPS design, for which the surface catalytic uncertainty associated with the jump in heating at the nose-IAD juncture is reduced by a minimum of 93%. In this paper, the aeroshell is assumed to be rigid and the inflatable fluid interaction effect is left for future investigations
Team Software Development for Aerothermodynamic and Aerodynamic Analysis and Design
NASA Technical Reports Server (NTRS)
Alexandrov, N.; Atkins, H. L.; Bibb, K. L.; Biedron, R. T.; Carpenter, M. H.; Gnoffo, P. A.; Hammond, D. P.; Jones, W. T.; Kleb, W. L.; Lee-Rausch, E. M.
2003-01-01
A collaborative approach to software development is described. The approach employs the agile development techniques: project retrospectives, Scrum status meetings, and elements of Extreme Programming to efficiently develop a cohesive and extensible software suite. The software product under development is a fluid dynamics simulator for performing aerodynamic and aerothermodynamic analysis and design. The functionality of the software product is achieved both through the merging, with substantial rewrite, of separate legacy codes and the authorship of new routines. Examples of rapid implementation of new functionality demonstrate the benefits obtained with this agile software development process. The appendix contains a discussion of coding issues encountered while porting legacy Fortran 77 code to Fortran 95, software design principles, and a Fortran 95 coding standard.
Aerothermodynamic Facilities And Measurement: Flow Characterization in Shock Tunnels
NASA Technical Reports Server (NTRS)
Cavolowsky, John A.; Edwards, Thomas A. (Technical Monitor)
1995-01-01
This presentation will examine the key performance aspects of shock tunnels as they relate to their use as aerothermodynamic flow simulation facilities. Assessment of shock tube reservoir conditions and flow contaminants generated in the shock tube will be presented along with their limiting impact on viable test envelopes, Facility nozzle performance as it pertains to test time assessment and nozzle exit flow quality (survey of pressure, temperature, and species) will be addressed. Also included will be a discussion of free stream flow diagnostics, both intrusive and nonintrusive, for measurement of critical flow properties not directly inferred from surface mounted transducers. The use of computational fluid dynamics for purposes of validating experimental measurements as well as predicting performance in regimes where measurements are not feasible or possible will be discussed. The use of CFD for facility research and design will also be presented.
Experimental and Computational Aerothermodynamics of a Mars Entry Vehicle
NASA Technical Reports Server (NTRS)
Hollis, Brian R.
1996-01-01
An aerothermodynamic database has been generated through both experimental testing and computational fluid dynamics simulations for a 70 deg sphere-cone configuration based on the NASA Mars Pathfinder entry vehicle. The aerothermodynamics of several related parametric configurations were also investigated. Experimental heat-transfer data were obtained at hypersonic test conditions in both a perfect gas air wind tunnel and in a hypervelocity, high-enthalpy expansion tube in which both air and carbon dioxide were employed as test gases. In these facilities, measurements were made with thin-film temperature-resistance gages on both the entry vehicle models and on the support stings of the models. Computational results for freestream conditions equivalent to those of the test facilities were generated using an axisymmetric/2D laminar Navier-Stokes solver with both perfect-gas and nonequilibrium thermochemical models. Forebody computational and experimental heating distributions agreed to within the experimental uncertainty for both the perfect-gas and high-enthalpy test conditions. In the wake, quantitative differences between experimental and computational heating distributions for the perfect-gas conditions indicated transition of the free shear layer near the reattachment point on the sting. For the high enthalpy cases, agreement to within, or slightly greater than, the experimental uncertainty was achieved in the wake except within the recirculation region, where further grid resolution appeared to be required. Comparisons between the perfect-gas and high-enthalpy results indicated that the wake remained laminar at the high-enthalpy test conditions, for which the Reynolds number was significantly lower than that of the perfect-gas conditions.
Accurate upwind-monotone (nonoscillatory) methods for conservation laws
NASA Technical Reports Server (NTRS)
Huynh, Hung T.
1992-01-01
The well known MUSCL scheme of Van Leer is constructed using a piecewise linear approximation. The MUSCL scheme is second order accurate at the smooth part of the solution except at extrema where the accuracy degenerates to first order due to the monotonicity constraint. To construct accurate schemes which are free from oscillations, the author introduces the concept of upwind monotonicity. Several classes of schemes, which are upwind monotone and of uniform second or third order accuracy are then presented. Results for advection with constant speed are shown. It is also shown that the new scheme compares favorably with state of the art methods.
The design and application of upwind schemes on unstructured meshes
NASA Technical Reports Server (NTRS)
Barth, Timothy J.; Jespersen, Dennis C.
1989-01-01
Solution and mesh generation algorithms for solving the Euler equations on unstructured meshes consisting of triangle and quadrilateral control volumes are presented. Cell-centered and mesh-vertex upwind finite-volume schemes are developed which utilize multi-dimensional monotone linear reconstruction procedures. These algorithms differ from existing algorithms (even on structured meshes). Numerical results in two dimensions are presented.
Compressible bubble dynamic simulations with central-upwind schemes
NASA Astrophysics Data System (ADS)
Koukouvinis, P.; Gavaises, M.; Georgoulas, A.; Marengo, M.
2015-12-01
This paper discusses the implementation of an explicit density-based solver, based on the central-upwind schemes originally suggested by Kurganov, for the simulation of cavitating bubble dynamic flows. Explicit density based solvers are suited for highly dynamic, violent flows, involving large density ratios, as is rather common in cavitating flows. Moreover, the central-upwind schemes have the advantage of avoiding direct evaluation of the Jacobian matrix or estimation of the wave pattern emerging from Euler equations. Second order accuracy can be achieved with TVD MUSCL schemes. Basic comparison with the predicted wave pattern of the central-upwind schemes is performed with the exact solution of the Riemann problem showing an excellent agreement. Then several different bubble configurations were tested, similar to the work of Lauer et al. (2012). The central-upwind schemes prove to be able to handle the large pressure and density ratios appearing in cavitating flows, giving similar predictions in the evolution of the bubble shape.
Experimental Aerothermodynamics In Support Of The Columbia Accident Investigation
NASA Technical Reports Server (NTRS)
Horvath, Thomas J.
2004-01-01
The technical foundation for the most probable damage scenario reported in the Columbia Accident Investigation Board's final report was largely derived from synergistic aerodynamic/aerothermodynamic wind tunnel measurements and inviscid predictions made at NASA Langley Research Center and later corroborated with engineering analysis, high fidelity numerical viscous simulations, and foam impact testing near the close of the investigation. This report provides an overview of the hypersonic aerothermodynamic wind tunnel program conducted at NASA Langley and illustrates how the ground-based heating measurements provided early insight that guided the direction and utilization of agency resources in support of the investigation. Global surface heat transfer mappings, surface streamline patterns, and shock shapes were measured on 0.0075 scale models of the Orbiter configuration with and without postulated damage to the thermal protection system. Test parametrics include angle of attack from 38 to 42 degs, sideslip angles of 38 to 42 degs, sideslip angles of plus or minus 1 deg, Reynolds numbers based upon model length from 0.05 x 10(exp 6) to 6.5 x 10(exp 6), and normal shock density ratios of 5 (Mach 6 Air) and 12 (Mach 6 CF4). The primary objective of the testing was to provide surface heating characteristics on scaled Orbiter models with outer mold line perturbations to simulate various forms of localized surface damage to the thermal protection system. Initial experimental testing conducted within two weeks of the accident simulated a broad spectrum of thermal protection system damage to the Orbiter windward surface and was used to refute several hypothesized forms of thermal protection system damage, which included gouges in the windward thermal protection system tiles, breaches through the wing new the main landing gear door, and protuberances along the wing leading edge that produced asymmetric boundary layer transition. As the forensic phase of the investigation
NASA Technical Reports Server (NTRS)
Arnold, James O.; Deiwert, G. S.
1997-01-01
The dream of producing an air-breathing, hydrogen fueled, hypervelocity aircraft has been before the aerospace community for decades. However, such a craft has not yet been realized, even in an experimental form. Despite the simplicity and beauty of the concept, many formidable problems must be overcome to make this dream a reality. This paper summarizes the aero/aerothermodynamic issues that must be addressed to make the dream a reality and discusses how aerothermodynamics facilities and their modem companion, real-gas computational fluid dynamics (CFD), can help solve the problems blocking the way to realizing the dream. The approach of the paper is first to outline the concept of an air-breathing hypersonic vehicle and then discuss the nose-to-tail aerothermodynamics issues and special aerodynamic problems that arise with such a craft. Then the utility of aerothermodynamic facilities and companion CFD analysis is illustrated by reviewing results from recent United States publications wherein these problems have been addressed. Papers selected for the discussion have k e n chosen such that the review will serve to survey important U.S. aero/aerothermodynamic real gas and conventional wind tunnel facilities that are useful in the study of hypersonic, hydrogen propelled hypervelocity vehicles.
Nonequilibrium effects on the aerothermodynamics of transatmospheric and aerobraking vehicles
NASA Technical Reports Server (NTRS)
Hassan, Basil; Candler, Graham V.
1993-01-01
A 3D CFD algorithm is used to study the effect of thermal and chemical nonequilibrium on slender and blunt body aerothermodynamics. Both perfect gas and reacting gas air models are used to compute the flow over a generic transatmospheric vehicle and a proposed lunar transfer vehicle. The reacting air is characterized by a translational-rotational temperature and a vibrational-electron-electronic temperature and includes eight chemical species. The effects of chemical reaction, vibrational excitation, and ionization on lift-to-drag ratio and trim angle are investigated. Results for the NASA Ames All-body Configuration show a significant difference in center of gravity location for a reacting gas flight case when compared to a perfect gas wind tunnel case at the same Mach number, Reynolds number, and angle of attack. For the same center of gravity location, the wind tunnel model trims at lower angle of attack than the full-scale flight case. Nonionized and ionized results for a proposed lunar transfer vehicle compare well to computational results obtained from a previously validated reacting gas algorithm. Under the conditions investigated, effects of weak ionization on the heat transfer and aerodynamic coefficients were minimal.
Aerothermodynamic Analysis of the Project FIRE II Afterbody Flow
NASA Technical Reports Server (NTRS)
Wright, Micheal J.; Loomis, Mark; Arnold, Jim (Technical Monitor)
2000-01-01
35 years later, the Project FIRE II ballistic reentry to Earth at a nominal velocity of 11.4 km/s remains one of the best sources of heating data for the design of sample return capsules. The data from this flight experiment encompass both the thermochemical non-equilibrium and equilibrium flow regimes and include measurements of both radiative and total heating on the forebody and afterbody. Because of this, a number of researchers have performed computational fluid dynamics (CFD) simulations of the forebody of the FIRE II entry vehicle, with generally good results. In particular, Olynick et. al. coupled a Navier-Stokes solver (GIANTS) with a radiation code (NOVAR) and showed excellent agreement in surface heat transfer over the FIRE II trajectory between 1634 and 1651 seconds (77 km to 37 km). However, in most cases the primary motivation of the previous work was to understand and model the coupling between shock layer radiation and aerothermodynamics, and thus the simulations concentrated on the forebody flow only. To our knowledge there have been no prior published attempts to reproduce the afterbody heating data. However, an understanding of this data is critical to our efforts to design the next generation of Earth and planetary entry vehicles and to assess our need for additional flight data.
Aerothermodynamics Feasibility Assessment of a Mars Atmoshperic Sample Return Mission
NASA Astrophysics Data System (ADS)
Ferracina, L.; Larranaga, J.; Falkner, P.
2011-02-01
ESA's optional Mars Robotic Exploration Preparation (MREP) programme is based on a long term collaboration with NASA, by taking Mars exploration as global objective, and Mars Sample Return (MSR) mission as long term goal to be achieved by the mid 2020's. Considering today's uncertainties, different missions are envisaged and prepared by ESA as possible alternative missions to MSR in the timeframe of 2020- 2026, in case the required technology readiness is not reached by 2015 or landed mass capabilities are exceeded for any of the MSR mission elements. One of the ESA considered missions within this framework is the Mars Atmospheric Sample Return Mission. This mission has been recently assessed by ESA using its Concurrent Design Facility (CDF), aiming to enter with a probe at Mars low altitudes (≈50 km), collect a sample of airborne atmosphere (gas and dust) and return the sample back to Earth. This paper aim at reporting the preliminary aerothermodynamic assessment of the design of the Martian entry probe conducted within the CDF study. Special attention has been paid to the selection of aerodynamically efficient vehicle concepts compare to blunt bodies and to the effect of the hot-temperature shock to the cavity placed at stagnation point and used in the atmospheric sampling system.
Aerothermodynamic environments of aerobraking vehicles for manned Mars missions
NASA Technical Reports Server (NTRS)
Ledoux, Stephen T.; Vas, Irwin E.
1991-01-01
The aerothermodynamic environments of manned spacraft aerobraking in the Martian and earth atmospheres are evaluated. Thermal performance of aerobrake concepts are examined for current cryogenic-aerobrake and advanced propulsion missions entailing three different modes of aerobraking: (1) aerocapture into an orbit about Mars, (2) descent and landing at Mars, and (3) Mars return direct entry at earth. Analyses for these vehicles and modes included both radiative and convective heating, where radiative heating is shown to be a significant portion of the total stagnation point heating induced on the vehicle. A comprehensive parametric study of the effects of ballistic coefficient, nose radius, entry velocity, and L/D on stagnation point heating is described. Optimal nose radii for ranges of ballistic coefficient and entry velocity are determined. The peak heating rates are shown to be 83 W/sq cm and 90 W/sq cm for a low and high L/D Mars transfer vehicle configuration, respectively. Heating profiles for these vehicles using boundary layer techniques show that a high L/D shape will result in a smaller high-temperature region provided the flow is laminar. An examination of a crew return vehicle for a Mars return direct entry trajectory shows that the thermal protection for this aerobrake will require an ablative material for heat rejection due to the large heating rates (about 1 kW/sq cm).
Mars Science Laboratory Entry Capsule Aerothermodynamics and Thermal Protection System
NASA Technical Reports Server (NTRS)
Edquist, Karl T.; Hollis, Brian R.; Dyakonov, Artem A.; Laub, Bernard; Wright, Michael J.; Rivellini, Tomasso P.; Slimko, Eric M.; Willcockson, William H.
2007-01-01
The Mars Science Laboratory (MSL) spacecraft is being designed to carry a large rover (greater than 800 kg) to the surface of Mars using a blunt-body entry capsule as the primary decelerator. The spacecraft is being designed for launch in 2009 and arrival at Mars in 2010. The combination of large mass and diameter with non-zero angle-of-attack for MSL will result in unprecedented convective heating environments caused by turbulence prior to peak heating. Navier-Stokes computations predict a large turbulent heating augmentation for which there are no supporting flight data1 and little ground data for validation. Consequently, an extensive experimental program has been established specifically for MSL to understand the level of turbulent augmentation expected in flight. The experimental data support the prediction of turbulent transition and have also uncovered phenomena that cannot be replicated with available computational methods. The result is that the flight aeroheating environments predictions must include larger uncertainties than are typically used for a Mars entry capsule. Finally, the thermal protection system (TPS) being used for MSL has not been flown at the heat flux, pressure, and shear stress combinations expected in flight, so a test program has been established to obtain conditions relevant to flight. This paper summarizes the aerothermodynamic definition analysis and TPS development, focusing on the challenges that are unique to MSL.
Aerothermodynamic Environment Definition for the Genesis Sample Return Capsule
NASA Technical Reports Server (NTRS)
Cheatwood, F. McNeil; Merski, N. Ronald, Jr.; Riley, Christopher J.; Mitcheltree, Robert A.
2001-01-01
NASA's Genesis sample return mission will be the first to return material from beyond the Earth-Moon system. NASA Langley Research Center supported this mission with aerothermodynamic analyses of the sample return capsule. This paper provides an overview of that effort. The capsule is attached through its forebody to the spacecraft bus. When the attachment is severed prior to Earth entry, forebody cavities remain. The presence of these cavities could dramatically increase the heating environment in their vicinity and downstream. A combination of computational fluid dynamics calculations and wind tunnel phosphor thermography tests were employed to address this issue. These results quantify the heating environment in and around the cavities, and were a factor in the decision to switch forebody heat shield materials. A transition map is developed which predicts that the flow aft of the penetrations will still be laminar at the peak heating point of the trajectory. As the vehicle continues along the trajectory to the peak dynamic pressure point, fully turbulent flow aft of the penetrations could occur. The integrated heat load calculations show that a heat shield sized to the stagnation point levels will be adequate for the predicted environment aft of the penetrations.
Second order upwind Lagrangian particle method for Euler equations
Samulyak, Roman; Chen, Hsin -Chiang; Yu, Kwangmin
2016-06-01
A new second order upwind Lagrangian particle method for solving Euler equations for compressible inviscid fluid or gas flows is proposed. Similar to smoothed particle hydrodynamics (SPH), the method represents fluid cells with Lagrangian particles and is suitable for the simulation of complex free surface / multiphase flows. The main contributions of our method, which is different from SPH in all other aspects, are (a) significant improvement of approximation of differential operators based on a polynomial fit via weighted least squares approximation and the convergence of prescribed order, (b) an upwind second-order particle-based algorithm with limiter, providing accuracy and longmore » term stability, and (c) accurate resolution of states at free interfaces. In conclusion, numerical verification tests demonstrating the convergence order for fixed domain and free surface problems are presented.« less
Are upwind techniques in multi-phase flow models necessary?
Park, C.-H.; Boettcher, N.; Wang, W.; Kolditz, O.
2011-09-10
Two alternatives of primary variables are compared for two-phase flow in heterogeneous media by solving fully established benchmarks. The first combination utilizes pressure of the wetting fluid and saturation of the non-wetting fluid as primary variables, while the second employs capillary pressure of the wetting fluid and pressure of the non-wetting fluid. While the standard Galerkin finite element method (SGFEM) is known to fail in the physical reproduction of two-phase flow in heterogeneous media (unless employing a fully upwind correction), the second scheme with capillary pressure as a primary variable without applying an upwind technique produces correct physical fluid behaviour in heterogeneous media, as observed from experiments.
NASA Astrophysics Data System (ADS)
Reynier, Philippe
2014-10-01
This contribution is a survey of aerodynamic and aerothermodynamics data related to Mars entry. The survey includes the studies carried out in the frame of projects aiming at preparing exploration missions involving entry probes into Mars atmosphere and the efforts have been concentrated on the aerothermodynamics developments. Russian (including former Soviet Union), European and NASA aerothermodynamics developments for preparing such missions have been accounted for. If a focus has been dedicated to the flight data gathered during Viking and Mars Pathfinder entries, the experimental and numerical activities carried out for the different projects have been also considered. The emphasis has been put on the post-flight analysis of flight experiments. The objective of the activity has been to develop a database of the developments performed for Mars entry that will be of interest for the preparation of future missions and for testing new models related to radiative transfer, and chemical kinetics schemes based on a state-to-state approach.
The use of the tethered satellite system to perform low density aerothermodynamics studies
NASA Technical Reports Server (NTRS)
Carlomagno, Giovanni M.; Deluca, Luigi; Siemers, Paul M.; Wood, George M., Jr.
1988-01-01
The Tethered Satellite System (TSS) is a cooperative space system development activity of the U.S.A. and Italy. It is comprised of the Tether Satellite (TS) and the deployer. Within TSS, the Shuttle Tethered Aerothermodynamic Research Facility (STARFAC) concept has the potential to provide access to vast portions of the upper atmosphere for atmospheric and aerothermodynamic research. The feasibility and capability of the TSS to operate as a continuous open wind tunnel and to perform low density aerothermodynamic studies are investigated. This is accomplished through a modified version of the TS simulation program (SKYHOOK). The results indicate that STARFAC concept is both feasible and practical. The TS can go below 100 km but, if thrust is used, large velocity variation (delta V) maneuvers and an attitude control are required; if a satellite lift is considered, large tether tension is produced and an attitude control is required.
In-flight measurement of upwind dynamic soaring in albatrosses
NASA Astrophysics Data System (ADS)
Sachs, Gottfried
2016-03-01
In-flight measurement results on upwind flight of albatrosses using dynamic soaring are presented. It is shown how the birds manage to make progress against the wind on the basis of small-scale dynamic soaring maneuvers. For this purpose, trajectory features, motion quantities and mechanical energy relationships as well as force characteristics are analyzed. The movement on a large-scale basis consists of a tacking type flight technique which is composed of dynamic soaring cycle sequences with alternating orientation to the left and right. It is shown how this is performed by the birds so that they can achieve a net upwind flight without a transversal large-scale movement and how this compares with downwind or across wind flight. Results on upwind dynamic soaring are presented for low and high wind speed cases. It is quantified how much the tacking trajectory length is increased when compared with the beeline distance. The presented results which are based on in-flight measurements of free flying albatrosses were achieved with an in-house developed GPS-signal tracking method yielding the required high precision for the small-scale dynamic soaring flight maneuvers.
X-38 NASA/DLR/ESA-Dassault Aviation Integrated Aerodynamic and Aerothermodynamic Activities
NASA Technical Reports Server (NTRS)
Labbe, Steve G.; Perez, Leo F.; Fitzgerald, Steve; Longo, Jose; Rapuc, Marc; Molina, Rafael; Nicholson, Leonard S. (Technical Monitor)
1999-01-01
The characterization of the aeroshape selected for the X-38 [Crew Return Vehicle (CRV) demonstrator] is presently being performed as a cooperative endeavour between NASA, DLR (through its TETRA Program), and European Space Agency (ESA) with Dassault Aviation integrating the aerodynamic and aerothermodynamic activities. The methodologies selected for characterizing the aerodynamic and aerothermodynamic environment of the X-38 are presented. Also, the implications for related disciplines such as Guidance Navigation and Control (GN&C) with its corresponding Flight Control System (FCS), Structural, and Thermal Protection System (TPS) design are discussed. An attempt is made at defining the additional activities required to support the design of a derived operational CRV.
Langley Aerothermodynamic Facilities Complex: Enhancements and Testing Capabilities
NASA Technical Reports Server (NTRS)
Micol, J. R.
1998-01-01
Description, capabilities, recent upgrades, and utilization of the NASA Langley Research Center (LaRC) Aerothermodynamic Facilities Complex (AFC) are presented. The AFC consists of five hypersonic, blow-down-to-vacuum wind tunnels that collectively provide a range of Mach number from 6 to 20, unit Reynolds number from 0.04 to 22 million per foot and, most importantly for blunt configurations, normal shock density ratio from 4 to 12. These wide ranges of hypersonic simulation parameters are due, in part, to the use of three different test gases (air, helium, and tetrafluoromethane), thereby making several of the facilities unique. The Complex represents nearly three-fourths of the conventional (as opposed to impulse)-type hypersonic wind tunnels operational in this country. AFC facilities are used to assess and optimize the hypersonic aerodynamic performance and aeroheating characteristics of aerospace vehicle concepts and to provide benchmark aerodynamic/aeroheating data fr generating the flight aerodynamic databook and final design of the thermal protection system (TPS) (e.g., establishment of flight limitations not to exceed TPS design limits). Modifications and enhancements of AFC hardware components and instrumentation have been pursued to increase capability, reliability, and productivity in support of programmatic goals. Examples illustrating facility utilization in recent years to generate essentially all of the experimental hypersonic aerodynamic and aeroheating information for high-priority, fast-paced Agency programs are presented. These programs include Phase I of the Reusable Launch Vehicle (RLV) Advanced Technology Demonstrator, X-33 program, PHase II of the X-33 program, X-34 program, the Hyper-X program ( a Mach 5,7, and 10 airbreathing propulsion flight experiment), and the X-38 program (Experimental Crew Return Vehicle, X-CRV). Current upgrades/enchancements and future plans for the AFC are discussed.
An efficient upwind/relaxation algorithm for the Euler and Navier-Stokes equations
NASA Astrophysics Data System (ADS)
Zha, Ge-Cheng; Liu, Dao-Zhi; Ma, Tie-You
1989-05-01
An efficient Euler and full Navier-Stokes solver based on a flux splitting scheme is presented. The original Van Leer flux vector splitting form is extended to arbitrary body-fitted co-ordinates in the physical domain so that it can be used with a finite volume scheme. The block matrix is inverted by Gauss-Seidel iteration. It is verified that the often used reflection boundary condition will produce incorrect flux crossing the wall and cause too large numerical dissipation if flux vector splitting is used. To remove such errors, an appropriate tretment of wall boundary conditions is suggested. Inviscid and viscous steady transonic internal flows are analyzed, including the case of shock-induced boundary layer separation.
Grid-independent upwind scheme for multidimensional flow
NASA Technical Reports Server (NTRS)
Parpia, Ijaz H.; Michalek, Donna J.
1993-01-01
Recent advances in the development of a grid-independent finite volume scheme for the Euler equations of gas dynamics are described. In the proposed method, flowfield gradient data are reconstructed locally (on a triangle) using five elementary planar waves, and an upwind numerical flux function for grid-oblique waves is developed to model the effect of the passage of these waves on the data in a cell. Numerical examples for several two-dimensional test problems are included. These results show high wave resolution and nearly monotone strong-wave transitions.
The upwind control volume scheme for unstructured triangular grids
NASA Technical Reports Server (NTRS)
Giles, Michael; Anderson, W. Kyle; Roberts, Thomas W.
1989-01-01
A new algorithm for the numerical solution of the Euler equations is presented. This algorithm is particularly suited to the use of unstructured triangular meshes, allowing geometric flexibility. Solutions are second-order accurate in the steady state. Implementation of the algorithm requires minimal grid connectivity information, resulting in modest storage requirements, and should enhance the implementation of the scheme on massively parallel computers. A novel form of upwind differencing is developed, and is shown to yield sharp resolution of shocks. Two new artificial viscosity models are introduced that enhance the performance of the new scheme. Numerical results for transonic airfoil flows are presented, which demonstrate the performance of the algorithm.
Multigrid properties of upwind-biased data reconstructions
NASA Technical Reports Server (NTRS)
Warren, Gary P.; Roberts, Thomas W.
1993-01-01
The multigrid properties of two data reconstruction methods used for achieving second-order spatial accuracy when solving the two-dimensional Euler equations are examined. The data reconstruction methods are used with an implicit upwind algorithm which uses linearized backward-Euler time-differencing. The solution of the resulting linear system is performed by an iterative procedure. In the present study only regular quadrilateral grids are considered, so a red-black Gauss-Seidel iteration is used. Although the Jacobian is approximated by first-order upwind extrapolation, two alternative data reconstruction techniques for the flux integral that yield higher-order spatial accuracy at steady state are examined. The first method, probably most popular for structured quadrilateral grids, is based on estimating the cell gradients using one-dimensional reconstruction along curvilinear coordinates. The second method is based on Green's theorem. Analysis and numerical results for the two dimensional Euler equations show that data reconstruction based on Green's theorem has superior multigrid properties as compared to the one-dimensional data reconstruction method.
Comments on the Diffusive Behavior of Two Upwind Schemes
NASA Technical Reports Server (NTRS)
Wood, William A.; Kleb, William L.
1998-01-01
The diffusive characteristics of two upwind schemes, multi-dimensional fluctuation splitting and locally one-dimensional finite volume, are compared for scalar advection-diffusion problems. Algorithms for the two schemes are developed for node-based data representation on median-dual meshes associated with unstructured triangulations in two spatial dimensions. Four model equations are considered: linear advection, non-linear advection, diffusion, and advection-diffusion. Modular coding is employed to isolate the effects of the two approaches for upwind flux evaluation, allowing for head-to-head accuracy and efficiency comparisons. Both the stability of compressive limiters and the amount of artificial diffusion generated by the schemes is found to be grid-orientation dependent, with the fluctuation splitting scheme producing less artificial diffusion than the finite volume scheme. Convergence rates are compared for the combined advection-diffusion problem, with a speedup of 2.5 seen for fluctuation splitting versus finite volume when solved on the same mesh. However, accurate solutions to problems with small diffusion coefficients can be achieved on coarser meshes using fluctuation splitting rather than finite volume, so that when comparing convergence rates to reach a given accuracy, fluctuation splitting shows a speedup of 29 over finite volume.
Diffusion Characteristics of Upwind Schemes on Unstructured Triangulations
NASA Technical Reports Server (NTRS)
Wood, William A.; Kleb, William L.
1998-01-01
The diffusive characteristics of two upwind schemes, multi-dimensional fluctuation splitting and dimensionally-split finite volume, are compared for scalar advection-diffusion problems. Algorithms for the two schemes are developed for node-based data representation on median-dual meshes associated with unstructured triangulations in two spatial dimensions. Four model equations are considered: linear advection, non-linear advection, diffusion, and advection-diffusion. Modular coding is employed to isolate the effects of the two approaches for upwind flux evaluation, allowing for head-to-head accuracy and efficiency comparisons. Both the stability of compressive limiters and the amount of artificial diffusion generated by the schemes is found to be grid-orientation dependent, with the fluctuation splitting scheme producing less artificial diffusion than the dimensionally-split finite volume scheme. Convergence rates are compared for the combined advection-diffusion problem, with a speedup of 2-3 seen for fluctuation splitting versus finite volume when solved on the same mesh. However, accurate solutions to problems with small diffusion coefficients can be achieved on coarser meshes using fluctuation splitting rather than finite volume, so that when comparing convergence rates to reach a given accuracy, fluctuation splitting shows a 20-25 speedup over finite volume.
Upwind schemes and bifurcating solutions in real gas computations
NASA Technical Reports Server (NTRS)
Suresh, Ambady; Liou, Meng-Sing
1992-01-01
The area of high speed flow is seeing a renewed interest due to advanced propulsion concepts such as the National Aerospace Plane (NASP), Space Shuttle, and future civil transport concepts. Upwind schemes to solve such flows have become increasingly popular in the last decade due to their excellent shock capturing properties. In the first part of this paper the authors present the extension of the Osher scheme to equilibrium and non-equilibrium gases. For simplicity, the source terms are treated explicitly. Computations based on the above scheme are presented to demonstrate the feasibility, accuracy and efficiency of the proposed scheme. One of the test problems is a Chapman-Jouguet detonation problem for which numerical solutions have been known to bifurcate into spurious weak detonation solutions on coarse grids. Results indicate that the numerical solution obtained depends both on the upwinding scheme used and the limiter employed to obtain second order accuracy. For example, the Osher scheme gives the correct CJ solution when the super-bee limiter is used, but gives the spurious solution when the Van Leer limiter is used. With the Roe scheme the spurious solution is obtained for all limiters.
An Implicit Upwind Algorithm for Computing Turbulent Flows on Unstructured Grids
NASA Technical Reports Server (NTRS)
Anerson, W. Kyle; Bonhaus, Daryl L.
1994-01-01
An implicit, Navier-Stokes solution algorithm is presented for the computation of turbulent flow on unstructured grids. The inviscid fluxes are computed using an upwind algorithm and the solution is advanced in time using a backward-Euler time-stepping scheme. At each time step, the linear system of equations is approximately solved with a point-implicit relaxation scheme. This methodology provides a viable and robust algorithm for computing turbulent flows on unstructured meshes. Results are shown for subsonic flow over a NACA 0012 airfoil and for transonic flow over a RAE 2822 airfoil exhibiting a strong upper-surface shock. In addition, results are shown for 3 element and 4 element airfoil configurations. For the calculations, two one equation turbulence models are utilized. For the NACA 0012 airfoil, a pressure distribution and force data are compared with other computational results as well as with experiment. Comparisons of computed pressure distributions and velocity profiles with experimental data are shown for the RAE airfoil and for the 3 element configuration. For the 4 element case, comparisons of surface pressure distributions with experiment are made. In general, the agreement between the computations and the experiment is good.
Development and application of computational aerothermodynamics flowfield computer codes
NASA Technical Reports Server (NTRS)
Venkatapathy, Ethiraj
1992-01-01
Presented is a collection of papers on research activities carried out during the funding period of October 1991 to March 1992. Topics covered include: blunt body flows in thermochemical equilibrium; thermochemical relaxation in high enthalpy nozzle flow; single expansion ramp nozzle simulations; lunar return aerobraking; line boundary problem for three dimensional grids; and unsteady shock induced combustion.
Development of X-33/X-34 Aerothermodynamic Data Bases: Lessons Learned and Future Enhancements
NASA Technical Reports Server (NTRS)
Miller, C. G.
2000-01-01
A synoptic of programmatic and technical lessons learned in the development of aerothermodynamic data bases for the X-33 and X-34 programs is presented in general terms and from the perspective of the NASA Langley Research Center Aerothermodynamics Branch. The format used is that of the "aerothermodynamic chain," the links of which are personnel, facilities, models/test articles, instrumentation, test techniques, and computational fluid dynamics (CFD). Because the aerodynamic data bases upon which the X-33 and X-34 vehicles will fly are almost exclusively from wind tunnel testing, as opposed to CFD, the primary focus of the lessons learned is on ground-based testing. The period corresponding to the development of X-33 and X-34 aerothermodynamic data bases was challenging, since a number of other such programs (e.g., X-38, X-43) competed for resources at a time of downsizing of personnel, facilities, etc., outsourcing, and role changes as NASA Centers served as subcontractors to industry. The impact of this changing environment is embedded in the lessons learned. From a technical perspective, the relatively long times to design and fabricate metallic force and moment models, delays in delivery of models, and a lack of quality assurance to determine the fidelity of model outer mold lines (OML) prior to wind tunnel testing had a major negative impact on the programs. On the positive side, the application of phosphor thermography to obtain global, quantitative heating distributions on rapidly fabricated ceramic models revolutionized the aerothermodynamic optimization of vehicle OMLs, control surfaces, etc. Vehicle designers were provided with aeroheating information prior to, or in conjunction with, aerodynamic information early in the program, thereby allowing trades to be made with both sets of input; in the past only aerodynamic data were available as input. Programmatically, failure to include transonic aerodynamic wind tunnel tests early in the assessment phase
An Upwind Solver for the National Combustion Code
NASA Technical Reports Server (NTRS)
Sockol, Peter M.
2011-01-01
An upwind solver is presented for the unstructured grid National Combustion Code (NCC). The compressible Navier-Stokes equations with time-derivative preconditioning and preconditioned flux-difference splitting of the inviscid terms are used. First order derivatives are computed on cell faces and used to evaluate the shear stresses and heat fluxes. A new flux limiter uses these same first order derivatives in the evaluation of left and right states used in the flux-difference splitting. The k-epsilon turbulence equations are solved with the same second-order method. The new solver has been installed in a recent version of NCC and the resulting code has been tested successfully in 2D on two laminar cases with known solutions and one turbulent case with experimental data.
Analysis of implicit second-order upwind-biased stencils
NASA Technical Reports Server (NTRS)
Roberts, Thomas W.; Warren, Gary P.
1993-01-01
Truncation error and stability properties of several implicit upwind schemes for the two-dimensional Euler equations are examined. The schemes use linear data reconstruction methods to achieve second-order flux integrations where the implicit Jacobian operators are first order. The stability properties of the schemes are examined by a Von Neumann analysis of the linearized, constant-coefficient Euler equations. The choice of the data reconstruction method used to evaluate the flux integral has a dramatic effect on the convergence properties of the implicit solution method. In particular, the typical one-dimensional data reconstruction methods used with structured grids exhibit poor convergence properties compared to the unstructured grid method considered. Of the schemes examined, the one with the superior convergence properties is well-suited for both unstructured and structured grids, which has important implications for the design of implicit methods.
Factorizable Upwind Schemes: The Triangular Unstructured Grid Formulation
NASA Technical Reports Server (NTRS)
Sidilkover, David; Nielsen, Eric J.
2001-01-01
The upwind factorizable schemes for the equations of fluid were introduced recently. They facilitate achieving the Textbook Multigrid Efficiency (TME) and are expected also to result in the solvers of unparalleled robustness. The approach itself is very general. Therefore, it may well become a general framework for the large-scale, Computational Fluid Dynamics. In this paper we outline the triangular grid formulation of the factorizable schemes. The derivation is based on the fact that the factorizable schemes can be expressed entirely using vector notation. without explicitly mentioning a particular coordinate frame. We, describe the resulting discrete scheme in detail and present some computational results verifying the basic properties of the scheme/solver.
Factorizable Upwind Schemes: the Triangular Unstructured Grid Formulation
NASA Technical Reports Server (NTRS)
Nielsen, Eric J.; Sidilkover, David
2001-01-01
The upwind factorizable schemes for the equations of fluid was introduced recently. They facilitate achieving the Textbook Multigrid Efficiency (TME) and are expected also to result in the solvers of unparalleled robustness. The approach itself is very general. Therefore, it may well become a general framework for the large-scale Computational Fluid Dynamics. In this paper we outline the triangular grid formulation of the factorizable schemes. The derivation is based on the fact that the factorizable schemes can be expressed entirely using vector notation, without explicitly mentioning a particular coordinate frame. We describe the resulting discrete scheme in detail and present some computational results verifying the basic properties of the scheme/solver.
High upwind concentrations observed during an upslope tracer event
Ciolek, J.T. Jr.
1993-10-01
In February of 1991 the Rocky Flats Plant conducted twelve tracer experiments to validate an emergency response dispersion model known as the Terrain-Responsive Atmospheric Code (TRAC) (Hodgin 1985). Experimenters released 140 to 260 kilograms of inert tracer gas (sulfur hexafloride) from the plant over an 11 hour period. During each release, one hundred and sixty-five samples, most of which formed concentric rings of 8 and 16 km radius from the plant, recorded cumulative hourly concentrations of the tracer at one meter above ground level (AGL). Figure 1 contains a depiction of the sampler location, the terrain, and the meteorological stations available within the tracer study area. Brown (1991) describes the experimental setup in more detail. The subject of this paper is an event that occurred early in the fifth experiment, on February 9, 1991. In this experiment, tracer material released from 13:00 to 17:00 LST appeared both downwind and upwind of the source, with the highest concentrations upwind. During the fifth experiment, high pressure in Utah produced mostly sunny skis around Rocky Flats. For most of the day, one could find moderate (5 to 10 ms{sup {minus}1}) northerly (from the North) flow within the 700 to 500 mb level of the atmosphere (approximately 3000 to 5500 meters above Mean Sea Level (MSL)). Synoptic scale motions were isolated enough from the surface layer and heating was great enough to produce a 1 km deep upslope flow (flow from the East to the West) by late afternoon. The winds reversed and became downslope at approximately 17:30 LST.
Comparative study of upwind schemes for transonic and supersonic internal flows
NASA Astrophysics Data System (ADS)
Niculescu, M. L.; Dǎnǎilǎ, S.
2013-10-01
A study of some popular upwind schemes applied to transonic internal flows using some well-known test cases is done in this paper. We focused on upwind schemes because the central space discretizations have symmetry with respect to a change in sign; therefore, the physical propagation of perturbations along characteristics, typical of hyperbolic equations is not considered in the definition of numerical model. In contrast to the central space discretizations, the upwind schemes whose origin may be due to Courant et al. [1] are directed towards the introduction of the physical properties of the flow equations into the discretized formulation that has leads to upwinding techniques such as flux vector splitting and flux difference splitting. In order to test the accuracy, robustness and efficiency of some popular upwind methods (van Leer scheme, Roe scheme and Liou's AUSM+ scheme); we used some well-known test cases.
Overview of X-38 Hypersonic Aerothermodynamic Wind Tunnel Data and Comparison with Numerical Results
NASA Technical Reports Server (NTRS)
Campbell, C.; Caram, J.; Berry, S.; Horvath, T.; Merski, N.; Loomis, M.; Venkatapathy, E.
2004-01-01
A NASA team of engineers has been organized to design a crew return vehicle for returning International Space Station crew members from orbit. The hypersonic aerothermodynamic characteristics of the X-23/X-24A derived X-38 crew return vehicle are being evaluated in various wind tunnels in support of this effort. Aerothermodynamic data from two NASA hypersonic tunnels at Mach 6 and Mach 10 has been obtained with cast ceramic models and a thermographic phosphorus digital imaging system. General windward surface heating features are described based on experimental surface heating images and surface oil flow patterns for the nominal hypersonic aerodynamic orientation. Body flap reattachment heating levels are examined. Computational Fluid Dynamics tools have been applied at the appropriate wind tunnel conditions to make comparisons with this data.
Numerical Simulations Of High-Altitude Aerothermodynamics Of A Prospective Spacecraft Model
NASA Astrophysics Data System (ADS)
Vashchenkov, P. V.; Kaskovsky, A. V.; Krylov, A. N.; Ivanov, M. S.
2011-05-01
The paper describes the computations of aerothermodynamic characteristics of a promising spacecraft (Prospective Piloted Transport System) along its de- scent trajectory at altitudes from 120 to 60 km. The computations are performed by the DSMC method with the use of the SMILE software system and by the engineering technique (local bridging method) with the use of the RuSat software system. The influence of real gas effects (excitation of rotational and vibrational energy modes and chemical reactions) on aerothermodynamic characteristics of the vehicle is studied. A comparison of results obtained by the approximate engineering method and the DSMC method allow the accuracy of prediction of aerodynamic characteristics by the local bridging method to be estimated.
Wind-tunnel based definition of the AFE aerothermodynamic environment. [Aeroassist Flight Experiment
NASA Technical Reports Server (NTRS)
Miller, Charles G.; Wells, W. L.
1992-01-01
The Aeroassist Flight Experiment (AFE), scheduled to be performed in 1994, will serve as a precursor for aeroassisted space transfer vehicles (ASTV's) and is representative of entry concepts being considered for missions to Mars. Rationale for the AFE is reviewed briefly as are the various experiments carried aboard the vehicle. The approach used to determine hypersonic aerodynamic and aerothermodynamic characteristics over a wide range of simulation parameters in ground-based facilities is presented. Facilities, instrumentation and test procedures employed in the establishment of the data base are discussed. Measurements illustrating the effects of hypersonic simulation parameters, particularly normal-shock density ratio (an important parameter for hypersonic blunt bodies), and attitude on aerodynamic and aerothermodynamic characteristics are presented, and predictions from computational fluid dynamic (CFD) computer codes are compared with measurement.
NASA Technical Reports Server (NTRS)
Hollis, Brian R.
1996-01-01
A computational algorithm has been developed which can be employed to determine the flow properties of an arbitrary real (virial) gas in a wind tunnel. A multiple-coefficient virial gas equation of state and the assumption of isentropic flow are used to model the gas and to compute flow properties throughout the wind tunnel. This algorithm has been used to calculate flow properties for the wind tunnels of the Aerothermodynamics Facilities Complex at the NASA Langley Research Center, in which air, CF4. He, and N2 are employed as test gases. The algorithm is detailed in this paper and sample results are presented for each of the Aerothermodynamic Facilities Complex wind tunnels.
The use of the Tethered Satellite System to perform low density aerothermodynamics studies
NASA Technical Reports Server (NTRS)
Carlomagno, Giovanni M.; De Luca, Luigi; Siemers, Paul M.; Wood, George M., Jr.
1988-01-01
The Tethered Satellite System (TSS) is a cooperative space system development activity being carried out by USA and Italy. Within TSS, the Shuttle Tethered Aerothermodynamic Research Facility (STARFAC) concept has the potential to provide access to vast portions of the upper atmosphere for the purpose of atmospheric and aerothermodynamic research. The implementation of this capability will push Tether System (TS) state of the art to its limits; the primary problems being tether/satellite drag, heating, tension control, deployment/retrieval control. In this paper parametric studies are accomplished to assess some of these problems and to delineate the tradeoffs available to missions design to meet the engineering constraints. The utilization of aerodynamic rather than spherical shapes - (TSS) - as well as elementary satellite thrusting and lift are included in the present study.
Aerothermodynamic flow phenomena of the airframe-integrated supersonic combustion ramjet
NASA Astrophysics Data System (ADS)
Walton, James T.
1992-11-01
The unique component flow phenomena is discussed of the airframe-integrated supersonic combustion ramjet (scramjet) in a format geared towards new players in the arena of hypersonic propulsion. After giving an overview of the scramjet aerothermodynamic cycle, the characteristics are then covered individually of the vehicle forebody, inlet, combustor, and vehicle afterbody/nozzle. Attention is given to phenomena such as inlet speeding, inlet starting, inlet spillage, fuel injection, thermal choking, and combustor-inlet interaction.
NASA Technical Reports Server (NTRS)
Andrews, E. H., Jr.; Mackley, E. A.
1976-01-01
The NASA Hypersonic Research Engine (HRE) Project was initiated for the purpose of advancing the technology of airbreathing propulsion for hypersonic flight. A large component (inlet, combustor, and nozzle) and structures development program was encompassed by the project. The tests of a full-scale (18 in. diameter cowl and 87 in. long) HRE concept, designated the Aerothermodynamic Integration Model (AIM), at Mach numbers of 5, 6, and 7. Computer program results for Mach 6 component integration tests are presented.
Aerothermodynamic flow phenomena of the airframe-integrated supersonic combustion ramjet
NASA Technical Reports Server (NTRS)
Walton, James T.
1992-01-01
The unique component flow phenomena is discussed of the airframe-integrated supersonic combustion ramjet (scramjet) in a format geared towards new players in the arena of hypersonic propulsion. After giving an overview of the scramjet aerothermodynamic cycle, the characteristics are then covered individually of the vehicle forebody, inlet, combustor, and vehicle afterbody/nozzle. Attention is given to phenomena such as inlet speeding, inlet starting, inlet spillage, fuel injection, thermal choking, and combustor-inlet interaction.
Hypersonic research engine project. Phase 2: Aerothermodynamic Integration Model (AIM) test report
NASA Technical Reports Server (NTRS)
Andersen, W. L.; Kado, L.
1975-01-01
The Hypersonic Research Engine-Aerothermodynamic Integration Model (HRE-AIM) was designed, fabricated, and tested in the Hypersonic Tunnel Facility. The HRE-AIM is described along with its installation in the wind tunnel facility. Test conditions to which the HRE-AIM was subjected and observations made during the tests are discussed. The overall engine performance, component interaction, and ignition limits for the design are evaluated.
Impact of ETO propellants on the aerothermodynamic analyses of propulsion components
NASA Technical Reports Server (NTRS)
Civinskas, K. C.; Boyle, R. J.; Mcconnaughey, H. V.
1988-01-01
The operating conditions and the propellant transport properties used in Earth-to-Orbit (ETO) applications affect the aerothermodynamic design of ETO turbomachinery in a number of ways. Some aerodynamic and heat transfer implications of the low molecular weight fluids and high Reynolds number operating conditions on future ETO turbomachinery are discussed. Using the current SSME high pressure fuel turbine as a baseline, the aerothermodynamic comparisons are made for two alternate fuel turbine geometries. The first is a revised first stage rotor blade designed to reduce peak heat transfer. This alternate design resulted in a 23 percent reduction in peak heat transfer. The second design concept was a single stage rotor to yield the same power output as the baseline two stage rotor. Since the rotor tip speed was held constant, the turbine work factor doubled. In this alternate design, the peak heat transfer remained the same as the baseline. While the efficiency of the single stage design was 3.1 points less than the baseline two stage turbine, the design was aerothermodynamically feasible, and may be structurally desirable.
NASA Astrophysics Data System (ADS)
Verand, Jean-Luc; Pelissier, Christian; Sourgen, Frederic; Fontaine, Joelle; Garcon, Francois; Spel, Martin; van Hauwaert, Pierre; Charbonnier, Dominique; Vos, Jan; Vallee, Jean-Jacques; Pibarot, Julien; Tribot, Jean-Pierre; Mareschi, Vincenzo; Ferrarella, Daniella; Rufolo, Giuseppe
2011-05-01
The Intermediate eXperimental Vehicle (IXV) project objectives are the design, development, manufacture and ground and flight verification of an autonomous European lifting and aerodynamically controlled re-entry system, which is highly flexible and manoeuvrable. The IXV vehicle is a flying test bed for securing the next step of operational space vehicle development by supporting technology demonstration and system concept through the following objectives: a) Aerothermodynamics b) Advanced In Flight Experiments c) Thermal Protection System d) Guidance Navigation and Control e) System design The assessment of the general aerothermodynamic environment of IXV vehicle is mainly performed considering a smooth simplified geometry. However, the thermal protection system of IXV includes a mono-block ceramic matrix composite nose and an assembly of shingles between which steps and gaps are generated. From an aerothermodynamic point of view, such a distributed roughness layout cannot be ignored in terms of modification of the interaction between the flow and the body. To assess this effect, dedicated Mach number 5.5 wind tunnel tests (ONERA, S3MA facility) and numerical simulations (RTECH and CFS Engineering) have been performed during the phase C2 of the project. The paper presents the general logic of the work, with emphasis on the wind tunnel model design, tests involving infrared thermal measurements as well as the CFD rebuilding of the flow in the wind tunnel and the extrapolation from ground-to-flight.
PRE_X Programme: Aerothermodynamic Objectives and Aeroshape Definition for in Flight Experiments
NASA Astrophysics Data System (ADS)
Lambert, O.; Tribot, J.-P.; Saint-Cloud, F.
2002-01-01
As the expendable launch vehicles (ELV) are limited in their trend to lower costs, the reusability (Reusable Launch Vehicle, RLV) could be the way to make drastic step. By the year 2001, CNES proposed through the ANGEL phase 1 programme to preprare the required technical maturity before that RLV's become alternatives to ELV's. In such way, system ,propulsion, ground based demonstrations, aero-thermo-dynamics as well as in flight experimentation are planned. This paper is focused on the aero-thermo-dynamics (ATD) and in flight demonstration activities with emphasis on the better understanding of ATD problems emerging from past programmes among them shock wave transitionnal boundary layer interaction on surface control, boundary layer transition, local aerothermodynamic effects, gas- surface interaction, catalycity, base flow prediction,...In order to minimize as small as possible the management risk a first generation of vehicle dubbed Pre_X is designed to validate technological choices and to have as soon as possible re-entry data to calibrate the various tools involved in the future RLV definition. In addition, the main requirement for PRE_X aeroshape definition and the two different design approaches considered by Dassault Aviation and EADS-LV are discussed. Then, the more promising concept for the PRE_X application is presented. Finally, the current status of the ATD activities is given as well as the perspectives.
ERIC Educational Resources Information Center
Ramey, Kyle
2004-01-01
Relaxed intensity refers to a professional philosophy, demeanor, and way of life. It is the key to being an effective educational leader. To be successful one must be relaxed, which means managing stress efficiently, having fun, and enjoying work. Intensity allows one to get the job done and accomplish certain tasks or goals. Educational leaders…
Numerical solution of 3D Navier-Stokes equations with upwind implicit schemes
NASA Technical Reports Server (NTRS)
Marx, Yves P.
1990-01-01
An upwind MUSCL type implicit scheme for the three-dimensional Navier-Stokes equations is presented. Comparison between different approximate Riemann solvers (Roe and Osher) are performed and the influence of the reconstructions schemes on the accuracy of the solution as well as on the convergence of the method is studied. A new limiter is introduced in order to remove the problems usually associated with non-linear upwind schemes. The implementation of a diagonal upwind implicit operator for the three-dimensional Navier-Stokes equations is also discussed. Finally the turbulence modeling is assessed. Good prediction of separated flows are demonstrated if a non-equilibrium turbulence model is used.
Applying Upwind Godunov Methods to Calculate Two—Phase Mixture Conservation Laws
NASA Astrophysics Data System (ADS)
Zeidan, D.
2010-09-01
This paper continues a previous work (ICNAAM 2009; AIP Conference Proceedings, 1168, 601-604) on solving a hyperbolic conservative model for compressible gas—solid mixture flow using upwind Godunov methods. The numerical resolution of the model from Godunov first—order upwind and MUSCL—Hancock methods are reported. Both methods are based on the HLL Riemann solver in the framework of finite volume techniques. Calculation results are presented for a series of one—dimensional test problems. The results show that upwind Godunov methods are accurate and robust enough for two—phase mixture conservation laws.
NASA Technical Reports Server (NTRS)
Bose, Deepak
2012-01-01
The design of entry vehicles requires predictions of aerothermal environment during the hypersonic phase of their flight trajectories. These predictions are made using computational fluid dynamics (CFD) codes that often rely on physics and chemistry models of nonequilibrium processes. The primary processes of interest are gas phase chemistry, internal energy relaxation, electronic excitation, nonequilibrium emission and absorption of radiation, and gas-surface interaction leading to surface recession and catalytic recombination. NASAs Hypersonics Project is advancing the state-of-the-art in modeling of nonequilibrium phenomena by making detailed spectroscopic measurements in shock tube and arcjets, using ab-initio quantum mechanical techniques develop fundamental chemistry and spectroscopic databases, making fundamental measurements of finite-rate gas surface interactions, implementing of detailed mechanisms in the state-of-the-art CFD codes, The development of new models is based on validation with relevant experiments. We will present the latest developments and a roadmap for the technical areas mentioned above
Development and application of computational aerothermodynamics flowfield computer codes
NASA Technical Reports Server (NTRS)
Venkatapathy, Ethiraj
1993-01-01
Computations are presented for one-dimensional, strong shock waves that are typical of those that form in front of a reentering spacecraft. The fluid mechanics and thermochemistry are modeled using two different approaches. The first employs traditional continuum techniques in solving the Navier-Stokes equations. The second-approach employs a particle simulation technique (the direct simulation Monte Carlo method, DSMC). The thermochemical models employed in these two techniques are quite different. The present investigation presents an evaluation of thermochemical models for nitrogen under hypersonic flow conditions. Four separate cases are considered. The cases are governed, respectively, by the following: vibrational relaxation; weak dissociation; strong dissociation; and weak ionization. In near-continuum, hypersonic flow, the nonequilibrium thermochemical models employed in continuum and particle simulations produce nearly identical solutions. Further, the two approaches are evaluated successfully against available experimental data for weakly and strongly dissociating flows.
The design of high-resolution upwind shock-capturing methods
Rider, W.J.
1992-06-01
The design and construction of high-resolution upwind shock-capturing methods is an effective means of solving conservation laws of physics numerically. In the past, the design of such methods was generally categorized into several distinct methods. This work shows how these methods can be viewed in a unified manner. A generalized flux-corrected transport (FCT) algorithm is shown to be total variation diminishing (TVD) under some conditions. The new algorithm has improved properties from the standpoint of use and analysis. This is shown to be especially important for systems of equations. The new formulation allows Riemann solvers to be used effectively with FCT methods. This directly leads to a geometric analog to symmetric TVD and FCT methods that is developed and expanded upon. This unifies these methods with high-order Godunov (HOG) methods. Two new variants of this are derived and shown to be uniformly non-oscillatory. Limiters are an effective means of designing these types of methods. More general classes of limiters are discussed with extensions to a wider class of schemes. Flux-corrected transport and total variation bounded (TVB) limiters are discussed, modified, and expanded. Two new classes of limiters are described: s-limiters and generalized average limiters. The recently defined ULTIMATE limiter is analyzed within the framework of the other limiters. The benefits of relaxing strict constraints on the limiters such as TVD requirements are also discussed. For coarse grids, limiters such as the TVB and the generalized average with bias improve resolution considerably. Lastly, the question of whether the polynomial reconstruction technique used in a HOG method should be based on cell-averages or point-values is studied.
Extension of a streamwise upwind algorithm to a moving grid system
NASA Technical Reports Server (NTRS)
Obayashi, Shigeru; Goorjian, Peter M.; Guruswamy, Guru P.
1990-01-01
A new streamwise upwind algorithm was derived to compute unsteady flow fields with the use of a moving-grid system. The temporally nonconservative LU-ADI (lower-upper-factored, alternating-direction-implicit) method was applied for time marching computations. A comparison of the temporally nonconservative method with a time-conservative implicit upwind method indicates that the solutions are insensitive to the conservative properties of the implicit solvers when practical time steps are used. Using this new method, computations were made for an oscillating wing at a transonic Mach number. The computed results confirm that the present upwind scheme captures the shock motion better than the central-difference scheme based on the beam-warming algorithm. The new upwind option of the code allows larger time-steps and thus is more efficient, even though it requires slightly more computational time per time step than the central-difference option.
Solving Upwind-Biased Discretizations: Defect-Correction Iterations
NASA Technical Reports Server (NTRS)
Diskin, Boris; Thomas, James L.
1999-01-01
This paper considers defect-correction solvers for a second order upwind-biased discretization of the 2D convection equation. The following important features are reported: (1) The asymptotic convergence rate is about 0.5 per defect-correction iteration. (2) If the operators involved in defect-correction iterations have different approximation order, then the initial convergence rates may be very slow. The number of iterations required to get into the asymptotic convergence regime might grow on fine grids as a negative power of h. In the case of a second order target operator and a first order driver operator, this number of iterations is roughly proportional to h-1/3. (3) If both the operators have the second approximation order, the defect-correction solver demonstrates the asymptotic convergence rate after three iterations at most. The same three iterations are required to converge algebraic error below the truncation error level. A novel comprehensive half-space Fourier mode analysis (which, by the way, can take into account the influence of discretized outflow boundary conditions as well) for the defect-correction method is developed. This analysis explains many phenomena observed in solving non-elliptic equations and provides a close prediction of the actual solution behavior. It predicts the convergence rate for each iteration and the asymptotic convergence rate. As a result of this analysis, a new very efficient adaptive multigrid algorithm solving the discrete problem to within a given accuracy is proposed. Numerical simulations confirm the accuracy of the analysis and the efficiency of the proposed algorithm. The results of the numerical tests are reported.
Simulation of upwind maneuvering of a sailing yacht
NASA Astrophysics Data System (ADS)
Harris, Daniel Hartrick
A time domain maneuvering simulation of an IACC class yacht suitable for the analysis of unsteady upwind sailing including tacking is presented. The simulation considers motions in six degrees of freedom. The hydrodynamic and aerodynamic loads are calculated primarily with unsteady potential theory supplemented by empirical viscous models. The hydrodynamic model includes the effects of incident waves. Control of the rudder is provided by a simple rate feedback autopilot which is augmented with open loop additions to mimic human steering. The hydrodynamic models are based on the superposition of force components. These components fall into two groups, those which the yacht will experience in calm water, and those due to incident waves. The calm water loads are further divided into zero Froude number, or "double body" maneuvering loads, hydrostatic loads, gravitational loads, free surface radiation loads, and viscous/residual loads. The maneuvering loads are calculated with an unsteady panel code which treats the instantaneous geometry of the yacht below the undisturbed free surface. The free surface radiation loads are calculated via convolution of impulse response functions derived from seakeeping strip theory. The viscous/residual loads are based upon empirical estimates. The aerodynamic model consists primarily of a database of steady state sail coefficients. These coefficients treat the individual contributions to the total sail force of a number of chordwise strips on both the main and jib. Dynamic effects are modeled by using the instantaneous incident wind velocity and direction as the independent variables for the sail load contribution of each strip. The sail coefficient database was calculated numerically with potential methods and simple empirical viscous corrections. Additional aerodynamic load calculations are made to determine the parasitic contributions of the rig and hull. Validation studies compare the steady sailing hydro and aerodynamic loads
A streamwise upwind algorithm applied to vortical flow over a delta wing
NASA Technical Reports Server (NTRS)
Goorjian, Peter M.; Obayashi, Shigeru
1989-01-01
Improvements were made to a streamwise upwind algorithm so that it can be used for calculating flows with vortices. A calculation is shown of flow over a delta wing at an angle of attack. The laminar, thin layer, Navier-Stokes equations are used for the calculation. The results are compared with another upwind method, a central differencing method, and experimental data. The present method shows improvements in accuracy and convergence properties.
NASA Technical Reports Server (NTRS)
Kamhawi, Hilmi N.
2011-01-01
This report documents the work performed during from March 2010 October 2011. The Integrated Design and Engineering Analysis (IDEA) environment is a collaborative environment based on an object-oriented, multidisciplinary, distributed environment using the Adaptive Modeling Language (AML) as the underlying framework. This report will focus on describing the work done in the area of extending the aerodynamics, and aerothermodynamics module using S/HABP, CBAERO, PREMIN and LANMIN. It will also detail the work done integrating EXITS as the TPS sizing tool.
NASA Technical Reports Server (NTRS)
Salas, Manuel D.
2007-01-01
The research program of the aerodynamics, aerothermodynamics and plasmadynamics discipline of NASA's Hypersonic Project is reviewed. Details are provided for each of its three components: 1) development of physics-based models of non-equilibrium chemistry, surface catalytic effects, turbulence, transition and radiation; 2) development of advanced simulation tools to enable increased spatial and time accuracy, increased geometrical complexity, grid adaptation, increased physical-processes complexity, uncertainty quantification and error control; and 3) establishment of experimental databases from ground and flight experiments to develop better understanding of high-speed flows and to provide data to validate and guide the development of simulation tools.
An inlet analysis for the NASA hypersonic research engine aerothermodynamic integration model
NASA Technical Reports Server (NTRS)
Andrews, E. H., Jr.; Russell, J. W.; Mackley, E. A.; Simmonds, A. L.
1974-01-01
A theoretical analysis for the inlet of the NASA Hypersonic Research Engine (HRE) Aerothermodynamic Integration Model (AIM) has been undertaken by use of a method-of-characteristics computer program. The purpose of the analysis was to obtain pretest information on the full-scale HRE inlet in support of the experimental AIM program (completed May 1974). Mass-flow-ratio and additive-drag-coefficient schedules were obtained that well defined the range effected in the AIM tests. Mass-weighted average inlet total-pressure recovery, kinetic energy efficiency, and throat Mach numbers were obtained.
NASA Technical Reports Server (NTRS)
1987-01-01
Environ Corporation's relaxation system is built around a body lounge, a kind of super easy chair that incorporates sensory devices. Computer controlled enclosure provides filtered ionized air to create a feeling of invigoration, enhanced by mood changing aromas. Occupant is also surrounded by multidimensional audio and the lighting is programmed to change colors, patterns, and intensity periodically. These and other sensory stimulators are designed to provide an environment in which the learning process is stimulated, because research has proven that while an individual is in a deep state of relaxation, the mind is more receptive to new information.
NASA Astrophysics Data System (ADS)
Zhou, Hu; Wan, Decheng
2015-03-01
Although the upwind configuration is more popular in the field of wind energy, the downwind one is a promising type for the offshore wind energy due to its special advantages. Different configurations have different aerodynamic performance and it is important to predict the performance of both downwind and upwind configurations accurately for designing and developing more reliable wind turbines. In this paper, a numerical investigation on the aerodynamic performance of National Renewable Energy Laboratory (NREL) phase VI wind turbine in downwind and upwind configurations is presented. The open source toolbox OpenFOAM coupled with arbitrary mesh interface (AMI) method is applied to tackle rotating problems of wind turbines. Two 3D numerical models of NREL phase VI wind turbine with downwind and upwind configurations under four typical working conditions of incoming wind velocities are set up for the study of different unsteady characteristics of the downwind and upwind configurations, respectively. Numerical results of wake vortex structure, time histories of thrust, pressure distribution on the blade and limiting streamlines which can be used to identify points of separation in a 3D flow are presented. It can be concluded that thrust reduction due to blade-tower interaction is small for upwind wind turbines but relatively large for downwind wind turbines and attention should be paid to the vibration at a certain frequency induced by the cyclic reduction for both configurations. The results and conclusions are helpful to analyze the different aerodynamic performance of wind turbines between downwind and upwind configurations, providing useful references for practical design of wind turbine.
NASA Technical Reports Server (NTRS)
Arnold, James O.; Deiwert, George S.
1997-01-01
This paper surveys the use of aerothermodynamic facilities which have been useful in the study of external flows and propulsion aspects of hypersonic, air-breathing vehicles. While the paper is not a survey of all facilities, it covers the utility of shock tunnels and conventional hypersonic blow-down facilities which have been used for hypersonic air-breather studies. The problems confronting researchers in the field of aerothermodynamics are outlined. Results from the T5 GALCIT tunnel for the shock-on lip problem are outlined. Experiments on combustors and short expansion nozzles using the semi-free jet method have been conducted in large shock tunnels. An example which employed the NASA Ames 16-Inch shock tunnel is outlined, and the philosophy of the test technique is described. Conventional blow-down hypersonic wind tunnels are quite useful in hypersonic air-breathing studies. Results from an expansion ramp experiment, simulating the nozzle on a hypersonic air-breather from the NASA Ames 3.5 Foot Hypersonic wind tunnel are summarized. Similar work on expansion nozzles conducted in the NASA Langley hypersonic wind tunnel complex is cited. Free-jet air-frame propulsion integration and configuration stability experiments conducted at Langley in the hypersonic wind tunnel complex on a small generic model are also summarized.
Bourgois, Jan G; Callewaert, Margot; Celie, Bert; De Clercq, Dirk; Boone, Jan
2016-01-01
This study investigates the physiological responses to upwind sailing on a laser emulation ergometer and analyses the components of the physical profile that determine the physiological responses related to sailing level. Ten male high-level laser sailors performed an upwind sailing test, incremental cycling test and quadriceps strength test. During the upwind sailing test, heart rate (HR), oxygen uptake, ventilation, respiratory exchange ratio, rating of perceived exertion (RPE) and lactate concentration were measured, combined with near-infrared spectroscopy (NIRS) and electromyography (EMG) registration of the M. Vastus lateralis. Repeated measures ANOVA showed for the cardio-respiratory, metabolic and muscles responses (mean power frequency [MPF], root mean square [RMS], deoxy[Hb+Mb]) during the upwind sailing test an initial significant increase followed by a stabilisation, despite a constant increase in RPE. Stepwise regression analysis showed that better sailing level was for 46.5% predicted by lower MPF decrease. Lower MPF decrease was for 57.8% predicted by a higher maximal isometric quadriceps strength. In conclusion, this study indicates that higher sailing level was mainly determined by a lower rate of neuromuscular fatigue during the upwind sailing test (as indicated by MPF decrease). Additionally, the level of neuromuscular fatigue was mainly determined by higher maximal isometric quadriceps strength stressing the importance of resistance training in the planning of training. PMID:26323461
NASA Technical Reports Server (NTRS)
Jensen, Douglas R.
1990-01-01
During the months of June and July 1987, the Marine Stratocumulus Intensive Field Observation Experiment of First ISCCP Regional Experiment (FIRE) was conducted in the Southern California offshore area in the vicinity of San Nicolas Island (SNI). The Naval Ocean Systems Center (NOSC) airborne platform was utilized during FIRE to investigate the upwind low level horizontal variability of the marine boundary layer structure to determine the representativeness of SNI-based measurements to upwind open ocean conditions. The NOSC airborne meteorological platform made three flights during FIRE, two during clear sky conditions (19 and 23 July), and one during two stratus conditions (15 July). The boundary layer structure variations associated with the stratus clouds of 15 July 1987 are discussed. Profiles of air temperature (AT) and relative humidity (RH) taken 'at' and 'upwind' of SNI do show differences between the so-called open ocean conditions and those taken near the island. However, the observed difference cannot be uniquely identified to island effects, especially since the upwind fluctuations of AT and RH bound the SNI measurements. Total optical depths measures at SNI do not appear to be greatly affected by any surface based aerosol effects created by the island and could therefore realistically represent open ocean conditions. However, if one were to use the SNI aerosol measurements to predict ship to ship EO propagation conditions, significant errors could be introduced due to the increased number of surface aerosols observed near SNI which may not be, and were not, characteristic of open ocean conditions. Sea surface temperature measurements taken at the island will not, in general, represent those upwind open ocean conditions. Also, since CTT's varied appreciably along the upwind radials, measurements of CTT over the island may not be representative of actual open ocean CTT's.
An Explicit Upwind Algorithm for Solving the Parabolized Navier-Stokes Equations
NASA Technical Reports Server (NTRS)
Korte, John J.
1991-01-01
An explicit, upwind algorithm was developed for the direct (noniterative) integration of the 3-D Parabolized Navier-Stokes (PNS) equations in a generalized coordinate system. The new algorithm uses upwind approximations of the numerical fluxes for the pressure and convection terms obtained by combining flux difference splittings (FDS) formed from the solution of an approximate Riemann (RP). The approximate RP is solved using an extension of the method developed by Roe for steady supersonic flow of an ideal gas. Roe's method is extended for use with the 3-D PNS equations expressed in generalized coordinates and to include Vigneron's technique of splitting the streamwise pressure gradient. The difficulty associated with applying Roe's scheme in the subsonic region is overcome. The second-order upwind differencing of the flux derivatives are obtained by adding FDS to either an original forward or backward differencing of the flux derivative. This approach is used to modify an explicit MacCormack differencing scheme into an upwind differencing scheme. The second order upwind flux approximations, applied with flux limiters, provide a method for numerically capturing shocks without the need for additional artificial damping terms which require adjustment by the user. In addition, a cubic equation is derived for determining Vegneron's pressure splitting coefficient using the updated streamwise flux vector. Decoding the streamwise flux vector with the updated value of Vigneron's pressure splitting improves the stability of the scheme. The new algorithm is applied to 2-D and 3-D supersonic and hypersonic laminar flow test cases. Results are presented for the experimental studies of Holden and of Tracy. In addition, a flow field solution is presented for a generic hypersonic aircraft at a Mach number of 24.5 and angle of attack of 1 degree. The computed results compare well to both experimental data and numerical results from other algorithms. Computational times required
Nguyen, B.T.; Hutchinson, S.A.
1995-07-01
The upwind leapfrog scheme for electromagnetic scattering is briefly described. Its application to the 3D Maxwell`s time domain equations is shown in detail. The scheme`s use of upwind characteristic variables and a narrow stencil result in a smaller demand in communication overhead, making it ideal for implementation on distributed memory parallel computers. The algorithm`s implementation on two message passing computers, a 1024-processor nCUBE 2 and a 1840-processor Intel Paragon, is described. Performance evaluation demonstrates that the scheme performs well with both good scaling qualities and high efficiencies on these machines.
2-D/Axisymmetric Formulation of Multi-dimensional Upwind Scheme
NASA Technical Reports Server (NTRS)
Wood, William A.; Kleb, William L.
2001-01-01
A multi-dimensional upwind discretization of the two-dimensional/axisymmetric Navier-Stokes equations is detailed for unstructured meshes. The algorithm is an extension of the fluctuation splitting scheme of Sidilkover. Boundary conditions are implemented weakly so that all nodes are updated using the base scheme, and eigen-value limiting is incorporated to suppress expansion shocks. Test cases for Mach numbers ranging from 0.1-17 are considered, with results compared against an unstructured upwind finite volume scheme. The fluctuation splitting inviscid distribution requires fewer operations than the finite volume routine, and is seen to produce less artificial dissipation, leading to generally improved solution accuracy.
An upwind, kinetic flux-vector splitting method for flows in chemical and thermal non-equilibrium
NASA Technical Reports Server (NTRS)
Eppard, W. M.; Grossman, B.
1993-01-01
We have developed new upwind kinetic difference schemes for flows with non-equilibrium thermodynamics and chemistry. These schemes are derived from the Boltzmann equation with the resulting Euler schemes developed as moments of the discretized Boltzmann scheme with a locally Maxwellian velocity distribution. Splitting the velocity distribution at the Boltzmann level is seen to result in a flux-split Euler scheme and is called Kinetic Flux Vector Splitting (KFVS). Extensions to flows with finite-rate chemistry and vibrational relaxation is accomplished utilizing nonequilibrium kinetic theory. Computational examples are presented comparing KFVS with the schemes of Van Leer and Roe for a quasi-one-dimensional flow through a supersonic diffuser, inviscid flow through two-dimensional inlet, and viscous flow over a cone at zero angle-of-attack. Calculations are also shown for the transonic flow over a bump in a channel and the transonic flow over an NACA 0012 airfoil. The results show that even though the KFVS scheme is a Riemann solver at the kinetic level, its behavior at the Euler level is more similar to the existing flux-vector splitting algorithms than to the flux-difference splitting scheme of Roe.
NASA Technical Reports Server (NTRS)
1991-01-01
Viewgraphs of briefings presented at the SSTAC/ARTS review of the draft Integrated Technology Plan (ITP) on aerothermodynamics, automation and robotics systems, sensors, and high-temperature superconductivity are included. Topics covered include: aerothermodynamics; aerobraking; aeroassist flight experiment; entry technology for probes and penetrators; automation and robotics; artificial intelligence; NASA telerobotics program; planetary rover program; science sensor technology; direct detector; submillimeter sensors; laser sensors; passive microwave sensing; active microwave sensing; sensor electronics; sensor optics; coolers and cryogenics; and high temperature superconductivity.
Aerothermodynamic heating and performance analysis of a high-lift aeromaneuvering AOTV concept
NASA Technical Reports Server (NTRS)
Menees, G. P.; Brown, K. G.; Wilson, J. F.; Davies, C. B.
1985-01-01
The thermal-control requirements for design-optimized aeromaneuvering performance are determined for space-based applications and low-earth orbit sorties involving large, multiple plane-inclination changes. The leading-edge heating analysis is the most advanced developed for hypersonic-rarefied flow over lifting surfaces at incidence. The effects of leading-edge bluntness, low-density viscous phenomena, and finite-rate flow-field chemistry and surface catalysis are accounted for. The predicted aerothermodynamic heating characteristics are correlated with thermal-control and flight-performance capabilities. The mission payload capability for delivery, retrieval, and combined operations is determined for round-trip sorties extending to polar orbits. Recommendations are given for future design refinements. The results help to identify technology issues required to develop prototype operational systems.
Legacy of the Space Shuttle from an Aerodynamic and Aerothermodynamic Perspective
NASA Technical Reports Server (NTRS)
Martin, Fred W.
2011-01-01
The development of the Space Shuttle Orbiter thermal protection system heating environment is described from a design stand point that began in the early 1970s. The desire for a light weight, reusable heat shield required the development of new technology, relative to previous manned spacecraft, and a systems approach to the design of the vehicle, entry guidance, and thermal protection system. Several unanticipated issues had to be resolved in both the entry and ascent phases of flight, which are discussed at a high level. During the life of the Program, significant improvements in computing power and numerical methods have been applied to Space Shuttle aerodynamic and aerothermodynamic issues, with the Shuttle Program often being the motivation, and or sponsor of the analysis development.
NASA Technical Reports Server (NTRS)
Sun, Y. H.; Gaede, A. E.; Sainio, W. C.
1975-01-01
Combustor test results of the NASA Aerothermodynamic Integration Model are presented of a ramjet engine developed for operation between Mach 3 and 8. Ground-based and flight experiments which provide the data required to advance the technology of hypersonic air-breathing propulsion systems as well as to evaluate facility and testing techniques are described. The engine was tested with synthetic air at Mach 5, 6, and 7. The hydrogen fuel was heated to 1500 R prior to injection to simulate a regeneratively cooled system. Combustor efficiencies up to 95 percent at Mach 6 were achieved. Combustor process in terms of effectiveness, pressure integral factor, total pressure recovery and Crocco's pressure-area relationship are presented and discussed. Interactions between inlet-combustor, combustor stages, combustor-nozzle, and the effects of altitude, combustor step, and struts are observed and analyzed.
NASA Astrophysics Data System (ADS)
Dutheil, Sylvain; Pibarot, Julien; Tran, Dac; Vallee, Jean-Jacques; Tribot, Jean-Pierre
2016-07-01
With the aim of placing Europe among the world's space players in the strategic area of atmospheric re-entry, several studies on experimental vehicle concepts and improvements of critical re-entry technologies have paved the way for the flight of an experimental space craft. The successful flight of the Intermediate eXperimental Vehicle (IXV), under ESA's Future Launchers Preparatory Programme (FLPP), is definitively a significant step forward from the Atmospheric Reentry Demonstrator flight (1998), establishing Europe as a key player in this field. The IXV project objectives were the design, development, manufacture and ground and flight verification of an autonomous European lifting and aerodynamically controlled reentry system, which is highly flexible and maneuverable. The paper presents, the role of aerodynamics aerothermodynamics as part of the key technologies for designing an atmospheric re-entry spacecraft and securing a successful flight.
Low density aerothermodynamics studies performed by means of the tethered satellite system
NASA Technical Reports Server (NTRS)
Carlomagno, Giovanni M.; De Luca, Luigi; Siemers, Paul M.; Wood, George M., Jr.
1986-01-01
Low density gas flow modeling and current ground wind-tunnel technologies are not presently able to produce fully reliable data concerning low density flow regimes. In order to answer some of these issues, the Shuttle Continuous Open Wind Tunnel (SCOWT) program has been proposed, which makes use of the tethered satellite system (TSS). SCOWT's objective is to investigate the energy and momentum transfer between the tethered satellite and its environmental medium within the range of the thermofluid-dynamic conditions experienced by TSS during its atmospheric flights. The feasibility and capability of SCOWT to perform low density aerothermodynamics studies are investigated. Some of the results, obtained by means of a tether simulation program, and the instrumentation and TSS design main requirements to meet SCOWT objectives are described.
NASA Technical Reports Server (NTRS)
Iliff, Kenneth W.; Shafer, Mary F.
1993-01-01
Aerodynamic and aerothermodynamic comparisons between flight and ground test for the Space Shuttle at hypersonic speeds are discussed. All of the comparisons are taken from papers published by researchers active in the Space Shuttle program. The aerodynamic comparisons include stability and control derivatives, center-of-pressure location, and reaction control jet interaction. Comparisons are also discussed for various forms of heating, including catalytic, boundary layer, top centerline, side fuselage, OMS pod, wing leading edge, and shock interaction. The jet interaction and center-of-pressure location flight values exceeded not only the predictions but also the uncertainties of the predictions. Predictions were significantly exceeded for the heating caused by the vortex impingement on the OMS pods and for heating caused by the wing leading-edge shock interaction.
Toro, P. G. P.; Minucci, M. A. S.; Chanes, J. B. Jr; Oliveira, A. C.; Gomes, F. A. A.; Myrabo, L. N.; Nagamatsu, Henry T.
2008-04-28
The new 0.60-m. nozzle exit diameter hypersonic shock tunnel was designed to study advanced air-breathing propulsion system such as supersonic combustion and/or laser technologies. In addition, it may be used for hypersonic flow studies and investigations of the electromagnetic (laser) energy addition for flow control. This new hypersonic shock tunnel was designed and installed at the Laboratory for of Aerothermodynamics and Hypersonics Prof. Henry T. Nagamatsu, IEAv-CTA, Brazil. The design of the tunnel enables relatively long test times, 2-10 milliseconds, suitable for the experiments performed at the laboratory. Free stream Mach numbers ranging from 6 to 25 can be produced and stagnation pressures and temperatures up to 360 atm. and up to 9,000 K, respectively, can be generated. Shadowgraph and schlieren optical techniques will be used for flow visualization.
Payne, J.L.; Hassan, B.
1998-09-01
Massively parallel computers have enabled the analyst to solve complicated flow fields (turbulent, chemically reacting) that were previously intractable. Calculations are presented using a massively parallel CFD code called SACCARA (Sandia Advanced Code for Compressible Aerothermodynamics Research and Analysis) currently under development at Sandia National Laboratories as part of the Department of Energy (DOE) Accelerated Strategic Computing Initiative (ASCI). Computations were made on a generic reentry vehicle in a hypersonic flowfield utilizing three different distributed parallel computers to assess the parallel efficiency of the code with increasing numbers of processors. The parallel efficiencies for the SACCARA code will be presented for cases using 1, 150, 100 and 500 processors. Computations were also made on a subsonic/transonic vehicle using both 236 and 521 processors on a grid containing approximately 14.7 million grid points. Ongoing and future plans to implement a parallel overset grid capability and couple SACCARA with other mechanics codes in a massively parallel environment are discussed.
NASA Astrophysics Data System (ADS)
Schmisseur, John D.; Erbland, Peter
2012-01-01
This article provides an introduction and overview to the efforts of NATO Research and Technology Organization Task Group AVT-136, Assessment of Aerothermodynamic Flight Prediction Tools through Ground and Flight Experimentation. During the period of 2006-2010, AVT-136 coordinated international contributions to assess the state-of-the-art and research challenges for the prediction of critical aerothermodynamic flight phenomena based on the extrapolation of ground test and numerical simulation. To achieve this goal, efforts were organized around six scientific topic areas: (1) Noses and leading edges, (2) Shock Interactions and Control Surfaces, (3) Shock Layers and Radiation, (4) Boundary Layer Transition, (5) Gas-Surface Interactions, and (6) Base and Afterbody Flows. A key component of the AVT-136 strategy was comparison of state-of-the-art numerical simulations with data to be acquired from planned flight research programs. Although it was recognized from the onset of AVT-136 activities that reliance on flight research data yet to be collected posed a significant risk, the group concluded the substantial benefit to be derived from comparison of computational simulations with flight data warranted pursuit of such a program of work. Unfortunately, program delays and failures in the flight programs contributing to the AVT-136 effort prevented timely access to flight research data. Despite this setback, most of the scientific topic areas developed by the Task Group made significant progress in the assessment of current capabilities. Additionally, the activities of AVT-136 generated substantial interest within the international scientific research community and the work of the Task Group was prominently featured in a total of six invited sessions in European and American technical conferences. In addition to this overview, reviews of the state-of-the-art and research challenges identified by the six research thrusts of AVT-136 are also included in this special
Odor-modulated upwind flight of the sphinx moth, Manduca sexta L.
Willis, M A; Arbas, E A
1991-10-01
1. Male and female Manduca sexta flew upwind in response to the odor of female sex-pheromone gland extract or fresh tobacco leaf respectively, and generated very similar zigzagging tracks along the odor plume. 2. After loss of odor during flight, males and females alike: (1) first flew slower and steered their flight more across the wind, then (2) stopped moving upwind, and finally (3) regressed downwind. 3. Males flying upwind in a pheromone plume in wind of different velocities maintained their ground speed near a relatively constant 'preferred' value by increasing their air speed as the velocity of the wind increased, and also maintained the average angle of their resultant flight tracks with respect to the wind at a preferred value by steering a course more precisely due upwind. 4. The inter-turn duration and turn rate, two measures of the temporal aspects of the flight track, were maintained, on average, with remarkable consistency across all wind velocities and in both sexes. The inter-turn durations also decreased significantly as moths approached the odor source, suggesting modulation of the temporal pattern of turning by some feature of the odor plume. This temporal regularity of turning appears to be one of the most stereotyped features of odor-modulated flight in M. sexta. PMID:1779417
NASA Technical Reports Server (NTRS)
Jilly, L. F. (Editor)
1975-01-01
The design and development of the Aerothermodynamic Integration Model (AIM) of the Hypersonic Research Engine (HRE) is described. The feasibility of integrating the various analytical and experimental data available for the design of the hypersonic ramjet engine was verified and the operational characteristic and the overall performance of the selected design was determined. The HRE-AIM was designed for operation at speeds of Mach 3 through Mach 8.
NASA Astrophysics Data System (ADS)
Finchenko, V. S.; Ivankov, A. A.; Shmatov, S. I.; Mordvinkin, A. S.
2015-12-01
The article presents the initial data for the ExoMars landing module aerothermodynamic calculations, used calculation methods, the calculation results of aerodynamic characteristics of the landing module shape and structural parameters of thermal protection selected during the conceptual design phase. Also, the test results of the destruction of the thermal protection material and comparison of the basic characteristics of the landing module with a front shield in the form of a cone and a spherical segment are presented.
Hybrid upwind discretization of nonlinear two-phase flow with gravity
NASA Astrophysics Data System (ADS)
Lee, S. H.; Efendiev, Y.; Tchelepi, H. A.
2015-08-01
Multiphase flow in porous media is described by coupled nonlinear mass conservation laws. For immiscible Darcy flow of multiple fluid phases, whereby capillary effects are negligible, the transport equations in the presence of viscous and buoyancy forces are highly nonlinear and hyperbolic. Numerical simulation of multiphase flow processes in heterogeneous formations requires the development of discretization and solution schemes that are able to handle the complex nonlinear dynamics, especially of the saturation evolution, in a reliable and computationally efficient manner. In reservoir simulation practice, single-point upwinding of the flux across an interface between two control volumes (cells) is performed for each fluid phase, whereby the upstream direction is based on the gradient of the phase-potential (pressure plus gravity head). This upwinding scheme, which we refer to as Phase-Potential Upwinding (PPU), is combined with implicit (backward-Euler) time discretization to obtain a Fully Implicit Method (FIM). Even though FIM suffers from numerical dispersion effects, it is widely used in practice. This is because of its unconditional stability and because it yields conservative, monotone numerical solutions. However, FIM is not unconditionally convergent. The convergence difficulties are particularly pronounced when the different immiscible fluid phases switch between co-current and counter-current states as a function of time, or (Newton) iteration. Whether the multiphase flow across an interface (between two control-volumes) is co-current, or counter-current, depends on the local balance between the viscous and buoyancy forces, and how the balance evolves in time. The sensitivity of PPU to small changes in the (local) pressure distribution exacerbates the problem. The common strategy to deal with these difficulties is to cut the timestep and try again. Here, we propose a Hybrid-Upwinding (HU) scheme for the phase fluxes, then HU is combined with implicit
... Top Doctors in the Nation Departments & Divisions Home Health Insights Stress & Relaxation Breathing and Relaxation Breathing and Relaxation Make ... Management Assess Your Stress Coping Strategies Identifying ... & Programs Health Insights Doctors & Departments Research & Science Education & Training Make ...
Relaxation Assessment with Varied Structured Milieu (RELAX).
ERIC Educational Resources Information Center
Cassel, Russell N.; Cassel, Susie L.
1983-01-01
Describes Relaxation Assessment with Varied Structured Milieu (RELAX), a clinical program designed to assess the degree to which an individual is able to demonstrate self-control for overall general relaxation. The program is designed for use with the Cassel Biosensors biofeedback equipment. (JAC)
Numerical methods for TVD transport and coupled relaxing processes in gases and plasmas
NASA Technical Reports Server (NTRS)
Cambier, Jean-Luc
1990-01-01
The construction of second-order upwind schemes for nonequilibrium plasmas, for both one- and two-fluid formulations is demonstrated. Coupled relaxation processes, including ionization kinetics and radiative processes and their algorithms for nonequilibrium, multiple temperature conditions are described as well. The paper applies the numerical techniques on some simple test cases, points out critical problems and their solutions, and makes qualitative comparisons with known results, whenever possible.
Second- and third-order upwind difference schemes for hyperbolic conservation laws
NASA Technical Reports Server (NTRS)
Yang, J. Y.
1984-01-01
Second- and third-order two time-level five-point explicit upwind-difference schemes are described for the numerical solution of hyperbolic systems of conservation laws and applied to the Euler equations of inviscid gas dynamics. Nonliner smoothing techniques are used to make the schemes total variation diminishing. In the method both hyperbolicity and conservation properties of the hyperbolic conservation laws are combined in a very natural way by introducing a normalized Jacobian matrix of the hyperbolic system. Entropy satisfying shock transition operators which are consistent with the upwind differencing are locally introduced when transonic shock transition is detected. Schemes thus constructed are suitable for shockcapturing calculations. The stability and the global order of accuracy of the proposed schemes are examined. Numerical experiments for the inviscid Burgers equation and the compressible Euler equations in one and two space dimensions involving various situations of aerodynamic interest are included and compared.
A third-order-accurate upwind scheme for Navier-Stokes solutions at high Reynolds numbers
NASA Astrophysics Data System (ADS)
Agarwal, R. K.
1981-01-01
A third-order-accurate upwind scheme is presented for solution of the steady two-dimensional Navier-Stokes equations in stream-function/vorticity form. The scheme is found to be accurate and stable at high Reynolds numbers. A series of test computations is performed on flows with large recirculating regions. In particular, highly accurate solutions are obtained for flow in a driven square cavity up to Reynolds numbers of 10,000. These computations are used to critically evaluate the accuracy of other existing first- and second-order-accurate upwind schemes. In addition, computations are carried out for flow in a channel with symmetric sudden expansion, flow in a channel with a symmetrically placed blunt base, and the flowfield of an impinging jet. Good agreement is obtained with the computations of other investigators as well as with the available experimental data.
NASA Technical Reports Server (NTRS)
Chen, S. C.; Liu, N. S.; Kim, H. D.
1992-01-01
An algorithm utilizing a first order upwind split flux technique and the diagonally dominant treatment is proposed to be the temporal operator for solving the Navier-Stokes equations. Given the limit of a five point stencil, the right hand side flux derivatives are formulated by several commonly used central and upwind schemes. Their performances are studied through a test case of free vortex convection in a uniform stream. From these results, a superior treatment for evaluating the flux term is proposed and compared with the rest. The application of the proposed algorithm to the full Navier-Stokes equations is demonstrated through a calculation of flow over a backward facing step. Results are compared against the calculation done by using the fourth order central differencing scheme with artificial damping.
Multidimensional single-step vector upwind schemes for highly convective transport problems
Felice, G. de; Denaro, F.M.; Meola, C. . Dept. di Energetica Termofluidodinamica e Condizionamenti Ambientali)
1993-06-01
After a synthesis of existing numerical approaches, a family of single-step time-marching upwind schemes is proposed for the convective-diffusive balance of a scalar in multidimensional incompressible laminar flows. A functional description of the new approach is given too. The schemes are adequate to the simulation of purely or highly convective transport phenomena. The proposed schemes were applied to a scalar test problem (rotating hill) and the results evaluated after an analysis of the intrinsic limitations of test itself. A third-order, fully upwind-biased scheme of the family is applied to the lid-driven and thermally driven square cavity problems. In the lid-driven case a steady-state solution is not achieved at Re = 10[sup 4].
Use of high-resolution upwind scheme for vortical flow simulations
NASA Technical Reports Server (NTRS)
Fujii, Kozo; Obayashi, Shigeru
1989-01-01
For vortical flow simulations at high Reynolds number, it is important to keep the artificial dissipation as small as possible since it induces unphysical decay of the vortex strength. One way to accomplish this is to decrease the grid spacing. Another way is to use computational schemes having little dissipation. In the present paper, one of the high-resolution upwind schemes called 'MUSCL with Roe's average'is applied to vortical flow simulations. Two examples are considered. One is the leading-edge separation-vortex flow over a strake-delta wing. The other is a high-angle of attack supersonic flow over a spaceplane-like configuration. The comparison with the central difference solutions indicates that the present upwind scheme is less dissipative and thus has better resolution for the vortical flows.
Hybrid Upwind Splitting (HUS) by a Field-by-Field Decomposition
NASA Technical Reports Server (NTRS)
Coquel, Frederic; Liou, Meng-Sing
1995-01-01
We introduce and develop a new approach for upwind biasing: the hybrid upwind splitting (HUS) method. This original procedure is based on a suitable hybridization of current prominent flux vector splitting (FVS) and flux difference splitting (FDS) methods. The HUS method is designed to naturally combine the respective strengths of the above methods while excluding their main deficiencies. Specifically, the HUS strategy yields a family of upwind methods that exhibit the robustness of FVS schemes in the capture of nonlinear waves and the accuracy of some FDS schemes in the resolution of linear waves. We give a detailed construction of the HUS methods following a general and systematic procedure directly performed at the basic level of the field by field (i.e. waves) decomposition involved in FDS methods. For such a given decomposition, each field is endowed either with FVS or FDS numerical fluxes, depending on the nonlinear nature of the field under consideration. Such a design principle is made possible thanks to the introduction of a convenient formalism that provides us with a unified framework for upwind methods. The HUS methods we propose bring significant improvements over current methods in terms of accuracy and robustness. They yield entropy-satisfying approximate solutions as they are strongly supported in numerical experiments. Field by field hybrid numerical fluxes also achieve fairly simple and explicit expressions and hence require a computational effort between that of the FVS and FDS. Several numerical experiments ranging from stiff 1D shock-tube to high speed viscous flows problems are displayed, intending to illustrate the benefits of the present approach. We assess in particular the relevance of our HUS schemes to viscous flow calculations.
Compressed Semi-Discrete Central-Upwind Schemes for Hamilton-Jacobi Equations
NASA Technical Reports Server (NTRS)
Bryson, Steve; Kurganov, Alexander; Levy, Doron; Petrova, Guergana
2003-01-01
We introduce a new family of Godunov-type semi-discrete central schemes for multidimensional Hamilton-Jacobi equations. These schemes are a less dissipative generalization of the central-upwind schemes that have been recently proposed in series of works. We provide the details of the new family of methods in one, two, and three space dimensions, and then verify their expected low-dissipative property in a variety of examples.
Bochev, Pavel Blagoveston
2011-06-01
We develop a new formulation of the Control Volume Finite Element Method (CVFEM) with a multidimensional Scharfetter-Gummel (SG) upwinding for the drift-diffusion equations. The formulation uses standard nodal elements for the concentrations and expands the flux in terms of the lowest-order Nedelec H(curl; {Omega})-compatible finite element basis. The SG formula is applied to the edges of the elements to express the Nedelec element degree of freedom on this edge in terms of the nodal degrees of freedom associated with the endpoints of the edge. The resulting upwind flux incorporates the upwind effects from all edges and is defined at the interior of the element. This allows for accurate evaluation of integrals on the boundaries of the control volumes for arbitrary quadrilateral elements. The new formulation admits efficient implementation through a standard loop over the elements in the mesh followed by loops over the element nodes (associated with control volume fractions in the element) and element edges (associated with flux degrees of freedom). The quantities required for the SG formula can be precomputed and stored for each edge in the mesh for additional efficiency gains. For clarity the details are presented for two-dimensional quadrilateral grids. Extension to other element shapes and three dimensions is straightforward.
Scanning of wind turbine upwind conditions: numerical algorithm and first applications
NASA Astrophysics Data System (ADS)
Calaf, Marc; Cortina, Gerard; Sharma, Varun; Parlange, Marc B.
2014-11-01
Wind turbines still obtain in-situ meteorological information by means of traditional wind vane and cup anemometers installed at the turbine's nacelle, right behind the blades. This has two important drawbacks: 1-turbine misalignment with the mean wind direction is common and energy losses are experienced; 2-the near-blade monitoring does not provide any time to readjust the profile of the wind turbine to incoming turbulence gusts. A solution is to install wind Lidar devices on the turbine's nacelle. This technique is currently under development as an alternative to traditional in-situ wind anemometry because it can measure the wind vector at substantial distances upwind. However, at what upwind distance should they interrogate the atmosphere? A new flexible wind turbine algorithm for large eddy simulations of wind farms that allows answering this question, will be presented. The new wind turbine algorithm timely corrects the turbines' yaw misalignment with the changing wind. The upwind scanning flexibility of the algorithm also allows to track the wind vector and turbulent kinetic energy as they approach the wind turbine's rotor blades. Results will illustrate the spatiotemporal evolution of the wind vector and the turbulent kinetic energy as the incoming flow approaches the wind turbine under different atmospheric stability conditions. Results will also show that the available atmospheric wind power is larger during daytime periods at the cost of an increased variance.
Central Upwind Scheme for a Compressible Two-Phase Flow Model
Ahmed, Munshoor; Saleem, M. Rehan; Zia, Saqib; Qamar, Shamsul
2015-01-01
In this article, a compressible two-phase reduced five-equation flow model is numerically investigated. The model is non-conservative and the governing equations consist of two equations describing the conservation of mass, one for overall momentum and one for total energy. The fifth equation is the energy equation for one of the two phases and it includes source term on the right-hand side which represents the energy exchange between two fluids in the form of mechanical and thermodynamical work. For the numerical approximation of the model a high resolution central upwind scheme is implemented. This is a non-oscillatory upwind biased finite volume scheme which does not require a Riemann solver at each time step. Few numerical case studies of two-phase flows are presented. For validation and comparison, the same model is also solved by using kinetic flux-vector splitting (KFVS) and staggered central schemes. It was found that central upwind scheme produces comparable results to the KFVS scheme. PMID:26039242
NASA Technical Reports Server (NTRS)
Arrington, J. P.; Woods, W. C.
1972-01-01
The longitudinal, directional, and lateral static stability and control characteristics of a delta lifting body and a delta-wing body were obtained at a Mach number of 20 in helium for operational Reynolds numbers over an angle-of-attack range of -4 deg to 55 deg. The aerodynamic characteristics of the wing body were then evaluated in an entry study to examine the effects of vehicle performance on the aerothermodynamic parameters associated with constant and variable angle-of-attack modes for a 1500-n. mi. cross range. The experimental results indicated that the vehicles were stable, except for neutral directional stability for the wing-body shape, and could be trimmed over the operational angle-of-attack range; however, the wing-body vehicle had adverse yaw due to roll control. This roll-yaw coupling was not examined for the lifting body. The trajectory analysis indicated that a 17-percent decrease in performance required little change in the constant angle-of-attack entry mode and, in turn, resulted in a small decrease in the total heat load. For the pitch-modulated entry, the performance decrease required the pitch maneuver to begin earlier during entry and to last longer in order to meet the 1500-n. mi. cross range without a major heating penalty. The performance reduction also had little effect on the maximum laminar radiation equilibrium temperature over a major portion of the lower surface of the wing-body vehicle regardless of the entry mode.
NASA Technical Reports Server (NTRS)
Gnoffo, Peter A.; Johnston, Christopher O.; Kleb, Bil
2010-01-01
Challenges to computational aerothermodynamic (CA) simulation and validation of hypersonic flow over planetary entry vehicles are discussed. Entry, descent, and landing (EDL) of high mass to Mars is a significant driver of new simulation requirements. These requirements include simulation of large deployable, flexible structures and interactions with reaction control system (RCS) and retro-thruster jets. Simulation of radiation and ablation coupled to the flow solver continues to be a high priority for planetary entry analyses, especially for return to Earth and outer planet missions. Three research areas addressing these challenges are emphasized. The first addresses the need to obtain accurate heating on unstructured tetrahedral grid systems to take advantage of flexibility in grid generation and grid adaptation. A multi-dimensional inviscid flux reconstruction algorithm is defined that is oriented with local flow topology as opposed to grid. The second addresses coupling of radiation and ablation to the hypersonic flow solver - flight- and ground-based data are used to provide limited validation of these multi-physics simulations. The third addresses the challenges of retro-propulsion simulation and the criticality of grid adaptation in this application. The evolution of CA to become a tool for innovation of EDL systems requires a successful resolution of these challenges.
Aerothermodynamics of compressible flow past a flat plate in the slip-flow regime
NASA Astrophysics Data System (ADS)
Cheng, Chi-Yang; Dai, Yi; Li, Genong; Hu, Yitao; Lai, Ming-Chia
2015-11-01
Compressible flow past a flat plate in the slip-flow regime features a very simple geometry and flow field, but it retains the most relevant and interesting physics in high-speed rarefied gas dynamics. In the slip-flow regime, the aerothermodynamic issues, especially the recovery factors and the convection heat transfer correlation, are the focus of this presentation. We first present the detailed similarity equations, especially the transformed Maxwell's slip and jump boundary conditions, and the equations for the Chapman-Rubesin parameter as well as how we incorporate the variable gas properties and the constitutive scaling model for the Knudsen layer in the similarity equations. The similarity solutions are compared with results published by E. R. van Driest [NACA Technical Note 2597, 1952]. We point out that van Driest's solutions were computed by using no-slip and no-jump boundary conditions. The recovery factor and Nusselt number of the plate are shown as functions of the Reynolds number and the Mach number. Finally, the similarity solutions are also compared with simulations of a two-dimensional computational fluid dynamics model solving the full Navier-Stokes-Fourier equations with slip and jump boundary conditions.
NASA Technical Reports Server (NTRS)
Campbell, Charles H.
2004-01-01
A graphic presentation of the aerothermodynamics analysis conducted in support of the STS-107 accident investigation. Investigation efforts were conducted as part of an integrated AATS team (Aero, Aerothermal, Thermal, Stress) directed by OVEWG. Graphics presented are: STS-107 Entry trajectory and timeline (1st off-nominal event to Post-LOS); Indications from OI telemetry data; Aero/aerothermo/thermal analysis process; Selected STS-107 side fuselage/OMS pod off-nominal temperatures; Leading edge structural subsystem; Relevant forensics evidence; External aerothermal environments; STS-107 Pre-entry EOM3 heating profile; Surface heating and temperatures; Orbiter wing leading edge damage survey; Internal aerothermal environments; Orbiter wing CAD model; Aerodynamic flight reconstruction; Chronology of aerodynamic/aerothermoydynamic contributions; Acreage TPS tile damage; Larger OML perturbations; Missing RCC panel(s); Localized damage to RCC panel/missing T-seal; RCC breach with flow ingestion; and Aero-aerothermal closure. NAIT served as the interface between the CAIB and NASA investigation teams; and CAIB requests for study were addressed.
NASA Technical Reports Server (NTRS)
Hollis, Brian R.
2010-01-01
Recent, current, and planned NASA missions that employ blunt-body entry vehicles pose aerothermodynamic problems that challenge the state-of-the art of experimental and computational methods. The issues of boundary-layer transition and turbulent heating on the heat shield have become important in the designs of both the Mars Science Laboratory and Crew Exploration Vehicle. While considerable experience in these general areas exists, that experience is mainly derived from simple geometries; e.g. sharp-cones and flat-plates, or from lifting bodies such as the Space Shuttle Orbiter. For blunt-body vehicles, application of existing data, correlations, and comparisons is questionable because an all, or mostly, subsonic flow field is produced behind the bow shock, as compared to the supersonic (or even hypersonic) flow of other configurations. Because of the need for design and validation data for projects such as MSL and CEV, many new experimental studies have been conducted in the last decade to obtain detailed boundary-layer transition and turbulent heating data on this class of vehicle. In this paper, details of several of the test programs are reviewed. The laminar and turbulent data from these various test are shown to correlate in terms of edge-based Stanton and Reynolds number functions. Correlations are developed from the data for transition onset and turbulent heating augmentation as functions of momentum thickness Reynolds number. These correlation can be employed as engineering-level design and analysis tools.
NASA Astrophysics Data System (ADS)
Chen, Xin; Liu, Li; Zhou, Sida; Yue, Zhenjiang
2016-04-01
Reduced order models(ROMs) based on the snapshots on the CFD high-fidelity simulations have been paid great attention recently due to their capability of capturing the features of the complex geometries and flow configurations. To improve the efficiency and precision of the ROMs, it is indispensable to add extra sampling points to the initial snapshots, since the number of sampling points to achieve an adequately accurate ROM is generally unknown in prior, but a large number of initial sampling points reduces the parsimony of the ROMs. A fuzzy-clustering-based adding-point strategy is proposed and the fuzzy clustering acts an indicator of the region in which the precision of ROMs is relatively low. The proposed method is applied to construct the ROMs for the benchmark mathematical examples and a numerical example of hypersonic aerothermodynamics prediction for a typical control surface. The proposed method can achieve a 34.5% improvement on the efficiency than the estimated mean squared error prediction algorithm and shows same-level prediction accuracy.
NASA Technical Reports Server (NTRS)
Wood, George M.; Wilmoth, Richard G.; Carlomagno, Giovanni M.; De Luca, Luigi
1990-01-01
Aerothermodynamic, aerodynamic, and atmospheric science data acquired between 55 and 150 km has been limited by the lack of vehicles or platforms capable of sustained operation at these altitudes. Tethered satellites, which have been under study for this purpose by NASA, the Italian Space Agency (ASI), and others for more than a decade, are expected to become a reality by mid-1991. This approach, in which an instrumented platform is maintained at a desired altitude by a tether attached to a host vehicle orbiting at higher altitudes, will provide the first opportunity to obtain steady state data over an extended period encompassing one or more orbital revolutions. This paper describes the objectives and measurement methods for the first of the facility-class satellites, the TSS-2, which is proposed for a 1995 deployment, and gives the status of the experiment definition. Monte Carlo modeling of the flow fields at 130 km around the baseline 1.6 m diameter sphere is discussed and illustrative results of the modeling given.
Aerothermodynamic Design of the Mars Science Laboratory Backshell and Parachute Cone
NASA Technical Reports Server (NTRS)
Edquist, Karl T.; Dyakonov, Artem A.; Wright, Michael J.; Tang, Chun Y.
2009-01-01
Aerothermodynamic design environments are presented for the Mars Science Laboratory entry capsule backshell and parachute cone. The design conditions are based on Navier-Stokes flowfield simulations on shallow (maximum total heat load) and steep (maximum heat flux) design entry trajectories from a 2009 launch. Transient interference effects from reaction control system thruster plumes were included in the design environments when necessary. The limiting backshell design heating conditions of 6.3 W/sq cm for heat flux and 377 J/sq cm for total heat load are not influenced by thruster firings. Similarly, the thrusters do not affect the parachute cover lid design environments (13 W/sq cm and 499 J/sq cm). If thruster jet firings occur near peak dynamic pressure, they will augment the design environments at the interface between the backshell and parachute cone (7 W/sq cm and 174 J/sq cm). Localized heat fluxes are higher near the thruster fairing during jet firings, but these areas did not require additional thermal protection material. Finally, heating bump factors were developed for antenna radomes on the parachute cone
NASA Technical Reports Server (NTRS)
Green, M. J.; Moss, J. N.; Wilson, J. F.
1984-01-01
This paper presents thermal protection system (TPS) requirements for a potential Titan aerocapture vehicle. Shock-layer solutions are obtained for a nominal trajectory through the current Titan model atmosphere. Fully laminar and fully turbulent solutions are presented along the blunted fore-cone in the windward symmetry plane of a bent-biconic vehicle. Using these solutions to define the aerothermodynamic environment, transient material-response solutions are obtained for a Galileo-type TPS with a carbon-phenolic ablator heat shield. Shock-layer results indicate that turbulent flow is the more realistic flow condition. They also show that the lengthy aerocapture heating pulse is dominated by convective heating. The TPS results show that the required insulation thickness is uniformly about 4 cm along the fore-cone because of the long heat-soak period. The total heat-shield thickness is 6.4 cm at the stagnation point, and 4.7 cm near the end of the fore-cone. These TPS requirements are greater than those presented in a previous Titan aerocapture study.
NASA Technical Reports Server (NTRS)
Drummond, J. Philip; Cockrell, Charles E., Jr.; Pellett, Gerald L.; Diskin, Glenn S.; Auslender, Aaron H.; Exton, Reginald J.; Guy, R. Wayne; Hoppe, John C.; Puster, Richard L.; Rogers, R. Clayton
2002-01-01
This White Paper examines the current state of Hypersonic Airbreathing Propulsion at the NASA Langley Research Center and the factors influencing this area of work and its personnel. Using this knowledge, the paper explores beyond the present day and suggests future directions and strategies for the field. Broad views are first taken regarding potential missions and applications of hypersonic propulsion. Then, candidate propulsion systems that may be applicable to these missions are suggested and discussed. Design tools and experimental techniques for developing these propulsion systems are then described, and approaches for applying them in the design process are considered. In each case, current strategies are reviewed and future approaches that may improve the techniques are considered. Finally, the paper concentrates on the needs to be addressed in each of these areas to take advantage of the opportunities that lay ahead for both the NASA Langley Research Center and the Aerodynamic Aerothermodynamic, and Aeroacoustics Competency. Recommendations are then provided so that the goals set forth in the paper may be achieved.
NASA Technical Reports Server (NTRS)
1976-01-01
Twelve aerothermodynamic space technology needs were identified to reduce the design uncertainties in aerodynamic heating and forces experienced by heavy lift launch vehicles, orbit transfer vehicles, and advanced single stage to orbit vehicles for the space transportation system, and for probes, planetary surface landers, and sample return vehicles for solar system exploration vehicles. Research and technology needs identified include: (1) increasing the fluid dynamics capability by at least two orders of magnitude by developing an advanced computer processor for the solution of fluid dynamic problems with improved software; (2) predicting multi-engine base flow fields for launch vehicles; and (3) developing methods to conserve energy in aerothermodynamic ground test facilities.
Computational Aerothermodynamic Assessment of Space Shuttle Orbiter Tile Damage: Open Cavities
NASA Technical Reports Server (NTRS)
Pulsonetti, Maria; Wood, William
2005-01-01
Computational aerothermodynamic simulations of Orbiter windside tile damage in flight were performed in support of the Space Shuttle Return-to-Flight effort. The simulations were performed for both hypervelocity flight and low-enthalpy wind tunnel conditions and contributed to the Return-to-Flight program by providing information to support a variety of damage scenario analyses. Computations at flight conditions were performed at or very near the peak heating trajectory point for multiple damage scenarios involving damage windside acreage reaction cured glass (RCG) coated silica tile(s). The cavities formed by the missing tile examined in this study were relatively short leading to flow features which indicated open cavity behavior. Results of the computations indicated elevated heating bump factor levels predicted for flight over the predictions for wind tunnel conditions. The peak heating bump factors, defined as the local heating to a reference value upstream of the cavity, on the cavity floor for flight simulation were 67% larger than the peak wind tunnel simulation value. On the downstream face of the cavity the flight simulation values were 60% larger than the wind tunnel simulation values. On the outer mold line (OML) downstream of the cavity, the flight values are about 20% larger than the wind tunnel simulation values. The higher heating bump factors observed in the flight simulations were due to the larger driving potential in terms of energy entering the cavity for the flight simulations. This is evidenced by the larger rate of increase in the total enthalpy through the boundary layer prior to the cavity for the flight simulation.
On the implementation of a class of upwind schemes for system of hyperbolic conservation laws
NASA Technical Reports Server (NTRS)
Yee, H. C.
1985-01-01
The relative computational effort among the spatially five point numerical flux functions of Harten, van Leer, and Osher and Chakravarthy is explored. These three methods typify the design principles most often used in constructing higher than first order upwind total variation diminishing (TVD) schemes. For the scalar case the difference in operation count between any two algorithms may be very small and yet the operation count for their system counterparts might be vastly different. The situation occurs even though one starts with two different yet equivalent representations for the scalar case.
On the relation between upwind-differencing schemes of Godunov, Engquist-Osher and Roe
NASA Technical Reports Server (NTRS)
Vanleer, B.
1981-01-01
The upwind differencing first order schemes of Godunov, Engquist-Osher and Roe are discussed on the basis of the inviscid Burgers equations. The differences between the schemes are interpreted as differences between the approximate Riemann solutions on which their numerical flux functions are based. Special attention is given to the proper formulation of these schemes when a source term is present. Second order two step schemes, based on the numerical flux functions of the first order schemes are also described. The schemes are compared in a numerical experiment, and recommendations on their use are included.
Upwind differencing and LU factorization for chemical non-equilibrium Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Shuen, Jian-Shun
1992-01-01
By means of either the Roe or the Van Leer flux-splittings for inviscid terms, in conjunction with central differencing for viscous terms in the explicit operator and the Steger-Warming splitting and lower-upper approximate factorization for the implicit operator, the present, robust upwind method for solving the chemical nonequilibrium Navier-Stokes equations yields formulas for finite-volume discretization in general coordinates. Numerical tests in the illustrative cases of a hypersonic blunt body, a ramped duct, divergent nozzle flows, and shock wave/boundary layer interactions, establish the method's efficiency.
Turbulence Intensity at Inlet of 80- by 120-Foot Wind Tunnel Caused by Upwind Blockage
NASA Technical Reports Server (NTRS)
Salazar, Denise; Yuricich, Jillian
2014-01-01
In order to estimate the magnitude of turbulence in the National Full-Scale Aerodynamics Complex (NFAC) 80- by 120-Foot Wind Tunnel (80 x 120) caused by buildings located upwind from the 80 x 120 inlet, a 150th-scale study was performed that utilized a nominal two-dimensional blockage placed ahead of the inlet. The distance of the blockage ahead of the inlet was varied. This report describes velocity measurements made in the plane of the 80 x 120 model inlet for the case of zero ambient (atmospheric) wind.
Parallel adaptive Cartesian upwind methods for shock-driven multiphysics simulation
Deiterding, Ralf
2011-01-01
The multiphysics fluid-structure interaction simulation of shock-loaded thin-walled structures requires the dynamic coupling of a shock-capturing flow solver to a solid mechanics solver for large deformations. By combining a Cartesian embedded boundary approach with dynamic mesh adaptation a generic software framework for such flow solvers has been constructed that allows easy exchange of the specific hydrodynamic finite volume upwind scheme and coupling to various explicit finite element solid dynamics solvers. The paper gives an overview of the computational approach and presents first simulations that couple the software to the general purpose solid dynamics code DYNA3D.
Computation of turbulent pipe and duct flow using third order upwind scheme
NASA Technical Reports Server (NTRS)
Kawamura, T.
1986-01-01
The fully developed turbulence in a circular pipe and in a square duct is simulated directly without using turbulence models in the Navier-Stokes equations. The utilized method employs a third-order upwind scheme for the approximation to the nonlinear term and the second-order Adams-Bashforth method for the time derivative in the Navier-Stokes equation. The computational results appear to capture the large-scale turbulent structures at least qualitatively. The significance of the artificial viscosity inherent in the present scheme is discussed.
NASA Astrophysics Data System (ADS)
Ozawa, Hiroshi
2016-04-01
This paper describes an experimental study that used a fast-response temperature-sensitive paint (TSP) to investigate the unsteady aerothermodynamic phenomena occurring on a shock-tube wall. To understand these phenomena in detail, a fast-response TSP with high temperature sensitivity developed for transient temperature measurement was applied to the wall. The shock-tube experiment was carried out under the over-tailored condition, with a pressure ratio of 110 for test gases of air in driver/driven tubes. The following aspects were clarified using the TSP: (a) the TSP could be used to visualize the unsteady aerothermodynamic phenomena and estimate the quantitative heat flux on the shock-tube wall; (b) an x-t diagram based on the TSP response showed shock-tube wall characteristics that included the incident/reflected shocks, laminar-to-turbulent boundary-layer transition, streaks in the turbulent boundary layer, reflected shock/turbulent boundary layer interaction, and waves reflected from a contact surface; (c) the TSP graphically showed that a transition front from the plate's leading edge and turbulent spots moved with 80% of the free-stream velocity behind the incident shock. In addition, the TSP could track the growth of the turbulent spots on the wall.
Candy, J.; Waltz, R.E.
2006-03-15
Equations which describe the evolution of volume-averaged gyrokinetic entropy are derived and added to GYRO [J. Candy and R.E. Waltz, J. Comput. Phys. 186, 545 (2003)], a Eulerian gyrokinetic turbulence simulation code. In particular, the creation of entropy through spatial upwind dissipation (there is zero velocity-space dissipation in GYRO) and the reduction of entropy via the production of fluctuations are monitored in detail. This new diagnostic has yielded several key confirmations of the validity of the GYRO simulations. First, fluctuations balance dissipation in the ensemble-averaged sense, thus demonstrating that turbulent GYRO simulations achieve a true statistical steady state. Second, at the standard spatial grid size, neither entropy nor energy flux is significantly changed by a 16-fold increase (from 32 to 512 grid points per cell) in the number of grid points in the two-dimensional velocity space. Third, the measured flux is invariant to an eightfold increase in the upwind dissipation coefficients. A notable conclusion is that the lack of change in entropy with grid refinement refutes the familiar but incorrect notion that Eulerian gyrokinetic codes miss important velocity-space structure. The issues of density and energy conservation and their relation to negligible second-order effects such as the parallel nonlinearity are also discussed.
Baker, T C; Hansson, B S; Löfstedt, C; Löfqvist, J
1988-01-01
A wind-borne plume of sex pheromone from a female moth or a synthetic source has a fine, filamentous structure that creates steep and rapid fluctuations in concentration for a male moth flying up the plume's axis. The firing rates from single antennal neurons on Agrotis segetum antennae decreased to nearly zero within seconds after the antennae were placed in a pheromone plume 70 cm downwind of a high-concentration source known from previous studies to cause in-flight arrestment of upwind progress. In a separate experiment, the fluctuating output from chilled neurons on Grapholita molesta antennae became attenuated in response to repetitive, experimentally delivered pheromone pulses. The attenuation was correlated with a previously reported higher percentage of in-flight arrestment exhibited by moths flying at cooler compared to warmer temperatures. These results indicate that two peripheral processes related to excessive concentration, complete adaptation of antennal neurons, or merely the attenuation of fluctuations in burst frequency, are important determinants of when upwind progress by a moth flying in a pheromone plume stops and changes to station keeping. Also, adaptation and attenuation may affect the sensation of blend quality by preferentially affecting cells sensitive to the most abundant components in airborne pheromone blends. PMID:3200859
Cardiorespiratory and muscular responses to simulated upwind sailing exercise in Optimist sailors.
Callewaert, Margot; Boone, Jan; Celie, Bert; De Clercq, Dirk; Bourgois, Jan G
2014-02-01
The aim of this work was to gain more insight into the cardiorespiratory and muscular (m. vastus lateralis) responses to simulated upwind sailing exercise in 10 high-level male and female Optimist sailors (10.8-14.4 years old). Hiking strap load (HSL) and cardiorespiratory variables were measured while exercising on a specially developed Optimist sailing ergometer. Electromyography (EMG) was used to determine mean power frequency (MPF) and root mean square (RMS). Near-infrared spectroscopy was used to measure deoxygenated Hemoglobin and Myoglobin concentration (deoxy[Hb+Mb]) and re-oxygenation. Results indicated that HSL and integrated EMG of the vastus lateralis muscle changed in accordance with the hiking intensity. Cardiorespiratory response demonstrated an initial significant increase and subsequently steady state in oxygen uptake (VO₂), ventilation (VE), and heart rate (HR) up to circa 40% VO₂peak, 30% VEpeak and 70% HRpeak respectively. At muscle level, results showed that highly trained Optimist sailors manage to stabilize the muscular demand and fatigue development during upwind sailing (after an initial increase). However, approaching the end of the hiking exercise, the MPF decrease, RMS increase, and deoxy[Hb+Mb] increase possibly indicate the onset of muscle fatigue. PMID:24018866
Baker, T C; Hansson, B S; Löfstedt, C; Löfqvist, J
1988-12-01
A wind-borne plume of sex pheromone from a female moth or a synthetic source has a fine, filamentous structure that creates steep and rapid fluctuations in concentration for a male moth flying up the plume's axis. The firing rates from single antennal neurons on Agrotis segetum antennae decreased to nearly zero within seconds after the antennae were placed in a pheromone plume 70 cm downwind of a high-concentration source known from previous studies to cause in-flight arrestment of upwind progress. In a separate experiment, the fluctuating output from chilled neurons on Grapholita molesta antennae became attenuated in response to repetitive, experimentally delivered pheromone pulses. The attenuation was correlated with a previously reported higher percentage of in-flight arrestment exhibited by moths flying at cooler compared to warmer temperatures. These results indicate that two peripheral processes related to excessive concentration, complete adaptation of antennal neurons, or merely the attenuation of fluctuations in burst frequency, are important determinants of when upwind progress by a moth flying in a pheromone plume stops and changes to station keeping. Also, adaptation and attenuation may affect the sensation of blend quality by preferentially affecting cells sensitive to the most abundant components in airborne pheromone blends. PMID:3200859
NASA Technical Reports Server (NTRS)
Narain, J. P.; Muramoto, K. K.; Lawrence, S. L.
1991-01-01
A three-dimensional parabolized Navier-Stokes computer code which employs an upwind algorithm is used to conduct a numerical study of an advanced maneuvering reentry vehicle configuration. Comparisons between numerical solutions and experimental data are presented for surface pressure, wall heat flux, and overall forces and moments. The effects of angle of attack, angle of yaw, and surface mass injection are investigated. Good agreement is observed between the calculated and measured data. The results of this investigation demonstrate the accuracy and efficiency of an upwind scheme in predicting the hypersonic flow field characteristics about a complex configuration.
NASA Astrophysics Data System (ADS)
Dobes, Jiri; Deconinck, Herman
2008-06-01
Multidimensional upwind residual distribution (RD) schemes have become an appealing alternative to more widespread finite volume and finite element methods (FEM) for solving compressible fluid flows. The RD approach allows to construct nonlinear second order and non-oscillatory methods at the same time. They are routinely used for steady state calculations of the complex flow problems, e.g., 3D turbulent transonic industrial-type simulations [H. Deconinck, K. Sermeus, R. Abgrall, Status of multidimensional upwind residual distribution schemes and applications in aeronautics, AAIA Paper 2000-2328, AIAA, 2000; K. Sermeus, H. Deconinck, Drag prediction validation of a multi-dimensional upwind solver, CFD-based aircraft drag prediction and reduction, VKI Lecture Series 2003-02, Von Karman Institute for Fluid Dynamics, Chausee do Waterloo 72, B-1640 Rhode Saint Genese, Belgium, 2003]. Despite its maturity, some problems are still present for the nonlinear schemes developed up to now: namely a poor iterative convergence for the transonic problems and a decrease of accuracy in smooth parts of the flow, caused by a weak L2 instability [M. Ricchiuto, Construction and analysis of compact residual discretizations for conservation laws on unstructured meshes. Ph.DE Thesis, Universite Libre de Bruxelles, Von Karman Institute for Fluid Dynamics, 2005]. We have developed a new formulation of a blended scheme between the second order linear LDA [R. Abgrall, M. Mezine, Residual distribution scheme for steady problems, 33rd Computational Fluid Dynamics course, VKI Lecture Series 2003-05, Von Karman Institute for Fluid Dynamics, Chausee do Waterloo 72, B-1640 Rhode Saint Genese, Belgium, 2003] scheme and the first order N scheme. The blending coefficient is based on a simple shock capturing operator and it is properly scaled such that second order accuracy is preserved. The approach is extended to unsteady flows problems using consistent formulation of the LDA scheme with the mass
Set up an Arc Welding Code with Enthalpy Method in Upwind Scheme
Ho, J.-E.
2010-05-21
In this study, a numerical code with enthalpy method in upwind scheme is proposed to estimate the distribution of thermal stress in the molten pool, which is primarily determined by the type of the input power and travel speed of heating source. To predict the cracker deficit inside the workpiece, a simulated program satisfying the diagonal domination and Scarborough criterion provides a stable iteration. Meantime, an experimental performance, operated by robot arm 'DR-400' to provide a steady and continuous arc welding, was also conducted to verify the simulated result. By surveying the consistence of molten pool bounded by contrast shade and simulated melting contour on the surface of workpiece, the validity of model proposed to predict the thermal cracker has been successfully identified.
Computation of viscous blast wave solutions with an upwind finite volume method
NASA Technical Reports Server (NTRS)
Molvik, Gregory A.
1987-01-01
A fully conservative, viscous, implicit, upwind, finite-volume scheme for the thin-layer Navier-Stokes equations is described with application to blast wave flow fields. In this scheme, shocks are captured without the oscillations typical of central differencing techniques and wave speeds are accurately predicted. The finite volume philosophy ensures conservation and since boundary conditions are also treated conservatively, accurate reflections of waves from surfaces are assured. Viscous terms in the governing equations are treated in a manner consistent with the finite volume philosophy, resulting in very accurate prediction of boundary layer quantities. Numerical results are presented for four viscous problems: a steady boundary layer, a shock-induced boundary layer, a blast wave/cylinder interaction and a blast wave/supersonic missile interaction. Comparisons of the results with an established boundary layer code, similarity solution, and experimental data show excellent agreement.
Upwind MacCormack Euler solver with non-equilibrium chemistry
NASA Technical Reports Server (NTRS)
Sherer, Scott E.; Scott, James N.
1993-01-01
A computer code, designated UMPIRE, is currently under development to solve the Euler equations in two dimensions with non-equilibrium chemistry. UMPIRE employs an explicit MacCormack algorithm with dissipation introduced via Roe's flux-difference split upwind method. The code also has the capability to employ a point-implicit methodology for flows where stiffness is introduced through the chemical source term. A technique consisting of diagonal sweeps across the computational domain from each corner is presented, which is used to reduce storage and execution requirements. Results depicting one dimensional shock tube flow for both calorically perfect gas and thermally perfect, dissociating nitrogen are presented to verify current capabilities of the program. Also, computational results from a chemical reactor vessel with no fluid dynamic effects are presented to check the chemistry capability and to verify the point implicit strategy.
Development of an upwind, finite-volume code with finite-rate chemistry
NASA Technical Reports Server (NTRS)
Molvik, Gregory A.
1994-01-01
Under this grant, two numerical algorithms were developed to predict the flow of viscous, hypersonic, chemically reacting gases over three-dimensional bodies. Both algorithms take advantage of the benefits of upwind differencing, total variation diminishing techniques, and a finite-volume framework, but obtain their solution in two separate manners. The first algorithm is a zonal, time-marching scheme, and is generally used to obtain solutions in the subsonic portions of the flow field. The second algorithm is a much less expensive, space-marching scheme and can be used for the computation of the larger, supersonic portion of the flow field. Both codes compute their interface fluxes with a temporal Riemann solver and the resulting schemes are made fully implicit including the chemical source terms and boundary conditions. Strong coupling is used between the fluid dynamic, chemical, and turbulence equations. These codes have been validated on numerous hypersonic test cases and have provided excellent comparison with existing data.
Keyes, D.E. . Dept. of Mechanical Engineering); Gropp, W.D. )
1990-01-01
Discrete systems arising in computational fluid dynamics applications often require wide stencils adapted to the local convective direction in order to accommodate higher-order upwind differencing, and involve multiple components perhaps coupling strongly at each point. Conventional exactly or approximately factored inverses of such operators are burdensome to apply globally, especially in problems complicated by non-tensor-product domain geometry or adaptive refinement, though their forward'' action is not. Such problems can be solved by iterative methods by using either point-block preconditioners or combination space-decoupled/component-decoupled preconditioners that are based on lower-order discretizations. Except for a global implicit solve on a coarse grid, each phase in the application of such preconditioners has simple locally exploitable structure. 16 refs., 2 figs., 3 tabs.
High resolution applications of the Osher upwind scheme for the Euler equations
NASA Technical Reports Server (NTRS)
Chakravarthy, S. R.; Osher, S.
1983-01-01
The 'Osher' scheme was introduced by Osher (1981). It represents an upwind finite-difference method for hyperbolic systems of conservation laws, including the Euler equations. In studies conducted by Osher (1981) and Osher and Solomon (1982), the method was applied to the nonisentropic form of the Euler equations in one dimension and the isentropically restricted form in two spatial dimensions, both in Cartesian coordinates. Chakravarthy and Osher (1982) have shown an approach for extending the Osher scheme to the Euler equations written for general geometries, taking into account the use of mappings to arbitrary curvilinear coordinate systems. The present investigation is concerned with the high resolution extension of the Osher scheme to second-order accuracy. Results are presented for several example problems, giving attention to quasi-one-dimensional Laval nozzle flow, a one-dimensional shock tube problem, and supersonic flow over a cylinder.
Multi-dimensional upwind fluctuation splitting scheme with mesh adaption for hypersonic viscous flow
NASA Astrophysics Data System (ADS)
Wood, William Alfred, III
A multi-dimensional upwind fluctuation splitting scheme is developed and implemented for two dimensional and axisymmetric formulations of the Navier-Stokes equations on unstructured meshes. Key features of the scheme are the compact stencil, full upwinding, and non-linear discretization which allow for second-order accuracy with enforced positivity. Throughout, the fluctuation splitting scheme is compared to a current state-of-the-art finite volume approach, a second-order, dual mesh upwind flux difference splitting scheme (DMFDSFV), and is shown to produce more accurate results using fewer computer resources for a wide range of test cases. The scalar test cases include advected shear, circular advection, non-linear advection with coalescing shock and expansion fans, and advection-diffusion. For all scalar cases the fluctuation splitting scheme is more accurate, and the primary mechanism for the improved fluctuation splitting performance is shown to be the reduced production of artificial dissipation relative to DMFDSFV. The most significant scalar result is for combined advection-diffusion, where the present fluctuation splitting scheme is able to resolve the physical dissipation from the artificial dissipation on a much coarser mesh than DMFDSFV is able to, allowing order-of-magnitude reductions in solution time. Among the inviscid test cases the converging supersonic streams problem is notable in that the fluctuation splitting scheme exhibits superconvergent third-order spatial accuracy. For the inviscid cases of a supersonic diamond airfoil, supersonic slender cone, and incompressible circular bump the fluctuation splitting drag coefficient errors are typically half the DMFDSFV drag errors. However, for the incompressible inviscid sphere the fluctuation splitting drag error is larger than for DMFDSFV. A Blasius flat plate viscous validation case reveals a more accurate v-velocity profile for fluctuation splitting, and the reduced artificial dissipation
Orographic effects of the subtropical and extratropical Andes on upwind precipitating clouds
NASA Astrophysics Data System (ADS)
Viale, Maximiliano; Garreaud, René
2015-05-01
The orographic effect of the Andes (30°S-55°S) on upwind precipitating clouds from midlatitude frontal systems is investigated using surface and satellite data. Rain gauges between 33°S and 44°S indicate that annual precipitation increases from the Pacific coast to the windward slopes by a factor of 1.8 ± 0.3. Hourly gauges and instantaneous satellite estimates reveal that the cross-barrier increase in annual precipitation responds to an increase in both the intensity and frequency of precipitation. CloudSat satellite data indicate that orographic effects of the Andes on precipitating ice clouds increase gradually from midlatitudes to subtropics, likely as a result of a reduction of synoptic forcing and an increase of the height of the Andes equatorward. To the south of 40°S, the thickness of clouds slightly decreases from offshore to the Andes. The total ice content increases substantially from the open ocean to the coastal zone (except to the south of 50°S, where there is no much variation over the ocean), and then experience little changes in the cross-mountain direction over the upstream and upslope sectors. Nevertheless, the maximum ice content over the upslope sector is larger and occurs at a lower level than their upwind counterparts. In the subtropics, the offshore clouds contain almost no ice, but the total and maximum ice content significantly increases toward the Andes, with values being much larger than their counterparts over the extratropical Andes. Further, the largest amounts of cloud ice are observed upstream of the tallest Andes, suggesting that upstream blocking dominates there.
NASA Astrophysics Data System (ADS)
Müller, K.; van Pinxteren, D.; Plewka, A.; Svrcina, B.; Kramberger, H.; Hofmann, D.; Bächmann, K.; Herrmann, H.
An extensive set of gaseous and particulate organic compounds was quantified before an orographic cloud passage at the upwind site of the research region in Thüringer Wald. Samples were collected with two different time resolutions, 2 h for gaseous species and spray absorber samples and the whole cloud event duration to determine the concentrations of ketones, aldehydes, monocarboxylic acids, dicarboxylic acids (DCA), hydrocarbons, biogenic sugars and alcohols in both the gas and particle phase. The measurement of different groups of organic compounds delivered size-segregated concentrations at the upwind site of a cloud experiment. The size distribution of DCA showed a peak in the mass-rich impactor stage 3 (0.42-1.2 μm). The concentrations of DCA from the filters, the impactor foils as well as the spray absorber samples decreased with increasing C-number. The time resolved measurements revealed an increasing mixing ratio from night time to midday for carboxylic and DCA, and related carbonyl compounds. The biogenic compounds xylitol (up to 103 ng m -3), levoglucosan (up to 62 ng m -3) and pinonaldehyde (up to 34 ng m -3) were the compounds found in highest concentrations in the particle phase beside the oxalate (up to 104 ng m -3). The organic trace gases with the highest mixing ratios identified were formaldehyde (up to 1.47 ppbv), acetaldehyde (up to 0.84 ppbv) and acetone (up to 0.65 ppbv), acetic acid (up to 0.43 ppbv) and formic acid (up to 0.41 ppbv).
NASA Technical Reports Server (NTRS)
Gaede, A. E.; Platte, W. (Editor)
1975-01-01
The data reduction program used to analyze the performance of the Aerothermodynamic Integration Model is described. Routines to acquire, calibrate, and interpolate the test data, to calculate the axial components of the pressure area integrals and the skin function coefficients, and to report the raw data in engineering units are included along with routines to calculate flow conditions in the wind tunnel, inlet, combustor, and nozzle, and the overall engine performance. Various subroutines were modified and used to obtain species concentrations and transport properties in chemical equilibrium at each of the internal and external engine stations. It is recommended that future test plans include the configuration, calibration, and channel assignment data on a magnetic tape generated at the test site immediately before or after a test, and that the data reduction program be designed to operate in a batch environment.
NASA Technical Reports Server (NTRS)
Loh, Ching Y.; Jorgenson, Philip C. E.
2007-01-01
A time-accurate, upwind, finite volume method for computing compressible flows on unstructured grids is presented. The method is second order accurate in space and time and yields high resolution in the presence of discontinuities. For efficiency, the Roe approximate Riemann solver with an entropy correction is employed. In the basic Euler/Navier-Stokes scheme, many concepts of high order upwind schemes are adopted: the surface flux integrals are carefully treated, a Cauchy-Kowalewski time-stepping scheme is used in the time-marching stage, and a multidimensional limiter is applied in the reconstruction stage. However even with these up-to-date improvements, the basic upwind scheme is still plagued by the so-called "pathological behaviors," e.g., the carbuncle phenomenon, the expansion shock, etc. A solution to these limitations is presented which uses a very simple dissipation model while still preserving second order accuracy. This scheme is referred to as the enhanced time-accurate upwind (ETAU) scheme in this paper. The unstructured grid capability renders flexibility for use in complex geometry; and the present ETAU Euler/Navier-Stokes scheme is capable of handling a broad spectrum of flow regimes from high supersonic to subsonic at very low Mach number, appropriate for both CFD (computational fluid dynamics) and CAA (computational aeroacoustics). Numerous examples are included to demonstrate the robustness of the methods.
A multiple-relaxation-time lattice Boltzmann method for high-speed compressible flows
NASA Astrophysics Data System (ADS)
Li, Kai; Zhong, Cheng-Wen
2015-05-01
This paper presents a coupling compressible model of the lattice Boltzmann method. In this model, the multiple-relaxation-time lattice Boltzmann scheme is used for the evolution of density distribution functions, whereas the modified single-relaxation-time (SRT) lattice Boltzmann scheme is applied for the evolution of potential energy distribution functions. The governing equations are discretized with the third-order Monotone Upwind Schemes for scalar conservation laws finite volume scheme. The choice of relaxation coefficients is discussed simply. Through the numerical simulations, it is found that compressible flows with strong shocks can be well simulated by present model. The numerical results agree well with the reference results and are better than that of the SRT version. Project supported by the Innovation Fund for Aerospace Science and Technology of China (Grant No. 2009200066) and the Aeronautical Science Fund of China (Grant No. 20111453012).
TEACHING NEUROMUSCULAR RELAXATION.
ERIC Educational Resources Information Center
NORRIS, JEANNE E.; STEINHAUS, ARTHUR H.
THIS STUDY ATTEMPTED TO FIND OUT WHETHER (1) THE METHODS FOR ATTAINING NEUROMUSCULAR RELAXATION THAT HAVE PROVED FRUITFUL IN THE ONE-TO-ONE RELATIONSHIP OF THE CLINIC CAN BE SUCCESSFULLY ADAPTED TO THE TEACHER-CLASS RELATIONSHIP OF THE CLASSROOM AND GYMNASIUM, AND (2) NEUROMUSCULAR RELAXATION CAN BE TAUGHT SUCCESSFULLY BY AN APPROPRIATELY TRAINED…
Mechanism of rotational relaxation.
NASA Technical Reports Server (NTRS)
Polanyi, J. C.; Woodall, K. B.
1972-01-01
A model is presented which describes the characteristic pattern of relaxation of a nonthermal rotational distribution of hydrogen halide, peaked initially at high rotational quantum number J, to a thermal distribution without generating a peak at intermediate J. A method for correcting infrared chemiluminiscence data for modest rotational relaxation is also suggested.
NASA Astrophysics Data System (ADS)
Jin, G.
2012-12-01
. The pressure field is then implicitly calculated from the pressure equation, which in turn results in the derived velocity field for directional flux calculation at each grid node. Directional flux at the center of each interaction surface is also calculated by interpolation from the element nodal fluxes using shape functions. The MPFA scheme is performed by a specific linear combination of all incoming fluxes into the upstream cell represented by either nodal fluxes or interpolated surface boundary fluxes to produce an upwind directional fluxed weighted relative mobility at the center of the interaction region boundary. Such an upwind weighted relative mobility is then used for calculating the saturations of each fluid phase explicitly. The proposed upwind weighting scheme has been implemented into a mixed finite element-finite volume (FE-FV) method, which allows for handling complex reservoir geometry with second-order accuracies in approximating primary variables. The numerical solver has been tested with several bench mark test problems. The application of the proposed scheme to migration path analysis of CO2 injected into deep saline reservoirs in 3-D has demonstrated its ability and robustness in handling multiphase flow with adverse mobility contrast in highly heterogeneous porous media.
Application of Central Upwind Scheme for Solving Special Relativistic Hydrodynamic Equations.
Yousaf, Muhammad; Ghaffar, Tayabia; Qamar, Shamsul
2015-01-01
The accurate modeling of various features in high energy astrophysical scenarios requires the solution of the Einstein equations together with those of special relativistic hydrodynamics (SRHD). Such models are more complicated than the non-relativistic ones due to the nonlinear relations between the conserved and state variables. A high-resolution shock-capturing central upwind scheme is implemented to solve the given set of equations. The proposed technique uses the precise information of local propagation speeds to avoid the excessive numerical diffusion. The second order accuracy of the scheme is obtained with the use of MUSCL-type initial reconstruction and Runge-Kutta time stepping method. After a discussion of the equations solved and of the techniques employed, a series of one and two-dimensional test problems are carried out. To validate the method and assess its accuracy, the staggered central and the kinetic flux-vector splitting schemes are also applied to the same model. The scheme is robust and efficient. Its results are comparable to those obtained from the sophisticated algorithms, even in the case of highly relativistic two-dimensional test problems. PMID:26070067
Bayesian Estimates of the Sensitivity of Orographic Precipitation to Upwind Sounding and Terrain
NASA Astrophysics Data System (ADS)
Delle Monache, L.; Rotunno, R.; Miglietta, M.; Tushaus, S.; Posselt, D. J.
2013-12-01
Orographic precipitation is a key component in the hydrologic cycle, and in many locations is associated with fresh water availability and/or incidence of flooding. The distribution, intensity, and amount of orographically generated rain and snowfall are closely related to the width and slope of the terrain, as well as the stability, temperature, and vapor content of the air upwind of the mountain. While the differences in precipitation associated with changes in any one or two of these variables are well known, the joint interaction between all of them is what ultimately results in the observed spatial and temporal distribution of rain and/or snow. Numerical models are useful tools in the examination of orographic precipitation; however, the number of simulations required to exhaustively explore all possible combinations of environment and terrain have rendered such an exercise computationally intractable. This presentation demonstrates how a Bayesian inverse model (a Markov chain Monte Carlo (MCMC) algorithm) can be used to explore sensitivity of orographic precipitation to variability in environment and topography. Specifically, the MCMC algorithm is used to generate those combinations of sounding and terrain consistent with a given precipitation distribution, without the need for an exhaustive perturbation analysis. The interactions between environment and topography are explored, as are the sensitivities inherent in the system.
Application of Central Upwind Scheme for Solving Special Relativistic Hydrodynamic Equations
Yousaf, Muhammad; Ghaffar, Tayabia; Qamar, Shamsul
2015-01-01
The accurate modeling of various features in high energy astrophysical scenarios requires the solution of the Einstein equations together with those of special relativistic hydrodynamics (SRHD). Such models are more complicated than the non-relativistic ones due to the nonlinear relations between the conserved and state variables. A high-resolution shock-capturing central upwind scheme is implemented to solve the given set of equations. The proposed technique uses the precise information of local propagation speeds to avoid the excessive numerical diffusion. The second order accuracy of the scheme is obtained with the use of MUSCL-type initial reconstruction and Runge-Kutta time stepping method. After a discussion of the equations solved and of the techniques employed, a series of one and two-dimensional test problems are carried out. To validate the method and assess its accuracy, the staggered central and the kinetic flux-vector splitting schemes are also applied to the same model. The scheme is robust and efficient. Its results are comparable to those obtained from the sophisticated algorithms, even in the case of highly relativistic two-dimensional test problems. PMID:26070067
A Five-Parameter Wind Field Estimation Method Based on Spherical Upwind Lidar Measurements
NASA Astrophysics Data System (ADS)
Kapp, S.; Kühn, M.
2014-12-01
Turbine mounted scanning lidar systems of focussed continuous-wave type are taken into consideration to sense approaching wind fields. The quality of wind information depends on the lidar technology itself but also substantially on the scanning technique and reconstruction algorithm. In this paper a five-parameter wind field model comprising mean wind speed, vertical and horizontal linear shear and homogeneous direction angles is introduced. A corresponding parameter estimation method is developed based on the assumption of upwind lidar measurements scanned over spherical segments. As a main advantage of this method all relevant parameters, in terms of wind turbine control, can be provided. Moreover, the ability to distinguish between shear and skew potentially increases the quality of the resulting feedforward pitch angles when compared to three-parameter methods. It is shown that minimal three measurements, each in turn from two independent directions are necessary for the application of the algorithm, whereas simpler measurements, each taken from only one direction, are not sufficient.
Development of an upwind, finite-volume code with finite-rate chemistry
NASA Technical Reports Server (NTRS)
Molvik, Gregory A.
1995-01-01
Under this grant, two numerical algorithms were developed to predict the flow of viscous, hypersonic, chemically reacting gases over three-dimensional bodies. Both algorithms take advantage of the benefits of upwind differencing, total variation diminishing techniques and of a finite-volume framework, but obtain their solution in two separate manners. The first algorithm is a zonal, time-marching scheme, and is generally used to obtain solutions in the subsonic portions of the flow field. The second algorithm is a much less expensive, space-marching scheme and can be used for the computation of the larger, supersonic portion of the flow field. Both codes compute their interface fluxes with a temporal Riemann solver and the resulting schemes are made fully implicit including the chemical source terms and boundary conditions. Strong coupling is used between the fluid dynamic, chemical and turbulence equations. These codes have been validated on numerous hypersonic test cases and have provided excellent comparison with existing data. This report summarizes the research that took place from August 1,1994 to January 1, 1995.
An upwind vertex centred Finite Volume solver for Lagrangian solid dynamics
NASA Astrophysics Data System (ADS)
Aguirre, Miquel; Gil, Antonio J.; Bonet, Javier; Lee, Chun Hean
2015-11-01
A vertex centred Jameson-Schmidt-Turkel (JST) finite volume algorithm was recently introduced by the authors (Aguirre et al., 2014 [1]) in the context of fast solid isothermal dynamics. The spatial discretisation scheme was constructed upon a Lagrangian two-field mixed (linear momentum and the deformation gradient) formulation presented as a system of conservation laws [2-4]. In this paper, the formulation is further enhanced by introducing a novel upwind vertex centred finite volume algorithm with three key novelties. First, a conservation law for the volume map is incorporated into the existing two-field system to extend the range of applications towards the incompressibility limit (Gil et al., 2014 [5]). Second, the use of a linearised Riemann solver and reconstruction limiters is derived for the stabilisation of the scheme together with an efficient edge-based implementation. Third, the treatment of thermo-mechanical processes through a Mie-Grüneisen equation of state is incorporated in the proposed formulation. For completeness, the study of the eigenvalue structure of the resulting system of conservation laws is carried out to demonstrate hyperbolicity and obtain the correct time step bounds for non-isothermal processes. A series of numerical examples are presented in order to assess the robustness of the proposed methodology. The overall scheme shows excellent behaviour in shock and bending dominated nearly incompressible scenarios without spurious pressure oscillations, yielding second order of convergence for both velocities and stresses.
High-Order Semi-Discrete Central-Upwind Schemes for Multi-Dimensional Hamilton-Jacobi Equations
NASA Technical Reports Server (NTRS)
Bryson, Steve; Levy, Doron; Biegel, Bryan (Technical Monitor)
2002-01-01
We present the first fifth order, semi-discrete central upwind method for approximating solutions of multi-dimensional Hamilton-Jacobi equations. Unlike most of the commonly used high order upwind schemes, our scheme is formulated as a Godunov-type scheme. The scheme is based on the fluxes of Kurganov-Tadmor and Kurganov-Tadmor-Petrova, and is derived for an arbitrary number of space dimensions. A theorem establishing the monotonicity of these fluxes is provided. The spacial discretization is based on a weighted essentially non-oscillatory reconstruction of the derivative. The accuracy and stability properties of our scheme are demonstrated in a variety of examples. A comparison between our method and other fifth-order schemes for Hamilton-Jacobi equations shows that our method exhibits smaller errors without any increase in the complexity of the computations.
High-Order Semi-Discrete Central-Upwind Schemes for Multi-Dimensional Hamilton-Jacobi Equations
NASA Technical Reports Server (NTRS)
Bryson, Steve; Levy, Doron; Biegel, Bran R. (Technical Monitor)
2002-01-01
We present high-order semi-discrete central-upwind numerical schemes for approximating solutions of multi-dimensional Hamilton-Jacobi (HJ) equations. This scheme is based on the use of fifth-order central interpolants like those developed in [1], in fluxes presented in [3]. These interpolants use the weighted essentially nonoscillatory (WENO) approach to avoid spurious oscillations near singularities, and become "central-upwind" in the semi-discrete limit. This scheme provides numerical approximations whose error is as much as an order of magnitude smaller than those in previous WENO-based fifth-order methods [2, 1]. Thee results are discussed via examples in one, two and three dimensions. We also pregnant explicit N-dimensional formulas for the fluxes, discuss their monotonicity and tl!e connection between this method and that in [2].
Hosangadi, A.; Sinha, N.; Dash, S.M. )
1992-01-01
A new Eulerian particulate solver whose numerical formulation is compatible with the numerics in state-of-the-art finite-volume upwind/implicit gas dynamic computer codes is presented. The heat transfer, drag, thermodynamic, and phase-change procedures in this code are derived from earlier, well established data fits and procedures. Performance for numerous flow problems with one- and two-way coupling is quite good. The solutions are nonoscillatory and robust and conserve flux balances very well. 18 refs.
Analytis, G.Th.
1995-09-01
As is well-known, both TRAC-BF1 and TRAC-PF are using the first upwind scheme when finite-differencing the phasic momentum equations. In contrast, RELAP5 uses the second upwind which is less diffusive. In this work, we shall assess the differences between the two schemes with our modified version of RELAP5/MOD3 by analyzing some transients of interest. These will include the LOFT LP-LB-1 and LOBI small break LOCA (SB-LOCA) BL34 tests, and a commercial PWR 200% hypothetical large break LOCA (LB-LOCA). In particular, we shall show that for some of these transients, the employment of the first upwind scheme results in significantly different code predictions than the ones obtained when the second upwind scheme is used.
Development of a 3-D upwind PNS code for chemically reacting hypersonic flowfields
NASA Technical Reports Server (NTRS)
Tannehill, J. C.; Wadawadigi, G.
1992-01-01
Two new parabolized Navier-Stokes (PNS) codes were developed to compute the three-dimensional, viscous, chemically reacting flow of air around hypersonic vehicles such as the National Aero-Space Plane (NASP). The first code (TONIC) solves the gas dynamic and species conservation equations in a fully coupled manner using an implicit, approximately-factored, central-difference algorithm. This code was upgraded to include shock fitting and the capability of computing the flow around complex body shapes. The revised TONIC code was validated by computing the chemically-reacting (M(sub infinity) = 25.3) flow around a 10 deg half-angle cone at various angles of attack and the Ames All-Body model at 0 deg angle of attack. The results of these calculations were in good agreement with the results from the UPS code. One of the major drawbacks of the TONIC code is that the central-differencing of fluxes across interior flowfield discontinuities tends to introduce errors into the solution in the form of local flow property oscillations. The second code (UPS), originally developed for a perfect gas, has been extended to permit either perfect gas, equilibrium air, or nonequilibrium air computations. The code solves the PNS equations using a finite-volume, upwind TVD method based on Roe's approximate Riemann solver that was modified to account for real gas effects. The dissipation term associated with this algorithm is sufficiently adaptive to flow conditions that, even when attempting to capture very strong shock waves, no additional smoothing is required. For nonequilibrium calculations, the code solves the fluid dynamic and species continuity equations in a loosely-coupled manner. This code was used to calculate the hypersonic, laminar flow of chemically reacting air over cones at various angles of attack. In addition, the flow around the McDonnel Douglas generic option blended-wing-body was computed and comparisons were made between the perfect gas, equilibrium air, and the
Volume 2: Explicit, multistage upwind schemes for Euler and Navier-Stokes equations
NASA Technical Reports Server (NTRS)
Elmiligui, Alaa; Ash, Robert L.
1992-01-01
The objective of this study was to develop a high-resolution-explicit-multi-block numerical algorithm, suitable for efficient computation of the three-dimensional, time-dependent Euler and Navier-Stokes equations. The resulting algorithm has employed a finite volume approach, using monotonic upstream schemes for conservation laws (MUSCL)-type differencing to obtain state variables at cell interface. Variable interpolations were written in the k-scheme formulation. Inviscid fluxes were calculated via Roe's flux-difference splitting, and van Leer's flux-vector splitting techniques, which are considered state of the art. The viscous terms were discretized using a second-order, central-difference operator. Two classes of explicit time integration has been investigated for solving the compressible inviscid/viscous flow problems--two-state predictor-corrector schemes, and multistage time-stepping schemes. The coefficients of the multistage time-stepping schemes have been modified successfully to achieve better performance with upwind differencing. A technique was developed to optimize the coefficients for good high-frequency damping at relatively high CFL numbers. Local time-stepping, implicit residual smoothing, and multigrid procedure were added to the explicit time stepping scheme to accelerate convergence to steady-state. The developed algorithm was implemented successfully in a multi-block code, which provides complete topological and geometric flexibility. The only requirement is C degree continuity of the grid across the block interface. The algorithm has been validated on a diverse set of three-dimensional test cases of increasing complexity. The cases studied were: (1) supersonic corner flow; (2) supersonic plume flow; (3) laminar and turbulent flow over a flat plate; (4) transonic flow over an ONERA M6 wing; and (5) unsteady flow of a compressible jet impinging on a ground plane (with and without cross flow). The emphasis of the test cases was validation of
Dielectric relaxation time spectroscopy.
Paulson, K S; Jouravleva, S; McLeod, C N
2000-11-01
A new mathematical method is developed to recover the permittivity relaxation spectrum of living tissue from measurements of the real and imaginary parts of the impedance. Aiming to derive information about electrical properties of living tissue without the prior selection of any impedance model, the procedure calculates the relaxation time distribution. It provides new characteristic independent parameters: time constants, their distribution, and the amplitudes of the associated dispersion. As the beta-dispersion is the most important in the area of electrical impedance spectroscopy of tissue, the paper gives an estimate of the essential frequency range to cover the whole relaxation spectrum in that area. Results are presented from both simulation and known lumped--constant element circuit. PMID:11077745
RELAX: detecting relaxed selection in a phylogenetic framework.
Wertheim, Joel O; Murrell, Ben; Smith, Martin D; Kosakovsky Pond, Sergei L; Scheffler, Konrad
2015-03-01
Relaxation of selective strength, manifested as a reduction in the efficiency or intensity of natural selection, can drive evolutionary innovation and presage lineage extinction or loss of function. Mechanisms through which selection can be relaxed range from the removal of an existing selective constraint to a reduction in effective population size. Standard methods for estimating the strength and extent of purifying or positive selection from molecular sequence data are not suitable for detecting relaxed selection, because they lack power and can mistake an increase in the intensity of positive selection for relaxation of both purifying and positive selection. Here, we present a general hypothesis testing framework (RELAX) for detecting relaxed selection in a codon-based phylogenetic framework. Given two subsets of branches in a phylogeny, RELAX can determine whether selective strength was relaxed or intensified in one of these subsets relative to the other. We establish the validity of our test via simulations and show that it can distinguish between increased positive selection and a relaxation of selective strength. We also demonstrate the power of RELAX in a variety of biological scenarios where relaxation of selection has been hypothesized or demonstrated previously. We find that obligate and facultative γ-proteobacteria endosymbionts of insects are under relaxed selection compared with their free-living relatives and obligate endosymbionts are under relaxed selection compared with facultative endosymbionts. Selective strength is also relaxed in asexual Daphnia pulex lineages, compared with sexual lineages. Endogenous, nonfunctional, bornavirus-like elements are found to be under relaxed selection compared with exogenous Borna viruses. Finally, selection on the short-wavelength sensitive, SWS1, opsin genes in echolocating and nonecholocating bats is relaxed only in lineages in which this gene underwent pseudogenization; however, selection on the functional
Relaxation techniques for stress
... Know. February 2013. Available at: nccih.nih.gov/health/stress/relaxation.htm . Accessed September 21, 2015. National Center ... A.D.A.M. Editorial team. Related MedlinePlus Health Topics Stress Browse the Encyclopedia A.D.A.M., Inc. ...
Application of program LAURA to perfect gas shock tube flows: A parametric study
NASA Technical Reports Server (NTRS)
Mitterer, K. F.; Mitcheltree, R. A.; Gnoffo, P. A.
1992-01-01
The Langley Aerothermodynamic Upwind Relaxation Algorithm (LAURA) was originally developed to solve steady-flow problems. The desire to validate the algorithm with shock tube experimental data motivated the development of a time-accurate version of the LAURA code. The current work presents a test of the Algorithm. Computational results are compared with the exact solution for a simple shock tube case. The parameters examined are Courant number, relaxation sweeps, grid spacing, and the inviscid relaxation factor. The results of the study indicate that LAURA is capable of producing accurate solutions when appropriate values are used for each parameter.
Hu, L H; Peng, W; Huo, R
2008-01-15
In case of a tunnel fire, toxic gas and smoke particles released are the most fatal contaminations. It is important to supply fresh air from the upwind side to provide a clean and safe environment upstream from the fire source for people evacuation. Thus, the critical longitudinal wind velocity for arresting fire induced upwind gas and smoke dispersion is a key criteria for tunnel safety design. Former studies and thus, the models built for estimating the critical wind velocity are all arbitrarily assuming that the fire takes place at the centre of the tunnel. However, in many real cases in road tunnels, the fire originates near the sidewall. The critical velocity of a near-wall fire should be different with that of a free-standing central fire due to their different plume entrainment process. Theoretical analysis and CFD simulation were performed in this paper to estimate the critical velocity for the fire near the sidewall. Results showed that when fire originates near the sidewall, it needs larger critical velocity to arrest the upwind gas and smoke dispersion than when fire at the centre. The ratio of critical velocity of a near-wall fire to that of a central fire was ideally estimated to be 1.26 by theoretical analysis. Results by CFD modelling showed that the ratio decreased with the increase of the fire size till near to unity. The ratio by CFD modelling was about 1.18 for a 500kW small fire, being near to and a bit lower than the theoretically estimated value of 1.26. However, the former models, including those of Thomas (1958, 1968), Dangizer and Kenndey (1982), Oka and Atkinson (1995), Wu and Barker (2000) and Kunsch (1999, 2002), underestimated the critical velocity needed for a fire near the tunnel sidewall. PMID:17544576
NASA Astrophysics Data System (ADS)
Sánchez Burillo, Guillermo; Beguería, Santiago; Latorre, Borja; Burguete, Javier
2014-05-01
Debris flows, snow and rock avalanches, mud and earth flows are often modeled by means of a particular realization of the so called shallow water equations (SWE). Indeed, a number of simulation models have been already developed [1], [2], [3], [4], [5], [6], [7]. Debris flow equations differ from shallow water equations in two main aspects. These are (a) strong bed gradient and (b) rheology friction terms that differ from the traditional SWE. A systematic analysis of the numerical solution of the hyperbolic system of equations rising from the shallow water equations with different rheological laws has not been done. Despite great efforts have been done to deal with friction expressions common in hydraulics (such as Manning friction), landslide rheologies are characterized by more complicated expressions that may deal to unphysical solutions if not treated carefully. In this work, a software that solves the time evolution of sliding masses over complex bed configurations is presented. The set of non- linear equations is treated by means of a first order upwind explicit scheme, and the friction contribution to the dynamics is treated with a suited numerical scheme [8]. In addition, the software incorporates various rheological models to accommodate for different flow types, such as the Voellmy frictional model [9] for rock and debris avalanches, or the Herschley-Bulkley model for debris and mud flows. The aim of this contribution is to release this code as a free, open source tool for the simulation of mass movements, and to encourage the scientific community to make use of it. The code uses as input data the friction coefficients and two input files: the topography of the bed and the initial (pre-failure) position of the sliding mass. In addition, another file with the final (post-event) position of the sliding mass, if desired, can be introduced to be compared with the simulation obtained result. If the deposited mass is given, an error estimation is computed by
NASA Astrophysics Data System (ADS)
Sánchez Burillo, Guillermo; Beguería, Santiago; Latorre, Borja; Burguete, Javier
2014-05-01
Debris flows, snow and rock avalanches, mud and earth flows are often modeled by means of a particular realization of the so called shallow water equations (SWE). Indeed, a number of simulation models have been already developed [1], [2], [3], [4], [5], [6], [7]. Debris flow equations differ from shallow water equations in two main aspects. These are (a) strong bed gradient and (b) rheology friction terms that differ from the traditional SWE. A systematic analysis of the numerical solution of the hyperbolic system of equations rising from the shallow water equations with different rheological laws has not been done. Despite great efforts have been done to deal with friction expressions common in hydraulics (such as Manning friction), landslide rheologies are characterized by more complicated expressions that may deal to unphysical solutions if not treated carefully. In this work, a software that solves the time evolution of sliding masses over complex bed configurations is presented. The set of non- linear equations is treated by means of a first order upwind explicit scheme, and the friction contribution to the dynamics is treated with a suited numerical scheme [8]. In addition, the software incorporates various rheological models to accommodate for different flow types, such as the Voellmy frictional model [9] for rock and debris avalanches, or the Herschley-Bulkley model for debris and mud flows. The aim of this contribution is to release this code as a free, open source tool for the simulation of mass movements, and to encourage the scientific community to make use of it. The code uses as input data the friction coefficients and two input files: the topography of the bed and the initial (pre-failure) position of the sliding mass. In addition, another file with the final (post-event) position of the sliding mass, if desired, can be introduced to be compared with the simulation obtained result. If the deposited mass is given, an error estimation is computed by
A Well-Balanced Central-Upwind Scheme for the 2D Shallow Water Equations on Triangular Meshes
NASA Technical Reports Server (NTRS)
Bryson, Steve; Levy, Doron
2004-01-01
We are interested in approximating solutions of the two-dimensional shallow water equations with a bottom topography on triangular meshes. We show that there is a certain flexibility in choosing the numerical fluxes in the design of semi-discrete Godunov-type central schemes. We take advantage of this fact to generate a new second-order, central-upwind method for the two-dimensional shallow water equations that is well-balanced. We demonstrate the accuracy of our method as well as its balance properties in a variety of examples.
... For Consumers Consumer Information by Audience For Women Hair Dye and Hair Relaxers Share Tweet Linkedin Pin it More sharing ... products. If you have a bad reaction to hair dyes and relaxers, you should: Stop using the ...
Dielectric Relaxation of Hexadeutero Dimethylsulfoxide
NASA Astrophysics Data System (ADS)
Betting, H.; Stockhausen, M.
1999-11-01
The dielectric relaxation parameters of the title substance (DMSO-d6) in its pure liquid state are determined from meas-urements up to 72 GHz at 20°C in comparison to protonated DMSO. While the relaxation strengths do not differ, the relax-ation time of DMSO-d 6 is significantly longer (21.3 ps) than that of DMSO (19.5 ps).
Relaxation in Physical Education Curricula.
ERIC Educational Resources Information Center
Coville, Claudia A.
1979-01-01
A theoretical framework for incorporating relaxation instruction in the physical education curriculum is presented based on the assumption that relaxation is a muscular-skeletal skill benefitting general motor skill acquisition. Theoretical principles, a definition of relaxation, and an analysis of stages of skill development are also used in the…
Relaxation phenomena in disordered systems
NASA Astrophysics Data System (ADS)
Sciortino, F.; Tartaglia, P.
1997-02-01
In this article we discuss how the assumptions of self-similarity imposed on the distribution of independently relaxing modes, as well as on their amplitude and characteristic times, manifest in the global relaxation phenomena. We also review recent applications of such approach to the description of relaxation phenomena in microemulsions and molecular glasses.
A Comparison of Relaxation Strategies.
ERIC Educational Resources Information Center
Matthews, Doris B.
Some researchers argue that all relaxation techniques produce a single relaxation response while others support a specific-effects hypothesis which suggests that progressive relaxation affects the musculoskeletal system and that guided imagery affects cognitive changes. Autogenics is considered a technique which is both somatic and cognitive. This…
NASA Astrophysics Data System (ADS)
Hamon, F. P.; Mallison, B.; Tchelepi, H.
2015-12-01
The systems of algebraic equations arising from implicit (backward-Euler) finite-volume discretization of the conservation laws governing multiphase flow in porous media are quite challenging for nonlinear solvers. In the presence of counter-current flow due to buoyancy, the coupling between flow (pressure) and transport (saturations) is often the cause of nonlinear problems when single-point Phase-Potential Upwinding (PPU) is used. To overcome such convergence problems in practice, the time step is reduced and Newton's method is restarted from the solution at the previous converged time step. Here, we generalize the work of Lee, Efendiev and Tchelepi [Advances in Water Resources, 2015] to propose an Implicit Hybrid Upwinding (IHU) scheme for coupled flow and transport. In the pure transport problem, we show that the numerical flux obtained with IHU is differentiable, monotone and consistent for two and three-phase flow. For coupled flow and transport, we prove saturation physical bounds as well as the existence of a solution to our scheme. Challenging two- and three-phase heterogeneous multi-dimensional numerical tests confirm that the new scheme is non-oscillatory and convergent, and illustrate the superior convergence rate of our IHU-based Newton solver for large time steps.
The Influence of the Upwind Sea Surface Temperature Gradient on a Heavy Rainfall Event During TiMREX
NASA Astrophysics Data System (ADS)
Toy, M. D.; Johnson, R. H.
2012-12-01
A long-lived precipitation area was observed off the southwest coast of Taiwan during the Terrain-influenced Monsoon Rainfall Experiment (TiMREX). Convective cells continually developed offshore and brought heavy amounts of rain to the coastal regions before dying out over land. Previous observational studies have attributed the development of the convective system to the interaction between a southwesterly low-level jet and a cold pool blocked by the downwind terrain. In this study we analyze the role of a strong sea surface temperature (SST) gradient that existed upwind of the island in initiating the offshore convection. Through model simulations with the Weather Research and Forecasting Model (WRF), we show that the airflow over the ocean front modifies the planetary boundary layer so as to enhance low-level convergence and initiate convection. Reducing the SST gradient in the model results in smaller simulated rainfall amounts. This indicates the importance of carefully accounting for the upwind SST pattern in forecasting heavy monsoon rainfall events.
NASA Astrophysics Data System (ADS)
McKenna Neuman, Cheryl
1998-10-01
Most natural surfaces containing non-erodible roughness elements are considerably more complex than those studied in modelling exercises and wind tunnel simulations. Unlike idealized roughness elements, which are uniform in size, shape (i.e., spheres or cylinders) and spacing, natural elements are challenging to measure in 3-dimensional space. Similarly, most deflation lag surfaces, such those as found on sandar and beaches, are spatially heterogeneous open systems in which sediment transport from an external supply is very likely. The development of irregular deflation lag surfaces, and the transport of sediment over these surfaces from an upwind source of sediment, was studied in a series of wind tunnel simulations. Surfaces prepared with crushed gravel and natural beach shingle respond conservatively in terms of the adjustment to the deflation and deposition of sediment. Deflation lag surfaces, prepared with no spacing between the roughness elements (i.e., close packed), demonstrate little to no change in coverage with the introduction of particles from an upwind source. Neither the element type nor the friction velocity affect this outcome. As the center-to-center element spacing increases to 60 mm, infilling of the lag surface eventually is observed, with the element coverage reduced by a factor between 2 and 4. For a given threshold ratio ( Rt), the roughness density ( λ) is smaller than observed in previous simulation studies based on idealized roughness elements.
NASA Astrophysics Data System (ADS)
Londrillo, P.; del Zanna, L.
2004-03-01
We present a general framework to design Godunov-type schemes for multidimensional ideal magnetohydrodynamic (MHD) systems, having the divergence-free relation and the related properties of the magnetic field B as built-in conditions. Our approach mostly relies on the constrained transport (CT) discretization technique for the magnetic field components, originally developed for the linear induction equation, which assures [∇.B]num=0 and its preservation in time to within machine accuracy in a finite-volume setting. We show that the CT formalism, when fully exploited, can be used as a general guideline to design the reconstruction procedures of the B vector field, to adapt standard upwind procedures for the momentum and energy equations, avoiding the onset of numerical monopoles of O(1) size, and to formulate approximate Riemann solvers for the induction equation. This general framework will be named here upwind constrained transport (UCT). To demonstrate the versatility of our method, we apply it to a variety of schemes, which are finally validated numerically and compared: a novel implementation for the MHD case of the second-order Roe-type positive scheme by Liu and Lax [J. Comput. Fluid Dyn. 5 (1996) 133], and both the second- and third-order versions of a central-type MHD scheme presented by Londrillo and Del Zanna [Astrophys. J. 530 (2000) 508], where the basic UCT strategies have been first outlined.
NASA Technical Reports Server (NTRS)
Grossman, B.; Garrett, J.; Cinnella, P.
1989-01-01
Several versions of flux-vector split and flux-difference split algorithms were compared with regard to general applicability and complexity. Test computations were performed using curve-fit equilibrium air chemistry for an M = 5 high-temperature inviscid flow over a wedge, and an M = 24.5 inviscid flow over a blunt cylinder for test computations; for these cases, little difference in accuracy was found among the versions of the same flux-split algorithm. For flows with nonequilibrium chemistry, the effects of the thermodynamic model on the development of flux-vector split and flux-difference split algorithms were investigated using an equilibrium model, a general nonequilibrium model, and a simplified model based on vibrational relaxation. Several numerical examples are presented, including nonequilibrium air chemistry in a high-temperature shock tube and nonequilibrium hydrogen-air chemistry in a supersonic diffuser.
Relaxing music for anxiety control.
Elliott, Dave; Polman, Remco; McGregor, Richard
2011-01-01
The purpose of this investigation was to determine the characteristics of relaxing music for anxiety control. Undergraduate students (N=84) were instructed to imagine themselves in an anxiety producing situation while listening to a selection of 30 music compositions. For each composition, level of relaxation, the factors that either enhanced or detracted from its relaxing potential and the emotional labels attached were assessed. Participants were also asked to state which music components (e.g., tempo, melody) were most conducive to relaxation. Additional information was obtained through the use of a focus group of 6 undergraduate music students. This paper presents details on the characteristics of relaxing-music for anxiety control and emotional labels attached to the relaxing compositions. Furthermore, an importance value has been attached to each of the music components under scrutiny, thus providing an indication of which music components should receive greatest attention when selecting music for anxiety control. PMID:22097099
NASA Astrophysics Data System (ADS)
Macák, Marek; Čunderlík, Róbert; Mikula, Karol; Minarechová, Zuzana
2016-02-01
The paper presents a novel original upwindbased approach for solving the oblique derivative boundary value problem by the finite volume method. In this approach, the oblique derivative boundary condition is interpreted as a stationary advection equation for the unknown disturbing potential. Its approximation is then performed by using the first order upwind scheme taking into account information from inflow parts of the finite volume boundary only. When the numerical scheme is derived, numerical simulations in 2D and 3D domains are performed and the experimental order of convergence of the proposed algorithm is studied. Moreover a comparison with a solution by the central scheme previously used for this kind of problem is performed. Finally we present numerical experiments dealing with the global and local gravity field modelling.
Kim, D.H. . Living System Research Lab.); Jeon, D.Y.; Hahn, S.Y. . Dept. of Electrical Engineering)
1999-05-01
In general, an electromagnetic apparatus such as linear induction motors, MAGLEV vehicles or electromagnetic launchers, involves conducting parts in motion. This paper presents a new algorithm based on the indirect boundary integral equation method to analyze the electromagnetic system with a moving conductor. The proposed algorithm yields relatively stable and accurate solutions because a fundamental Green's function of diffusion type is used which is valid for any value of the Peclet number. In addition, computer memory and computing time for 3D computation can be saved considerably by using the boundary integral equations of minimum order and the singular property of the Green's function. In order to prove these, numerical results obtained by the proposed algorithm and the upwind finite element method are compared with their analytic solutions.
Renormalized reaction and relaxation rates
NASA Astrophysics Data System (ADS)
Gorbachev, Yuriy E.
2016-06-01
Impact of the non-equilibrium on the reaction and relaxation rates (called as generalized relaxation rates - GRR), for the spatially inhomogeneous gas mixture is considered. Discarding the assumption that the 'chemical' part of the collisional integral is a small correction to non-reactive part, the expression for the zero-order GRR is derived. They are represented as a renormalization of the traditional reaction and relaxation rates, which means mixing of all corresponding processes. Thus all reactions and relaxation processes are entangled.
Comet Bursting Through Relaxation
NASA Astrophysics Data System (ADS)
Jacobson, Seth A.; Scheeres, D. J.
2012-10-01
Comets may be excited and occupy non-principal axis (complex) rotation states for a large fraction of their lifetimes. Many comet nuclei have been identified or are suspected to occupy non-principal axis (complex) rotation [Belton 2005, etc.] as well as have evolving rotation rates [Belton 2011, etc.]. Comet orbits drive these rotation states through cycles of excitation due to surface jets and relaxation due to time variable internal stresses that dissipate energy in the anelastic comet interior. Furthermore, relaxation from complex rotation can increase the loads along the symmetry axis of prolate comets. These loads stretch the body along the symmetry axis and may be the cause of the characteristic ``bowling pin’’ shape and eventually may lead to failure. This is an alternative model for comet bursting. Each cycle deposits only a small amount of energy and stress along the axis, but this process is repeated every orbit during which jets are activated. Our model for the evolution of comet nuclei includes torques due to a number of discrete jets located on the surface based on Neishtadt et al. [2002]. The model also includes internal dissipation using an approach developed by Sharma et al. [2005] and Vokrouhlicky et al. [2009]. These equations are averaged over the instantaneous spin state and the heliocentric orbit so the long-term evolution of the comet can be determined. We determine that even after the inclusion of internal dissipation there still exist non-principal axis equilibrium states for certain jet geometries. For ranges of dissipation factors and jet geometries, prolate comets are found to occupy states that have time variable internal loads over long time periods. These periodic loadings along the symmetry axis may lead to ``necking’’ as the body extends along the axis to release the stress and eventually disruption.
NASA Technical Reports Server (NTRS)
Wood, William A., III
2002-01-01
A multi-dimensional upwind fluctuation splitting scheme is developed and implemented for two-dimensional and axisymmetric formulations of the Navier-Stokes equations on unstructured meshes. Key features of the scheme are the compact stencil, full upwinding, and non-linear discretization which allow for second-order accuracy with enforced positivity. Throughout, the fluctuation splitting scheme is compared to a current state-of-the-art finite volume approach, a second-order, dual mesh upwind flux difference splitting scheme (DMFDSFV), and is shown to produce more accurate results using fewer computer resources for a wide range of test cases. A Blasius flat plate viscous validation case reveals a more accurate upsilon-velocity profile for fluctuation splitting, and the reduced artificial dissipation production is shown relative to DMFDSFV. Remarkably, the fluctuation splitting scheme shows grid converged skin friction coefficients with only five points in the boundary layer for this case. The second half of the report develops a local, compact, anisotropic unstructured mesh adaptation scheme in conjunction with the multi-dimensional upwind solver, exhibiting a characteristic alignment behavior for scalar problems. The adaptation strategy is extended to the two-dimensional and axisymmetric Navier-Stokes equations of motion through the concept of fluctuation minimization.
Hassan, Y.A.
1985-02-01
The upwind finite difference scheme used in the COMMIX-1A code to approximate the convective transport term was shown to contribute more artificial diffusivity to the computed temperature field than the total physical diffusivity that could be inferred from experimental observations of the thermal mixing processes. With aim to reduce such numerical error, a new finite difference scheme called the mass-flow-weighted skew-upwind approach was developed, incorporated in COMMIX-1A, and evaluated by comparison with known analytical solutions and experimental data. This report describes the follow-on activities pursued to demonstrate the applicability of the new finite difference scheme to simulations of pressurized thermal shock events. These activities include benchmarking the code against the Creare one-fifth-scale test 51 and simulating mixing in the full-scale cold leg and downcomer of a Babcock and Wilcox plant using both the upwind and the mass-flow-weighted skew-upwind finite difference scehmes. The new scheme is shown to considerably reduce the amount of artifical diffusivity contributed to the computed temperature field by the convective term approximation.
NASA Astrophysics Data System (ADS)
Hansen, D. Flemming; Led, Jens J.
2001-08-01
A new experiment for selective determination of the relaxation rates of fast relaxing NMR signals is presented. The experiment is derived from the conventional inversion recovery experiment by substituting the 180° inversion pulse of this experiment with a signal eliminating relaxation filter (SERF) consisting of three 180° pulses separated by two variable delays, Δ1 and Δ2. The SERF experiment allows a selective suppression of signals with relaxation rates below a given limit while monitoring the relaxation of faster relaxing signals. The experiment was tested on a sample of 20% oxidized plastocyanin from Anabaena variabilis, where the fast exchange of an electron between the reduced (diamagnetic) and the oxidized (paramagnetic) form results in a series of average signals with widely different relaxation rates. To ensure an optimum extraction of information from the experimental data, the relaxation rates were obtained from the SERF experiment by a simultaneous analysis of all the FIDs of the experiment using a fast linear prediction model method developed previously. The reliability of the relaxation rates obtained from the SERF experiment was confirmed by a comparison of the rates with the corresponding rates obtained from a conventional inversion recovery experiment.
Aerothermodynamic radiation studies
NASA Technical Reports Server (NTRS)
Donohue, K.; Reinecke, W. G.; Rossi, D.; Marinelli, W. J.; Krech, R. H.; Caledonia, G. E.
1991-01-01
We have built and made operational a 6 in. electric arc driven shock tube which alloys us to study the non-equilibrium radiation and kinetics of low pressure (0.1 to 1 torr) gases processed by 6 to 12 km/s shock waves. The diagnostic system allows simultaneous monitoring of shock radiation temporal histories by a bank of up to six radiometers, and spectral histories with two optical multi-channel analyzers. A data set of eight shots was assembled, comprising shocks in N2 and air at pressures between 0.1 and 1 torr and velocities of 6 to 12 km/s. Spectrally resolved data was taken in both the non-equilibrium and equilibrium shock regions on all shots. The present data appear to be the first spectrally resolved shock radiation measurements in N2 performed at 12 km/s. The data base was partially analyzed with salient features identified.