Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Convergence analysis”

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 1,423 records · Page 79

Mechanical verification of a schematic Byzantine clock synchronization algorithm

Schneider generalizes a number of protocols for Byzantine fault tolerant clock synchronization and presents a uniform proof for their correctness. The authors present a machine checked proof of this schematic protocol that revises some of the details in Schneider's original analysis. The verification was carried out with the EHDM system developed at the SRI Computer Science Laboratory. The mechanically checked proofs include the verification that the egocentric mean function used in Lamport and Melliar-Smith's Interactive Convergence Algorithm satisfies the requirements of Schneider's protocol.

Shankar, Natarajan↗

Numerical study of shock-wave/boundary layer interactions in premixed hydrogen-air hypersonic flows

A computational study of shock wave/boundary layer interactions involving premixed combustible gases, and the resulting combustion processes is presented. The analysis is carried out using a new fully implicit, total variation diminishing (TVD) code developed for solving the fully coupled Reynolds-averaged Navier-Stokes equations and species continuity equations in an efficient manner. To accelerate the convergence of the basic iterative procedure, this code is combined with vector extrapolation methods. The chemical nonequilibrium processes are simulated by means of a finite-rate chemistry model for hydrogen-air combustion. Several validation test cases are presented and the results compared with experimental data or with other computational results. The code is then applied to study shock wave/boundary layer interactions in a ram accelerator configuration. Results indicate a new combustion mechanism in which a shock wave induces combustion in the boundary layer, which then propagates outwards and downstream. At higher Mach numbers, spontaneous ignition in part of the boundary layer is observed, which eventually extends along the entire boundary layer at still higher values of the Mach number.

Yungster, Shaye↗

Constrained optimization using design of experiment surfaces

An algorithm for solving constrained optimization problems is presented. First, design of experiment techniques are used to survey the design space. After evaluating the objective and constraint functions, as specified by Taguchi orthogonal arrays, analytical models of these functions are generated using a least-squares regression analysis. Next, a nonlinear programming package is used to optimize the analytical model. Based on the optimization information, the design space is reduced so as to close in around the minimum, and the entire procedure is repeated until convergence. An important feature of the algorithm is that function gradients are not required; therefore, for problems in which gradients would have to be estimated using finite-differences the number of function evaluations required for the optimization is significantly reduced, when compared with traditional nonlinear programming techniques. In addition, there is no requirement that the gradients must be smooth and continuous.

Bolt, Marvin Vance↗

Modeling of Rolling Element Bearing Mechanics: Computer Program Updates

The Rolling Element Bearing Analysis System (REBANS) extends the capability available with traditional quasi-static bearing analysis programs by including the effects of bearing race and support flexibility. This tool was developed under contract for NASA-MSFC. The initial version delivered at the close of the contract contained several errors and exhibited numerous convergence difficulties. The program has been modified in-house at MSFC to correct the errors and greatly improve the convergence. The modifications consist of significant changes in the problem formulation and nonlinear convergence procedures. The original approach utilized sequential convergence for nested loops to achieve final convergence. This approach proved to be seriously deficient in robustness. Convergence was more the exception than the rule. The approach was changed to iterate all variables simultaneously. This approach has the advantage of using knowledge of the effect of each variable on each other variable (via the system Jacobian) when determining the incremental changes. This method has proved to be quite robust in its convergence. This technical memorandum documents the changes required for the original Theoretical Manual and User's Manual due to the new approach.

Ryan, S. G.↗

Uncertainty Analysis of Decomposing Polyurethane Foam

