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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

PDF turbulence modeling and DNS

The problem of time discontinuity (or jump condition) in the coalescence/dispersion (C/D) mixing model is addressed in probability density function (pdf). A C/D mixing model continuous in time is introduced. With the continuous mixing model, the process of chemical reaction can be fully coupled with mixing. In the case of homogeneous turbulence decay, the new model predicts a pdf very close to a Gaussian distribution, with finite higher moments also close to that of a Gaussian distribution. Results from the continuous mixing model are compared with both experimental data and numerical results from conventional C/D models. The effect of Coriolis forces on compressible homogeneous turbulence is studied using direct numerical simulation (DNS). The numerical method used in this study is an eight order compact difference scheme. Contrary to the conclusions reached by previous DNS studies on incompressible isotropic turbulence, the present results show that the Coriolis force increases the dissipation rate of turbulent kinetic energy, and that anisotropy develops as the Coriolis force increases. The Taylor-Proudman theory does apply since the derivatives in the direction of the rotation axis vanishes rapidly. A closer analysis reveals that the dissipation rate of the incompressible component of the turbulent kinetic energy indeed decreases with a higher rotation rate, consistent with incompressible flow simulations (Bardina), while the dissipation rate of the compressible part increases; the net gain is positive. Inertial waves are observed in the simulation results.

Hsu, A. T.↗
NASA NTRS1991

Retrieval of wind temperature and pressure from single Doppler radar and a numerical model

A 4D data assimilation algorithm to obtain 3D wind and thermodynamic fields from radial velocities and a numerical model is derived. It is a blend of the continuous updating technique of Charney et al. (1969) and a diagnostic pressure and temperature retrieval technique suggested by Gal-Chen (1978). This model is tested only against model-generated data which are viewed as 'data observed from a real atmosphere'. A thermal bubble developing in a dry, neutral environment is simulated. The thermal recovery technique is found to be a necessary procedure to obtain successful data assimilation results. The present assimilation method is found to be capable of reducing the 'observational error' and making the assimilation run converge toward the control run. Without data insertion, errors in the nonobserved wind component are large.

Liou, Yu-Chieng↗
NASA NTRS1993

Multi Wave Link: An interactive database for co-ordination of multiwavelength programs

Multi Wave Link is an operational database used by the international astronomical community to plan and coordinate multiwavelength observational campaigns on various celestial objects. This facility has been and continues to be an essential tool for the organization of numerous simultaneous ground-based and satellite observations. The Multi Wave Link system has both an interactive component and an electronic mail alert notice component. The interactive portion contains a database of schedules from several of the world's larger observatories and orbiting satellites. Lists of observers, observatories, and their available instrumentation are additionally contained in the system. Another feature added to the software allows the user to search the Royal Greenwich Observatory e-mail listing of international astronomers. This capability is often used by astronomers to find the electronic address of colleagues for improved communication. However, the most productive aspect of the Multi Wave Link system appears to be one of the simplest parts of the system: the electronic mail alert notices for observing campaigns. Information useful in planning the structure of future astronomical databases can be obtained by studying patterns of use of the Multi Wave Link system. This paper analyzes 22 months of Multi Wave Link usage from Sep. 1991 - Jul. 1993.

Cordova, France↗
NASA NTRS1997

Wave Number Selection for Incompressible Parallel Jet Flows Periodic in Space

The temporal instability of a spatially periodic parallel flow of an incompressible inviscid fluid for various jet velocity profiles is studied numerically using Floquet Analysis. The transition matrix at the end of a period is evaluated by direct numerical integration. For verification, a method based on approximating a continuous function by a series of step functions was used. Unstable solutions were found only over a limited range of wave numbers and have a band type structure. The results obtained are analogous to the behavior observed in systems exhibiting complexity at the edge of order and chaos.

Miles, Jeffrey Hilton↗
NASA NTRS1961

Calculation of Flow Fields from Bow-Wave Profiles for the Downstream Region of Blunt-Nosed Circular Cylinders in Axial Hypersonic Flight

A method by which known bow-wave profiles may be analyzed to give the flow fields around blunt-nosed cylinders in axial hypersonic flow is presented. In the method, the assumption is made that the pressure distribution curve in a transverse plane is similar to that given by blast- wave theory. Numerical analysis based on the one-dimensional energy and continuity equations then leads to distributions of all the flow variables in the cross section, for either a perfect gas or a real gas. The entire flow field need not be solved. Attention can be confined to any desired station. The critical question is the validity of the above assumption. It is tested for the case of a hemisphere cylinder in flight at 20,000 ft/sec. The flow is analyzed for three stations along the cylindrical afterbody, and found to compare very closely with the results of an exact (inviscid) solution. The assumed form of the pressure distribution occurs at stations as close as 1.2 diameters to the body nose. However, it is suggested that the assumption may not apply this far forward in general, particularly when bodies of nonsmooth contour are considered.

