Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Incremental Computing”

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 199 records · Page 11

Applying Incremental, Inductive Model Checking to the Modal Mu Calculus

This work aims to demonstrate that an incremental inductive model checking algorithm built on top of Boolean satisfiability (SAT) solvers can be extended to support modal mu calculus (MMC) formulas. The resulting algorithm, called modal mu calculus model checking using myopic constraints (MC3), solves MMC model checking problems over Boolean labeled transition systems (LTSs). MMC subsumes simple invariance/reachability (as solved by the IC3 algorithm), linear temporal logic (LTL, as solved by the fair algorithm), computation tree logic (CTL, as solved by IICTL), and CTL* (which in turn subsumes LTL and CTL, but was not previously supported by any incremental inductive algorithm). The algorithm is implemented in a prototype solver, mc3.

97 MATHEMATICS AND COMPUTING↗

Then and Now: Improving Software Portability, Productivity, and 100× Performance

The US Exascale Computing Project (ECP) has succeeded in preparing applications to run efficiently on the first reported Exascale supercomputers in the world. To achieve this, it modernized the whole leadership software stack, from libraries to simulation codes. In this article, we contrast selected leadership software before and after ECP. We discuss how sustainable research software development for leadership computing can embrace the conversation with the hardware vendors, the leadership computing facilities, the software community, and the domain scientists who are the application developers and integrators of software products. We elaborate on how software needs to take portability as a central design principle and to benefit from interdependent teams; we also demonstrate how moving to programming languages with high momentum, like modern C++, can help improve the sustainability, interoperability, and performance of research software. Finally, we showcase how cross-institutional efforts can enable algorithm advances that are beyond incremental performance optimization.

97 MATHEMATICS AND COMPUTING↗

Transonic Drag Prediction on a DLR-F6 Transport Configuration Using Unstructured Grid Solvers

A second international AIAA Drag Prediction Workshop (DPW-II) was organized and held in Orlando Florida on June 21-22, 2003. The primary purpose was to inves- tigate the code-to-code uncertainty. address the sensitivity of the drag prediction to grid size and quantify the uncertainty in predicting nacelle/pylon drag increments at a transonic cruise condition. This paper presents an in-depth analysis of the DPW-II computational results from three state-of-the-art unstructured grid Navier-Stokes flow solvers exercised on similar families of tetrahedral grids. The flow solvers are USM3D - a tetrahedral cell-centered upwind solver. FUN3D - a tetrahedral node-centered upwind solver, and NSU3D - a general element node-centered central-differenced solver. For the wingbody, the total drag predicted for a constant-lift transonic cruise condition showed a decrease in code-to-code variation with grid refinement as expected. For the same flight condition, the wing/body/nacelle/pylon total drag and the nacelle/pylon drag increment predicted showed an increase in code-to-code variation with grid refinement. Although the range in total drag for the wingbody fine grids was only 5 counts, a code-to-code comparison of surface pressures and surface restricted streamlines indicated that the three solvers were not all converging to the same flow solutions- different shock locations and separation patterns were evident. Similarly, the wing/body/nacelle/pylon solutions did not appear to be converging to the same flow solutions. Overall, grid refinement did not consistently improve the correlation with experimental data for either the wingbody or the wing/body/nacelle pylon configuration. Although the absolute values of total drag predicted by two of the solvers for the medium and fine grids did not compare well with the experiment, the incremental drag predictions were within plus or minus 3 counts of the experimental data. The correlation with experimental incremental drag was not significantly changed by specifying transition. Although the sources of code-to-code variation in force and moment predictions for the three unstructured grid codes have not yet been identified, the current study reinforces the necessity of applying multiple codes to the same application to assess uncertainty.

Lee-Rausch, E. M.↗

Materials constitutive models for nonlinear analysis of thermally cycled structures

Effects of inelastic materials models on computed stress-strain solutions for thermally loaded structures were studied by performing nonlinear (elastoplastic creep) and elastic structural analyses on a prismatic, double edge wedge specimen of IN 100 alloy that was subjected to thermal cycling in fluidized beds. Four incremental plasticity creep models (isotropic, kinematic, combined isotropic kinematic, and combined plus transient creep) were exercised for the problem by using the MARC nonlinear, finite element computer program. Maximum total strain ranges computed from the elastic and nonlinear analyses agreed within 5 percent. Mean cyclic stresses, inelastic strain ranges, and inelastic work were significantly affected by the choice of inelastic constitutive model. The computing time per cycle for the nonlinear analyses was more than five times that required for the elastic analysis.

Kaufman, A.↗

Static geometric and material nonlinear analysis.

