Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “solve time”

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 199 records · Page 11

An upwind differencing scheme for the time-accurate incompressible Navier-Stokes equations

The two-dimensional incompressible Navier-Stokes equations are solved in a time-accurate manner in using the method of pseudocompressibility. Using this method, subiterations in pseudotime are required to satisfy the continuity equation at each time step. An upwind differencing scheme based on flux-difference splitting is used to compute the convective terms. The upwind differencing is biased based on the sign of the local eigenvalue of the Jacobian matrix. Third-order or fifth-order spatial accuracy is maintained throughout the interior grid points. The equations are solved using an implicit line-relaxation scheme. This solution scheme is stable and is capable of running at large time steps in pseudotime, leading to fast convergence for each physical time step. A variety of computed results are presented to validate the present scheme. Results for the flow over an oscillating plate are compared with the exact analytic solution, good agreement is seen. Excellent comparison is obtained between the computed solution and the analytical results for inviscid channel flow with an oscillating back pressure. Flow solutions over a circular cylinder with vortex shedding are also presented. Finally, the flow past an airfoil at -90 deg angle-of-attack is also computed.

Rogers, Stuart E.↗

Boson Star from Repulsive Light Scalars and Gravitational Waves

We study properties of boson stars consisting of ultra-light scalar dark matter with repulsive self-interactions. We investigate the origin of the maximum mass of spherically symmetric stable stars which emerges only when solving the full equations of motion in curved space-time, but not when solving the approximated Schrödinger-Newton equations. When the repulsion is weak, the backreaction of the curvature on the scalars acts as an additional source of attraction and can overcome the repulsion, resulting in a maximum star mass and compactness. We also point out that the potential in a UV completed particle physics model of light scalar dark matter is generally more complicated than the widely used φ 4 interaction. Additional interactions beyond φ 4 in the potential can dramatically change the properties of boson stars as well as modify the prospect of LIGO gravitational wave detection for binary mergers of boson stars.

Djuna Croon↗

Real-time operational planning for the U.S. air traffic system

This paper describes an integrated planning model for the U.S. air traffic system. The approach incorporates the dual objectives of monitoring collision risk while minimizing transportation costs. Specialized solution algorithms exploit the underlying structure of the model - especially for large-scale examples. The proposed formulation is tested with real-world data for the Indianapolis control sector. Additional experiments with a CRAY X-MP/24 supercomputer show that a full-scale model can be solved under real time conditions. Despite these advances, additional work is required in developing a practical system. Suggestions are made for combining advances in computer graphics and mathematical modeling.

Mulvey, John M.↗

Modeling of the WSTF frictional heating apparatus in high pressure systems

In order to develop a computer program able to model the frictional heating of metals in high pressure oxygen or nitrogen a number of additions have been made to the frictional heating model originally developed for tests in low pressure helium. These additions include: (1) a physical property package for the gases to account for departures from the ideal gas state; (2) two methods for spatial discretization (finite differences with quadratic interpolation or orthogonal collocation on finite elements) which substantially reduce the computer time required to solve the transient heat balance; (3) more efficient programs for the integration of the ordinary differential equations resulting from the discretization of the partial differential equations; and (4) two methods for determining the best-fit parameters via minimization of the mean square error (either a direct search multivariable simplex method or a modified Levenburg-Marquardt algorithm). The resulting computer program has been shown to be accurate, efficient and robust for determining the heat flux or friction coefficient vs. time at the interface of the stationary and rotating samples.

Skowlund, Christopher T.↗

Bubble Combustion

A method of energy production that is capable of low pollutant emissions is fundamental to one of the four pillars of NASA s Aeronautics Blueprint: Revolutionary Vehicles. Bubble combustion, a new engine technology currently being developed at Glenn Research Center promises to provide low emissions combustion in support of NASA s vision under the Emissions Element because it generates power, while minimizing the production of carbon dioxide (CO2) and nitrous oxides (NOx), both known to be Greenhouse gases. and allows the use of alternative fuels such as corn oil, low-grade fuels, and even used motor oil. Bubble combustion is analogous to the inverse of spray combustion: the difference between bubble and spray combustion is that spray combustion is spraying a liquid in to a gas to form droplets, whereas bubble combustion involves injecting a gas into a liquid to form gaseous bubbles. In bubble combustion, the process for the ignition of the bubbles takes place on a time scale of less than a nanosecond and begins with acoustic waves perturbing each bubble. This perturbation causes the local pressure to drop below the vapor pressure of the liquid thus producing cavitation in which the bubble diameter grows, and upon reversal of the oscillating pressure field, the bubble then collapses rapidly with the aid of the high surface tension forces acting on the wall of the bubble. The rapid and violent collapse causes the temperatures inside the bubbles to soar as a result of adiabatic heating. As the temperatures rise, the gaseous contents of the bubble ignite with the bubble itself serving as its own combustion chamber. After ignition, this is the time in the bubble s life cycle where power is generated, and CO2, and NOx among other species, are produced. However, the pollutants CO2 and NOx are absorbed into the surrounding liquid. The importance of bubble combustion is that it generates power using a simple and compact device. We conducted a parametric study using CAVCHEM, a computational model developed at Glenn, that simulates the cavitational collapse of a single bubble in a liquid (water) and the subsequent combustion of the gaseous contents inside the bubble. The model solves the time-dependent, compressible Navier-Stokes equations in one-dimension with finite-rate chemical kinetics using the CHEMKIN package. Specifically, parameters such as frequency, pressure, bubble radius, and the equivalence ratio were varied while examining their effect on the maximum temperature, radius, and chemical species. These studies indicate that the radius of the bubble is perhaps the most critical parameter governing bubble combustion dynamics and its efficiency. Based on the results of the parametric studies, we plan on conducting experiments to study the effect of ultrasonic perturbations on the bubble generation process with respect to the bubble radius and size distribution.

