Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 523 records · Page 29

Simulation of Wake Vortex Radiometric Detection via Jet Exhaust Proxy

This paper describes an analysis of the potential of an airborne hyperspectral imaging IR instrument to infer wake vortices via turbine jet exhaust as a proxy. The goal was to determine the requirements for an imaging spectrometer or radiometer to effectively detect the exhaust plume, and by inference, the location of the wake vortices. The effort examines the gas spectroscopy of the various major constituents of turbine jet exhaust and their contributions to the modeled detectable radiance. Initially, a theoretical analysis of wake vortex proxy detection by thermal radiation was realized in a series of simulations. The first stage used the SLAB plume model to simulate turbine jet exhaust plume characteristics, including exhaust gas transport dynamics and concentrations. The second stage used these plume characteristics as input to the Line By Line Radiative Transfer Model (LBLRTM) to simulate responses from both an imaging IR hyperspectral spectrometer or radiometer. These numerical simulations generated thermal imagery that was compared with previously reported wake vortex temperature data. This research is a continuation of an effort to specify the requirements for an imaging IR spectrometer or radiometer to make wake vortex measurements. Results of the two-stage simulation will be reported, including instrument specifications for wake vortex thermal detection. These results will be compared with previously reported results for IR imaging spectrometer performance.

Daniels, Taumi S.↗

History and Status of the CDISC Aerodynamic Design Method

This paper will review the development and application of the CDISC aerodynamic design method that began in the early 1980s and continues in active use today. The method uses an iterative, knowledge-based approach that has been coupled with numerous flow solvers utilizing a variety of grid types and levels of flow physics. It has been applied to configurations varying in complexity from 2-D airfoils to full 3-D aircraft at flow conditions ranging from low-speed, high-lift to supersonic cruise. The knowledge-based approach provides a very rapid design process, with designs typically completed in approximately 1-2x the time required for an initial converged analysis of the baseline. The evolution of the method as it tracked the expansive growth of grid and flow solver technology, along with some major applications over 5 decades, are given, followed by a description of the philosophy and process used in design. Finally, an example of its use in a preliminary design environment is given.

CDISC↗

History and Status of the CDISC Aerodynamic Design Method

This paper will review the development and application of the CDISC aerodynamic design method that began in the early 1980s and continues in active use today. The method uses an iterative, knowledge-based approach that has been coupled with numerous flow solvers utilizing a variety of grid types and levels of flow physics. It has been applied to configurations varying in complexity from 2-D airfoils to full 3-D aircraft at flow conditions ranging from low-speed, high-lift to supersonic cruise. The knowledge-based approach provides a very rapid design process, with designs typically completed in approximately 1-2x the time required for an initial converged analysis of the baseline. The evolution of the method as it tracked the expansive growth of grid and flow solver technology, along with some major applications over 5 decades, are given, followed by a description of the philosophy and process used in design. Finally, an example of its use in a preliminary design environment is given.

CDISC↗

Three Dimensional Constraint Effects on the Estimated (Delta)CTOD during the Numerical Simulation of Different Fatigue Threshold Testing Techniques

Variation in constraint through the thickness of a specimen effects the cyclic crack-tip-opening displacement (DELTA CTOD). DELTA CTOD is a valuable measure of crack growth behavior, indicating closure development, constraint variations and load history effects. Fatigue loading with a continual load reduction was used to simulate the load history associated with fatigue crack growth threshold measurements. The constraint effect on the estimated DELTA CTOD is studied by carrying out three-dimensional elastic-plastic finite element simulations. The analysis involves numerical simulation of different standard fatigue threshold test schemes to determine how each test scheme affects DELTA CTOD. The American Society for Testing and Materials (ASTM) prescribes standard load reduction procedures for threshold testing using either the constant stress ratio (R) or constant maximum stress intensity (K(sub max)) methods. Different specimen types defined in the standard, namely the compact tension, C(T), and middle cracked tension, M(T), specimens were used in this simulation. The threshold simulations were conducted with different initial K(sub max) values to study its effect on estimated DELTA CTOD. During each simulation, the DELTA CTOD was estimated at every load increment during the load reduction procedure. Previous numerical simulation results indicate that the constant R load reduction method generates a plastic wake resulting in remote crack closure during unloading. Upon reloading, this remote contact location was observed to remain in contact well after the crack tip was fully open. The final region to open is located at the point at which the load reduction was initiated and at the free surface of the specimen. However, simulations carried out using the constant Kmax load reduction procedure did not indicate remote crack closure. Previous analysis results using various starting K(sub max) values and different load reduction rates have indicated DELTA CTOD is independent of specimen size. A study of the effect of specimen thickness and geometry on the measured DELTA CTOD for various load reduction procedures and its implication in the estimation of fatigue crack growth threshold values is discussed.