This paper presents a survey of four aspects of nonlinear analysis - basic formulation, plasticity relations, computational procedures, and methods of solution. This survey is for the most part limited to static geometric and material nonlinear analysis under the assumption of small strains. The two formulations covered are the total Lagrangian and the incremental moving coordinate formulations. A comparison of the two formulations is presented and an attempt is made to clarify the underlying basis of each formulation. The survey of computational procedures reveals that most researchers are now using some form of numerical integration to evaluate nonlinear contributions.

Stricklin, J. A.↗

Experimental and theoretical analysis of carbon driven detonation waves in a heterogeneously premixed Rotating Detonation Engine

Coal dust explosions can be hazardous; however, they can also generate a significant rise in stagnation pressure if adequately harnessed. Rotating detonation combustors seek to take advantage of the stagnation pressure rise phenomenon in a more sustained and controlled manner via confinement to a physical annulus, leading to increased overall thermodynamic efficiency. Here this investigation presents an analysis of detonations fueled by Carbon Black, a solid particulate consisting of virtually pure carbon molecules and lean Hydrogen-Air mixtures. It is realized that with the addition of Carbon Black, an increase of lean mixture detonability and detonation velocities extending the operating limit over that of a pure hydrogen-air mixture is experienced. For all testing conditions, the total equivalence ratio is held at φ = 1, while the fuel mixture's carbon mass fraction is increased from 0 to 0.7 while the hydrogen is decreased. Detonation wave velocities are extracted from high-speed imaging through applying a Discrete Fourier Transform algorithm to determine changes to the wave speed as Carbon Black particles are introduced. As a result, due to the addition of Carbon Black as an auxiliary fuel source, detonations were formed instead of deflagrations in operating conditions where one would expect deflagrations at the same hydrogen-air equivalence ratios without Carbon Black addition. The detonation formation provides evidence that the coal particles are reacting within the detonation wave in a large enough capacity to support a detonation wave within the annulus. Furthermore, the wave speed is shown to increase with the additional of carbon particles. At a constant global equivalence ratio, the detonation wave velocities were found to decrease with hydrogen's incremental replacement with coal particles. Whereby, through a theoretical comparison of the heat of combustion as computed from the experimentally derived detonation wave velocities, a linear relationship of the two was shown to exist. Therefore, the heat of combustion has the potential to describe an operational limit to sustaining a detonation wave.

42 ENGINEERING↗

An overview of the value of parabolic dish solar thermal systems in industrial cogeneration applications

The essential elements of the cogeneration system configuration to be captured were the displacement of thermal energy by collection and use of the Brayton exhaust stream, and the sale back to the utility of any electricity production in excess of on-site requirements. In contrast to simply dumping these energy flows, their use or sale obviously serves, by itself, to increase gross value of the solar thermal energy system. Net allowable cost of the parabolic dish modules may or may not be increased, however. A consideration is that the waste heat capture and delivery subsystems are not free. This study does not address the incremental cost of adding waste heat capture, transport, and conversion (to steam, if necessary). It does compute a value for the thermal energy thereby displaced. This value can serve as a first-round input to any detailed economic evaluation of waste heat recovery.

Source record↗

An in-flight investigation of ground effect on a forward-swept wing airplane

A limited flight experiment was conducted to document the ground-effect characteristics of the X-29A research airplane. This vehicle has an aerodynamic platform which includes a forward-swept wing and close-coupled, variable incidence canard. The flight-test program obtained results for errors in the airdata measurement and for incremental normal force and pitching moment caused by ground effect. Correlations with wind-tunnel and computational analyses were made. The results are discussed with respect to the dynamic nature of the flight measurements, similar data from other configurations, and pilot comments. The ground-effect results are necessary to obtain an accurate interpretation of the vehicle's landing characteristics. The flight data can also be used in the development of many modern aircraft systems such as autoland and piloted simulations.

Curry, Robert E.↗

An in-flight investigation of ground effect on a forward-swept wing airplane

A limited flight experiment was conducted to document the ground effect characteristics of the X-29A research aircraft. This vehicle has an aerodynamic platform which includes a forward-swept wing and close-coupled, variable incidence canard. The flight-test program obtained results for errors in the air data measurement and for incremental normal force and pitching moment caused by ground effect. Correlations with wind-tunnel and computational analyses were made. The results are discussed with respect to the dynamic nature of the flight measurements, similar data from other configurations, and pilot comments. The ground effect results are necessary to obtain an accurate interpretation of the vehicle's landing characteristics. The flight data can also be used in the development of many modern aircraft systems such as autoland and piloted simulation.

Curry, Robert E.↗

Validation of Long Bone Mechanical Properties from Densitometry

