Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “computational costs”

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 127 records · Page 7

Rotor transient analysis

Undamped modes approximate dynamic behavior of rotors and bearings. Application of modal analysis to uncouple equations of motion simplifies stability, steady-state unbalance response, and transient response analysis of system; nonlinear stability is predicted from calculated frequency spectra. Analysis provides designers with complete information without involving large-scale computational costs. Programs are written in FORTRAN IV for use on CDC 6600 computer.

Allaire, P. E.↗

High-Fidelity Multidisciplinary Sensitivity Analysis Framework for Multipoint Rotorcraft Optimization

A multidisciplinary, gradient-based sensitivity-analysis methodology is evaluated for optimization of rotorcraft configurations. The tightly coupled discipline models include physics-based fluid dynamics and rotorcraft comprehensive analysis. A discretely-consistent adjoint method accounts for sensitivities of the unsteady flow and unstructured, dynamic, over set grids, while sensitivities of structural responses to aerodynamic loads are computed using a complex-variable method. The methodology is applied to optimize the shape of UH-60A Blackhawk helicopter blades for hover and forward flight conditions. The objective of the multipoint design is to simultaneously increase the rotorcraft figure of merit in a hover flight and reduce the rotor power in a forward flight. Trimmed loose-coupling solutions for the baseline configuration are used to initiate the tight-coupling multidisciplinary analysis. Target thrust and rolling and pitching moments are enforced as optimization constraints. The optimized configuration improves the optimization metrics at both design points. The improved performance and all constraints are maintained over many revolutions beyond the optimization interval, satisfying the required flight conditions. Computational cost of the optimization cycle is assessed in a high-performance computing environment and found affordable for design of rotorcraft in general level-flight conditions.

Wang, Li↗

RNS Applications for Interacting Sub- and Supersonic Flows

A solution based grid adaptation method that combines elements of the multigrid method for solution acceleration and the domain decomposition philosophy for grid optimization is described. Unlike other solution based adaptive gridding schemes, wherein the overhead of recomputing the grid and re-evaluating the solution on the adapted grid leads to higher computational costs compared to a non-adapted calculation, the present methodology reduces the computational time required to obtain the solution. The computational effort involved in the present calculation is significantly lower than a non-adapted calculation that utilizes the multigrid method purely as a convergence acceleration tool. In addition to convergence acceleration, the multigrid framework provides a mechanism of information transfer from regions wherein grid refinement is specified to unrefined coarse grid regions. The basis for domain decomposition in the current procedure is the variation in grid refinement requirements for each coordinate direction in different portions of the flow field. The method is demonstrated herein on an efficient set of governing equations termed the reduced Navier Stokes equations, applied in conjunction with a set of physical boundary conditions. The governing equations are discretized through a pressure based flux splitting procedure that is uniformly applicable from incompressible to supersonic Mach numbers.

Rubin, Stanley G.↗

Effects of Spatial Resolution on Retropropulsion Aerodynamics in an Atmospheric Environment

Development of a powered descent capability for atmospheric environments is heavily reliant on computational simulation. The prohibitive computational cost of such simulations motivates an improvement in the understanding of the minimum computational fidelity re-quired to accurately characterize aerodynamic-propulsive interference for such applications. This work examines the applicability of detached eddy simulation methods for retropropulsion in atmospheric environments through utilization of a GPU-accelerated computational framework, yielding data that are largely unachievable with conventional high-performance computing resources. This effort was specifically designed to quantitatively assess the effects of spatial resolution on vehicle aerodynamics for nominal operation of a low lift-to-drag ratio, human-scale Mars lander concept. The test matrix and scaling approach span relevant nozzle expansion conditions as well as mid-supersonic to high-subsonic operating conditions. Solutions were generated using computational grids ranging from 143 million to 1.14 billion grid points (degrees of freedom). This paper will provide an overview of the computational campaign, approach, and discussion of preliminary results focused on a range of operating conditions for a conceptual low lift-to-drag, human-scale Mars lander.

Ashley M Korzun↗

Modeling Deformable Linear Objects for Autonomous Robotic Outfitting of Lunar Surface Systems