Seshadri, Banavara R.↗

Unified approach for incompressible flows

An unified approach for solving both compressible and incompressible flows was investigated in this study. The difference in CFD code development between incompressible and compressible flows is due to the mathematical characteristics. However, if one can modify the continuity equation for incompressible flows by introducing pseudocompressibility, the governing equations for incompressible flows would have the same mathematical characters as compressible flows. The application of a compressible flow code to solve incompressible flows becomes feasible. Among numerical algorithms developed for compressible flows, the Centered Total Variation Diminishing (CTVD) schemes possess better mathematical properties to damp out the spurious oscillations while providing high-order accuracy for high speed flows. It leads us to believe that CTVD schemes can equally well solve incompressible flows. In this study, the governing equations for incompressible flows include the continuity equation and momentum equations. The continuity equation is modified by adding a time-derivative of the pressure term containing the artificial compressibility. The modified continuity equation together with the unsteady momentum equations forms a hyperbolic-parabolic type of time-dependent system of equations. The continuity equation is modified by adding a time-derivative of the pressure term containing the artificial compressibility. The modified continuity equation together with the unsteady momentum equations forms a hyperbolic-parabolic type of time-dependent system of equations. Thus, the CTVD schemes can be implemented. In addition, the boundary conditions including physical and numerical boundary conditions must be properly specified to obtain accurate solution. The CFD code for this research is currently in progress. Flow past a circular cylinder will be used for numerical experiments to determine the accuracy and efficiency of the code before applying this code to more specific applications.

Chang, Tyne-Hsien↗

User's Manual for LINER: FORTRAN Code for the Numerical Simulation of Plane Wave Propagation in a Lined Two-Dimensional Channel

LINER is a system of Fortran 77 codes which performs a 2D analysis of acoustic wave propagation and noise suppression in a rectangular channel with a continuous liner at the top wall. This new implementation is designed to streamline the usage of the several codes making up LINER, resulting in a useful design tool. Major input parameters are placed in two main data files, input.inc and nurn.prm. Output data appear in the form of ASCII files as well as a choice of GNUPLOT graphs. Section 2 briefly describes the physical model. Section 3 discusses the numerical methods; Section 4 gives a detailed account of program usage, including input formats and graphical options. A sample run is also provided. Finally, Section 5 briefly describes the individual program files.

Reichert, R, S.↗

Some Recent Observations on the Burning of Isolated N-Heptane and Alcohol Droplets

In a joint program involving Prof. F.A. Williams of the University of California, San Diego and Dr. V. Nayagam of the National Center for Microgravity Research on Combustion and Fluid Dynamics, the combustion of liquid fuel droplets of n-heptane, n-decane, methanol, methanol-water, ethanol and ethanol-water having initial diameters between about 1 mm and 6 mm continues to be studied. The objectives of the work are to improve fundamental knowledge of droplet combustion dynamics for pure fuels and fuel-water mixtures through microgravity experiments and theoretical analyses. The Princeton contributions to the collaborative program supports the engineering design, data analysis, and data interpretation requirements for the study of initially single component, spherically symmetric, isolated droplet combustion studies through experiments and numerical modeling. The complementary UCSD contributions apply asymptotic theoretical analyses and are described in the published literature and in a companion communication in this conference. The combined program continues to focus on analyses of results obtained from Fiber Supported Droplet Combustion (FSDC) experiments (FSDC-2, STS- 94) conducted with the above fuels in shuttle cabin air and Droplet Combustion Experiment (DCE) data obtained for unsupported and fiber supported droplets of n-heptane in Helium-Oxygen mixtures and cabin air (STS-83, STS-94). The program is preparing for a second DCE experimental mission using methanol/methanol-water as fuels and helium-oxygen-nitrogen environments. DCE-2 is to be conducted aboard the International Space Station. Emphases of recent Princeton work are on the study of simple alcohols (methanol, ethanol) and alcohol/water mixtures as fuels, with time-dependent measurements of drop size, flame-standoff, liquid-phase composition, and finally, extinction. Ground based experiments have included bench-scale studies at Princeton and collaborative experimental studies in the 2.2 second drop tower at NASA-Glenn Research Center.

