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 595 records · Page 33

Numerical investigations in three-dimensional internal flows

In the present reporting period, the 3D STUFF code was used to solve the underbody flow for the waverider configuration vehicle. In order to start the space-marched version of the code, the time-marched version (TUFF) was used to solve the forward portion of the underside of the forebody flowfield. A grid was generated which went from the tip of the nose to the location of the cowl lip. This includes all of the inlet ramp system. Previous indications were that three-dimensional effects could be expected on the ramps of such an aircraft. For purposes of the present study, no sidewalls were assumed. The sidewalls were eliminated to simplify the calculations and to show the potential effects of three-dimensional flow in the absence of a full sidewall. Further, this flow was also analyzed using the newly released OVERFLOW code and comparisons between the two codes were made. In addition to these 3D calculations, 2D calculations using the OVERFLOW code were also obtained for the Mach 5 inlet model in this reporting period. Comparisons between the experimental data, previous computational fluid dynamics (CFD) results and those from OVERFLOW were made.

Rose, William C.↗

Anticipation of the landing shock phenomenon in flight simulation

An aircraft landing may be described as a controlled crash because a runway surface is intercepted. In a simulation model the transition from aerodynamic flight to weight on wheels involves a single computational cycle during which stiff differential equations are activated; with a significant probability these initial conditions are unrealistic. This occurs because of the finite cycle time, during which large restorative forces will accompany unrealistic initial oleo compressions. This problem was recognized a few years ago at Ames Research Center during simulation studies of a supersonic transport. The mathematical model of this vehicle severely taxed computational resources, and required a large cycle time. The ground strike problem was solved by a described technique called anticipation equations. This extensively used technique has not been previously reported. The technique of anticipating a significant event is a useful tool in the general field of discrete flight simulation. For the differential equations representing a landing gear model stiffness, rate of interception and cycle time may combine to produce an unrealistic simulation of the continuum.

Mcfarland, Richard E.↗

Full potential unsteady computations including aeroelastic effects

A unified formulation is presented based on the full potential framework coupled with an appropriate structural model to compute steady and unsteady flows over rigid and flexible configurations across the Mach number range. The unsteady form of the full potential equation in conservation form is solved using an implicit scheme maintaining time accuracy through internal Newton iterations. A flux biasing procedure based on the unsteady sonic reference conditions is implemented to compute hyperbolic regions with moving sonic and shock surfaces. The wake behind a trailing edge is modeled using a mathematical cut across which the pressure is satisfied to be continuous by solving an appropriate vorticity convection equation. An aeroelastic model based on the generalized modal deflection approach interacts with the nonlinear aerodynamics and includes both static as well as dynamic structural analyses capability. Results are presented for rigid and flexible configurations at different Mach numbers ranging from subsonic to supersonic conditions. The dynamic response of a flexible wing below and above its flutter point is demonstrated.

Shankar, Vijaya↗

Crew Performance Support System to Aid in Anomaly Resolution: Concept of Operations