This paper presents structural models of deformable linear objects (DLOs). DLOs are a subclass of deformable objects that encompasses common outfitting elements such as cables and ropes. Models are validated through hardware experiments, and integration in a robotic autonomy architecture for space environments is discussed. A persistent human presence on the lunar surface is one of the next major milestones in space exploration. This requires the development of robust extraplanetary construction technologies including structures and materials modeling and robotic systems. Previous robotic construction technology development has primarily focused on structural assembly, with significantly less focus on robotically performed outfitting tasks to instantiate subsystems providing power, data, life support, etc. These tasks involve manipulation of highly flexible elements, which are difficult to model, such as cable harnesses, ropes, and hoses. Robotic manipulation of DLOs, especially cable harnesses, is an active area of research as cable harnesses are essential for providing power and data to space assets. DLO models that can be used for robot manipulator trajectory generation are necessary for autonomous operation of lunar infrastructure. There are many proposed methods for modeling DLOs, and they primarily fall into three types: 1) discrete model-based, 2) continuum model-based, and 3) Neural Network-based. These types each have pros and cons, and the tradeoff between model accuracy and computational speed informs which type should be used. An understanding of this trade-off is imperative for real-time control of autonomous systems. High computational requirements reduce the speed of the model, making real-time control difficult, while accuracy is critical to preventing collisions. Discrete models, such as a mass-spring multibody representation, require relatively few calculations, and accuracy is directly tied to the step size of the discretization. Continuum models, such as a B-spline representation or a Cosserat rod model (a mix of continuous and discrete), are more informed of the structural properties of the cable and are much more accurate than a rigid body mass-spring model, but at significant computational cost. A Neural Network approach can provide an online solution with very few computational steps, but properly generating training data can be difficult and validation for an in-space application is not trivial. This paper explores the trade-off between different modeling approaches and compares accuracy and computational speed/complexity of the three types mentioned above. Model accuracy is evaluated using a cable in a static configuration. True cable shape is obtained using a depth camera for RGB images and point-cloud segmentation. The purpose of this experiment is to evaluate the trade-offs of different approaches to the DLO modeling problem. Understanding the tradeoffs between different cable modeling techniques paves the way for developing robotic control and planning architectures necessary for real-time manipulation of DLOs for lunar infrastructure outfitting. Real-time control is required for robotic systems to be able to actively manipulate a cable in a harsh environment where model and sensor errors compound, and environmental conditions can cause significant disturbances. Cable routing must be performed in areas with high density of objects/obstacles: through truss structures, near solar panels or mirror arrays, next to bundles of electrical equipment. Understanding the best way to plan and manipulate a cable without disrupting the environment or damaging the cable is imperative to robotic outfitting operations on the lunar surface.

Amy M Quartaro↗

Simulation of the Rotorwash Induced by a Quadrotor Urban Air Taxi in Ground Effect

A rotorcraft hovering near the ground causes downwash and outwash. The impact of these induced velocities on urban air mobility operations has not been extensively quantified. This paper explores the opportunity of using CFD to perform dedicated studies on the aerodynamics of rotorcraft in ground effect (IGE). For this purpose, we compare load and flow predictions obtained with the high-fidelity CFD solver OVERFLOW and with a medium-fidelity vortex particle-mesh (VPM) method, for a single rotor IGE. We show that the computational cost of high-fidelity simulation makes it impractical to sweep through a large number of configurations or operating conditions. However, a small set of cases can be used to verify a medium-fidelity tool. The latter provides a better trade-off between “accuracy” and computational intensity, which enables numerous, longer simulations at a more affordable cost. As an example, the outwash flow is computed for two different designs of a quadrotor air taxi in hover and reveals the existence of increased velocities between the rotors. This work illustrates how CFD can help identify dangerous areas for passengers and ground personnel when they approach the vehicle.

HECC↗

Uncertainty Quantification and Certification Prediction of Low-Boom Supersonic Aircraft Configurations