Corrigan, Jackie↗

A General-Purpose Optimization Engine for Multi-Disciplinary Design Applications

A general purpose optimization tool for multidisciplinary applications, which in the literature is known as COMETBOARDS, is being developed at NASA Lewis Research Center. The modular organization of COMETBOARDS includes several analyzers and state-of-the-art optimization algorithms along with their cascading strategy. The code structure allows quick integration of new analyzers and optimizers. The COMETBOARDS code reads input information from a number of data files, formulates a design as a set of multidisciplinary nonlinear programming problems, and then solves the resulting problems. COMETBOARDS can be used to solve a large problem which can be defined through multiple disciplines, each of which can be further broken down into several subproblems. Alternatively, a small portion of a large problem can be optimized in an effort to improve an existing system. Some of the other unique features of COMETBOARDS include design variable formulation, constraint formulation, subproblem coupling strategy, global scaling technique, analysis approximation, use of either sequential or parallel computational modes, and so forth. The special features and unique strengths of COMETBOARDS assist convergence and reduce the amount of CPU time used to solve the difficult optimization problems of aerospace industries. COMETBOARDS has been successfully used to solve a number of problems, including structural design of space station components, design of nozzle components of an air-breathing engine, configuration design of subsonic and supersonic aircraft, mixed flow turbofan engines, wave rotor topped engines, and so forth. This paper introduces the COMETBOARDS design tool and its versatility, which is illustrated by citing examples from structures, aircraft design, and air-breathing propulsion engine design.

Patnaik, Surya N.↗

Solving Nonlinear Euler Equations with Arbitrary Accuracy

A computer program that efficiently solves the time-dependent, nonlinear Euler equations in two dimensions to an arbitrarily high order of accuracy has been developed. The program implements a modified form of a prior arbitrary- accuracy simulation algorithm that is a member of the class of algorithms known in the art as modified expansion solution approximation (MESA) schemes. Whereas millions of lines of code were needed to implement the prior MESA algorithm, it is possible to implement the present MESA algorithm by use of one or a few pages of Fortran code, the exact amount depending on the specific application. The ability to solve the Euler equations to arbitrarily high accuracy is especially beneficial in simulations of aeroacoustic effects in settings in which fully nonlinear behavior is expected - for example, at stagnation points of fan blades, where linearizing assumptions break down. At these locations, it is necessary to solve the full nonlinear Euler equations, and inasmuch as the acoustical energy is of the order of 4 to 5 orders of magnitude below that of the mean flow, it is necessary to achieve an overall fractional error of less than 10-6 in order to faithfully simulate entropy, vortical, and acoustical waves.

Dyson, Rodger W.↗

Analysis of turbulent free jet hydrogen-air diffusion flames with finite chemical reaction rates

The nonequilibrium flow field resulting from the turbulent mixing and combustion of a supersonic axisymmetric hydrogen jet in a supersonic parallel coflowing air stream is analyzed. Effective turbulent transport properties are determined using the (K-epsilon) model. The finite-rate chemistry model considers eight reactions between six chemical species, H, O, H2O, OH, O2, and H2. The governing set of nonlinear partial differential equations is solved by an implicit finite-difference procedure. Radial distributions are obtained at two downstream locations of variables such as turbulent kinetic energy, turbulent dissipation rate, turbulent scale length, and viscosity. The results show that these variables attain peak values at the axis of symmetry. Computed distributions of velocity, temperature, and mass fraction are also given. A direct analytical approach to account for the effect of species concentration fluctuations on the mean production rate of species (the phenomenon of unmixedness) is also presented. However, the use of the method does not seem justified in view of the excessive computer time required to solve the resulting system of equations.

Sislian, J. P.↗

Computational methods to obtain time optimal jet engine control

