Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iteration method”

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 343 records · Page 19

Collisional-radiative switching - A powerful technique for converging non-LTE calculations

A very simple technique has been developed to converge statistical equilibrium and model atmospheric calculations in extreme non-LTE conditions when the usual iterative methods fail to converge from an LTE starting model. The proposed technique is based on a smooth transition from a collision-dominated LTE situation to the desired non-LTE conditions in which radiation dominates, at least in the most important transitions. The proposed approach was used to successfully compute stellar models with He abundances of 0.20, 0.30, and 0.50; Teff = 30,000 K, and log g = 2.9.

Hummer, D. G.↗

Cosmic-ray-modified stellar winds. III - A numerical iterative approach

A numerical iterative method is used to determine the modification of a stellar wind flow with a termination shock by the galactic cosmic rays. A two-fluid model consisting of cosmic rays and thermal stellar wind gas is used in which the cosmic rays are coupled to the background flow via scattering with magnetohydrodynamic waves or irregularities. A polytropic model is used to describe the thermal stellar wind gas, and the cosmic-rays are modeled as a hot, low-density gas with negligible mass flux. The positive galactic cosmic-ray pressure gradient serves to brake the outflowing stellar wind gas, and the cosmic rays modify the location of the critical point of the wind, the location of the shock, the wind fluid velocity profile, and the thermal gas entropy constants on both sides of the shock. The transfer of energy to the cosmic rays results in an outward radial flux of cosmic-ray energy.

Ko, C. M.↗

A three dimensional calculation of elastic equilibrium for composite materials

A compact scheme is applied to three-dimensional elasticity problems for composite materials, involving simple geometries. The mathematical aspects of this approach are discussed, in particular the iteration method. A vector processor code implementing the compact scheme is presented, and several numerical experiments are summarized.

Lustman, Liviu R.↗

High-resolution spectroscopy of central stars of planetary nebulae - LSS 1362

High resolution spectra of the hot sdO star LSS 1362, a central star of a planetary nebula, are analyzed by means of model atmosphere techniques. Improved NLTE model atmospheres based on the newly developed 'accelerated lambda iteration' method are used. The very high effective temperature (100000 K) derived previously from the UV energy distribution is confirmed. The surface gravity (log g = 5.3) and a normal helium abundance are derived. A mass of 0.65 solar mass and a luminosity of 10 exp 3.9 solar luminosities result from a comparison with evolutionary tracks. The distance is 1.1 kpc and the nebula has a diameter of about 0.25 pc to 0.5 pc. Weak photospheric emission lines of He II, CIV, Nv, and O v are detected. These lines are slightly blue-shifted with respect to the absorption line spectrum which is regarded as evidence for an atmospheric velocity field. It is conjectured that the observed broadening of the metal lines (25 km/s) may be caused by the velocity field rather than by rotation.

Heber, U.↗

Numerical simulation of hypersonic flow over an aeroassist flight experiment vehicle

Solutions of wind-tunnel and entry-flight flow around the vehicle are obtained from the Navier-Stokes equations coupled with the chemical species continuity equations if needed. The time-iterative method employs several techniques: shock fitting, chemistry-split ADI and an algebraic grid in conformal spherical-polar space. Sensitivities of the results to numerical parameters and to frozen, equilibrium and finite rate reactions are investigated in the forebody computation. Quantitative results are obtained for the shock layer and the near wake for the entire vehicle corresponding to both ground test and flight conditions. Complex flow characteristics are analyzed on the basis of the complete flowfield over the aerobrake and simplified afterbodies. The method is stable and cost effective, and has yielded shock locations and wall pressure distributions which are in good agreement with wind-tunnel data.

Li, C. P.↗

Paradigm For Statistical Analysis Of Threshold Detection

Relative values assigned to double mistakes and changed iteratively. Method for statistical analysis of threshold detection saves experimental time by enabling use of same set of measurements with respect to two thresholds. Accounts for double mistakes, unresolvable by threshold-detection technique.

Diner, Daniel B.↗

Variable-Metric Algorithm For Constrained Optimization

Variable Metric Algorithm for Constrained Optimization (VMACO) is nonlinear computer program developed to calculate least value of function of n variables subject to general constraints, both equality and inequality. First set of constraints equality and remaining constraints inequalities. Program utilizes iterative method in seeking optimal solution. Written in ANSI Standard FORTRAN 77.

Frick, James D.↗

On the constrained Chebyshev approximation problem on ellipses