The objective of this study was to assess whether cross-sectional areal properties, calculated from densitometry, correlate to the true flexural properties. Right and left male embalmed tibiae were used in the study. Prior to scanning, the proximal end of each tibia was potted in a fixture with registration pins, flushed thoroughly with water under pressure to remove trapped air, and then placed in a constant thickness water bath attached to a precision indexer. Two sets of three scans of the entire tibia were taken with an Hologic QDR 1000/W densitometer at rotations of 0, 45, and 90 degrees about the tibia long axis. An aluminum step phantom and a bone step phantom, machined from bovine cortical bone, were also in the bath and scanned separately. Pixel attenuation data from the two sets of scans were averaged to reduce noise. Pixel data from the high energy beam were then converted to equivalent thicknesses using calibration equations. Cross-sectional areal properties (centroid, principal area moments and principal angle) along the length were computed from the three registered scans using methods developed in our laboratory. Flexural rigidities. Four strain gages were bonded around the circumference of each of 5 cross-sections encompassing the entire diaphysis. A known transverse load was then applied to the distal end and the bone was rotated 360 degrees in eight increments of 45 degrees each. Strains from the eight orientations were analyzed along with the known applied bending moments at each section to compute section centroids, curvatures, principal flexural rigidities and principal angle. Reference axes between the two methods were maintained within +/- 0.5 degrees using an electronic inclinometer. Principal angles (flexural - areal) differed by -2.0 +/- 4.0 degrees, and 1.0 +/- 2.5 degrees for the right and left tibia, respectively. Section principal flexural rigidities were highly correlated to principal areal moments (right: r(sup 2)= 0.997; left: r(sup 2)= 0.978) indicating a nearly constant effective flexural modulus. Right and left tibia exhibited a very high degree of symmetry when comparing either flexural or areal properties. To our knowledge this is the first study to validate the use of densitometry (DXA) to predict three dimensional structural properties of long bones. Our initial results support the conclusion that bone mineral and its distribution are the primary determinants of flexural modulus and rigidity.

Whalen, R.↗

Methodology for Sensitivity Analysis, Approximate Analysis, and Design Optimization in CFD for Multidisciplinary Applications

An incremental iterative formulation together with the well-known spatially split approximate-factorization algorithm, is presented for solving the large, sparse systems of linear equations that are associated with aerodynamic sensitivity analysis. This formulation is also known as the 'delta' or 'correction' form. For the smaller two dimensional problems, a direct method can be applied to solve these linear equations in either the standard or the incremental form, in which case the two are equivalent. However, iterative methods are needed for larger two-dimensional and three dimensional applications because direct methods require more computer memory than is currently available. Iterative methods for solving these equations in the standard form are generally unsatisfactory due to an ill-conditioned coefficient matrix; this problem is overcome when these equations are cast in the incremental form. The methodology is successfully implemented and tested using an upwind cell-centered finite-volume formulation applied in two dimensions to the thin-layer Navier-Stokes equations for external flow over an airfoil. In three dimensions this methodology is demonstrated with a marching-solution algorithm for the Euler equations to calculate supersonic flow over the High-Speed Civil Transport configuration (HSCT 24E). The sensitivity derivatives obtained with the incremental iterative method from a marching Euler code are used in a design-improvement study of the HSCT configuration that involves thickness. camber, and planform design variables.

Taylor, Arthur C., III↗

Methodology for sensitivity analysis, approximate analysis, and design optimization in CFD for multidisciplinary applications

In this study involving advanced fluid flow codes, an incremental iterative formulation (also known as the delta or correction form) together with the well-known spatially-split approximate factorization algorithm, is presented for solving the very large sparse systems of linear equations which are associated with aerodynamic sensitivity analysis. For smaller 2D problems, a direct method can be applied to solve these linear equations in either the standard or the incremental form, in which case the two are equivalent. Iterative methods are needed for larger 2D and future 3D applications, however, because direct methods require much more computer memory than is currently available. Iterative methods for solving these equations in the standard form are generally unsatisfactory due to an ill-conditioning of the coefficient matrix; this problem can be overcome when these equations are cast in the incremental form. These and other benefits are discussed. The methodology is successfully implemented and tested in 2D using an upwind, cell-centered, finite volume formulation applied to the thin-layer Navier-Stokes equations. Results are presented for two sample airfoil problems: (1) subsonic low Reynolds number laminar flow; and (2) transonic high Reynolds number turbulent flow.

Taylor, Arthur C., III↗

Disentangling the physics of the attractive Hubbard model as a fully interacting model of fermions via the accessible and symmetry-resolved entanglement entropies