Dynamic Programming and the Fletcher-Reeves Conjugate Gradient Method are two existing methods which can be applied to solve a general class of unconstrained fixed time, free right end optimal control problems. New techniques are developed to adapt these methods to solve a time optimal control problem with state variable and control constraints. Specifically, they are applied to compute a time optimal control for a jet engine control problem.

Basso, R. J.↗

Time-dependent aerodynamic analysis of ducted and unducted propfans at angle of attack

A three-dimensional unsteady aerodynamic analysis is described for predicting the time-dependent flow about ducted and unducted propfans operating at angle of attack. Although the freestream is assumed to be uniform, the flow relative to the rotating blades varies with circumferential position, resulting in an inherent unsteadiness due to the nonaxial inflow. The time-dependent Euler equations are solved utilizing a Runge-Kutta time-stepping scheme. The analysis is based on a finite-volume discretization employing a multiple-block grid network. To permit the use of large calculation time steps, an implicit residual smoothing scheme previously tested for unsteady flow calculations in two dimensions is extended to three spatial dimensions. For unducted propfans, a single H-type grid block is used for each blade passage to determine the time-periodic flowfield. For ducted propfans (ultra-high bypass fans) a body-centered C-type grid is wrapped about the cowl to improve the accuracy of the analysis in the high gradient flow region near the cowl leading edge. Numerical results are compared with available data for both ducted and unducted propfans operating at angle of attack.

Hall, Edward J.↗

The Thin Oil Film Equation

A thin film of oil on a surface responds primarily to the wall shear stress generated on that surface by a three-dimensional flow. The oil film is also subject to wall pressure gradients, surface tension effects and gravity. The partial differential equation governing the oil film flow is shown to be related to Burgers' equation. Analytical and numerical methods for solving the thin oil film equation are presented. A direct numerical solver is developed where the wall shear stress variation on the surface is known and which solves for the oil film thickness spatial and time variation on the surface. An inverse numerical solver is also developed where the oil film thickness spatial variation over the surface at two discrete times is known and which solves for the wall shear stress variation over the test surface. A One-Time-Level inverse solver is also demonstrated. The inverse numerical solver provides a mathematically rigorous basis for an improved form of a wall shear stress instrument suitable for application to complex three-dimensional flows. To demonstrate the complexity of flows for which these oil film methods are now suitable, extensive examination is accomplished for these analytical and numerical methods as applied to a thin oil film in the vicinity of a three-dimensional saddle of separation.

Brown, James L.↗

Development of a fractional-step method for the unsteady incompressible Navier-Stokes equations in generalized coordinate systems

A fractional step method is developed for solving the time-dependent three-dimensional incompressible Navier-Stokes equations in generalized coordinate systems. The primitive variable formulation uses the pressure, defined at the center of the computational cell, and the volume fluxes across the faces of the cells as the dependent variables, instead of the Cartesian components of the velocity. This choice is equivalent to using the contravariant velocity components in a staggered grid multiplied by the volume of the computational cell. The governing equations are discretized by finite volumes using a staggered mesh system. The solution of the continuity equation is decoupled from the momentum equations by a fractional step method which enforces mass conservation by solving a Poisson equation. This procedure, combined with the consistent approximations of the geometric quantities, is done to satisfy the discretized mass conservation equation to machine accuracy, as well as to gain the favorable convergence properties of the Poisson solver. The momentum equations are solved by an approximate factorization method, and a novel ZEBRA scheme with four-color ordering is devised for the efficient solution of the Poisson equation. Several two- and three-dimensional laminar test cases are computed and compared with other numerical and experimental results to validate the solution method. Good agreement is obtained in all cases.

Rosenfeld, Moshe↗

Viscosity Relaxation in Molten HgZnTe

Rotating cup measurements of the viscosity of the pseudo-binary melt, HgZnTe have shown that the isothermal liquid with zinc mole fraction 0.16 requires tens of hours of equilibration time before a steady viscous state can be achieved. Over this relaxation period, the viscosity at 790 C increases by a factor of two, while the viscosity at 810 C increases by 40%. Noting that the Group VI elements tend to polymerize when molten, we suggest that the viscosity of the melt is enhanced by the slow formation of Te atom chains. To explain the build-up of linear Te n-mers, we propose a scheme, which contains formation reactions with second order kinetics that increase the molecular weight, and decomposition reactions with first order kinetics that inactivate the chains. The resulting rate equations can be solved for the time dependence of each molecular weight fraction. Using these molecular weight fractions, we calculate the time dependence of the average molecular weight. Using the standard semi-empirical relation between polymer average molecular weight and viscosity, we then calculate the viscosity relaxation curve. By curve fitting, we find that the data imply that the rate constant for n-mer formation is much smaller than the rate constant for n-mer deactivation, suggesting that Te atoms only weakly polymerize in molten HgZnTe. The steady state toward which the melt relaxes occurs as the rate of formation of an n-mer becomes exactly balanced by the sum of the rate for its deactivation and the rate for its polymerization to form an (n+1)-mer.