Dryer, F. L.↗

Efficient numerical techniques for complex fluid flows

The aim of the present research program is the development of more efficient and reliable calculation schemes for the coupled momentum and continuity equations. The resulting schemes would significantly reduce the expense of computing complex flows such as those in combustion chambers, gas turbines, and heat exchangers.

Patankar, Suhas V.↗

Numerical Study of Steady and Unsteady Canard-Wing-Body Aerodynamics

The use of canards in advanced aircraft for control and improved aerodynamic performance is a topic of continued interest and research. In addition to providing maneuver control and trim, the influence of canards on wing aerodynamics can often result in increased maximum lift and decreased trim drag. In many canard-configured aircraft, the main benefits of canards are realized during maneuver or other dynamic conditions. Therefore, the detailed study and understanding of canards requires the accurate prediction of the non-linear unsteady aerodynamics of such configurations. For close-coupled canards, the unsteady aerodynamic performance associated with the canard-wing interaction is of particular interest. The presence of a canard in close proximity to the wing results in a highly coupled canard-wing aerodynamic flowfield which can include downwash/upwash effects, vortex-vortex interactions and vortex-surface interactions. For unsteady conditions, these complexities of the canard-wing flowfield are further increased. The development and integration of advanced computational technologies provide for the time-accurate Navier-Stokes simulations of the steady and unsteady canard-wing-body flox,fields. Simulation, are performed for non-linear flight regimes at transonic Mach numbers and for a wide range of angles of attack. For the static configurations, the effects of canard positioning and fixed deflection angles on aerodynamic performance and canard-wing vortex interaction are considered. For non-static configurations, the analyses of the canard-wing body flowfield includes the unsteady aerodynamics associated with pitch-up ramp and pitch oscillatory motions of the entire geometry. The unsteady flowfield associated with moving canards which are typically used as primary control surfaces are considered as well. The steady and unsteady effects of the canard on surface pressure integrated forces and moments, and canard-wing vortex interaction are presented in detail including the effects of the canard on the static and dynamic stability characteristics. The current study provides an understanding of the steady and unsteady canard-wing-body flowfield. Emphasis is placed on the effects of the canard on aerodynamic performance as well as the detailed flow physics of the canard-wing flowfield interactions. The computational tools developed to accurately predict the time-accurate flowfield of moving canards provides for the capability of coupled fluids-controls simulations desired in the detailed design and analysis of advanced aircraft.

Eugene, L. Tu↗

Numerical Simulations of a Co-Axial Supersonic-Combusting Free-Jet Experiment

CFD calculations using the Reynolds-averaged Navier-Stokes equations coupled with species continuity equations have been made for a supersonic coaxial-jet CFD-validation experiment to determine the sensitivity of the external flowfield to the main-nozzle exit profile. Four different nozzle exit profiles were used in the study: a uniform profile, one computed using only the nozzle geometry, one computed using the nozzle geometry and part of the upstream facility combustor, and one using the nozzle and the full facility combustor. Two cases were examined using the four profiles: a non-reacting case without coflow and a reacting case with hydrogen coflow. Results show that the nozzle exit profile has a significant effect on the external flowfield. The uniform profile produced the longest jet while the profile created with the full combustor produced the shortest jet. The nozzle-only and part-combustor profiles fell between the other two profiles. The reacting flow was found to be more sensitive to the nozzle exit profile since it affects the downstream mixing and combustion. These calculations indicate the importance of properly setting the nozzle-exit profile for this type of calculation.

Gaffney, Richard L.↗

The Wallops Flight Facility Rapid Response Range Operations Initiative

