Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical techniques”

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 73 records · Page 4

Spectrum formation in supernovae - Numerical techniques

The study combines several novel techniques for spectrum simulation in the Eddington computer program which solves the comoving frame equation of transfer coupled with the statistical and radiative equilibrium equations. One of these is a generalization of the accelerated lambda iteration (ALI) scheme to include an approximate frequency-derivative operator. This greatly enhances the convergence rate of ALI in optically thick, high-velocity shear flows. Another is a partial linearization technique which is capable of efficiently solving a very large number of rate equations on a moderately sized computer. An expansion opacity and emissivity approximation is derived which makes it possible to determine the effect on the transfer and statistical equilibrium of a very large number of lines not explicitly represented in the frequency grid and additionally to treat line-blanketing from species not explicitly included in the rate equations. The utility of these techniques is illustrated with models of two supernovae.

Eastman, Ronald G.↗

Comparisons among second-order, fourth-order, and pseudospectral advection techniques in a mesoscale numerical model

The reported investigation is concerned with mesoscale squall-line simulation utilizing numerical techniques. Particular attention is given to fourth-order results and their implications for the dynamical evolution of mesoscale squall-line systems which contain severe local storms. An adiabatic inviscid set of prognostic equations is utilized in a z coordinate system for the fundamental experiment. The fourth-order advection is employed for all time-dependent equations. The Euler-backward time marching scheme is utilized with a 60 second time step. The horizontal mesh length is 42 km. Horizontal diffusion is accomplished by utilizing a smoother-desmoother. Lateral boundary conditions are designed to reduce the development of strong gradients in dependent variables near the boundaries by bringing interior values from the grid to the boundaries.

Kaplan, M.↗

Mean Element Propagations Using Numerical Averaging

The long-term evolution characteristics (and stability) of an orbit are best characterized using a mean element propagation of the perturbed two body variational equations of motion. The averaging process eliminates short period terms leaving only secular and long period effects. In this study, a non-traditional approach is taken that averages the variational equations using adaptive numerical techniques and then numerically integrating the resulting EOMs. Doing this avoids the Fourier series expansions and truncations required by the traditional analytic methods. The resultant numerical techniques can be easily adapted to propagations at most solar system bodies.

orbit propagations↗

Numerical derivative techniques for trajectory optimization

The adoption of robust numerical optimization techniques in trajectory simulation programs has resulted in powerful design and analysis tools. These trajectory simulation/optimization programs are widely used, and a representative list includes the GTS system, the POST program, and newer collocation methods such as OTIS and FONPAC. All of these programs rely on optimization algorithms which require objective function and constraint gradient data during the iteration process. However, most trajectory optimization problems lack simple analytical expressions for these derivatives. In the general case a function evaluation involves integrating aerodynamic, propulsive, and gravity forces over multiple trajectory phases with complex control models. With the newer collocation methods, the integration is replaced by defect constraints and cubic approximations for the state. While analytic gradient expressions can sometimes be derived for trajectory optimization problems, the derivation is cumbersome, time consuming, and prone to mistakes. Fortunately, an alternate method exists for the gradient evaluation, namely finite difference approximations. In this paper some finite difference gradient techniques developed for use with the GTS system are presented. These techniques include methods for computing first and second partial derivatives of single and multiple sets of functions. A key feature of these methods is an error control mechanism which automatically adjusts the perturbation size to obtain accurate derivative values.

Hallman, Wayne P.↗

A Comparison of Metamodeling Techniques via Numerical Experiments

This paper presents a comparative analysis of a few metamodeling techniques using numerical experiments for the single input-single output case. These experiments enable comparing the models' predictions with the phenomenon they are aiming to describe as more data is made available. These techniques include (i) prediction intervals associated with a least squares parameter estimate, (ii) Bayesian credible intervals, (iii) Gaussian process models, and (iv) interval predictor models. Aspects being compared are computational complexity, accuracy (i.e., the degree to which the resulting prediction conforms to the actual Data Generating Mechanism), reliability (i.e., the probability that new observations will fall inside the predicted interval), sensitivity to outliers, extrapolation properties, ease of use, and asymptotic behavior. The numerical experiments describe typical application scenarios that challenge the underlying assumptions supporting most metamodeling techniques.