Sensitivity/uncertainty analyses are necessary to determine where to allocate resources for improved predictions in support of our nation's nuclear safety mission. Yet, sensitivity/uncertainty analyses are not commonly performed on complex combustion models because the calculations are time consuming, CPU intensive, nontrivial exercises that can lead to deceptive results. To illustrate these ideas, a variety of sensitivity/uncertainty analyses were used to determine the uncertainty associated with thermal decomposition of polyurethane foam exposed to high radiative flux boundary conditions. The polyurethane used in this study is a rigid closed-cell foam used as an encapsulant. Related polyurethane binders such as Estane are used in many energetic materials of interest to the JANNAF community. The complex, finite element foam decomposition model used in this study has 25 input parameters that include chemistry, polymer structure, and thermophysical properties. The response variable was selected as the steady-state decomposition front velocity calculated as the derivative of the decomposition front location versus time. An analytical mean value sensitivity/uncertainty (MV) analysis was used to determine the standard deviation by taking numerical derivatives of the response variable with respect to each of the 25 input parameters. Since the response variable is also a derivative, the standard deviation was essentially determined from a second derivative that was extremely sensitive to numerical noise. To minimize the numerical noise, 50-micrometer element dimensions and approximately 1-msec time steps were required to obtain stable uncertainty results. As an alternative method to determine the uncertainty and sensitivity in the decomposition front velocity, surrogate response surfaces were generated for use with a constrained Latin Hypercube Sampling (LHS) technique. Two surrogate response surfaces were investigated: 1) a linear surrogate response surface (LIN) and 2) a quadratic response surface (QUAD). The LHS techniques do not require derivatives of the response variable and are subsequently relatively insensitive to numerical noise. To compare the LIN and QUAD methods to the MV method, a direct LHS analysis (DLHS) was performed using the full grid and timestep resolved finite element model. The surrogate response models (LIN and QUAD) are shown to give acceptable values of the mean and standard deviation when compared to the fully converged DLHS model.

Hobbs, Michael L.↗

Modeling and Analysis of the Reverse Water Gas Shift Process for In-Situ Propellant Production

This report focuses on the development of mathematical models and simulation tools developed for the Reverse Water Gas Shift (RWGS) process. This process is a candidate technology for oxygen production on Mars under the In-Situ Propellant Production (ISPP) project. An analysis of the RWGS process was performed using a material balance for the system. The material balance is very complex due to the downstream separations and subsequent recycle inherent with the process. A numerical simulation was developed for the RWGS process to provide a tool for analysis and optimization of experimental hardware, which will be constructed later this year at Kennedy Space Center (KSC). Attempts to solve the material balance for the system, which can be defined by 27 nonlinear equations, initially failed. A convergence scheme was developed which led to successful solution of the material balance, however the simplified equations used for the gas separation membrane were found insufficient. Additional more rigorous models were successfully developed and solved for the membrane separation. Sample results from these models are included in this report, with recommendations for experimental work needed for model validation.

Whitlow, Jonathan E.↗

Vestibular adaptation to space in monkeys

Otolith-induced eye movements of rhesus monkeys were studied before and after the 1989 COSMOS 2044 and the 1992 to 1993 COSMOS 2229 flights. Two animals flew in each mission for approximately 2 weeks. After flight, spatial orientation of the angular vestibulo-ocular reflex was altered. In one animal the time constant of postrotatory nystagmus, which had been shortened by head tilts with regard to gravity before flight, was unaffected by the same head tilts after flight. In another animal, eye velocity, which tended to align with a gravitational axis before flight, moved toward a body axis after flight. This shift of orientation disappeared by 7 days after landing. After flight, the magnitude of compensatory ocular counter-rolling was reduced by about 70% in both dynamic and static tilts. Modulation in vergence in response to naso-occipital linear acceleration during off-vertical axis rotation was reduced by more than 50%. These changes persisted for 11 days after recovery. An up and down asymmetry of vertical nystagmus was diminished for 7 days. Gains of the semicircular canal-induced horizontal and vertical angular vestibulo-ocular reflexes were unaffected in both flights, but the gain of the roll angular vestibulo-ocular reflex was decreased. These data indicate that there are short- and long-term changes in otolith-induced eye movements after adaptation to microgravity. These experiments also demonstrate the unique value of the monkey as a model for studying effects of vestibular adaptation in space. Eye movements can be measured in three dimensions in response to controlled vestibular and visual stimulation, and the results are directly applicable to human beings. Studies in monkeys to determine how otolith afferent input and central processing is altered by adaptation to microgravity should be an essential component of future space-related research.

Cosmos 2229 Project↗

Analytical and Numerical Results for an Adhesively Bonded Joint Subjected to Pure Bending

