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 541 records · Page 30

Full numerical simulation of coflowing, axisymmetric jet diffusion flames

The near field of a non-premixed flame in a low speed, coflowing axisymmetric jet is investigated numerically using full simulation. The time-dependent governing equations are solved by a second-order, explicit finite difference scheme and a single-step, finite rate model is used to represent the chemistry. Steady laminar flame results show the correct dependence of flame height on Peclet number and reaction zone thickness on Damkoehler number. Forced simulations reveal a large difference in the instantaneous structure of scalar dissipation fields between nonbuoyant and buoyant cases. In the former, the scalar dissipation marks intense reaction zones, supporting the flamelet concept; however, results suggest that flamelet modeling assumptions need to be reexamined. In the latter, this correspondence breaks down, suggesting that modifications to the flamelet modeling approach are needed in buoyant turbulent diffusion flames.

Mahalingam, S.↗

Calculation of unsteady linearized Euler flows in cascades using harmonically deforming grids

A method for calculating unsteady, inviscid, compressible flows in cascades is presented. Using the linearized Euler technique, the flow is decomposed into a steady or mean flow plus a harmonically varying small disturbance flow. The equations that describe the small disturbance flow are linear variable coefficient equations, and are solved using a pseudo-time time marching Lax-Wendroff technique. Unlike previous linearized methods, however, the solution is computed on a harmonically deforming computational grid that conforms to the motion of the vibrating airfoils. The mean flow and perturbation flow solutions are defined in the deforming coordinate system rather than in a coordinate system fixed in space. Hence, no extrapolation terms are required to implement the upwash boundary conditions at the airfoil surfaces, significantly improving the accuracy of the method. For transonic flow calculations, unsteady shock motions are modelled using shock capturing. The unsteady loads due to the shock motion are then seen as pressure impulses. Representative computational results are presented for transonic channel flows and subsonic and transonic cascade flows.

Hall, Kenneth C.↗

Adjoint-operators and non-adiabatic learning algorithms in neural networks

Adjoint sensitivity equations are presented, which can be solved simultaneously (i.e., forward in time) with the dynamics of a nonlinear neural network. These equations provide the foundations for a new methodology which enables the implementation of temporal learning algorithms in a highly efficient manner.

Toomarian, N.↗

The effects of self-generated and applied magnetic fields on the computation of flow over a Mars return aerobrake

A CFD technique is developed to calculate the electromagnetic phenomena simultaneously with the fluid flow in the shock layer over an axisymmetric blunt body in a thermal-equilibrium chemical-nonequilibrium environment. The flowfield is solved using an explicit time-marching, first-order spatially accurate scheme. The electromagnetic phenomena are coupled to the real-gas flow solver through an iterative procedure. The electromagnetic terms introduce a strong stiffness, which was overcome by using significantly smaller time steps for the electromagnetic conservation equation. The technique is applied in calculating the flow over a Mars return aerobrake vehicle entering the Earth's atmosphere. For the case where no external field is applied, the electromagnetic effects have little impact on the flowfield.

Palmer, Grant↗

Prediction of unsteady aerodynamic loads in cascades 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. Using the linearized Euler technique, one decomposes the flow into a mean (or steady) flow plus an unsteady, harmonically varying, small disturbance flow. Linear variable coefficient equations describe the small disturbance behavior of the flow, and are solved using a pseudo-time marching Lax-Wendroff scheme. For the blade motion problem, a harmonically deforming computational rid that conforms to the motion of vibrating blades eliminates large error producing mean flow gradient terms that would otherwise appear in the unsteady flow tangency boundary condition. The paper also presents a new, numerically exact, nonreflecting far-field boundary condition based on an eigenanalysis of the discretized equations. Computed flow solutions demonstrate the computational accuracy and efficiency of the present method. The solution of the linearized Euler equations requires one to two orders of magnitude less computer time than solution of the nonlinear Euler equations using traditional time-accurate time-marching techniques. In addition, the deformable grid significantly improves the accuracy of the solution.

Hall, Kenneth C.↗

Computation of incompressible viscous flows through artificial heart devices with moving boundaries