As missions progress into deep space, communication delays and disruptions will disenable the crew’s reliance on Earth experts. There are also limitations in the amount of data that can be downlinked to the ground. It is prudent to assume that critical, complex vehicle or habitat sub-systems will malfunction at a time when a lunar or Mars’ crew cannot rely on the Earth-Support team to detect, diagnose and resolve the problem and it is impractical to expect a small crew to step-in with the same level of expertise as 50+ authorities. The crew will need novel processes and advanced technological support to independently identify and resolve safety- and time-critical anomalies. That a self-reliant crew is unable to respond appropriately to time-critical anomalies is a significant risk to crew safety and mission success. This risk is driven by several factors; novel and unanticipated anomalies would not have been trained pre-flight, the crew could forget their pre-flight training or spaceflight stressors could impair the crew’s problem-solving ability. At last year’s IWS, Beard reported that a single spaceflight stressor (elevated CO2) could undermine the crew’s ability to independently respond to emergencies. Concept of Operations (ConOps) provide a common view of future system functions to all stakeholders. For the current project, a ConOps was developed that describes the operational processes, practices and capabilities needed by a crew of astronauts on deep space missions to autonomously respond to anticipated and unanticipated anomalies. It is crucial to recognize that, as of August 2018 existing technologies are unable to effectively support crew anomaly response to unanticipated events. “Intelligent technology” has not reached a maturity level that permits generalizing a solution to novel situations. For example, to train intelligent technology requires volumes of data that do not exist. The complexities involved in a manned mission to Mars cannot be compared to sending rovers to Mars using scripted software. This ConOps proposes a Crew Performance Support System (CPSS) that will push NASA and its industry partners toward what will be required for a safe and successful manned mission to Mars. Anomaly resolution during a deep space mission will take place within a dynamic, or changing, context. The figure to the left shows five broad contextual variables: the organizational culture, mission context, system characteristics, team characteristics and individual characteristics. The yellow arrow indicates that spaceflight and task-related stressors can affect system, team and individual crewmember characteristics and therefore anomaly response potential. The figure depicts a protective umbrella of Human-System Integration (HSI) principles that should be instituted during CPSS development including a balanced workload, shared situation awareness and building an appropriate level of trust in the automation. The figure also depicts two interrelated and cooperative components, an HSI Data System and other Enabling Capabilities will be required to support crew anomaly response and Earth-Support situation awareness. As we journey from ISS to Gateway to Mars, multiple, simultaneous and integrated research and development efforts (i.e., support systems co-evolution) must be implemented to meet the problem-solving challenges a self-reliant crew will face on a Mars’ mission. The crossovers between the capabilities are just as important as the discrete capabilities themselves. As the capabilities mature, the lines between the support subdomains will blur and an integrated system will emerge. The ConOps summarizes current knowledge about how highly trained people solve anomalies in safety- and time-critical situations, describes a group of capabilities that could help to reduce the extant risk and documents requirements levied on additional systems that provides critical inputs to the CPSS. Scenarios are used to promote a shared understanding of processes, practices and technological goals needed for safe and productive manned missions beyond LEO.

HSIA risk↗

Explicit large time-step schemes for the shallow water equations

Modifications to explicit finite difference schemes for solving the shallow water equations for meteorological applications by increasing the time step for the fast gravity waves are analyzed. Terms associated with the gravity waves in the shallow water equations are treated on a coarser grid than those associated with the slow Rossby waves, which contain much more of the available energy and must be treated with higher accuracy, enabling a several-fold increase in time step without degrading the accuracy of the solution. The method is presented in Cartesian and spherical coordinates for a rotating earth, using generalized leapfrog, frozen coefficient, and Fourier filtering finite difference schemes. Computational results verify the numerical stability of the approach.

Turkel, E.↗

Leading-edge receptivity to a vortical freestream disturbance: A numerical analysis

The receptivity to freestream vorticity of the boundary layer over a flat plate with an elliptic leading edge is investigated numerically. The flow is simulated by solving the incompressible Navier-Stokes system in general curvilinear coordinates with the vorticity and stream function as dependent variables. A finite-difference scheme which is second-order accurate in both space and time is used. As a first step, the steady basic-state solution is computed. Then a small amplitude vortical disturbance is introduced at the upstream boundary and the governing equations are solved time-accurately to evaluate the spatial and temporal growth of the perturbations leading to instability waves (Tollmien-Schlichting waves) inside the boundary layer. Preliminary results for a symmetric, 2-D disturbance reveal the presence of Tollmien-Schlichting waves aft of the flat-plate/ellipse juncture.

Buter, Thomas A.↗

An efficient parallel algorithm for the solution of a tridiagonal linear system of equations

Tridiagonal linear systems of equations are solved on conventional serial machines in a time proportional to N, where N is the number of equations. The conventional algorithms do not lend themselves directly to parallel computations on computers of the ILLIAC IV class, in the sense that they appear to be inherently serial. An efficient parallel algorithm is presented in which computation time grows as log sub 2 N. The algorithm is based on recursive doubling solutions of linear recurrence relations, and can be used to solve recurrence relations of all orders.

Stone, H. S.↗

An efficient parallel algorithm for the solution of a tridiagonal linear system of equations.

Tridiagonal linear systems of equations can be solved on conventional serial machines in a time proportional to N, where N is the number of equations. The conventional algorithms do not lend themselves directly to parallel computation on computers of the Illiac IV class, in the sense that they appear to be inherently serial. An efficient parallel algorithm is presented in which computation time grows as log(sub-2) N. The algorithm is based on recursive doubling solutions of linear recurrence relations, and can be used to solve recurrence relations of all orders.

Stone, H. S.↗

User's guide for the computer code COLTS for calculating the coupled laminar and turbulent flow over a Jovian entry probe