A one-dimensional, semi-analytical methodology that was previously developed for evaluating adhesively bonded joints composed of anisotropic adherends and adhesives that exhibit inelastic material behavior is further verified in the present paper. A summary of the first-order differential equations and applied joint loading used to determine the adhesive response from the methodology are also presented. The method was previously verified against a variety of single-lap joint configurations from the literature that subjected the joints to cases of axial tension and pure bending. Using the same joint configuration and applied bending load presented in a study by Yang, the finite element analysis software ABAQUS was used to further verify the semi-analytical method. Linear static ABAQUS results are presented for two models, one with a coarse and one with a fine element meshing, that were used to verify convergence of the finite element analyses. Close agreement between the finite element results and the semi-analytical methodology were determined for both the shear and normal stress responses of the adhesive bondline. Thus, the semi-analytical methodology was successfully verified using the ABAQUS finite element software and a single-lap joint configuration subjected to pure bending.

Smeltzer, Stanley S., III↗

Design, Analysis, and On-Sun Evaluation of Reflective Strips Under Controlled Buckling

Solar concentrators are envisioned for use in a variety of space-based applications, including applications involving in situ resource utilization. Identifying solar concentrators that minimize mass and cost are of great interest, especially since launch cost is driven in part by the mass of the payload. Concentrators must also be able to survive the wide temperature excursions on the lunar surface. Identifying smart structures which compensate for changes in concentrator geometry brought about by temperature extremes are of interest. Some applications may benefit from the ability to change the concentrators focal pattern at will. This paper addresses a method of designing a single reflective strip to produce a desired focal pattern through the use of controlled buckling. Small variations in the cross section over the length of the reflective strip influence the distribution of light in the focal region. A finite element method of analysis is utilized here which calculates the curve produced for a given strip cross section and axial load. Varying axial force and strip cross section over the length of the reflective strip provide a means of optimizing ray convergence in the focal region. Careful selection of a tapered cross section yields a reflective strip that approximates a parabola. An array of reflective strips under controlled buckling produces a light weight concentrator and adjustments in the compression of individual strips provide a means of compensating for temperature excursions or changing the focal pattern at will.

solar collectors↗

Design, Analysis, and On-Sun Evaluation of Reflective Strips Under Controlled Buckling

Solar concentrators are envisioned for use in a variety of space-based applications, including applications involving in situ resource utilization. Identifying solar concentrators that minimize mass and cost are of great interest, especially since launch cost is driven in part by the mass of the payload. Concentrators must also be able to survive the wide temperature excursions on the lunar surface. Identifying smart structures which compensate for changes in concentrator geometry brought about by temperature extremes are of interest. Some applications may benefit from the ability to change the concentrator's focal pattern at will. This paper addresses a method of designing a single reflective strip to produce a desired focal pattern through the use of controlled buckling. Small variations in the cross section over the length of the reflective strip influence the distribution of light in the focal region. A finite element method of analysis is utilized here which calculates the curve produced for a given strip cross section and axial load. Varying axial force and strip cross section over the length of the reflective strip provide a means of optimizing ray convergence in the focal region. Careful selection of a tapered cross section yields a reflective strip that approximates a parabola. An array of reflective strips under controlled buckling produces a light weight concentrator and adjustments in the compression of individual strips provide a means of compensating for temperature excursions or changing the focal pattern at will.

buckling↗

Analysis of the Capacity Potential of Current Day and Novel Configurations for New York's John F. Kennedy Airport

In 2015, a series of systems analysis studies were conducted on John F. Kennedy Airport in New York (NY) in a collaborative effort between NASA and the Port Authority of New York and New Jersey (PANYNJ). This work was performed to build a deeper understanding of NY airspace and operations to determine the improvements possible through operational changes with tools currently available, and where new technology is required for additional improvement. The analysis was conducted using tool-based mathematical analyses, video inspection and evaluation using recorded arrival/departure/surface traffic captured by the Aerobahn tool (used by Kennedy Airport for surface metering), and aural data archives available publically through the web to inform the video segments. A discussion of impacts of trajectory and operational choices on capacity is presented, including runway configuration and usage (parallel, converging, crossing, shared, independent, staggered), arrival and departure route characteristics (fix sharing, merges, splits), and how compression of traffic is staged. The authorization in March of 2015 for New York to use reduced spacing under the Federal Aviation Administration (FAA) Wake Turbulence Recategorization (RECAT) also offers significant capacity benefit for New York airports when fully transitioned to the new spacing requirements, and the impact of these changes for New York is discussed. Arrival and departure capacity results are presented for each of the current day Kennedy Airport configurations. While the tools allow many variations of user-selected conditions, the analysis for these studies used arrival-priority, no-winds, additional safety buffer of 5% to the required minimum spacing, and a mix of traffic typical for Kennedy. Two additional "novel" configurations were evaluated. These configurations are of interest to Port Authority and to their airline customers, and are believed to offer near-term capacity benefit with minimal operational and equipage changes. One of these is the addition of an Optimized Profile Descent (OPD) route to runways 22L and 22R, and the other is the simultaneous use of 4 runways, which is not currently done at Kennedy. The background and configuration for each of these is described, and the capacity results are presented along with a discussion of drawbacks and enablers for each.