The complicated ways in which electrons interact in many-body systems such as molecules and materials have long been viewed through the lens of local electron correlation and associated correlation functions. However, quantum information science has demonstrated that more global diagnostics of quantum states like the entanglement entropy can provide a complementary and clarifying lens on electronic behavior. One particularly useful measure that can be used to distinguish between quantum and classical sources of entanglement is the accessible entanglement, the entanglement available as a quantum resource for systems subject to conservation laws, such as fixed particle number, due to superselection rules. In this work, we introduce an algorithm and demonstrate how to compute accessible and symmetry-resolved entanglements for interacting fermion systems. This is accomplished by combining an incremental version of the swap algorithm with a recursive auxiliary field quantum Monte Carlo algorithm recently developed by the authors. We apply these tools to study the pairing and charge density waves exhibited in the paradigmatic attractive Hubbard model via entanglement. We find that the particle and spin symmetry-resolved entanglements and their related full probability distribution functions show very clear—and unique—signatures of the underlying electronic behavior even when those features are less pronounced in conventional correlation functions. Altogether, this work provides a systematic means of characterizing the entanglement within quantum systems that can grant a deeper understanding of the complicated electronic behavior that underlies quantum phase transitions and crossovers in many-body systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Experimental aerodynamic coefficients on a Shuttle-like vehicle at Mach 6 and 10 and comparison to prediction

Aerodynamic coefficients measured for an Orbiter-like configuration are presented for a range of angle of attack from 15 deg to 45 deg at Mach 6 and 10. The low-gamma aspect of a real gas on aerodynamic characteristics and shock detachment distance were simulated by testing this configuration in Mach 6 air (freestream gamma equal to 1.4) and Mach 6 CF4 (freestream gamma equal to 1.17). The effects of Mach number, Reynolds number, and gamma on these measurements are examined and comparisons made to an inviscid flowfield computer code known as HALIS. Pitching-moment coefficients measured in CF4 revealed a nose-up increment when compared with measurements in perfect air, indicative of real-gas effects. In general, the HALIS code accurately predicted the measured aerodynamic coefficients in air and CF4.

Micol, J. R.↗

HOLEGAGE 1.0 - Strain-Gauge Drilling Analysis Program

Interior stresses inferred from changes in surface strains as hole is drilled. Computes stresses using strain data from each drilled-hole depth layer. Planar stresses computed in three ways: least-squares fit for linear variation with depth, integral method to give incremental stress data for each layer, and/or linear fit to integral data. Written in FORTRAN 77.

Hampton, Roy V.↗

The Application of Modeling and Simulation to the Behavioral Deficit of Autism

This abstract describes a research effort to apply technological advances in virtual reality simulation and computer-based games to create behavioral modification programs for individuals with Autism Spectrum Disorder (ASD). The research investigates virtual social skills training within a 3D game environment to diminish the impact of ASD social impairments and to increase learning capacity for optimal intellectual capability. Individuals with autism will encounter prototypical social contexts via computer interface and will interact with 3D avatars with predefined roles within a game-like environment. Incremental learning objectives will combine to form a collaborative social environment. A secondary goal of the effort is to begin the research and development of virtual reality exercises aimed at triggering the release of neurotransmitters to promote critical aspects of synaptic maturation at an early age to change the course of the disease.

Anton, John J.↗

Reduction of computer usage costs in predicting unsteady aerodynamic loadings caused by control surface motions: Analysis and results

Results of theoretical and numerical investigations conducted to develop economical computing procedures were applied to an existing computer program that predicts unsteady aerodynamic loadings caused by leading and trailing edge control surface motions in subsonic compressible flow. Large reductions in computing costs were achieved by removing the spanwise singularity of the downwash integrand and evaluating its effect separately in closed form. Additional reductions were obtained by modifying the incremental pressure term that account for downwash singularities at control surface edges. Accuracy of theoretical predictions of unsteady loading at high reduced frequencies was increased by applying new pressure expressions that exactly satisified the high frequency boundary conditions of an oscillating control surface. Comparative computer result indicated that the revised procedures provide more accurate predictions of unsteady loadings as well as providing reduction of 50 to 80 percent in computer usage costs.

Rowe, W. S.↗

Computational methods for efficient structural reliability and reliability sensitivity analysis

This paper presents recent developments in efficient structural reliability analysis methods. The paper proposes an efficient, adaptive importance sampling (AIS) method that can be used to compute reliability and reliability sensitivities. The AIS approach uses a sampling density that is proportional to the joint PDF of the random variables. Starting from an initial approximate failure domain, sampling proceeds adaptively and incrementally with the goal of reaching a sampling domain that is slightly greater than the failure domain to minimize over-sampling in the safe region. Several reliability sensitivity coefficients are proposed that can be computed directly and easily from the above AIS-based failure points. These probability sensitivities can be used for identifying key random variables and for adjusting design to achieve reliability-based objectives. The proposed AIS methodology is demonstrated using a turbine blade reliability analysis problem.

Wu, Y.-T.↗