Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical continuation”

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 325 records · Page 18

Multigrid Methods for Fully Implicit Oil Reservoir Simulation

In this paper we consider the simultaneous flow of oil and water in reservoir rock. This displacement process is modeled by two basic equations: the material balance or continuity equations and the equation of motion (Darcy's law). For the numerical solution of this system of nonlinear partial differential equations there are two approaches: the fully implicit or simultaneous solution method and the sequential solution method. In the sequential solution method the system of partial differential equations is manipulated to give an elliptic pressure equation and a hyperbolic (or parabolic) saturation equation. In the IMPES approach the pressure equation is first solved, using values for the saturation from the previous time level. Next the saturations are updated by some explicit time stepping method; this implies that the method is only conditionally stable. For the numerical solution of the linear, elliptic pressure equation multigrid methods have become an accepted technique. On the other hand, the fully implicit method is unconditionally stable, but it has the disadvantage that in every time step a large system of nonlinear algebraic equations has to be solved. The most time-consuming part of any fully implicit reservoir simulator is the solution of this large system of equations. Usually this is done by Newton's method. The resulting systems of linear equations are then either solved by a direct method or by some conjugate gradient type method. In this paper we consider the possibility of applying multigrid methods for the iterative solution of the systems of nonlinear equations. There are two ways of using multigrid for this job: either we use a nonlinear multigrid method or we use a linear multigrid method to deal with the linear systems that arise in Newton's method. So far only a few authors have reported on the use of multigrid methods for fully implicit simulations. Two-level FAS algorithm is presented for the black-oil equations, and linear multigrid for two-phase flow problems with strong heterogeneities and anisotropies is studied. Here we consider both possibilities. Moreover we present a novel way for constructing the coarse grid correction operator in linear multigrid algorithms. This approach has the advantage in that it preserves the sparsity pattern of the fine grid matrix and it can be extended to systems of equations in a straightforward manner. We compare the linear and nonlinear multigrid algorithms by means of a numerical experiment.

Molenaar, J.↗

Alternatives to Re-Planning: Methods for Plan Re-Evaluation at Runtime

Current planning algorithms have difficulty handling the complexity that is due to an increase in domain uncertainty, and especially in the case of multi-dimensional continuous spaces. Therefore, they produce plans that do not take into account numerous situations that can occur at runtime, such as faults or other changes in the planning domain itself. Thus there is a gap between the plan generation and the reality experienced at runtime. Here we present two methods that allow the plan conditionals to be revised w.r.t. uncertainty on the system as estimated at runtime.

Benazera, Emmanuel↗

Self-Assembling, Flexible, Pre-Ceramic Composite Preforms

In this innovation, light weight, high temperature, compact aerospace structures with increased design options are made possible by using self-assembling, flexible, pre-ceramic composite materials. These materials are comprised of either ceramic or carbon fiber performs, which are infiltrated with polymer precursors that convert to ceramics upon thermal exposure. The preform architecture can vary from chopped fibers formed into blankets or felt, to continuous fibers formed into a variety of 2D or 3D weaves or braids. The matrix material can also vary considerably. For demonstration purposes, a 2D carbon weave was infiltrated with a SiC polymer precursor. The green or unfired material is fabricated into its final shape while it is still pliable. It is then folded or rolled into a much more compact shape, which will occupy a smaller space. With this approach, the part remains as one continuous piece, rather than being fabricated as multiple sections, which would require numerous seals for eventual component use. The infiltrated preform can then be deployed in-situ. The component can be assembled into its final shape by taking advantage of the elasticity of the material, which permits the structure to unfold and spring into its final form under its own stored energy. The pre-ceramic composites are converted to ceramics and rigidized immediately after deployment. The final ceramic composite yields a high-temperature, high-strength material suitable for a variety of aerospace structures. The flexibility of the material, combined with its high-temperature structural capacity after rigidization, leads to a less complex component design with an increased temperature range. The collapsibility of these structures allows for larger components to be designed and used, and also offers the potential for increased vehicle performance. For the case of collapsible nozzle extensions, a larger nozzle, and thus a larger nozzle exit plane, is possible because interference with surrounding structures can be avoided in the collapsed state. The larger exit plane leads to an increase in expansion area ratio, which has the potential to increase thrust and overall rocket performance. In general, the use of advanced ceramic materials can lead to improved engine and vehicle performance. The ceramics can run hotter, so less cooling is required. Fuel to coolant ratios can be balanced more readily to reduce weight. Engine efficiency can also be increased with hotter combustion and exhaust temperatures. In addition, the ceramic composites themselves can reduce the component weight by as much as 50 percent, which can translate into greater payload for the vehicle

Jaskowiak, Martha H.↗

Streamlining Pre-Launch Support for Missions with Highly Variable Trajectories

As the Goddard Space Flight Center's (GSFC) Flight Dynamics Facility (FDF) continues to ensure mission success in support of Human Space Flight (HSF) and Launch Vehicle (LV) missions, it continues to improve on different areas of mission support. The FDF processes numerous different types of trajectory profiles, and, for each trajectory, the FDF provides multiple acquisition data vectors to the Space Network (SN) White Sands Complex (WSC) for further processing and uplinking to the SN's Tracking and Data Relay Satellites (TDRSs), tracking the vehicle and maintaining communications throughout required support. Some trajectories target a specific orbit plane through Right Ascension of the Ascending Node (RAAN), or yaw, steering, which can yield to high variability across the launch window. Because of this, there can be a large number of trajectory profiles to process and generate acquisition data for that ensure TDRS is accurately pointing to the launch vehicle. Processing many trajectory profiles by FDF and WSC increases the complexity and level of effort associated with support. FDF employs a legacy boundary on the range difference between trajectories, which determines the number of discrete trajectory cases required to maintain communication with TDRS over the full launch window. An analysis was performed to determine a new boundary, taking into consideration current support requirements. The results show the boundary can be expanded beyond the currently employed legacy value without an impact on communication. In doing so, the number of trajectory profiles processed can be reduced, reducing the complexity and level of effort required, with no impact to mission safety.

Flight Operations↗

Ground-Based Vision Tracker for Advanced Air Mobility and Urban Air Mobility

Advanced Air Mobility (AAM) Air Mobility and Urban Air Mobility (UAM) require aircraft surveillance and monitoring for safety and security. Persistent tracking of flying objects provides Air Traffic Control (ATC) and Air Traffic Management (ATM) continuous coverage and knowledge of the national airspace (NAS). Since there are numerous more AAM and UAM aircraft than commercial aircraft, it will be challenging to utilize the same ATC/ATM architectures. A first step in creating a similar ATC/ATM architecture for AAM/UAM will require ground-based and airborne-based sensors to provide monitoring, which will be difficult in urban environments due to GPS degradation. This paper proposes a vision-based tracking method with static cameras by utilizing image subtraction and blob detection, which avoids adding additional electromagnetic interferences in the environment with sensors such as radar. The ground-based vision tracker (GBVT) outputs the detected objects' azimuth and elevation angles from unmanned aerial system (UAS) flight tests. Future and ongoing work includes sending the detected objects' azimuth and elevation angles as inputs for an extended Kalman filter (EKF) to estimate the position and velocity of the detected object.

distributed sensing↗

A numerical study of mixing in supersonic combustors with hypermixing injectors

A numerical study was conducted to evaluate the performance of wall mounted fuel-injectors designed for potential Supersonic Combustion Ramjet (SCRAM-jet) engine applications. The focus of this investigation was to numerically simulate existing combustor designs for the purpose of validating the numerical technique and the physical models developed. Three different injector designs of varying complexity were studied to fully understand the computational implications involved in accurate predictions. A dual transverse injection system and two streamwise injector designs were studied. The streamwise injectors were designed with swept ramps to enhance fuel-air mixing and combustion characteristics at supersonic speeds without the large flow blockage and drag contribution of the transverse injection system. For this study, the Mass-Averaged Navier-Stokes equations and the chemical species continuity equations were solved. The computations were performed using a finite-volume implicit numerical technique and multiple block structured grid system. The interfaces of the multiple block structured grid systems were numerically resolved using the flux-conservative technique. Detailed comparisons between the computations and existing experimental data are presented. These comparisons show that numerical predictions are in agreement with the experimental data. These comparisons also show that a number of turbulence model improvements are needed for accurate combustor flowfield predictions.

Lee, J.↗

A numerical study of mixing in supersonic combustors with hypermixing injectors

A numerical study was conducted to evaluate the performance of wall mounted fuel-injectors designed for potential Supersonic Combustion Ramjet (SCRAM-jet) engine applications. The focus of this investigation was to numerically simulate existing combustor designs for the purpose of validating the numerical technique and the physical models developed. Three different injector designs of varying complexity were studied to fully understand the computational implications involved in accurate predictions. A dual transverse injection system and two streamwise injector designs were studied. The streamwise injectors were designed with swept ramps to enhance fuel-air mixing and combustion characteristics at supersonic speeds without the large flow blockage and drag contribution of the transverse injection system. For this study, the Mass-Average Navier-Stokes equations and the chemical species continuity equations were solved. The computations were performed using a finite-volume implicit numerical technique and multiple block structured grid system. The interfaces of the multiple block structured grid systems were numerically resolved using the flux-conservative technique. Detailed comparisons between the computations and existing experimental data are presented. These comparisons show that numerical predictions are in agreement with the experimental data. These comparisons also show that a number of turbulence model improvements are needed for accurate combustor flowfield predictions.

Lee, J.↗

Cloud-Based Numerical Weather Prediction for Near Real-Time Forecasting and Disaster Response

The use of cloud computing resources continues to grow within the public and private sector components of the weather enterprise as users become more familiar with cloud‐computing concepts, and competition among service providers continues to reduce costs and other barriers to entry. Cloud resources can also provide capabilities similar to high‐performance computing environments, supporting multi‐node systems required for near real‐time, regional weather predictions. Referred to as "Infrastructure as a Service", or IaaS, the use of cloud-based computing hardware in an on‐demand payment system allows for rapid deployment of a modeling system in environments lacking access to a large, supercomputing infrastructure. Use of IaaS capabilities to support regional weather prediction may be of particular interest to developing countries that have not yet established large supercomputing resources, but would otherwise benefit from a regional weather forecasting capability. Recently, collaborators from NASA Marshall Space Flight Center and Ames Research Center have developed a scripted, on‐demand capability for launching the NOAA/NWS Science and Training Resource Center (STRC) Environmental Modeling System (EMS), which includes pre‐compiled binaries of the latest version of the Weather Research and Forecasting (WRF) model. The WRF‐EMS provides scripting for downloading appropriate initial and boundary conditions from global models, along with higher‐resolution vegetation, land surface, and sea surface temperature data sets provided by the NASA Short‐term Prediction Research and Transition (SPoRT) Center. This presentation will provide an overview of the modeling system capabilities and benchmarks performed on the Amazon Elastic Compute Cloud (EC2) environment. In addition, the presentation will discuss future opportunities to deploy the system in support of weather prediction in developing countries supported by NASA's SERVIR Project, which provides capacity building activities in environmental monitoring and prediction across a growing number of regional hubs throughout the world. Capacity‐building applications that extend numerical weather prediction to developing countries are intended to provide near real‐time applications to benefit public health, safety, and economic interests, but may have a greater impact during disaster events by providing a source for local predictions of weather‐related hazards, or impacts that local weather events may have during the recovery phase.

Molthan, Andrew↗

Vortex sheet modeling with higher order curved panels

A numerical technique is presented for modeling the vortex sheet with a deformable surface definition, along which a continuous vortex strength distribution in the spanwise direction is applied, so that by repeatedly modifying its shape, its true configuration is approached, in the proximity of its generating wing. Design problems requiring the inclusion of a realistic configuration of the vortex sheet are numerous. Examples discussed include: control effectiveness and stability derivatives, longitudinal stability, lateral stability, canards, propellers and helicopter rotors, and trailing vortex hazards.

Nagati, M. G.↗

Evaluation of a Drag-Free Control Concept for Missions in Low Earth Orbit

Atmospheric drag causes the greatest uncertainty in the equations of motion for spacecraft in Low Earth Orbit (LEO). If atmospheric drag eflects can be continuously and autonomously counteracted through the use of a drag-fee control system, drag may essentially be eliminated from the equations of motion for the spacecraft. The main perturbations on the spacecraft will then be those due to the gravitational field, which are much more easily predicted Through dynamical analysis and numerical simulation, this paper presents some potential costs and benefits associated with the fuel used during continuous drag compensation. In light of this cost-benefit analysis, simulation results are used to validate the concept of drag-free control for LEO spacecraft missions having certain characteristics.

Fleck, Melissa E.↗

Final Report on NASA Portable Laser Coating Removal Systems Field Demonstrations and Testing

Processes currently used throughout the National Aeronautics and Space Administration (NASA) to remove corrosion and coatings from structures, ground service equipment, small parts and flight components result in waste streams consisting of toxic chemicals, spent media blast materials, and waste water. When chemicals are used in these processes they are typically high in volatile organic compounds (VOC) and are considered hazardous air pollutants (HAP). When blast media is used, the volume of hazardous waste generated is increased significantly. Many of the coatings historically used within NASA contain toxic metals such as hexavalent chromium, and lead. These materials are highly regulated and restrictions on worker exposure continue to increase. Most recently the Occupational Safety and Health Administration (OSHA) reduced the permissible exposure limit (PEL) for hexavalent chromium (CrVI) from 52 to 5 micrograms per cubic meter of air as an 8-hour time-weighted average. Hexavalent chromium is found in numerous pretreatment and primer coatings used within the Space Shuttle Program. In response to the need to continue to protect assets within the agency and the growing concern over these new regulations, NASA is researching different ways to continue the required maintenance of both facility and flight equipment in a safe, efficient, and environmentally preferable manner. The use of laser energy to prepare surfaces for a variety of processes, such as corrosion and coating removal, weld preparation, and non destructive evaluation (NDE) is a relatively new application of the technology that has been proven to be environmentally preferable and in many cases less labor intensive than currently used removal methods. The novel process eliminates VOCs and blast media and captures the removed coatings with an integrated vacuum system. This means that the only waste generated are the coatings that are removed, resulting in an overall cleaner process. The development of a Portable Laser Coating Removal System (PLCRS) started as the goal of a Joint Group on Pollution Prevention (JG-PP) project, led by the Air Force, where several types of lasers in several configurations were thoroughly evaluated. Following this project, NASA decided to evaluate the best performers on processes and coatings specific to the agency. Laser systems used during this project were all of a similar design, between 40 and 500 Watts, most of which had integrated vacuum systems in order to collect materials removed from substrate surfaces during operation.

Rothgeb, Matthew J↗

Band models and correlations for infrared radiation

Absorption of infrared radiation by various line and band models are briefly reviewed. Narrow band model relations for absorptance are used to develop 'exact' formulations for total absorption by four wide band models. Application of a wide band model to a particular gas largely depends upon the spectroscopic characteristic of the absorbing-emitting molecule. Seven continuous correlations for the absorption of a wide band model are presented and each one of these is compared with the exact (numerical) solutions of the wide band models. Comparison of these results indicate the validity of a correlation for a particular radiative transfer application. In radiative transfer analyses, use of continuous correlations for total band absorptance provides flexibilities in various mathematical operations.

Tiwari, S. N.↗

Accurate computation and continuation of homoclinic and heteroclinic orbits for singular perturbation problems

In earlier papers, Doedel and the authors have developed a numerical method and derived error estimates for the computation of branches of heteroclinic orbits for a system of autonomous ordinary differential equations in R(exp n). The idea of the method is to reduce a boundary value problem on the real line to a boundary value problem on a finite interval by using a local (linear or higher order) approximation of the stable and unstable manifolds. A practical limitation for the computation of homoclinic and heteroclinic orbits has been the difficulty in obtaining starting orbits. Typically these were obtained from a closed form solution or via a homotopy from a known solution. Here we consider extensions of our algorithm which allow us to obtain starting orbits on the continuation branch in a more systematic way as well as make the continuation algorithm more flexible. In applications, we use the continuation software package AUTO in combination with some initial value software. The examples considered include computation of homoclinic orbits in a singular perturbation problem and in a turbulent fluid boundary layer in the wall region problem.

Vaughan, William W.↗

An implicit time-marching method for studying unsteady flow with massive separation

A fully implicit time-marching method is developed such that all spatial derivatives are approximated using central differences, but no use is made of any artificial dissipation. The numerical method solves the discretized equations using Alternating Direction Implicit-Block Gaussian Elimination technique. The method is implemented in the unsteady analysis, which solves the incompressible Navier-Stokes equations in terms of vorticity and stream function in generalized orthogonal coordinates. A clustered conformal C-grid is employed, and every effort is made to resolve the various length scales in the flow problem. The metric discontinuity at the branch-cut is treated appropriately using analytic continuation. Introduction of the BGE reordering permits implicit treatment of the branch cut in the numerical method. The vorticity singularity at the cusped trailing edge is also appropriately treated. This accurate and efficient implicit method is used to study flow at Re = 1000, past a 12-percent thick symmetric Joukowski airfoil at high angle of attack 30 and 53 deg.

Osswald, G. A.↗

Report of the Panel on Computer and Information Technology

Aircraft have become more and more dependent on computers (information processing) for improved performance and safety. It is clear that this activity will grow, since information processing technology has advanced by a factor of 10 every 5 years for the past 35 years and will continue to do so. Breakthroughs in device technology, from vacuum tubes through transistors to integrated circuits, contribute to this rapid pace. This progress is nearly matched by similar, though not as dramatic, advances in numerical software and algorithms. Progress has not been easy. Many technical and nontechnical challenges were surmounted. The outlook is for continued growth in capability but will require surmounting new challenges. The technology forecast presented in this report has been developed by extrapolating current trends and assessing the possibilities of several high-risk research topics. In the process, critical problem areas that require research and development emphasis have been identified. The outlook assumes a positive perspective; the projected capabilities are possible by the year 2000, and adequate resources will be made available to achieve them. Computer and information technology forecasts and the potential impacts of this technology on aeronautics are identified. Critical issues and technical challenges underlying the achievement of forecasted performance and benefits are addressed.

Lundstrom, Stephen F.↗

Numerical simulations of positively-biased probes and dielectric-conductor disks in a plasma

Plasma densities in the low-earth orbit range may be sufficient to cause difficulties for spacecraft operating at high voltages in this environment. The present investigation is concerned with the results of a continuing effort to develop a particle-in-cell (PIC), cylindrically-symmetric, 2-1/2 dimensional, self-consistent numerical simulation code. The simulation has the objective to explore the interactions of an ambient plasma with a conducting disk, which may be partially covered by a dielectric material. The disk and the surrounding dielectric material represent a hole in an insulator covering a conductor. Attention is given to a review of the simulation model, plain disk calculations, the disk and dielectric configuration, and the 'pinhole' effect.

Brandon, S. T.↗

A theory of supercritical wing sections, with computer programs and examples.

Mathematical methods for the design of supercritical wings, which depend on the numerical solution of the partial differential equations of two-dimensional gas dynamics, are developed. The main contribution is a computer program for the design of shockless transonic airfoils using the hodograph transformation and analytic continuation into the complex domain. The mathematical theory is described, and a manual for users of the programs is provided. Numerical examples are given and computational results are discussed, and the computer programs themselves are listed. The analysis routine can be used to ascertain whether the profiles behave well at off-design conditions, or to smooth coordinates and obtain a desirable shape more quickly when perfectly shockless flow is not essential.

Bauer, F.↗

Numerical Simulation of Combustion and Extinction of a Solid Cylinder in Low-Speed Cross Flow

The combustion and extinction behavior of a diffusion flame around a solid fuel cylinder (PMMA) in low-speed forced flow in zero gravity was studied numerically using a quasi-steady gas phase model. This model includes two-dimensional continuity, full Navier Stokes' momentum, energy, and species equations with a one-step overall chemical reaction and second-order finite-rate Arrhenius kinetics. Surface radiation and Arrhenius pyrolysis kinetics are included on the solid fuel surface description and a parameter Phi, representing the percentage of gas-phase conductive heat flux going into the solid, is introduced into the interfacial energy balance boundary condition to complete the description for the quasi-steady gas-phase system. The model was solved numerically using a body-fitted coordinate transformation and the SIMPLE algorithm. The effects of varying freestream velocity and Phi were studied. These parameters have a significant effect on the flame structure and extinction limits. Two flame modes were identified: envelope flame and wake flame. Two kinds of flammability limits were found: quenching at low-flow speeds due to radiative loss and blow-off at high flow speeds due to insufficient gas residence time. A flammability map was constructed showing the existence of maximum Phi above which the solid is not flammable at any freestream velocity.

Tien, J. S.↗