Seiff, Alvin↗
NASA NTRS1999

In Situ Monitoring of Crystal Growth Using MEPHISTO, Mission STS 87-Program USMP-4: Experimental Analysis

This report summarizes the results of the In situ Monitoring of Crystal Growth Using MEPHISTO (Material por l'Etude des Phenomenes Interessant de la Solidification sur Terre et en Orbite) experiment on USMP-4. The report includes microstructural and compositional data obtained during the first year of the post flight analysis, as well as numerical simulation of the flight experiment. Additional analyses are being continued and will be reported in the near future. The experiments utilized MEPHISTO hardware to study the solidification and melting behavior of bismuth alloyed with 1 at% tin. The experiments involved repeated melting and solidification of three samples, each approximately 90 cm long and 6mm in diameter. Half of each sample also included a 2 mm. diameter growth capillary, to assist in the formation of single grain inside. One sample provided the Seebeck voltage generated during melting and freezing processes. Another one provided temperature data and Peltier pulsed demarcation of the interface shape for post flight analysis. The third sample provided resistance and velocity measurements, as well as additional thermal data. The third sample was also quenched at the end of the mission to preserve the interface composition for post flight determination. A total of more than 45cm of directionally solidified alloy were directionally solidified at the end of the flight for post mission structural and compositional characterization. Metallurgical analysis of the samples has shown that the interfacial kinetics play a key role in controlling the morphological stability of faceted alloys. Substantial differences were observed in the Seebeck signal between the ground-based experiments and the space-based experiments. The temperature gradient in the liquid for the ground-based experiments was also significantly lower than the temperature gradient in the liquid for the space-based experiments. Both of these observations indicate significant influence of liquid convection for the ground-based experiments.

Abbaschian, Reza↗
NASA NTRS1999

Bubble Dynamics on a Heated Surface

In this work, we study steady and oscillatory thermocapillary and natural convective flows generated by a bubble on a heated solid surface. The interaction between gas and vapor bubbles with the surrounding fluid is of interest for both space and ground-based processing. A combined numerical-experimental approach is adopted here. The temperature field is visualized using Mach-Zehnder and/or Wollaston Prism Interferometry and the flow field is observed by a laser sheet flow visualization technique. A finite element numerical model is developed which solves the transient two-dimensional continuity, momentum, and energy equations and includes the effects of temperature-dependent surface tension and bubble surface deformation. Below the critical Marangoni number, the steady state low-g and 1-g temperature and velocity fields predicted by the finite element model are in excellent agreement with both the visualization experiments in our laboratory and recently published experimental results in the literature. Above the critical Marangoni number, the model predicts an oscillatory flow which is also closely confirmed by experiments. It is shown that the dynamics of the oscillatory flow are directly controlled by the thermal and hydrodynamic interactions brought about by combined natural and thermocapillary convection. Therefore, as numerical simulations show, there are considerable differences between the 1-g and low-g temperature and flow fields at both low and high Marangoni numbers. This has serious implications for both materials processing and fluid management in space.

Kassemi, M.↗
NASA NTRS2008

Ocean Carbon and Biogeochemistry Scoping Workshop on Terrestrial and Coastal Carbon Fluxes in the Gulf of Mexico, St. Petersburg, FL

Despite their relatively small surface area, ocean margins may have a significant impact on global biogeochemical cycles and, potentially, the global air-sea fluxes of carbon dioxide. Margins are characterized by intense geochemical and biological processing of carbon and other elements and exchange large amounts of matter and energy with the open ocean. The area-specific rates of productivity, biogeochemical cycling, and organic/inorganic matter sequestration are high in coastal margins, with as much as half of the global integrated new production occurring over the continental shelves and slopes (Walsh, 1991; Doney and Hood, 2002; Jahnke, in press). However, the current lack of knowledge and understanding of biogeochemical processes occurring at the ocean margins has left them largely ignored in most of the previous global assessments of the oceanic carbon cycle (Doney and Hood, 2002). A major source of North American and global uncertainty is the Gulf of Mexico, a large semi-enclosed subtropical basin bordered by the United States, Mexico, and Cuba. Like many of the marginal oceans worldwide, the Gulf of Mexico remains largely unsampled and poorly characterized in terms of its air-sea exchange of carbon dioxide and other carbon fluxes. The goal of the workshop was to bring together researchers from multiple disciplines studying terrestrial, aquatic, and marine ecosystems to discuss the state of knowledge in carbon fluxes in the Gulf of Mexico, data gaps, and overarching questions in the Gulf of Mexico system. The discussions at the workshop were intended to stimulate integrated studies of marine and terrestrial biogeochemical cycles and associated ecosystems that will help to establish the role of the Gulf of Mexico in the carbon cycle and how it might evolve in the face of environmental change.

Air sea exchanges↗
NASA NTRS2011

Impacts of Aerosol Direct Effects on the South Asian Climate: Assessment of Radiative Feedback Processes Using Model Simulations and Satellite/Surface Measurements

Current assessment of aerosol radiative effect is hindered by our incomplete knowledge of aerosol optical properties, especially absorption, and our current inability to quantify physical and microphysical processes. In this research, we investigate direct aerosol radiative effect over heavy aerosol loading areas (e.g., Indo-Gangetic Plains, South/East Asia) and its feedbacks on the South Asian climate during the pre-monsoon season (March-June) using the Purdue Regional Climate Model (PRCM) with prescribed aerosol data derived by the NASA Goddard Earth Observing System Model (GEOS-5). Our modeling domain covers South and East Asia (60-140E and 0-50N) with spatial resolutions of 45 km in horizontal and 28 layers in vertical. The model is integrated from 15 February to 30 June 2008 continuously without nudging (i.e., only forced by initial/boundary conditions). Two numerical experiments are conducted with and without the aerosol-radiation effects. Both simulations are successful in reproducing the synoptic patterns on seasonal-to-interannual time scales and capturing a pre-monsoon feature of the northward rainfall propagation over Indian region in early June which shown in Tropical Rainfall Measuring Mission (TRMM) observation. Preliminary result suggests aerosol-radiation interactions mainly alter surface-atmosphere energetics and further result in an adjustment of the vertical temperature distribution in lower atmosphere (below 700 hPa). The modifications of temperature and associated rainfall and circulation feedbacks on the regional climate will be discussed in the presentation.

Wang, Sheng-Hsiang↗
NASA NTRS1995

The method of space-time conservation element and solution element-applications to one-dimensional and two-dimensional time-marching flow problems

A nontraditional numerical method for solving conservation laws is being developed. The new method is designed from a physicist's perspective, i.e., its development is based more on physics than numerics. Even though it uses only the simplest approximation techniques, a 2D time-marching Euler solver developed recently using the new method is capable of generating nearly perfect solutions for a 2D shock reflection problem used by Helen Yee and others. Moreover, a recent application of this solver to computational aeroacoustics (CAA) problems reveals that: (1) accuracy of its results is comparable to that of a 6th order compact difference scheme even though nominally the current solver is only of 2nd-order accuracy; (2) generally, the non-reflecting boundary condition can be implemented in a simple way without involving characteristic variables; and (3) most importantly, the current solver is capable of handling both continuous and discontinuous flows very well and thus provides a unique numerical tool for solving those flow problems where the interactions between sound waves and shocks are important, such as the noise field around a supersonic over- or under-expansion jet.

Chang, Sin-Chung↗
NASA NTRS1990

Multi-blocking strategies for the INS3D incompressible Navier-Stokes code

With the continuing development of bigger and faster supercomputers, computational fluid dynamics (CFD) has become a useful tool for real-world engineering design and analysis. However, the number of grid points necessary to resolve realistic flow fields numerically can easily exceed the memory capacity of available computers. In addition, geometric shapes of flow fields, such as those in the Space Shuttle Main Engine (SSME) power head, may be impossible to fill with continuous grids upon which to obtain numerical solutions to the equations of fluid motion. The solution to this dilemma is simply to decompose the computational domain into subblocks of manageable size. Computer codes that are single-block by construction can be modified to handle multiple blocks, but ad-hoc changes in the FORTRAN have to be made for each geometry treated. For engineering design and analysis, what is needed is generalization so that the blocking arrangement can be specified by the user. INS3D is a computer program for the solution of steady, incompressible flow problems. It is used frequently to solve engineering problems in the CFD Branch at Marshall Space Flight Center. INS3D uses an implicit solution algorithm and the concept of artificial compressibility to provide the necessary coupling between the pressure field and the velocity field. The development of generalized multi-block capability in INS3D is described.

Gatlin, Boyd↗
NASA NTRS1949

Response of a Rotating Propeller to Aerodynamic Excitation

The flexural vibration of a rotating propeller blade with clamped shank is analyzed with the object of presenting, in matrix form, equations for the elastic bending moments in forced vibration resulting from aerodynamic forces applied at a fixed multiple of rotational speed. Matrix equations are also derived which define the critical speeds end mode shapes for any excitation order and the relation between critical speed and blade angle. Reference is given to standard works on the numerical solution of matrix equations of the forms derived. The use of a segmented blade as an approximation to a continuous blade provides a simple means for obtaining the matrix solution from the integral equation of equilibrium, so that, in the numerical application of the method presented, the several matrix arrays of the basic physical characteristics of the propeller blade are of simple form, end their simplicity is preserved until, with the solution in sight, numerical manipulations well-known in matrix algebra yield the desired critical speeds and mode shapes frame which the vibration at any operating condition may be synthesized. A close correspondence between the familiar Stodola method and the matrix method is pointed out, indicating that any features of novelty are characteristic not of the analytical procedure but only of the abbreviation, condensation, and efficient organization of the numerical procedure made possible by the use of classical matrix theory.

Arnoldi, Walter E.↗
NASA NTRS1978

A composition and thickness model for lunar impact crater and basin deposits

A model was developed for determining the composition and thicknesses of continuous deposits of large lunar impact craters and basins. Results of a photogeologic study and topographic analysis of continuous deposits of the lunar crater Delisle, which support the model, show that numerous secondary craters and crater chains with concentric dunes on their uprange rims occur well within the mapped deposits. At any given radius from Delisle, the secondary craters are of equal freshness and not vastly different in size, but with increasing radius they become better defined. The upper surface of the continuous deposits reaches the level of the pre-existing mare surface well within the mapped boundary of the deposits.

Morrison, R. H.↗
NASA NTRS1998

Studying Turbulence Using Numerical Simulation Databases: Proceedings of the Summer Program - No. 7

The Seventh Summer Program of the Center for Turbulence Research took place in the four-week period, July 5 to July 31, 1998. This was the largest CTR Summer Program to date, involving thirty-six participants from the U. S. and nine other countries. Thirty-one Stanford and NASA-Ames staff members facilitated and contributed to most of the Summer projects. A new feature, and perhaps a preview of the future programs, was that many of the projects were executed on non-NASA computers. These included supercomputers located in Europe as well as those operated by the Departments of Defense and Energy in the United States. In addition, several simulation programs developed by the visiting participants at their home institutions were used. Another new feature was the prevalence of lap-top personal computers which were used by several participants to carry out some of the work that in the past were performed on desk-top workstations. We expect these trends to continue as computing power is enhanced and as more researchers (many of whom CTR alumni) use numerical simulations to study turbulent flows. CTR's main role continues to be in providing a forum for the study of turbulence for engineering analysis and in facilitating intellectual exchange among the leading researchers in the field. Once again the combustion group was the largest. Turbulent combustion has enjoyed remarkable progress in using simulations to address increasingly complex and practically more relevant questions. The combustion group's studies included such challenging topics as fuel evaporation, soot chemistry, and thermonuclear reactions. The latter study was one of three projects related to the Department of Energy's ASCI Program (www.llnl.gov/asci); the other two (rocket propulsion and fire safety) were carried out in the turbulence modeling group. The flow control and acoustics group demonstrated a successful application of the so-called evolution algorithms which actually led to a previously unknown forcing strategy for jets yielding increased spreading rate. A very efficient algorithm for flow in complex geometries with moving boundaries based on the immersed boundary forcing technique was tested with very encouraging results. Also a new strategy for the destruction of aircraft trailing vortices was introduced and tested. The Reynolds Averaged Modeling (RANS) group demonstrated that the elliptic relaxation concept for RANS calculations is also applicable to transonic flows with shocks; however, prediction of laminar/turbulent transition remains an important pacing item. A large fraction of the LES effort was devoted to the development and testing of a new algorithmic procedure (as opposed to phenomenological model) for subgrid scale modeling based on regularized de-filtering of the flow variables. This appears to be a very promising approach, and a significant effort is currently underway to assess its robustness in high Reynolds number flows and in conjunction with numerical methods for complex flows. As part of the Summer Program two review tutorials were given on Turbulent structures in hydrocarbon pool fires (Sheldon Tieszen), and Turbulent combustion modeling: from RANS to LES via DNS (Luc Vervisch); and two seminars entitled Assessment of turbulence models for engineering applications (Paul Durbin) and Subgrid-scale modeling for non-premixed, turbulent reacting flows (James Riley) were presented. A number of colleagues from universities, government agencies, and industry attended the final presentations of the participants on July 31 and participated in the discussions. There are twenty-six papers in this volume grouped in five areas. Each group is preceded with an overview by its coordinator.

Source record↗
NASA NTRS1989

Parallel and vector computation for stochastic optimal control applications

A general method for parallel and vector numerical solutions of stochastic dynamic programming problems is described for optimal control of general nonlinear, continuous time, multibody dynamical systems, perturbed by Poisson as well as Gaussian random white noise. Possible applications include lumped flight dynamics models for uncertain environments, such as large scale and background random atmospheric fluctuations. The numerical formulation is highly suitable for a vector multiprocessor or vectorizing supercomputer, and results exhibit high processor efficiency and numerical stability. Advanced computing techniques, data structures, and hardware help alleviate Bellman's curse of dimensionality in dynamic programming computations.

Hanson, F. B.↗
NASA NTRS1994

Prediction of Undsteady Flows in Turbomachinery Using the Linearized Euler Equations on Deforming Grids

A linearized Euler solver for calculating unsteady flows in turbomachinery blade rows due to both incident gusts and blade motion is presented. The model accounts for blade loading, blade geometry, shock motion, and wake motion. Assuming that the unsteadiness in the flow is small relative to the nonlinear mean solution, the unsteady Euler equations can be linearized about the mean flow. This yields a set of linear variable coefficient equations that describe the small amplitude harmonic motion of the fluid. These linear equations are then discretized on a computational grid and solved using standard numerical techniques. For transonic flows, however, one must use a linear discretization which is a conservative linearization of the non-linear discretized Euler equations to ensure that shock impulse loads are accurately captured. Other important features of this analysis include a continuously deforming grid which eliminates extrapolation errors and hence, increases accuracy, and a new numerically exact, nonreflecting far-field boundary condition treatment based on an eigenanalysis of the discretized equations. Computational results are presented which demonstrate the computational accuracy and efficiency of the method and demonstrate the effectiveness of the deforming grid, far-field nonreflecting boundary conditions, and shock capturing techniques. A comparison of the present unsteady flow predictions to other numerical, semi-analytical, and experimental methods shows excellent agreement. In addition, the linearized Euler method presented requires one or two orders-of-magnitude less computational time than traditional time marching techniques making the present method a viable design tool for aeroelastic analyses.

Clark, William S.↗
NASA NTRS1993

Numerical prediction of transition of the F-16 wing at supersonic speeds

This work is part of the high speed research program currently underway at NASA. This project has the goal of gaining understanding of the technical requirements for supersonic-hypersonic flight. Specifically, this research is part of a continuing project to study the laminar flow over swept wings at high speeds and involves the numerical prediction of the flow about the F-16XL wing. The research uses the CNS/ARC3D codes and the resulting crossflow velocity components in order to estimate transition locations on the wing. Effects of angle of attack on the extent of laminar flows was found to be minimal. This result can be attributed to the fact that a laminar flow airfoil was used in this study, which has a continuous favorable pressure gradient over approximately the first 20 percent of the chord for angles of attacks up to 10 degrees. It should also be noted that even after 20 percent chord the pressure gradient either slowly continued to increase, but never decreased before 90 percent chord, except for the higher swept cases when separation occurs. Angles of attack greater than 10 degrees were not considered since this study assumes natural laminar flow for normal supersonic cruise flight conditions.

Cummings, Russell M.↗
NASA NTRS1987

Application of a flux-split algorithm to chemically relaxing, hypervelocity blunt-body flows

Viscous, nonequilibrium, hypervelocity flow fields over two axisymmetric configurations are numerically simulated using a factored, implicit, flux-split algorithm. The governing gas-dynamic and species-continuity equations for laminar flow are presented. The gas-dynamics/nonequilibrium-chemistry coupling procedure is developed as part of the solution procedure and is described in detail. Numerical solutions are presented for hypervelocity flows over a hemisphere and over an axisymmetric aeroassisted orbital transfer vehicle using three different chemistry models. The gas models considered are those for an ideal gas, for a frozen gas, and for chemically relaxing air consisting of five species. The calculated results are compared with existing numerical solutions in the literature along the stagnation line of the hemisphere. The effects of free-stream Reynolds number on the nonequilibrium flow field are discussed.

Balakrishnan, A.↗