Su, Ching-Hua↗

Calculation of fully three-dimensional separated flow with an unsteady viscous-inviscid interaction method

Three previous papers have shown that the viscous inviscid numerical methods were capable of calculating separated turbulent flows. The 'Semi-inverse method' and the models of the first author allow the computation of massive separations, stalled flows, and shock wave-boundary layer interactions, in two-dimensional or quasi-three-dimensional conditions, from low speeds to supersonic speeds. The 'Semi-implicit' method for time-consistent coupling allows for the computation of time-accurate transonic separated flow and buffer computations in two-dimensions. The present paper shows that the viscous-inviscid interaction approach is also able to compute the fully three-dimensional flow separation phenomena. The method is based on a thin-layer approximation of the theory of 'Defect-Formulation' that provides the viscous-inviscid splitting of the Navier-Stokes equations. A parametric analytical modelling of the 3D-turbulent velocity profiles is involved. Numerically, the 3D-velocity profiles are discretized in the normal z-direction and driven by parametric integral 3D-equations in direct or inverse modes in the x-direction. The viscous-inviscid coupling is fully 3D and solved the time-consistency problem with an extension of the 'Semi-implicit' method previously suggested in two-dimensions. A 3D inviscid subroutine with TSP approximation is used. Results are obtained for transonic steady flows over wings with shock-induced transonic separation. The method provides results for highly three-dimensional flow separations, such as induced by a three-dimensional through at the wall of a flat plate. The 3D viscous-inviscid coupling and the 3D model of the velocity field provide three-dimensional instantaneous skin-friction lines whose patterns exhibit the same complex topology as Navier-Stokes solvers, with foci, nodes, and saddle-points.

Leballeur, J. C.↗

Optimal guidance law development for an advanced launch system

The objective of this research effort was to develop a real-time guidance approach for launch vehicles ascent to orbit injection. Various analytical approaches combined with a variety of model order and model complexity reduction have been investigated. Singular perturbation methods were first attempted and found to be unsatisfactory. The second approach based on regular perturbation analysis was subsequently investigated. It also fails because the aerodynamic effects (ignored in the zero order solution) are too large to be treated as perturbations. Therefore, the study demonstrates that perturbation methods alone (both regular and singular perturbations) are inadequate for use in developing a guidance algorithm for the atmospheric flight phase of a launch vehicle. During a second phase of the research effort, a hybrid analytic/numerical approach was developed and evaluated. The approach combines the numerical methods of collocation and the analytical method of regular perturbations. The concept of choosing intelligent interpolating functions is also introduced. Regular perturbation analysis allows the use of a crude representation for the collocation solution, and intelligent interpolating functions further reduce the number of elements without sacrificing the approximation accuracy. As a result, the combined method forms a powerful tool for solving real-time optimal control problems. Details of the approach are illustrated in a fourth order nonlinear example. The hybrid approach is then applied to the launch vehicle problem. The collocation solution is derived from a bilinear tangent steering law, and results in a guidance solution for the entire flight regime that includes both atmospheric and exoatmospheric flight phases.

Calise, Anthony J.↗

Unsteady transonic aerodynamics and aeroelastic calculations at low-supersonic freestreams

A computational procedure is presented to simulate transonic unsteady flows and corresponding aeroelasticity of wings at low-supersonic freestreams. The flow is modeled by using the transonic small-perturbation theory. The structural equations of motions are modeled using modal equations of motion directly coupled with aerodynamics. Supersonic freestreams are simulated by properly accounting for the boundary conditions based on pressure waves along the flow characteristics in streamwise planes. The flow equations are solved using the time-accurate, alternating-direction implicit finite-difference scheme. The coupled aeroelastic equations of motion are solved by an integration procedure based on the time-accurate, linear-acceleration method. The flow modeling is verified by comparing calculations with experiments for both steady and unsteady flows at supersonic freestreams. The unsteady computations are made for oscillating wings. Comparisons of computed results with experiments show good agreement. Aeroelastic responses are computed for a rectangular wing at Mach numbers ranging from subtransonic to upper-transonic (supersonic) freestreams. The extension of the transonic dip into the upper transonic regime is illustrated.

Guruswamy, Guru P.↗

Exact solutions to Brans-Dicke cosmologies in flat Friedmann universes.

The Brans-Dicke cosmological equations for flat Friedmann-type expanding universes are solved parametrically for time, density, expansion parameter, and scalar field. These results reduce to a previously obtained exact solution to the radiation cosmology. Although the scalar field may be undetectable at the present epoch, it is felt that, if it exists, it must play an important role as one approaches the initial singularity of the cosmology.

Morganstern, R. E.↗