The primary objective of this work was to develop and demonstrate a process for accurate and efficient uncertainty quantification and certification prediction of low-boom, supersonic, transport aircraft. High-fidelity computational fluid dynamics models of multiple low-boom configurations were investigated including the Lockheed Martin SEEB-ALR body of revolution, the NASA 69 Delta Wing, and the Lockheed Martin 1021-01 configuration. A nonintrusive polynomial chaos surrogate modeling approach was used for reduced computational cost of propagating mixed, inherent (aleatory) and model-form (epistemic) uncertainty from both the computation fluid dynamics model and the near-field to ground level propagation model. A methodology has also been introduced to quantify the plausibility of a design to pass a certification under uncertainty. Results of this study include the analysis of each of the three configurations of interest under inviscid and fully turbulent flow assumptions. A comparison of the uncertainty outputs and sensitivity analyses between the configurations is also given. The results of this study illustrate the flexibility and robustness of the developed framework as a tool for uncertainty quantification and certification prediction of low-boom, supersonic aircraft.

West, Thomas K., IV↗

FPGA Implementation of Stereo Disparity with High Throughput for Mobility Applications

High speed stereo vision can allow unmanned robotic systems to navigate safely in unstructured terrain, but the computational cost can exceed the capacity of typical embedded CPUs. In this paper, we describe an end-to-end stereo computation co-processing system optimized for fast throughput that has been implemented on a single Virtex 4 LX160 FPGA. This system is capable of operating on images from a 1024 x 768 3CCD (true RGB) camera pair at 15 Hz. Data enters the FPGA directly from the cameras via Camera Link and is rectified, pre-filtered and converted into a disparity image all within the FPGA, incurring no CPU load. Once complete, a rectified image and the final disparity image are read out over the PCI bus, for a bandwidth cost of 68 MB/sec. Within the FPGA there are 4 distinct algorithms: Camera Link capture, Bilinear rectification, Bilateral subtraction pre-filtering and the Sum of Absolute Difference (SAD) disparity. Each module will be described in brief along with the data flow and control logic for the system. The system has been successfully fielded upon the Carnegie Mellon University's National Robotics Engineering Center (NREC) Crusher system during extensive field trials in 2007 and 2008 and is being implemented for other surface mobility systems at JPL.

Random access memory↗

Modeling Low-G Slosh Using Negative Mass

The modeling of the behavior of fluids in low- or microgravity is a field with little validated data but of great interest to space agencies and companies in the face of humanity’s expanding presence beyond low Earth orbit. Efforts to refuel and repair satellites in geostationary orbit require an understanding of the behavior of sloshing fuel in microgravity. The stability of remote sensing equipment and satellites depends on robust fluid slosh dampening or very accurate modeling, and the Artemis program includes the refueling of large landing systems in its architecture. However, current methodologies are quite binary: either use a relatively low-accuracy equivalent mechanical model or use computational fluid dynamics which are not practical for onboard systems due to their relatively high computational cost. In this work, a novel method of tracking the bubble rather than tracking the fuel itself is presented. Computational fluid dynamics simulations are used to validate the model based on available experimental data, and the dynamics of a refueler spacecraft based on Gateway’s Logistics Module are simulated both for the sloshing case and without sloshing considered.

Matthew M. Wittal↗

Explicit modeling and concurrent processing in the simulation of multibody dynamic systems

The objective is to present the activities at TRW in developing the capability to simulate the behavior of large flexible multibody space structures. The features of the simulation tools are: (1) to accommodate all rigid/flexible body degrees-of-freedom which incorporate the control system models and external forces, (2) to provide the flexibility to incorporate engineering-defined models and to retain parameters of significance to the engineer, (3) to reduce the computation cost by one order of magnitude (two orders of magnitude compared to a CRAY 1S), and (4) to keep it versatile so that radical variations in anticipated space structures can be accommodated. The current computer tools to simulate multibody systems appear not only to be very costly and time consuming, but also do not produce the desired fidelity of the mathematical models. In summary, a multibody simulation tool will be developed in the near future which will allow solution of the dynamics and controls of the deployment of the LDR backup structure, or the problem associated with the robotic assembly of the structure. The tools will allow the engineer to define the modeling technique and solve problems in less time and at reduced cost.