Crespo, Luis G.↗

Interactive phenomena in supersonic jet mixing problems. I Phenomenology and numerical modeling techniques

The interactive phenomena that occur in supersonic jet mixing flowfields, and numerical modeling techniques developed to analyze such phenomena are discussed. A spatial marching procedure based on solving the parabolized Navier-Stokes jet mixing equations is presented. This procedure combines shock-capturing methodology for the analysis of supersonic mixing regions with pressure-split methodology for the analysis of subsonic mixing regions. The two regions are coupled at viscous sonic lines utilizing a viscous-characteristic coupling procedure. Specialized techniques for the treatment of jet boundary growth, strong discontinuties (Mach disks), and small embedded subsonic zones (behind Mach disks) are presented. Turbulent processes are represented by two-equation turbulence model formulations. In Part II of this article, numerical studies are presented for a variety of supersonic jet interactive phenomena.

Dash, S. M.↗

Sound diffraction around screens and wedges for arbitrary point source locations

A development is presented for the Green's function for a point source in the vicinity of a rigid wedge. The diffraction contributions to the Green's function for arbitrary source and listener location is expressed in a form which can be readily evaluated using the Laguerre technique for numerical integration. The present approach offers the advantages of efficient numerical evaluation and of relatively straightforward reduction to well-known analytical approximations in limiting cases. Comparisons with previously obtained experimental and numerical results obtained by Ambaud and Bergassoli (1972) are presented. The comparison with the experimental results is excellent; the advantages of the present numerical technique, vis a vis that of Ambaud and Bergassoli, are pointed out.

Hadden, W. J., Jr.↗

Triple flame structure and diffusion flame stabilization

The stabilization of diffusion flames is studied using asymptotic techniques and numerical tools. The configuration studied corresponds to parallel streams of cold oxidizer and fuel initially separated by a splitter plate. It is shown that stabilization of a diffusion flame may only occur in this situation by two processes. First, the flame may be stabilized behind the flame holder in the wake of the splitter plate. For this case, numerical simulations confirm scalings previously predicted by asymptotic analysis. Second, the flame may be lifted. In this case a triple flame is found at longer distances downstream of the flame holder. The structure and propagation speed of this flame are studied by using an actively controlled numerical technique in which the triple flame is tracked in its own reference frame. It is then possible to investigate the triple flame structure and velocity. It is shown, as suggested from asymptotic analysis, that heat release may induce displacement speeds of the triple flame larger than the laminar flame speed corresponding to the stoichiometric conditions prevailing in the mixture approaching the triple flame. In addition to studying the characteristics of triple flames in a uniform flow, their resistance to turbulence is investigated by subjecting triple flames to different vortical configurations.

Veynante, D.↗

Three-dimensional magnetotail equilibria by numerical relaxation techniques

Consideration is given to a numerical method to iteratively approach 3D magnetostatic force equilibria, with modeling emphasis on the development of a suitable model of the earth's magnetotail, including a portion of the inner magnetosphere, i.e., on models which violate the 'tail approximation' commonly employed in analytical models. The numerical approach is compared to methods developed for laboratory plasma physics. The method is applied to the magnetotail outside of 10 R(E), using Tsyganenko's (1987) model for the quiet magnetosphere as an initial condition. The changes of the magnetic field necessary to yield an equilibrium configuration and the resulting distribution of the self-consistently derived pressure are discussed. It is shown that a self-consistent magnetotail equilibrium based on a close approximation to Tsyganenko's requires a region 1 type current system, which is not present in the initial configuration.

Hesse, Michael↗

Hourglass control in staggered-grid hydrodynamics using virtual element stabilization techniques