While the dominant focus on short response missions has appropriately centered on the launch vehicle and spacecraft, often overlooked or afterthought phases of these missions have been launch site operations and the activities of launch range organizations. Throughout the history of organized spaceflight, launch ranges have been the bane of flight programs as the source of expense, schedule delays, and seemingly endless requirements. Launch Ranges provide three basic functions: (1) provide an appropriate geographical location to meet orbital other mission trajectory requirements, (2) provide project services such as processing facilities, launch complexes, tracking and data services, and expendable products, and (3) assure safety and property protection to participating personnel and third-parties. The challenge with which launch site authorities continuously struggle, is the inherent conflict arising from projects whose singular concern is execution of their mission, and the range s need to support numerous simultaneous customers. So, while tasks carried out by a launch range committed to a single mission pale in comparison to efforts of a launch vehicle or spacecraft provider and could normally be carried out in a matter of weeks, major launch sites have dozens of active projects separate sponsoring organizations. Accommodating the numerous tasks associated with each mission, when hardware failures, weather, maintenance requirements, and other factors constantly conspire against the range resource schedulers, make the launch range as significant an impediment to responsive missions as launch vehicles and their cargo. The obvious solution to the launch site challenge was implemented years ago when the Department of Defense simply established dedicated infrastructure and personnel to dedicated missions, namely an Inter Continental Ballistic Missile. This however proves to be prohibitively expensive for all but the most urgent of applications. So the challenge becomes how can a launch site provide acceptably responsive mission services to a particular customer without dedicating extensive resources and while continuing to serve other projects? NASA's Wallops Flight Facility (WFF) is pursuing solutions to exactly this challenge. NASA, in partnership with the Virginia Commercial Space Flight Authority, has initiated the Rapid Response Range Operations Initiative (R3Ops). R3Ops is a multi-phased effort to incrementally establish and demonstrate increasingly responsive launch operations, with an ultimate goal of providing ELV-class services in a maximum of 7-10 days from initial notification routinely, and shorter schedules possible with committed resources. This target will be pursued within the reality of simultaneous concurrent programs, and ideally, largely independent of specialized flight system configurations. WFF has recently completed Phase 1 of R3Ops, an in-depth collection (through extensive expert interviews) and software modeling of individual steps by various range disciplines. This modeling is now being used to identify existing inefficiencies in current procedures, to identify bottlenecks, and show interdependencies. Existing practices are being tracked to provide a baseline to benchmark against as new procedures are implemented. This paper will describe in detail the philosophies behind WFF's R3Ops, the data collected and modeled in Phase 1, and strategies for meeting responsive launch requirements in a multi-user range environment planned for subsequent phases of this initiative.

Underwood, Bruce E.↗

Pressure Drop Across Woven Screens Under Uniform and Nonuniform Flow Conditions

Tests were conducted to determine the experimental pressure drop and velocity data for water flowing through woven screens. The types of materials used are dutch twill and square weave fabrics. Pressure drop measures were made at four locations in a rectangular channel. The data are presented as change in pressure compared with the average entry velocity and the numerical relationship is determined by dividing the volumetric flow rate by the screen area open to flow. The equations of continuity and momentum are presented. A computer program listing an extension of a theoretical model and data from that computer program are included.

Ludewig, M.↗

Interaction of NaCl/g/ and HCl/g/ with condensed Na2SO4

Na2SO4(l)-NaCl(g) interactions were studied at a total pressure of one atmosphere of air or oxygen for various temperatures of Na2SO4(l) and for various partial pressures of NaCl(g) and H2O(g). Mass spectrometric sampling techniques were used to identify and monitor gas phase species. Continuous recording thermomicrogravimetric measurements were conducted to determine condensed phase weight change rates. Experimental measurements were supplemented with thermodynamic calculations. Numerous experiments were performed at sample temperatures of 900 and 1000 C with 300 ppm NaCl(g). In these experiments, the reproducibility of the Na2SO4 vaporization weight loss rate and initial weight gain upon addition of NaCl(g) were found to be satisfactory. It was found that the addition of NaCl(g) to air flowing over Na2SO4(l) at 900 and 1000 C enhances the rate of weight loss of the Na2SO4(l). This enhancement increases when H2O(g) is also added to the air flow.

Stearns, C. A.↗

The interplanetary electric field, cleft currents and plasma convection in the polar caps