The extension of computational fluid dynamics techniques to artificial heart flow simulations is illustrated. Unsteady incompressible Navier-Stokes equations written in 3-D generalized curvilinear coordinates are solved iteratively at each physical time step until the incompressibility condition is satisfied. The solution method is based on the pseudo compressibility approach and uses an implicit upwind differencing scheme together with the Gauss-Seidel line relaxation method. The efficiency and robustness of the time accurate formulation of the algorithm are tested by computing the flow through model geometries. A channel flow with a moving indentation is computed and validated with experimental measurements and other numerical solutions. In order to handle the geometric complexity and the moving boundary problems, a zonal method and an overlapping grid embedding scheme are used, respectively. Steady state solutions for the flow through a tilting disk heart valve was compared against experimental measurements. Good agreement was obtained. The flow computation during the valve opening and closing is carried out to illustrate the moving boundary capability.

Kiris, Cetin↗

Simulations of surface winds at the Viking Lander sites using a one-level model

The one-level model developed by Mass and Dempsey for use in predicting surface flows in regions of complex terrain was adapted to simulate surface flows at the Viking lander sites on Mars. In the one-level model, prediction equations for surface winds and temperatures are formulated and solved. Surface temperatures change with time in response to diabatic heating, horizontal advection, adiabatic heating and cooling effects, and horizontal diffusion. Surface winds can change in response to horizontal advection, pressure gradient forces, Coriolis forces, surface drag, and horizontal diffusion. Surface pressures are determined by integration of the hydrostatic equation from the surface to some reference level. The model has successfully simulated surface flows under a variety of conditions in complex-terrain regions on Earth.

Bridger, Alison F. C.↗

Application of a monotonic upstream-biased transport scheme to three-dimensional constituent transport calculations

The application of van Leer's scheme, a monotonic, upstream-biased differencing scheme, to three-dimensional constituent transport calculations is shown. The major disadvantage of the scheme is shown to be a self-limiting diffusion. A major advantage of the scheme is shown to be its ability to maintain constituent correlations. The scheme is adapted for a spherical coordinate system with a hybrid sigma-pressure coordinate in the vertical. Special consideration is given to cross-polar flow. The vertical wind calculation is shown to be extremely sensitive to the method of calculating the divergence. This sensitivity implies that a vertical wind formulation consistent with the transport scheme is essential for accurate transport calculations. The computational savings of the time-splitting method used to solve this equation are shown. Finally, the capabilities of this scheme are illustrated by an ozone transport and chemistry model simulation.

Allen, Dale J.↗

Numerical studies of natural transition in a decelerating boundary layer

The laminar-turbulent transition in a decelerating flat-plate boundary layer is investigated by solving numerically the incompressible time-dependent three-dimensional Navier-Stokes equations. In the early stages of the transition, two-dimensional and slightly oblique waves are found to grow because of an inflexional instability. Their subsequent nonlinear interactions trigger the breakdown and determine the pattern of the Lambda vortices. The tips of the Lambda vortices are rarely aligned in the flow direction and appear in local patches, often near the dislocation in the initial wave pattern. The observed features are accounted for by a simple wave-interference model.

Yang, Kyung-Soo↗

Three-dimensional calculations of supersonic reacting flows using an LU scheme

An implicit finite volume lower-upper time-marching method which efficiently solves the complete Navier-Stokes and specied equations in a fully coupled fashion is the basis of the present 3D numerical program for simulating the supersonic reacting flows of H2 in air. The chemistry model incorporated has nine species and 18 reaction steps. Calculations are presented for flowfields of underexpanded hydrogen jets that are transversely injected into the supersonic airstream within scramjet combustors; the shock structure, separated flow regions around the injector, and combustion-product distributions are clearly represented.

Yu, Sheng-Tao↗

Analog Processor To Solve Optimization Problems

Proposed analog processor solves "traveling-salesman" problem, considered paradigm of global-optimization problems involving routing or allocation of resources. Includes electronic neural network and auxiliary circuitry based partly on concepts described in "Neural-Network Processor Would Allocate Resources" (NPO-17781) and "Neural Network Solves 'Traveling-Salesman' Problem" (NPO-17807). Processor based on highly parallel computing solves problem in significantly less time.

Duong, Tuan A.↗

General-Purpose Time-Accurate Adaptive-Grid Method

Adaptive-grid method provides for generation of time-accurate grids used to solve wide variety of problems numerically. Intended use in computation of compressible flows. Part of solution method including algebraic grid-movement scheme; algorithm solving, by use of explicit numerical-integration scheme, Euler equations of unsteady flow; and temporal coupling routine utilizing grid-prediction/correction method to link dynamic grid with Euler-equation-solving algorithm. Improves ability to adapt grid to multiple features of solution and involves use of coarse grid in prediction stage of temporal coupling reducing central-processing-unit cost of obtaining new grid.

Bockelie, Michael J.↗

Finite difference solution for transient radiative cooling of a conducting semitransparent square region

