A second-order feedback method for optimal control computations.
Second order feedback method for computing optimal control corrections from state perturbations, noting results for nonlinear examples based of van der Pol equation
Bullock, T. E.↗
SEARCH · Search NASA
Results for “computational methods”
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.
Numerical simulation of jet aerodynamics using the three-dimensional Navier-Stokes code PAB3D
This report presents a unified method for subsonic and supersonic jet analysis using the three-dimensional Navier-Stokes code PAB3D. The Navier-Stokes code was used to obtain solutions for axisymmetric jets with on-design operating conditions at Mach numbers ranging from 0.6 to 3.0, supersonic jets containing weak shocks and Mach disks, and supersonic jets with nonaxisymmetric nozzle exit geometries. This report discusses computational methods, code implementation, computed results, and comparisons with available experimental data. Very good agreement is shown between the numerical solutions and available experimental data over a wide range of operating conditions. The Navier-Stokes method using the standard Jones-Launder two-equation kappa-epsilon turbulence model can accurately predict jet flow, and such predictions are made without any modification to the published constants for the turbulence model.
Pao, S. Paul↗
A finite element method for the computation of transonic flow past airfoils
A finite element method for the computation of the transonic flow with shocks past airfoils is presented using the artificial viscosity concept for the local supersonic regime. Generally, the classic element types do not meet the accuracy requirements of advanced numerical aerodynamics requiring special attention to the choice of an appropriate element. A series of computed pressure distributions exhibits the usefulness of the method.
Eberle, A.↗
Preconditioned upwind methods to solve 3-D incompressible Navier-Stokes equations for viscous flows
A computational method for calculating low-speed viscous flowfields is developed. The method uses the implicit upwind-relaxation finite-difference algorithm with a nonsingular eigensystem to solve the preconditioned, three-dimensional, incompressible Navier-Stokes equations in curvilinear coordinates. The technique of local time stepping is incorporated to accelerate the rate of convergence to a steady-state solution. An extensive study of optimizing the preconditioned system is carried out for two viscous flow problems. Computed results are compared with analytical solutions and experimental data.
Hsu, C.-H.↗
2025 Workshop on Envisioning Frontiers in AI and Computing for Biological Research: Position Papers
This workshop aims to identify key research directions for transforming biology using artificial intelligence (AI), machine learning (ML) and computational methods to facilitate the discovery of new behaviors, mechanisms, and designs of biological processes relevant to DOE missions, underpinning a broader U.S. bioeconomy. By developing novel AI/ML technologies to analyze and interpret complex biological data, researchers can organize and simulate biological processes at various scales as well as advance predictive understanding and manipulation of biological systems. This integration of computation, experimentation, and next-generation experimental technologies can lead to discoveries in new biological behaviors and mechanisms relevant to DOE missions. The focus is on how advanced computational and mathematical methods can impact this mission by exploring digital twins, foundation models, automated laboratory experiments, modeling of complex living systems, and data-driven approaches for the biodesign of plants and microbial systems. While data management is important, it is not the primary focus of this workshop, which will assess the current state, trends, and AI/ML challenges at the interface between biology and computational science to identify opportunities for high-impact research at their intersection. The goal is to define research needs and opportunities that align with biological sciences, computational sciences, and applied mathematics research.
59 BASIC BIOLOGICAL SCIENCES↗
Implementation of an ADI method on parallel computers
The implementation of an ADI method for solving the diffusion equation on three parallel/vector computers is discussed. The computers were chosen so as to encompass a variety of architectures. They are: the MPP, an SIMD machine with 16K bit serial processors; FLEX/32, an MIMD machine with 20 processors; and CRAY/2, an MIMD machine with four vector processors. The Gaussian elimination algorithm is used to solve a set of tridiagonal systems on the FLEX/32 and CRAY/2 while the cyclic elimination algorithm is used to solve these systems on the MPP. The implementation of the method is discussed in relation to these architectures and measures of the performance on each machine are given. Simple performance models are used to describe the performance. These models highlight the bottlenecks and limiting factors for this algorithm on these architectures. Finally, conclusions are presented.
Fatoohi, Raad A.↗
Local and overall aerodynamic coefficients for bodies in hypersonic, rarefied flow
A computational method is given for the prediction of local pressure and viscous shear stress on windward surfaces of bluff, convex, axisymmetric or quasi-axisymmetric, hypersonic bodies in the transitional, rarefied flow regime. Overall aerodynamic forces and moments are computed by integration of the local quantities. The method is based upon a correlation of local pressure and shear stress computed by the direct simulation Monte Carlo (DSMC) numerical technique for cold wall, real gas conditions and some supplemental data from low-density, hypersonic wind tunnels. The relative simplicity of the method makes it feasible to do the necessary calculations with a personal computer. Two-dimensional shapes and leeward surfaces are not included in the scope of the method as it is presented here. Results are compared with DSMC computations for both local and overall coefficients. The latter includes sphere and blunt cone drag as well as lift and pitching moment coefficients for the NASA AFE vehicle at various angles of attack. Very satisfactory agreement is shown.
Potter, J. Leith↗
Numerical simulation of the transonic flowfield for wing/nacelle configurations
An efficient grid-interfacing zonal algorithm has been developed for computing the three-dimensional transonic flow field about wing/nacelle multicomponent configurations. The algorithm uses the full-potential formulation and the AF2 fully-implicit approximate factorization scheme. The flow field position is computed using a component-adaptive grid approach in which separate grids are employed for the individual components in the multicomponent configuration, where each component grid is optimized for a particular geometry such as the wing or nacelle. The wing and nacelle component grids are allowed to overlap, and flow field information is transmitted from one grid to another through the overlap region using trivariate interpolation. This paper presents a discussion of the computational methods used to generate both the wing and nacelle component grids, the technique used to interface the component grids, and the method used to obtain the inviscid multicomponent flow field solution. Computed results and correlations with experiment are presented to illustrate application of the analysis.
Atta, E. H.↗
Algorithms and software for nonlinear structural dynamics
The objective of this research is to develop efficient methods for explicit time integration in nonlinear structural dynamics for computers which utilize both concurrency and vectorization. As a framework for these studies, the program WHAMS, which is described in Explicit Algorithms for the Nonlinear Dynamics of Shells (T. Belytschko, J. I. Lin, and C.-S. Tsay, Computer Methods in Applied Mechanics and Engineering, Vol. 42, 1984, pp 225 to 251), is used. There are two factors which make the development of efficient concurrent explicit time integration programs a challenge in a structural dynamics program: (1) the need for a variety of element types, which complicates the scheduling-allocation problem; and (2) the need for different time steps in different parts of the mesh, which is here called mixed delta t integration, so that a few stiff elements do not reduce the time steps throughout the mesh.
Belytschko, Ted↗
Low speed aerodynamics CFD code validation
The objective of the research project is to develop and validate analytical methods for low-speed aerodynamics. The experimental needs for computational methods are presented. All data and results are presented in viewgraph format.
South, Jerry C.↗
Increasing Precision Of Temperature Sensors Of Liquid H2
Commercial silicon-diode temperature sensors intended for use in boiling or nearly boiling liquid hydrogen at temperatures near 37 degrees R recalibrated to greater precision by method involving careful attention to details of design, operation, and computation. Method based on fundamental electrical and thermodynamic principles and good engineering practice, also applicable to recalibration of other temperature sensors intended for use in other boiling or nearly boiling liquids.
Dempsey, Paula J.↗
DSMC Simulations of Ballute Aerothermodynamics Under Hypersonic Rarefied Conditions
This paper presents computational results obtained with the direct simulation Monte Carlo (DSMC) method for towed ballute applications. A ballute is an inflatable drag device that can be used to create a large amount of drag at high altitudes. Consequently, ballutes provide a potential technology for achieving aerocapture when the primary spacecraft velocity reduction (ΔV) is achieved at much higher altitudes than with the conventional rigid aeroshell. Since the ΔV is achieved at relatively high altitudes, rarefaction can be significant and is the motivation for the current study with the DSMC method. Computed surface and flow-field results are presented for a toroidal ballute, isolated tethers when exposed to free-stream flow conditions, and the flow interactions resulting from a toroidal ballute when towed by a six meter diameter Mars Pathfinder shaped (without tethers) spacecraft. All results presented are for Earth entry at velocities of 14 to 7 km/s (primary focus is at 8.55 km/s, same as some previous Titan aerocapture studies) and altitudes of 200 to 100 km. Variations of drag and heating coefficients as a function of rarefaction are presented. A description of the flow structure is provided and also an explanation of how it is affected by shock interactions produced solely by the ballute and those resulting from the two body combination of towed ballute and spacecraft is also given.
James N Moss↗
Aeroheating Test of CEV Entry Vehicle at Turbulent Conditions
An investigation of the aeroheating environment of the Project Orion Crew Entry Vehicle has been performed in the Arnold Engineering Development Center Tunnel 9. Data were measured on a approx. 3.5% scale model (0.1778m/7-inch diam.) of the vehicle using coaxial thermocouples in the Mach 8 and Mach 10 nozzles of Tunnel 9. Runs were performed at free stream Reynolds numbers of 1 106/ft to 20 10(exp 6)/ft in the Mach 10 nozzle and 8 10(exp 6)/ft to 48 10(exp 6)/ft in the Mach 8 nozzle. The test gas in Tunnel 9 is pure N2, which at these operating conditions remains un-dissociated and may be treated as a perfect gas. At these conditions, laminar, transitional, and turbulent flow was produced on the model at Mach 10, and transitional and turbulent conditions were produced on the model at Mach 8. The majority of runs were made on a clean, smooth-surface model configuration and a limited number of runs were made in which inserts with varying boundary-layer trips configurations were used to force the occurrence of transition. Laminar and turbulent predictions were generated for all wind tunnel test conditions and comparisons were performed with the data for the purpose of helping to define uncertainty margins for the computational method. Data from both the wind tunnel test and the computations are presented herein. Figure 1 shows a schematic of the thermocouple locations on the model and figures 2 and 3 show a photo and schematic of the AEDC Hypervelocity Tunnel 9. Figure 4 shows a typical grid used in the computations. From the comparisons shown in figures 5 through 8 it was concluded that for perfect-gas conditions, the computations could predict either fully-laminar or full-turbulent flow to within +/-10% of the experimental data. The experimental data showed that transition began on the leeside of the heatshield at a free stream Reynolds number of 9 10(exp 6)/ft in the Mach 10 nozzle and fully-developed turbulent flow was produced at 20 10(exp 6)/ft. In the Mach 8 nozzle, transition on the leeside of the heat-shield was observed for all test conditions, and full-developed turbulent flow occurred at a free stream Reynolds number of 18 10(exp 6)/ft. On the aftbody of the vehicle no evidence of turbulence was detected at Mach 10 conditions, and at Mach 8 conditions, transition appeared to begin on the windside of the aftbody at free stream Reynolds number of 18x10(exp 6)/ft with fully-developed turbulent flow occurring only at the highest test condition of 48x10(exp 6)/ft.
Hollis, Brian R.↗
Evolutionary Computation for the Identification of Emergent Behavior in Autonomous Systems
Over the past several years the Center for Evolutionary Computation and Automated Design at the Jet Propulsion Laboratory has developed a technique based on Evolutionary Computational Methods (ECM) that allows for the automated optimization of complex computationally modeled systems. An important application of this technique is for the identification of emergent behaviors in autonomous systems. Mobility platforms such as rovers or airborne vehicles are now being designed with autonomous mission controllers that can find trajectories over a solution space that is larger than can reasonably be tested. It is critical to identify control behaviors that are not predicted and can have surprising results (both good and bad). These emergent behaviors need to be identified, characterized and either incorporated into or isolated from the acceptable range of control characteristics. We use cluster analysis of automatically retrieved solutions to identify isolated populations of solutions with divergent behaviors.
Terrile, Richard J.↗
Model-Invariant Hybrid Computations of Separated Flows for RCA Standard Test Cases
NASA's Revolutionary Computational Aerosciences (RCA) subproject has identified several smooth-body separated flows as standard test cases to emphasize the challenge these flows present for computational methods and their importance to the aerospace community. Results of computations of two of these test cases, the NASA hump and the FAITH experiment, are presented. The computations were performed with the model-invariant hybrid LES-RANS formulation, implemented in the NASA code VULCAN-CFD. The model- invariant formulation employs gradual LES-RANS transitions and compensation for model variation to provide more accurate and efficient hybrid computations. Comparisons revealed that the LES-RANS transitions employed in these computations were sufficiently gradual that the compensating terms were unnecessary. Agreement with experiment was achieved only after reducing the turbulent viscosity to mitigate the effect of numerical dissipation. The stream-wise evolution of peak Reynolds shear stress was employed as a measure of turbulence dynamics in separated flows useful for evaluating computations.
Woodruff, Stephen↗
Magnetic Field Suppression of Flow in Semiconductor Melt
One of the most promising approaches for the reduction of convection during the crystal growth of conductive melts (semiconductor crystals) is the application of magnetic fields. Current technology allows the experimentation with very intense static fields (up to 80 KGauss) for which nearly convection free results are expected from simple scaling analysis in stabilized systems (vertical Bridgman method with axial magnetic field). However, controversial experimental results were obtained. The computational methods are, therefore, a fundamental tool in the understanding of the phenomena accounting during the solidification of semiconductor materials. Moreover, effects like the bending of the isomagnetic lines, different aspect ratios and misalignments between the direction of the gravity and magnetic field vectors can not be analyzed with analytical methods. The earliest numerical results showed controversial conclusions and are not able to explain the experimental results. Although the generated flows are extremely low, the computational task is a complicated because of the thin boundary layers. That is one of the reasons for the discrepancy in the results that numerical studies reported. Modeling of these magnetically damped crystal growth experiments requires advanced numerical methods. We used, for comparison, three different approaches to obtain the solution of the problem of thermal convection flows: (1) Spectral method in spectral superelement implementation, (2) Finite element method with regularization for boundary layers, (3) Multiquadric method, a novel method with global radial basis functions, that is proven to have exponential convergence. The results obtained by these three methods are presented for a wide region of Rayleigh and Hartman numbers. Comparison and discussion of accuracy, efficiency, reliability and agreement with experimental results will be presented as well.
Fedoseyev, A. I.↗
Overview of CFD methods and comparison with flight aerothermal data
An overview of previously published aerothermal investigations which demonstrate the capabilities of detailed computational fluid dynamics and engineering codes to predict the aerothermal environment about an entry vehicle is presented. The overview consists of a brief discussion of the computational methods and experimental data and includes comparisons between the computed results and data. The overview focuses primarily on analyses of flight data since these data provide the unique capability to assess the real-gas chemistry options in the codes. The computed results are based on a series of codes which are employed by the Aerothermodynamics Branch of the Space System Division at the Langley Research Center. The flight data, which were measured on the Reentry F, the Space Shuttle, and the Fire II vehicles, represent a wide range of vehicle configurations and freestream conditions. Also, results of one recent set of ground tests are included since the tests provide data on a model of a pending flight project. The comparisons of the predicted results and data demonstrate the adequacy of the present computational fluid dynamics capabilities and indicate the potential to predict the aerothermal environment about future flight vehicles.
Sutton, Kenneth↗
Analysis of steranes and triterpanes in geolipid extracts by automatic classification of mass spectra
A computer method is described for the automatic classification of triterpanes and steranes into gross structural type from their mass spectral characteristics. The method has been applied to the spectra obtained by gas-chromatographic/mass-spectroscopic analysis of two mixtures of standards and of hydrocarbon fractions isolated from Green River and Messel oil shales. Almost all of the steranes and triterpanes identified previously in both shales were classified, in addition to a number of new components. The results indicate that classification of such alkanes is possible with a laboratory computer system. The method has application to diagenesis and maturation studies as well as to oil/oil and oil/source rock correlations in which rapid screening of large numbers of samples is required.
Wardroper, A. M. K.↗