Glaab, Patricia↗

Order-by-order uncertainties of nucleon-nucleon Wolfenstein amplitudes in chiral effective field theory

Quantum mechanical invariance principles dictate the most general operator structure that can be present in the nucleon-nucleon (NN) interaction. Five independent operators appear in the on-shell NN amplitude together with five corresponding coefficient functions. The usual choice for these coefficient functions is known as the NN Wolfenstein amplitudes. We analyze the order-by-order convergence of each of the five NN Wolfenstein amplitudes predicted by a semilocal coordinate space potential implementation of chiral effective field theory (𝜒⁢EFT). We do this at laboratory kinetic energies between 25 and 200 MeV for both neutron-proton and proton-proton scattering. Our analysis uses the Gaussian-process methods developed by the BUQEYE Collaboration to describe the contributions of each 𝜒⁢EFT order, and so yields truncation uncertainties for each Wolfenstein amplitude that are correlated across scattering angles. We combine information on the size of different orders in the EFT to infer the 𝜒⁢EFT breakdown scale for each amplitude, finding, on average, Λ 𝑏 between 750 and 800 MeV. Furthermore, with this choice of Λ 𝑏 , the EFT truncation uncertainties cover both higher-order results and empirical Wolfenstein amplitudes well for all orders other than the leading order.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Analysis and design of steady transonic flow over airfoils by the method of parametric differentiation

The results reported here are based on applying the method of parametric differentiation (MPD) to transform the nonlinear differential equation governing small-disturbance transonic flow to a linear equation. Implicit approximate factorization and monotone methods were used to accelerate convergence of the linear problem by an order of magnitude over successive line over-relaxation. The relative merits of using MPD are discussed in comparison to conventional small-disturbance applications. Several MPD analyses are performed on an array of airfoils. A design procedure utilizing MPD is discussed and demonstrated for two nonlifting cases.

Halt, D. W.↗

Chemolithoautotroph distributions across the subsurface of a convergent margin

Subducting oceanic crusts release fluids rich in biologically relevant compounds into the overriding plate, fueling subsurface chemolithoautotrophic ecosystems. To understand the impact of subsurface geochemistry on microbial communities, we collected fluid and sediments from 14 natural springs across a ~200 km transect across the Costa Rican convergent margin and performed shotgun metagenomics. The resulting 404 metagenome-assembled genomes (MAGs) cluster into geologically distinct regions based on MAG abundance patterns: outer forearc-only (25% of total relative abundance), forearc/arc-only (38% of total relative abundance), and delocalized (37% of total relative abundance) clusters. In the outer forearc, Thermodesulfovibrionia, Candidatus Bipolaricaulia, and Firmicutes have hydrogenotrophic sulfate reduction and Wood-Ljungdahl (WL) carbon fixation pathways. In the forearc/arc, Anaerolineae, Ca. Bipolaricaulia, and Thermodesulfovibrionia have sulfur oxidation, nitrogen cycling, microaerophilic respiration, and WL, while Aquificae have aerobic sulfur oxidation and reverse tricarboxylic acid carbon fixation pathway. Transformation-based canonical correspondence analysis shows that MAG distribution corresponds to concentrations of aluminum, iron, nickel, dissolved inorganic carbon, and phosphate. While delocalized MAGs appear surface-derived, the subsurface chemolithoautotrophic, metabolic, and taxonomic landscape varies by the availability of minerals/metals and volcanically derived inorganic carbon. Furthermore, the WL pathway persists across all samples, suggesting that this versatile, energy-efficient carbon fixation pathway helps shape convergent margin subsurface ecosystems.

59 BASIC BIOLOGICAL SCIENCES↗

Free-Surface Liquid Lithium Flow Modeling and Stability Analysis for Fusion Applications