Numerical simulations using the staggered-grid hydrodynamics (SGH) discretization suffer from hourglass instabilities. In this work, we develop a stabilization method to suppress the hourglass instabilities using techniques from the virtual element method (VEM). The stiffness matrix of the VEM consists of two terms: the consistency matrix which is rank deficient and the stability matrix. Here, we first show that in two dimensions and on general polygons, the stiffness matrix of the SGH is identical to the consistency matrix of the linear VEM for both the diffusion equation and the linear elasticity equation. These analyses explain the origin of the hourglass instabilities of the SGH discretization method, and establish a theoretical foundation for our proposed stabilization method by augmenting the stiffness matrix of the SGH discretization using the VEM stability matrix. Then, we present numerical examples using Lagrangian SGH simulations. The numerical experiments demonstrate that the proposed VEM stabilization method is effective at eliminating hourglass modes in the SGH discretization.

97 MATHEMATICS AND COMPUTING↗

Inlet Flow Distortion and Unsteady Blade Response in a Transonic Axial-Compressor Rotor

This paper describes the unsteady blade surface pressures on the first-stage rotor blades of a two-stage transonic axial flow compressor experiencing inlet flow distortion. This study was conducted to demonstrate the ability of a full annulus unsteady Reynolds-averaged Navier-Stokes numerical technique to predict unsteady pressures on the rotor blades operating in a distorted inflow. A total pressure distortion produced by a variable mesh screen mounted near the inlet was used to excite the unsteady blade loading on the rotor. On-blade pressure transducers were used to measure the unsteady blade surface pressure. These pressures and the resulting transient load on the rotor blades were compared to the numerical prediction. It is important to develop numerical techniques to predict these transient loads to better understand the response of compressor blades to forcing functions. With this enhanced understanding and ability to predict these transient forces, more robust compressors can be developed. In the study, a high definition of the inlet flow distortion was achieved by rotating the distortion screens. In this manner the inlet flow distortion and the distortion at the first stage stator leading edge were measured at approximately every 0.7 degrees. This full annulus high definition of the inlet flow distortion was used as the inlet boundary condition for the numerical technique. The experimental measurements and numerical analyses are highly complementary in this study. Detailed comparisons between the measurements and the numerical analyses indicate that the current numerical procedure calculates the unsteady aerodynamic pressure on the blade surfaces reasonably well. Further, the agreement of the measured and predicted rotor exit flow distortion at the first stage stator leading edge provides verification of the numerical technique.

Rabe, D. C.↗

The Analysis and Design of Low Boom Configurations Using CFD and Numerical Optimization Techniques