A user's guide for a computer code 'COLTS' (Coupled Laminar and Turbulent Solutions) is provided which calculates the laminar and turbulent hypersonic flows with radiation and coupled ablation injection past a Jovian entry probe. Time-dependent viscous-shock-layer equations are used to describe the flow field. These equations are solved by an explicit, two-step, time-asymptotic finite-difference method. Eddy viscosity in the turbulent flow is approximated by a two-layer model. In all, 19 chemical species are used to describe the injection of carbon-phenolic ablator in the hydrogen-helium gas mixture. The equilibrium composition of the mixture is determined by a free-energy minimization technique. A detailed frequency dependence of the absorption coefficient for various species is considered to obtain the radiative flux. The code is written for a CDC-CYBER-203 computer and is capable of providing solutions for ablated probe shapes also.

Kumar, A.↗

On the Implementation and Further Validation of a Time Domain Boundary Element Method Broadband Impedance Boundary Condition

A time domain boundary integral equation with Burton-Miller reformulation is presented for acoustic scattering by surfaces with liners in a uniform mean flow. The Ingard-Myers impedance boundary condition is implemented using a broadband multipole impedance model and converted into time domain differential equations to augment the boundary integral equation. The coupled integral-differential equations are solved numerically by a March-On-in-Time (MOT) scheme. While the Ingard-Myers condition is known to support Kelvin-Helmholtz instability due to its use of a vortex sheet interface between the flow and the liner surface, it is found that by neglecting a second derivative term in the current time domain impedance boundary condition formulation, the instability an be effectively suppressed in computation. The proposed formulation and implementation are validated using a NASA Langley Research Center Grazing Flow Impedance Tube (GFIT) experimental dataset with satisfactory results. Moreover, a minimization procedure for finding the poles and coefficients of the broadband multiple impedance model is formulated in this paper by which, unlike the commonly used vector-fitting method, passivity of the model is ensured. Numerical tests show the proposed minimization approach is effective for modeling liners that are commonly used in aeroacoustic applications.

Time Domain Impedance Boundary Condition↗

Acoustic eigenmodes of corrugated ducts

Acoustic propagation in two-dimensional corrugated ducts is analyzed. A transform is made to a coordinate system that matches the duct walls. The resulting wave equation separates in space and time. The spatial differential equation is then solved by expanding in the eigenfunctions of a flat-walled duct. The prescription for this expansion is given by the appropriate application of time-independent perturbation theory. Calculations of first-order corrections for a periodic duct of both finite and infinite length are presented.

Succi, G. P.↗

Numerical Solution of the Incompressible Navier-Stokes Equations in Three-dimensional Generalized Curvilinear Coordinates

Numerically solving the incompressible Navier-Stokes equations is known to be time consuming and expensive. Testing of the INS3D computers code, which solves these equations with the use of the pseudocompressibility method, shows this method to be an efficient way to obtain the steady state solution. The effects of the waves introduced by the pseudocompressibility method are analyzed and criteria are set and tested for the choice of the pseudocompressibility parameter which governs the artificial sound speed. The code is tested using laminar flow over a two dimensional backward-facing step, and laminar flow over a two dimensional circular cylinder. The results of the computations over the backward-facing step are in excellent agreement with experimental results. The transient solution of the flow over the cylinder impulsively started from rest is in good agreement with experimental results. However, the computed frequency of periodic shedding of vortices behind the cylinder is not in agreement with the experimental value. For a three dimensional test case, computations were conducted for a cylinder end wall junction. The saddle point separation and horseshoe vortex system appear in the computed field. The solution also shows secondary vortex filaments which wrap around the cylinder and spiral up in the wake.

Rogers, S. E.↗

Simulations of Solar Wind Turbulence

Recently we have restructured our approach to simulating magnetohydrodynamic (MHD) turbulence in the solar wind. Previously, we had defined a 'virtual' heliosphere that contained, for example, a tilted rotating current sheet, microstreams, quasi-two-dimensional fluctuations as well as Alfven waves. In this new version of the code, we use the global, time-stationary, WKB Alfven wave-driven solar wind model developed by Usmanov and described in Usmanov and Goldstein [2003] to define the initial state of the system. Consequently, current sheets, and fast and slow streams are computed self-consistently from an inner, photospheric, boundary. To this steady-state configuration, we add fluctuations close to, but above, the surface where the flow become super-Alfvenic. The time-dependent MHD equations are then solved using a semi-discrete third-order Central Weighted Essentially Non-Oscillatory (CWENO) numerical scheme. The computational domain now includes the entire sphere; the geometrical singularity at the poles is removed using the multiple grid approach described in Usmanov [1996]. Wave packets are introduced at the inner boundary such as to satisfy Faraday's Law [Yeh and Dryer, 1985] and their nonlinear evolution are followed in time.