Liquid metal plasma facing components are considered an attractive design choice for fusion devices including pilot plants. Virtual prototyping of such devices includes modeling of free-surface flow of the electrically conductive liquid, which requires computational fluid dynamics (CFD) and magnetohydrodynamics (MHD) simulations. Numerical tools capable of simulating flows and heat transfer in the free-surface MHD flow were developed at PPPL based on the customized ANSYS CFX. Here, MHD is introduced using a magnetic vector potential approach. Free-surface flow capabilities are available in the code and were tested. Special stabilization procedures were derived and applied to improve convergence of the momentum equations with the source terms due to the Lorentz force and surface tension. Important characteristics of the fusion-relevant liquid metal flow is free surface smoothness and stability. Heat flux from the plasma impacts the liquid surface at a very acute angle, so any change of the free surface from axisymmetry can dramatically increase the local heat flux density and thus create excessive evaporation of liquid lithium into the plasma, which is detrimental to operations. Stability analysis of the liquid metal film flow was performed to determine applicable flow regimes. Thin film flow along horizontal wall is considered. Effects of gravity, magnetic field, and surface tension are included in the analysis.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Comparison of Photovoltaic Module Luminescence Imaging Techniques: Assessing the Influence of Lateral Currents in High-Efficiency Device Structures

We investigate the implications of using partial or patterned illumination for luminescence imaging of photovoltaic modules. Partial illumination induces local photovoltage variations that drive lateral current flow into non-illuminated cell regions, causing the average injection level to vary over the course of luminescence measurement. The extent of lateral current depends on material properties and is especially significant for high-efficiency devices such as the silicon heterojunction with intrinsic thin layer (HIT) modules that we investigate. We compare four module imaging methods including large-area illumination, scanning-laser illumination, contactless electroluminescence, and pattern-illuminated photoluminescence. By evaluating the statistical correlation among the different methods under varying laser fluences, we conclude that most methods converge to similar results at matched photon dose per cell (not laser power density) for the modules that we study, while pattern-illuminated photoluminescence yields distinct images from the other techniques regardless of the illumination power. Our combined analysis using multiple luminescence imaging techniques gives specific insight into the causes of HIT module degradation, including suspected degradation of the transparent conductive oxide and passivation layer.

14 SOLAR ENERGY↗

The many-body expansion for aqueous systems revisited: III. Hofmeister ion – water interactions

We report a Many Body Energy (MBE) analysis of aqueous ionic clusters containing anions and cations at the two opposite ends of the Hofmeister series, viz. the kosmotropes Ca2+, SO42- and chaotropes NH4+ and ClO4- with 9 water molecules to quantify the how these ions in altering the interaction between the water molecules in their immediate surrounding. The current results are contrasted to the ones reported earlier for water clusters as well as for alkali metal and halide ion aqueous clusters of the same size, which lie in the middle of the Hofmeister series. Through this analysis, noteworthy differences between the MBE of kosmotropes and chaotropes were identified. The MBE of kosmotropes is dominated by ion-water interactions that extends beyond the 4-body term, the point at which the MBE of pure water converges. The percentage contribution of the 2- B to the total cluster binding energy is noticeably larger. The disruption due to the dominant ion results in weak, unfavorable water-water interactions. The MBE for chaotropes, on the other hand, was found to converge more quickly as it more closely resembles that of pure water clusters. Chaotropes exhibit weaker overall binding energies and ion-water interactions with more favorable water-water interactions, somewhat recovering the pattern of the 2-4 body terms exemplified by pure water clusters. More importantly, both kosmotropic and chaotropic ions exhibit an anticorrelation between the 2-B ion-water (I-W) and water-water (W-W) interactions as well as between the 3-B (I-W-W) and (I-W) interactions. The consideration of two different structural arrangements (ion inside and outside of a water cluster) suggests that fully solvated (ion inside) chaotropes disrupt the hydrogen bonding network in a similar manner as partially solvated (ion outside) kosmotropes and offer useful insights into the modeling requirements of bulk vs. an interface. Finally, the 2-B contribution to the total Basis Set Superposition Error (BSSE) correction for the kosmotropic and chaotropic ions follows the previously reported erf profile vs. intermolecular distance. When scaled for the corresponding dimer energies and distances, a single profile fits the current results together with all previously reported ones for the pure water and halide water clusters.

Herman, Kristina M.↗