The use of computational fluid dynamics (CFD) for the analysis of sonic booms generated by aircraft has been shown to increase the accuracy and reliability of predictions. CFD takes into account important three-dimensional and nonlinear effects that are generally neglected by modified linear theory (MLT) methods. Up to the present time, CFD methods have been primarily used for analysis or prediction. Some investigators have used CFD to impact the design of low boom configurations using trial and error methods. One investigator developed a hybrid design method using a combination of Modified Linear Theory (e.g. F-functions) and CFD to provide equivalent area due to lift driven by a numerical optimizer to redesign or modify an existing configuration to achieve a shaped sonic boom signature. A three-dimensional design methodology has not yet been developed that completely uses nonlinear methods or CFD. Constrained numerical optimization techniques have existed for some time. Many of these methods use gradients to search for the minimum of a specified objective function subject to a variety of design variable bounds, linear and nonlinear constraints. Gradient based design optimization methods require the determination of the objective function gradients with respect to each of the design variables. These optimization methods are efficient and work well if the gradients can be obtained analytically. If analytical gradients are not available, the objective gradients or derivatives with respect to the design variables must be obtained numerically. To obtain numerical gradients, say, for 10 design variables, might require anywhere from 10 to 20 objective function evaluations. Typically, 5-10 global iterations of the optimizer are required to minimize the objective function. In terms of using CFD as a design optimization tool, the numerical evaluation of gradients can require anywhere from 100 to 200 CFD computations per design for only 10 design variables. If one CFD computation requires an hour of computational time on a Cray computer, one can see that the use of constrained numerical optimization quickly becomes impractical.Hence, in order to practically couple a numerical design optimization technique with a CFD method, the CFD method must be extremely efficient with running times on the order of only minutes. The CFD Euler code developed under NASA sponsorship and referred to as MIM3D-SB for the most part fulfills these efficiency requirements. Analysis of wing- body configurations can be computed in a matter of a few minutes. The present study will concentrate on the feasibility of the use of this CFD code in conjunction with a numerical design optimization technique for the sonic boom reduction of candidate HSCT configurations. A preliminary supersonic aircraft design system has been established that utilizes the numerical design optimization code NPSOL developed at Stanford University coupled with the supersonic NUM3D-SB CFD code. Many questions still need to be answered in regard to using CFD and numerical optimizers as design tools. There are difficulties related to both the CFD codes and the numerical optimizers. Numerical optimizers can converge to a local minima rather than a global minima. This behavior is largely a function of the initial guess in the design space. The optimizer also is searching for a minimum of the function in terms of its derivative without any regard to the actual function value. Numerically (i.e. CFD) determined gradients can also generate spurious numerical local minima. In addition, for the sonic boom problem, grid fineness will also determine the accuracy of the final design solution. Design optimization methods work well on problems defined by continuous objective functions. The sonic boom signature design problem is not necessarily defined by a continuous objective function. The signature can have a variety of shapes; i.e. from N-wave to multiple shocks. The far-field or ground signature may not transition continuously from one shape to another and hence, may exhibit discontinuous behavior. This is also a source of difficulty in using design optimization methods.In the following sections, several low boom and one reference aircraft configuration will be analyzed to predict their sonic boom signature characteristics. Modifications to some of these designs will also be presented to demonstrate the feasibility of using CFD as a design tool and to demonstrate the feasibility of designing shaped sonic boom signatures. Design modifications to some configurations will be presented to demonstrate the feasibility of achieving shaped signatures with reduced levels and not necessarily to represent realistic or aerodynamically efficient design modifications. Fuselage volume or camber are used as design variables in order to have a minimal effect on the primary wing aerodynamics. The paper will also seek to demonstrate whether a hybrid or ramped signature is feasible to achieve. For the low-boom configurations, the CFD predicted signatures will be compared qualitatively to their MLT design signatures.

Siclari, Michael J.↗

Liquid rocket boosters for the Space Shuttle - An outline of trade issues and analysis techniques

Numerous liquid rocket booster configurations were studied to replace the Solid Rocket Boosters on the current Space Shuttle. While improved safety, reliability, performance and economy of operation were sought, the designs also sought to minimize volumetric growth to reduce aerodynamic loads and structural modifications. Aero loads could be significant as a function of booster diameter if extrapolated directly from solid rocket characteristics. Procedures to meet these goals with analytical booster sizing and simple dynamic simulations are shown. Engine chamber pressure proved to be a good independent variable for implications of pressure-fed and pump-fed engines, fuel and material selection. For some configurations, ascent trajectory modifications were suggested: longer vertical rise, lofting, heads-up orientation, and higher first and second stage target altitudes. Since engine manufacture would dominate booster costs, booster simplification of reuse are discussed based on cost estimates.

Kelly, W. D.↗

Stability of numerical integration techniques for transient rotor dynamics

A finite element model of a rotor bearing system was analyzed to determine the stability limits of the forward, backward, and centered Euler; Runge-Kutta; Milne; and Adams numerical integration techniques. The analysis concludes that the highest frequency mode determines the maximum time step for a stable solution. Thus, the number of mass elements should be minimized. Increasing the damping can sometimes cause numerical instability. For a uniform shaft, with 10 mass elements, operating at approximately the first critical speed, the maximum time step for the Runge-Kutta, Milne, and Adams methods is that which corresponds to approximately 1 degree of shaft movement. This is independent of rotor dimensions.

Kascak, A. F.↗

Improved numerical methods for turbulent viscous recirculating flows

The objective of the present study is to improve both the accuracy and computational efficiency of existing numerical techniques used to predict viscous recirculating flows in combustors. A review of the status of the study is presented along with some illustrative results. The effort to improve the numerical techniques consists of the following technical tasks: (1) selection of numerical techniques to be evaluated; (2) two dimensional evaluation of selected techniques; and (3) three dimensional evaluation of technique(s) recommended in Task 2.

Vandoormaal, J. P.↗