Constrained Chebyshev approximation problems of the type with minimum (p is an element of Pi(sub n):p(c)=1) and maximum (z is an element of E) with /p(z)/ are considered. Here Pi(sub n) denotes the set of all complex polynomials of degree at most n, E is any ellipse in the complex plane, and c is an element of C/E. Such approximation problems arise in the context of optimizing semi-iterative methods for the solution of large, sparse systems of linear equations Ax=b with complex non-Hermitian coefficient matrices A. The problem of obtaining optimal polynomial preconditioners for conjugate gradient type methods for Ax=b also leads to problems of this type. A new family of polynomials -- q(sub n)(z;c), n is an element of N, and c is an element of C/E -- are introduced as the polynomials which are optimal for a modified version of the Chebyshev approximation problem with Pi(sub n) replaced by a certain subfamily. Some simple properties of q(sub n) are also listed. A necessary and sufficient condition for q(sub n) to be the extremal polynomial for the approximation problem is then derived. Finally, it is shown that q(sub n) is indeed optimal for the problem for all fixed n whenever the distance between c and E is sufficiently large. Results of some numerical tests are presented.

Freund, Roland↗

Finite elements: Theory and application

Recent advances in FEM techniques and applications are discussed in reviews and reports presented at the ICASE/LaRC workshop held in Hampton, VA in July 1986. Topics addressed include FEM approaches for partial differential equations, mixed FEMs, singular FEMs, FEMs for hyperbolic systems, iterative methods for elliptic finite-element equations on general meshes, mathematical aspects of FEMS for incompressible viscous flows, and gradient weighted moving finite elements in two dimensions. Consideration is given to adaptive flux-corrected FEM transport techniques for CFD, mixed and singular finite elements and the field BEM, p and h-p versions of the FEM, transient analysis methods in computational dynamics, and FEMs for integrated flow/thermal/structural analysis.

Dwoyer, D. L.↗

A computational efficient modelling of laminar separation bubbles

In predicting the aerodynamic characteristics of airfoils operating at low Reynolds numbers, it is often important to account for the effects of laminar (transitional) separation bubbles. Previous approaches to the modelling of this viscous phenomenon range from fast but sometimes unreliable empirical correlations for the length of the bubble and the associated increase in momentum thickness, to more accurate but significantly slower displacement-thickness iteration methods employing inverse boundary-layer formulations in the separated regions. Since the penalty in computational time associated with the more general methods is unacceptable for airfoil design applications, use of an accurate yet computationally efficient model is highly desirable. To this end, a semi-empirical bubble model was developed and incorporated into the Eppler and Somers airfoil design and analysis program. The generality and the efficiency was achieved by successfully approximating the local viscous/inviscid interaction, the transition location, and the turbulent reattachment process within the framework of an integral boundary-layer method. Comparisons of the predicted aerodynamic characteristics with experimental measurements for several airfoils show excellent and consistent agreement for Reynolds numbers from 2,000,000 down to 100,000.

Dini, Paolo↗

Parallelization of implicit finite difference schemes in computational fluid dynamics

Implicit finite difference schemes are often the preferred numerical schemes in computational fluid dynamics, requiring less stringent stability bounds than the explicit schemes. Each iteration in an implicit scheme involves global data dependencies in the form of second and higher order recurrences. Efficient parallel implementations of such iterative methods are considerably more difficult and non-intuitive. The parallelization of the implicit schemes that are used for solving the Euler and the thin layer Navier-Stokes equations and that require inversions of large linear systems in the form of block tri-diagonal and/or block penta-diagonal matrices is discussed. Three-dimensional cases are emphasized and schemes that minimize the total execution time are presented. Partitioning and scheduling schemes for alleviating the effects of the global data dependencies are described. An analysis of the communication and the computation aspects of these methods is presented. The effect of the boundary conditions on the parallel schemes is also discussed.

Decker, Naomi H.↗

Investigation of ultra wideband multi-channel dichroic beamsplitters from 0.3 to 52 microns

The development of a set of multi-channel dichroics which includes a 6 channel dichroic operating over the wavelength region from 0.3 to 52 microns is described. In order to achieve the optimum performance, the optical constraints of PbTe, Ge, and CdTe coatings in the strongly absorptive region have been determined by use of a new iterative method using normal incidence reflectance measurement of the multilayer together with initial values of energy gap and infinite refractive index for the semiconductor model. The design and manufacture of the dichroics is discussed and the final results are presented.

Zhang, K. Q.↗

On the joint inversion of geophysical data for models of the coupled core-mantle system

Joint inversion of magnetic, earth rotation, geoid, and seismic data for a unified model of the coupled core-mantle system is proposed and shown to be possible. A sample objective function is offered and simplified by targeting results from independent inversions and summary travel time residuals instead of original observations. These data are parameterized in terms of a very simple, closed model of the topographically coupled core-mantle system. Minimization of the simplified objective function leads to a nonlinear inverse problem; an iterative method for solution is presented. Parameterization and method are emphasized; numerical results are not presented.