Goldstein, Melvyn L.↗

Numerical simulation of axisymmetric turbulent jet flow

The transport equations for the turbulent Reynolds stresses and energy dissipation rate in conjunction with the governing equations for the vorticity and stream function are numerically solved for in the case of an axisymmetric jet into stagnant surroundings. The time dependence of the equations is retained allowing for any transient results to be interpreted in terms of conditionally averaged flow quantities. The general form of the governing differential equations is presented as well as the numerical procedure used in the solution. Comments concerning the inherent limitations in solving the time-dependent set of equations are made and the versatility of the approach examined. Comparisons of the various components of the stationary Reynolds stress tensor are made with experimental results.

Gatski, T. B.↗

Multigrid, Fractional-Step Computation Of Flow

Speed of computer code solving three-dimensional Navier-Stokes equations of flow of incompressible fluid by fractional-step method increased significantly by use of multigrid procedures. In method, equations solved on general nonorthogonal curvilinear coordinate grid, using volume fluxes. At each time step, computations performed in two fractional steps. In first step, equations of conservation of momentum solved by use of gradient of pressure from previous time step via explicit approximate-factorization method, yielding approximate flow field that does not satisfy equation of conservation of mass. In second step, discrete Poisson-like equation with Neumann-type boundary conditions, formed by combining equations of conservation of momentum and mass, solved iteratively.

Kwak, Dochan↗

Time-split finite-volume method for three-dimensional blunt-body flow

An efficient numerical method for calculating plane, axisymmetric, and fully three-dimensional blunt-body flow is presented. It is a second-order-accurate, time-dependent finite-volume procedure that solves the Euler equations in integral conservation-law form. These equations are written with respect to a Cartesian coordinate system in which an embedded mesh adjusts in time to the motion of the bow shock that is automatically captured as part of the weak solution. With such an adjusting mesh, oscillations in flow properties near the shock are shown to be virtually eliminated. The scheme uses a time-splitting concept that accelerates the convergence appreciably. Comparisons are made between computed and experimental results.

Rizzi, A. W.↗

Analysis of turbulent flow about an isolated airfoil using a time dependent Navier-Stokes procedure

The procedure solves the Navier-Stokes equations by the consistently split linearized block implicit method of Briley and McDonald in a body fitted coordinate system. The procedure is described and results are presented for flow about an airfoil whose incidence changes from 6 degrees to 19 degrees at a Reynolds number of one million and Mach number of 0.2. In addition, the unsteady flow about an airfoil held at a constant 19 degree incidence is examined and compared to data.

Shamroth, S. J.↗

High-Resolution Genuinely Multidimensional Solution of Conservation Laws by the Space-Time Conservation Element and Solution Element Method

In this overview paper, we review the basic principles of the method of space-time conservation element and solution element for solving the conservation laws in one and two spatial dimensions. The present method is developed on the basis of local and global flux conservation in a space-time domain, in which space and time are treated in a unified manner. In contrast to the modern upwind schemes, the approach here does not use the Riemann solver and the reconstruction procedure as the building blocks. The drawbacks of the upwind approach, such as the difficulty of rationally extending the 1D scalar approach to systems of equations and particularly to multiple dimensions is here contrasted with the uniformity and ease of generalization of the Conservation Element and Solution Element (CE/SE) 1D scalar schemes to systems of equations and to multiple spatial dimensions. The assured compatibility with the simplest type of unstructured meshes, and the uniquely simple nonreflecting boundary conditions of the present method are also discussed. The present approach has yielded high-resolution shocks, rarefaction waves, acoustic waves, vortices, ZND detonation waves, and shock/acoustic waves/vortices interactions. Moreover, since no directional splitting is employed, numerical resolution of two-dimensional calculations is comparable to that of the one-dimensional calculations. Some sample applications displaying the strengths and broad applicability of the CE/SE method are reviewed.

Himansu, Ananda↗