Gluck, R.↗

Rapid Aero Modeling for Urban Air Mobility Aircraft in Computational Experiments

Rapid Aero Modeling (RAM) applied to computational testing, RAM-C, is an approach to efficiently and automatically obtain aerodynamic models during computational investigations. RAM-C is designed to estimate models appropriate for flight dynamics studies and simulations. The approach responds to a demand for experimental efficiency and model fidelity that has increased with growing aircraft complexity and aerodynamic nonlinearities associated with hybrid and electric vertical takeoff and landing (eVTOL) aircraft. In an Urban Air Mobility (UAM) transportation system, it is expected that aircraft will embrace many features from both airplanes and rotorcraft. These vehicles present many more factors than conventional aircraft which can lead to increased computational costs and missed key factor interactions when applying traditional testing and modeling methods. RAM-C provides feedback loops around computational codes to rapidly guide testing toward aerodynamic models meeting user-defined fidelity goals. It combines and extends concepts from design of experiment theory and aircraft system identification theory that allow the user the freedom to choose, in advance of the test, a specific level of fidelity in terms of prediction error. RAM-C only collects enough data required to meet the user-specified prediction error requirements thus saving computational time and resources. The overall achievable fidelity of the final model also depends on the accuracy of the test facility, or in this case, the computational modeling approach. Previous studies to support development of the RAM-T process were conducted in wind tunnel tests to assess potential metrics, algorithms, and procedures. This paper presents results from the next steps taken and tests conducted for the development of RAM-C technology and highlights some of the unique features of RAM applied eVTOL configurations in a computational study.

Aerodynamics↗

Verifying X-Radiographs With Computed Tomographs

Nondestructive technique gives added confidence in inspection. Ambiguous indications in radiographic inspections of metal castings checked by computed tomography. Fast and inexpensive conventional x-ray inspection used to make film image of overall casting, and slower, more costly computed tomography used to reinspect relatively few parts of casting presenting possible diffraction patterns or other difficult-to-interpret features. Method effective in resolving ambiguities in radiographs of turbine blades. Provides same information as metallurgical sectioning.

Roy, Jagatjit↗

Aerodynamic design optimization using sensitivity analysis and computational fluid dynamics

A new and efficient method is presented for aerodynamic design optimization, which is based on a computational fluid dynamics (CFD)-sensitivity analysis algorithm. The method is applied to design a scramjet-afterbody configuration for an optimized axial thrust. The Euler equations are solved for the inviscid analysis of the flow, which in turn provides the objective function and the constraints. The CFD analysis is then coupled with the optimization procedure that uses a constrained minimization method. The sensitivity coefficients, i.e. gradients of the objective function and the constraints, needed for the optimization are obtained using a quasi-analytical method rather than the traditional brute force method of finite difference approximations. During the one-dimensional search of the optimization procedure, an approximate flow analysis (predicted flow) based on a first-order Taylor series expansion is used to reduce the computational cost. Finally, the sensitivity of the optimum objective function to various design parameters, which are kept constant during the optimization, is computed to predict new optimum solutions. The flow analysis of the demonstrative example are compared with the experimental data. It is shown that the method is more efficient than the traditional methods.

Baysal, Oktay↗

PCLAM: a Python Module for Computing Surface Lineloads and Moments

Lineloads serve a unique and important role in aerodynamic database development as well as configuration design and analysis. A new software suite was constructed which can compute lineloads rapidly enough that the calculations can run in tandem with high-fidelity fluid flow solvers. This enables the calculation of iteration- or time-dependent lineloads, which have thus far been too computationally costly to create for complicated systems like launch vehicles. The suite is organized into a standalone Python module named PCLAM (PCLAM Computes Lineloads And Moments) which can be imported into other software with minimal restructuring by the user or developer. The computed lineloads are integrated with a $C^0$ numerical quality and exhibit the expected sensitivity to underlying grid resolution. Even at low grid and lineload resolutions, the computed lineloads were found to be in strong agreement with several analytical test cases.

sectional loads↗

Near-Wall Models in Large Eddy Simulations of Flow Behind a Backward-Facing Step