Voorhies, Coerte V.↗

Combined AIE/EBE/GMRES approach to incompressible flows

Adaptive implicit-explicit (AIE), grouped element-by-element (GEBE), and generalized minimum residuals (GMRES) solution techniques for incompressible flows are combined. In this approach, the GEBE and GMRES iteration methods are employed to solve the equation systems resulting from the implicitly treated elements, and therefore no direct solution effort is involved. The benchmarking results demonstrate that this approach can substantially reduce the CPU time and memory requirements in large-scale flow problems. Although the description of the concepts and the numerical demonstration are based on the incompressible flows, the approach presented here is applicable to larger class of problems in computational mechanics.

Liou, J.↗

Study of shock-induced combustion using an implicit TVD scheme

The supersonic combustion flowfields associated with various hypersonic propulsion systems, such as the ram accelerator, the oblique detonation wave engine, and the scramjet, are being investigated using a new computational fluid dynamics (CFD) code. The code solves the fully coupled Reynolds-averaged Navier-Stokes equations and species continuity equations in an efficient manner. It employs an iterative method and a second order differencing scheme to improve computational efficiency. The code is currently being applied to study shock wave/boundary layer interactions in premixed combustible gases, and to investigate the ram accelerator concept. Results obtained for a ram accelerator configuration indicate a new combustion mechanism in which a shock wave induces combustion in the boundary layer, which then propagates outward and downstream. The combustion process creates a high pressure region over the back of the projectile resulting in a net positive thrust forward.

Yungster, Shayne↗

Acoustic modes in fluid networks

Pressure and flow rate eigenvalue problems for one-dimensional flow of a fluid in a network of pipes are derived from the familiar transmission line equations. These equations are linearized by assuming small velocity and pressure oscillations about mean flow conditions. It is shown that the flow rate eigenvalues are the same as the pressure eigenvalues and the relationship between line pressure modes and flow rate modes is established. A volume at the end of each branch is employed which allows any combination of boundary conditions, from open to closed, to be used. The Jacobi iterative method is used to compute undamped natural frequencies and associated pressure/flow modes. Several numerical examples are presented which include acoustic modes for the Helium Supply System of the Space Shuttle Orbiter Main Propulsion System. It should be noted that the method presented herein can be applied to any one-dimensional acoustic system involving an arbitrary number of branches.

Michalopoulos, C. D.↗

Stability investigations of airfoil flow by global analysis

As the result of global, non-parallel flow stability analysis the single value of the disturbance growth-rate and respective frequency is obtained. This complex value characterizes the stability of the whole flow configuration and is not referred to any particular flow pattern. The global analysis assures that all the flow elements (wake, boundary and shear layer) are taken into account. The physical phenomena connected with the wake instability are properly reproduced by the global analysis. This enhances the investigations of instability of any 2-D flows, including ones in which the boundary layer instability effects are known to be of dominating importance. Assuming fully 2-D disturbance form, the global linear stability problem is formulated. The system of partial differential equations is solved for the eigenvalues and eigenvectors. The equations, written in the pure stream function formulation, are discretized via FDM using a curvilinear coordinate system. The complex eigenvalues and corresponding eigenvectors are evaluated by an iterative method. The investigations performed for various Reynolds numbers emphasize that the wake instability develops into the Karman vortex street. This phenomenon is shown to be connected with the first mode obtained from the non-parallel flow stability analysis. The higher modes are reflecting different physical phenomena as for example Tollmien-Schlichting waves, originating in the boundary layer and having the tendency to emerge as instabilities for the growing Reynolds number. The investigations are carried out for a circular cylinder, oblong ellipsis and airfoil. It is shown that the onset of the wake instability, the waves in the boundary layer, the shear layer instability are different solutions of the same eigenvalue problem, formulated using the non-parallel theory. The analysis offers large potential possibilities as the generalization of methods used till now for the stability analysis.

Morzynski, Marek↗

Radiative equilibrium in rapidly expanding shells

The equation of radiative equilibrium for rapidly expanding spherical shells including the special relativistic terms is discussed, and a temperature-correction procedure for the solution of this equation is given. The temperature-correction scheme combines a partial linearization approach with the approximate Lambda-operator iteration method. For some test problems, the temperature structures of atmospheres in radiative equilibrium in the Lagrangian frame are calculated. The effects of the velocity field on the temperature stratification and the luminosity are also discussed. We find that the temperature structure of a moving photosphere in radiative equilibrium in the Lagrangian frame differs significantly from that of a static photosphere and that, at least, all first-order terms, in particular the advection and aberration terms, have to be included consistently both in the radiative transfer and the energy equation in studies of supernovae atmospheres and other high-speed radiating flows.

Hauschildt, P. H.↗