Transient solutions were obtained for a square region of heat conducting semitransparent material cooling by thermal radiation. The region is in a vacuum environment, so energy is dissipated only by radiation from within the medium leaving through its boundaries. The effect of heat conduction during the transient is to partially equalize the internal temperature distribution. As the optical thickness of the region is increased, the temperature gradients increase near the boundaries and corners, unless heat conduction is large. The solution procedure must provide accurate temperature distributions in these regions to prevent error in the calculated radiation losses. Two-dimensional numerical Gaussian integration is used to obtain the local radiative source term. A finite difference procedure with variable space and time increments is used to solve the transient energy equation. Variable spacing was used to concentrate grid points in regions with large temperature gradients.

Siegel, R.↗

A finite element study of incompressible flows past oscillating cylinders and aerofoils

We present our numerical results for certain unsteady flows past oscillating cylinders and aerofoils. The computations are based on the stabilized space-time finite element formulation. The implicit equation systems resulting from the space-time finite element discretizations are solved using iterative solution techniques. One of the problems studied is flow past a cylinder which is forced to oscillate in the horizontal direction. In this case we observe a change from an unsymmetric mode of vortex shedding to a symmetric one. An extensive study was carried out for the case in which a cylinder is mounted on lightly damped springs and allowed to oscillate in the vertical direction. In this case the motion of the cylinder needs to be determined as part of the solution, and under certain conditions this motion changes the vortex-shedding pattern of the flow field significantly. This nonlinear fluid-structure interaction exhibits certain interesting behavior such as 'lock-in' and 'hysteresis', which are in good agreement with the laboratory experiments carried out by other researchers in the past. Preliminary results for flow past a pitching aerofoil are also presented.

Mittal, S.↗

Operability engineering in the Deep Space Network

Many operability problems exist at the three Deep Space Communications Complexes (DSCC's) of the Deep Space Network (DSN). Four years ago, the position of DSN Operability Engineer was created to provide the opportunity for someone to take a system-level approach to solving these problems. Since that time, a process has been developed for personnel and development engineers and for enforcing user interface standards in software designed for the DSCC's. Plans are for the participation of operations personnel in the product life-cycle to expand in the future.

Wilkinson, Belinda↗

Transient Sub-Poissonian Distribution for Single-Mode Lasers

In this paper, the transient photon statistics for single-mode lasers is investigated by making use of the theory of quantum electrodynamics. By taking into account of the transitive time l,we obtain the master equation for Jaynes-Cummings model. The relation between the Mandel factor and the time is obtained by directly solving the master equation. The result shows that a transient phenomenon from the transient super-Poissonian distribution to the transient sub-Poissonian distribution occurs for single-mode lasers. In addition, the influences of the thermal light field and the cavity loss on the transient sub-Poissonian distribution are also studied.

Zang, J. Y.↗

Time Dependent Models of Grain Formation Around Carbon Stars

Carbon-rich Asymptotic Giant Branch stars are sites of dust formation and undergo mass loss at rates ranging from 10(exp -7) to 10(exp -4) solar mass/yr. The state-of-the-art in modeling these processes is time-dependent models which simultaneously solve the grain formation and gas dynamics problem. We present results from such a model, which also includes an exact solution of the radiative transfer within the system.

Egan, M. P.↗

Instability Mechanisms of Thermally-Driven Interfacial Flows in Liquid-Encapsulated Crystal Growth

During the past year, a great deal of effort was focused on the enhancement and refinement of the computational tools developed as part of our previous NASA grant. In particular, the interface mollification algorithm developed earlier was extended to incorporate the effects of surface-rheological properties in order to allow the study of thermocapillary flows in the presence of surface contamination. These tools will be used in the computational component of the proposed research in the remaining years of this grant. A detailed description of the progress made in this area is provided elsewhere. Briefly, the method developed allows for the convection and diffusion of bulk-insoluble surfactants on a moving and deforming interface. The novelty of the method is its grid independence: there is no need for front tracking, surface reconstruction, body-fitted grid generation, or metric evaluations; these are all very expensive computational tasks in three dimensions. For small local radii of curvature there is a need for local grid adaption so that the smearing thickness remains a small fraction of the radius of curvature. A special Neumann boundary condition was devised and applied so that the calculated surfactant concentration has no variations normal to the interface, and it is hence truly a surface-defined quantity. The discretized governing equations are solved subsequently using a time-split integration scheme which updates the concentration and the shape successively. Results demonstrate excellent agreement between the computed and exact solutions.

Haj-Hariri, Hossein↗