Accurate large eddy simulation (LES) of a wall-bounded flow generally requires a near-wall resolution comparable to that in direct numerical simulation (DNS). As much as 50% of the total grid points and computational costs are expended in the near-wall regions in a typical simulation. This limits LES to fairly low Reynolds numbers on current computers. To perform practical flow applications at realistically high Reynolds numbers, such as flow over an airfoil, it is desirable to replace very thin, near-wall regions in the LES with easily and inexpensively computed wall models to specify the near-wall boundary conditions.

Cabot, W.↗

Spiking Neurons for Analysis of Patterns

Artificial neural networks comprising spiking neurons of a novel type have been conceived as improved pattern-analysis and pattern-recognition computational systems. These neurons are represented by a mathematical model denoted the state-variable model (SVM), which among other things, exploits a computational parallelism inherent in spiking-neuron geometry. Networks of SVM neurons offer advantages of speed and computational efficiency, relative to traditional artificial neural networks. The SVM also overcomes some of the limitations of prior spiking-neuron models. There are numerous potential pattern-recognition, tracking, and data-reduction (data preprocessing) applications for these SVM neural networks on Earth and in exploration of remote planets. Spiking neurons imitate biological neurons more closely than do the neurons of traditional artificial neural networks. A spiking neuron includes a central cell body (soma) surrounded by a tree-like interconnection network (dendrites). Spiking neurons are so named because they generate trains of output pulses (spikes) in response to inputs received from sensors or from other neurons. They gain their speed advantage over traditional neural networks by using the timing of individual spikes for computation, whereas traditional artificial neurons use averages of activity levels over time. Moreover, spiking neurons use the delays inherent in dendritic processing in order to efficiently encode the information content of incoming signals. Because traditional artificial neurons fail to capture this encoding, they have less processing capability, and so it is necessary to use more gates when implementing traditional artificial neurons in electronic circuitry. Such higher-order functions as dynamic tasking are effected by use of pools (collections) of spiking neurons interconnected by spike-transmitting fibers. The SVM includes adaptive thresholds and submodels of transport of ions (in imitation of such transport in biological neurons). These features enable the neurons to adapt their responses to high-rate inputs from sensors, and to adapt their firing thresholds to mitigate noise or effects of potential sensor failure. The mathematical derivation of the SVM starts from a prior model, known in the art as the point soma model, which captures all of the salient properties of neuronal response while keeping the computational cost low. The point-soma latency time is modified to be an exponentially decaying function of the strength of the applied potential. Choosing computational efficiency over biological fidelity, the dendrites surrounding a neuron are represented by simplified compartmental submodels and there are no dendritic spines. Updates to the dendritic potential, calcium-ion concentrations and conductances, and potassium-ion conductances are done by use of equations similar to those of the point soma. Diffusion processes in dendrites are modeled by averaging among nearest-neighbor compartments. Inputs to each of the dendritic compartments come from sensors. Alternatively or in addition, when an affected neuron is part of a pool, inputs can come from other spiking neurons. At present, SVM neural networks are implemented by computational simulation, using algorithms that encode the SVM and its submodels. However, it should be possible to implement these neural networks in hardware: The differential equations for the dendritic and cellular processes in the SVM model of spiking neurons map to equivalent circuits that can be implemented directly in analog very-large-scale integrated (VLSI) circuits.

Huntsberger, Terrance↗

Economical Third-Order Methods for Accurate Surface Heating Predictions on Simplex Element Meshes

A node-centered edge-based finite volume discretization of the compressible Navier-Stokes equations is presented with the heat flux reformulated as a first order system. A dissipation vector is derived for the reformulated system, such that the heat flux can be upgraded to O(h^3) on simplex element meshes in the same fashion as the inviscid fluxes. The method of manufactured solutions is used to demonstrate this uniform order property in heat flux. This new system is shown to efficiently produce accurate surface heating predictions on hypersonic hemisphere flow using an anisotropic simplex element mesh, achieving O(h^3) accuracy at relatively low computational cost compared to similar methods.

Computational Fluid Dynamics↗