The relationship between the pattern of plasma convection in the polar cleft and the dynamics of the interplanetary electric field (IEF) is examined theoretically. It is shown that owing to the geometrical properties of the magnetosphere, the East-West component of the IEF will drive field-aligned currents which connect to the ionosphere at points lying on either side of noon, while currents associated with the North-South component of the IEF will connect the two polar caps as sheet currents, also centered at 12 MLT. In order to describe the consequences of the Interplanetary Magnetic Field (IMF) effects upon high-latitude electric fields and convection patterns, a series of numerical simulations was carried out. The simulations were based on a solution to the steady-state equation of current continuity in a height-integrated ionospheric current. The simulations demonstrate that a simple hydrodynamical model can account for the narrow 'throats' of strong dayside antisunward convection observed during periods of southward interplanetary IMF drift, as well as the sunward convection observed during periods of strongly northward IMF drift.

Banks, P. M.↗

Bulk Formation of Metallic Glasses and Amorphous Silicon from the Melt

By using metallic glass compositions with a high relative glass transition temperature, such as Pd40Ni40P20, homogeneous nucleation also becomes negligible. Large (5g) masses of this alloys were obtained using a molten B2O3 flux. Presently, bulk glass formation in iron based glasses is being investigated. It is expected that if an undercooling of about 250K can be achieved in a Ge or Si melt, formation of the amorphous semiconductor phase (rather than the crystal) may be kinetically favored. The volumetric behavior of undercooled liquid Ga droplet dispersion is investigated by dilatometry. A theoretical model (both analytical and numerical) was developed for transient nucleation in glass forming melts. The model, originally designed for isothermal conditions, was extended to continuous quenching. It is being applied to glass formation in various metallic and oxide systems. A further refinement will be the inclusion of diffusion controlled interfacial rearrangements governing the growth of the crystal embryos.

Spaepen, F.↗

Viscous shock-layer analysis of hypersonic flows over long slender bodies

Numerical solutions from the time-steady viscous shock-layer equations are presented for the hypersonic laminar and turbulent flow of a perfect gas over long slender bodies. These results are obtained from a spatial-marching implicit finite-difference technique. Detailed comparisons have been made with other predictions and experimental data to assess the accuracy of the present numerical technique, especially for slender-body flows. Results from the present method show that coupling the normal momentum and continuity equations and the use of the Vigneron pressure condition in the subsonic nose region give quite accurate and stable results. These results compare (with some exceptions) favorably with those obtained using PNS and other VSL methods.

Gupta, R. N.↗

Emittance bounds for transient radiative cooling of a scattering rectangular region

Transient cooling was analyzed for a two-dimensional gray rectangular region that emits, absorbs, and isotropically scatters radiation. The region, initially at uniform temperature, is placed in surroundings at a much lower temperature. The cooling analysis yields two simultaneous equations for the transient temperature and scattering source function distributions. The region starts to cool with the medium at uniform temperature. For this condition, the local emittance around the perimeter was obtained by a numerical solution of the radiative equations, and the rate of heat loss from the entire medium was evaluated. As cooling continues, the emissive ability decreases because of the lower temperatures in the outer regions of the rectangle. Based on some previous work, a lower limit for the transient overall emittance was found by obtaining a similarity solution. For some ranges of optical thickness and scattering albedo, the lower limit is only a small amount below the initial emittance. For these conditions, the initial emittance can be used to compute the entire transient resulting in a considerable simplification.

Siegel, Robert↗

Magmatic intrusions and hydrothermal systems on Mars

We are continuing our investigation of Martian hydrothermal systems and the formation of fluvial valleys on Mars. Here we present our initial numerical modeling results of hydrothermal systems associated with magmatic intrusions on Mars. To model such hydrothermal systems, we consider single, cylindrical intrusions of 4 km height 2 km below the surface. Our preliminary results of modeling hydrothermal systems associated with magmatic intrusions in the Martian environment suggest that such systems, if associated with intrusions of several 10(exp 2) km(exp 3) or larger, in the presence of permeable, water-rich subsurface should be able to provide adequate discharges of water over the time periods needed to form fluvial valleys. However, it is important to note that the amount of water ultimately available for erosion also depends on the surface and subsurface lithology.

Gulick, V. C.↗