Search NASA⌕ Search

SEARCH · Search NASA

Results for “iterative 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 487 records · Page 27

Methodology to determine printability criteria of highly concentrated pastes through rheological characterization

Material extrusion is an additive manufacturing technique that enables the creation of reproducible and complex hardware by depositing a viscous, shear-thinning ink onto a substrate in a custom-pattern via extrusion through a syringe. Here, the ability of an ink to be extruded onto a substrate in many layers, and maintain the desired shape is what defines the printability. Printability is often investigated by formulating, printing, and postmortem analysis of final parts in an iterative manner. Investigations of printability through rheological characterization have often been concerned with inks that straddle the line between printable and too thin, leaving out an entire class of inks that are highly-filled pastes, where extrudability is the limiting factor. Highly-filled pastes continue to pose issues for researchers as the effect of filler morphology, size, loading, and packing fraction on the ink rheology and corresponding printability is not understood. While traditional rheological characterizations may be useful for some inks, we show that protocols utilizing steady-shear, or large-amplitude oscillatory shear are difficult and unreliable for highly-filled pastes. Through transient rheology paired with real-time images we show that each traditional protocol produces inhomogeneous deformations that violate the assumptions that underly common rheological definitions. Instead, we demonstrate metrics measured with small-amplitude oscillatory shear that are correlated to the printability of various ink formulations ranging in loading. The rheological measures that accurately predict the printability of the inks are the axial stress measured at small amplitudes, and the critical stress amplitude above which rheological characterizations become impossible. In addition, we estimate the maximum packing fraction for each filler, based on the exponent common to hard sphere models, and show that the printability of each ink can be predicted by the ratio of the packing fraction to the theoretical maximum. We show how small-amplitude oscillatory shear allows users to develop printability criteria for any ink to enhance the workflow in the development of new inks, increase the performance of material extrusion printing, and improve the stability of printed parts, with less wasted time and materials.

36 MATERIALS SCIENCE↗

TOUGH-UDEC: A simulator for coupled multiphase fluid flows, heat transfers and discontinuous deformations in fractured porous media

A numerical simulator entitled TOUGH-UDEC is presented for the assessment of coupled thermal-hydraulic-mechanical processes in fractured porous media. Two existing well-established codes, TOUGH2 and UDEC, are coupled to model multiphase fluid flows, heat transfers, and discontinuous deformations in fractured porous media by means of discrete fracture representation. TOUGH2 is widely used for the modeling of heat transfers and multiphase multicomponent fluid flows, and UDEC is a well-known distinct element code for rock mechanics. The two codes are solved sequentially, with coupling parameters passed to each equation at specific intervals. After solving thermal-hydraulic equations within the TOUGH2 code, pressure and temperature information is imported into the UDEC model. Upon solving the mechanical equation within the UDEC code the calculated fracture apertures are converted to the equivalent permeability and porosity values for a TOUGH2 flow analysis. The solution is calculated by iteratively following an explicit sequence for numerical efficiency. Verifications are presented to demonstrate the capabilities of the coupled TOUGH-UDEC simulator. Three application examples of (1) shear dilation due to increased pore pressure, (2) thermal stress and (3) CO 2 injection, show that the new simulator can be an effective tool for geoengineering applications involving shear activation of fractures and faults.

54 ENVIRONMENTAL SCIENCES↗

Study of plasma start-up initiated by second harmonic electron cyclotron resonance heating on WEGA experiment

Although both 1st harmonic ordinary mode (O1) and 2nd harmonic extra-ordinary mode (X2) have been successfully used to initiate pre-ionization and breakdown in many devices, a complete theoretical model is still missing to explain the success of this method. Moreover, some experimental observations are not completely understood, such as what occurs during the delay time between the turn-on of ECRH power and first signals of density or light measurements. Since during this free period the ECRH power has to be absorbed by in-vessel components, it is of prime importance to know what governs this delay time. Recently, dedicated start-up experiments have been performed on WEGA, using a 28 GHz ECRH system in X2-mode. This machine has the interesting capability to be run also as a tokamak allowing comparative experiments between stellarator (ι/2π > 0) and tokamak (ι/2π = 0) configurations. Different scans in heating power, neutral gas pressure, and rotational transform (ι) show clearly that the start-up is a two step process. A first step following the turn-on of the ECRH power during which no measurable electron density (or just above the noise level in some cases), ECE and radiated power is detected. Its duration depends strongly on the level of injected power. The second step corresponds to the gas ionization and plasma expansion phase, with a velocity of density build-up and filling-up of the vessel volume depending mainly on pressure, gas and rotational transform. Moreover, an interesting scenario of ECRH pre-ionization without loop voltage in tokamak configuration by applying a small optimal vertical field is relevant for start-up assistance on future experiments like ITER. The results from this experimental parametric study are useful for the modeling of the start-up assisted by the second harmonic electron cyclotron resonance heating. The aim of this work is to establish predictive scenarios for both ITER and W7-X operation.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Modeling L 2,3 -edge X-ray absorption spectroscopy with linear response exact two-component relativistic time-dependent density functional theory

X-ray absorption spectroscopy (XAS) is a powerful tool that can provide physical insights into element-specific chemical processes and reactivities. Although relativistic time-dependent density functional theory (TDDFT) has been previously applied to model the L-edge region in XAS, there has not been a more comprehensive study of the choices of basis sets and density functional kernels available for variational relativistic excited state methods. In this work, we introduce the implementation of the generalized preconditioned locally harmonic residual algorithm to solve the complex-valued relativistic TDDFT for modeling the L-edge X-ray absorption spectra. Further, we investigate the L 2,3 -edge spectra of a series of molecular complexes using relativistic linear response TDDFT with a hybrid iterative diagonalization algorithm. A systematic error analysis was carried out with a focus on the energetics, intensities, and magnitude of L 2 –L 3 splitting compared to experiments. Additionally, the results from relativistic TDDFT calculations are compared to those computed using other theoretical methods, and the multideterminantal effects on the L-edge XAS were investigated.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Analysis of drift effects on the tokamak power scrape-off width using SOLPS-ITER

SOLPS-ITER, a comprehensive 2D scrape-off layer modeling package, is used to examine the physical mechanisms that set the scrape-off width (${{\lambda}_{q}}$ ) for inter-ELM power exhaust. Guided by Goldston's heuristic drift (HD) model, which shows remarkable quantitative agreement with experimental data, this research examines drift effects on ${{\lambda}_{q}}$ in a DIII-D H-mode magnetic equilibrium. As a numerical expedient, a low target recycling coefficient of 0.9 is used in the simulations, resulting in outer target plasma that is sheath limited instead of conduction limited as in the experiment. Scrape-off layer (SOL) particle diffusivity (D SOL ) is scanned from 1 to 0.1 m2 s –1 . Across this diffusivity range, outer divertor heat flux is dominated by a narrow (~3–4mm when mapped to the outer midplane) electron convection channel associated with thermoelectric current through the SOL from outer to inner divertor. An order-unity up–down ion pressure asymmetry allows net ion drift flux across the separatrix, facilitated by an artificial mechanism that mimics the anomalous electron transport required for overall ambipolarity in the HD model. At ${{D}_{\text{SOL}}}=0.1$ m2 s –1 , the density fall-off length is similar to the electron temperature fall-off length, as predicted by the HD model and as seen experimentally. Furthermore, this research represents a step toward a deeper understanding of the power scrape-off width, and serves as a basis for extending fluid modeling to more experimentally relevant, high-collisionality regimes.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Metallographic autopsies of full-scale ITER prototype cable-in-conduit conductors after full cyclic testing in SULTAN: III. The importance of strand surface roughness in long twist pitch conductors

As part of the ITER conductor qualification process, 3 m long Cable-in-Conduit Conductors (CICCs) were tested at the SULTAN facility under conditions simulating ITER operation so as to establish the current sharing temperature, T cs , as a function of multiple full Lorentz force loading cycles. After a comprehensive evaluation of both the Toroidal Field (TF) and the Central Solenoid (CS) conductors, it was found that T cs degradation was common in long twist pitch TF conductors while short twist pitch CS conductors showed some T cs increase. However, one kind of TF conductors containing superconducting strand fabricated by the Bochvar Institute of Inorganic Materials (VNIINM) avoided T cs degradation despite having long twist pitch. In our earlier metallographic autopsies of long and short twist pitch CS conductors, we observed a substantially greater transverse strand movement under Lorentz force loading for long twist pitch conductors, while short twist pitch conductors had negligible transverse movement. With help from the literature, we concluded that the transverse movement was not the source of T cs degradation but rather an increase of the compressive strain in the Nb 3 Sn filaments possibly induced by longitudinal movement of the wires. Like all TF conductors this TF VNIINM conductor showed large transverse motions under Lorentz force loading, but Tcs actually increased, as in all short twist pitch CS conductors. We here propose that the high surface roughness of the VNIINM strand may be responsible for the suppression of the compressive strain enhancement (characteristic of long twist pitch conductors). Furthermore, it appears that increasing strand surface roughness could improve the performance of long twist pitch CICCs.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

Comparative Convergence Analysis of Nonlinear AMLI-Cycle Multigrid

The purpose of our paper is to provide a comprehensive convergence analysis of the nonlinear algebraic multilevel iteration (AMLI)-cycle multigrid (MG) method for symmetric positive definite problems. We show that the nonlinear AMLI-cycle MG method is uniformly convergent, based on classical assumptions for approximation and smoothing properties. Furthermore, under only the assumption that the smoother is convergent, we show that the nonlinear AMLI-cycle method is always better (or not worse) than the respective V-cycle MG method. Finally, numerical experiments are presented to illustrate the theoretical results.

97 MATHEMATICS AND COMPUTING↗

Conditioning Analysis of Incomplete Cholesky Factorizations with Orthogonal Dropping

The analysis of preconditioners based on incomplete Cholesky factorization in which the neglected (dropped) components are orthogonal to the approximations being kept is presented. General estimate for the condition number of the preconditioned system is given which only depends on the accuracy of individual approximations. The estimate is further improved if, for instance, only the newly computed rows of the factor are modified during each approximation step. In this latter case it is further shown to be sharp. The analysis is illustrated with some existing factorizations in the context of discretized elliptic partial differential equations.

97 MATHEMATICS AND COMPUTING↗

Improved tomographic reconstruction of large-scale real-world data by filter optimization

In advanced tomographic experiments, large detector sizes and large numbers of acquired datasets can make it difficult to process the data in a reasonable time. At the same time, the acquired projections are often limited in some way, for example having a low number of projections or a low signal-to-noise ratio. Direct analytical reconstruction methods are able to produce reconstructions in very little time, even for large-scale data, but the quality of these reconstructions can be insufficient for further analysis in cases with limited data. Iterative reconstruction methods typically produce more accurate reconstructions, but take significantly more time to compute, which limits their usefulness in practice. In this paper, we present the application of the SIRT-FBP method to large-scale real-world tomographic data. The SIRT-FBP method is able to accurately approximate the simultaneous iterative reconstruction technique (SIRT) method by the computationally efficient filtered backprojection (FBP) method, using precomputed experiment-specific filters. We specifically focus on the many implementation details that are important for application on large-scale real-world data, and give solutions to common problems that occur with experimental data. We show that SIRT-FBP filters can be computed in reasonable time, even for large problem sizes, and that precomputed filters can be reused for future experiments. Reconstruction results are given for three different experiments, and are compared with results of popular existing methods. The results show that the SIRT-FBP method is able to accurately approximate iterative reconstructions of experimental data. Furthermore, they show that, in practice, the SIRT-FBP method can produce more accurate reconstructions than standard direct analytical reconstructions with popular filters, without increasing the required computation time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Upper-surface blowing nacelle design study for a swept wing airplane at cruise conditions

A study was made to design two types of overwing nacelles for an existing wing-body at a design condition of Mach = 0.8 and C sub L = 0.2. Internal and external surface contours were developed for nacelles having either a D-shaped nozzle or a high-aspect-ratio nozzle for upper-surface blowing in the powered-lift mode of operation. The goal of the design was the development of external nacelle lines that would minimize high-speed aerodynamic interference effects. Each nacelle type was designed for both two- and four-engine airplanes using an iterative process of aerodynamic potential flow analysis. Incremental nacelle drag estimates were made for flow-through wind tunnel models of each configuration.

Gillette, W. B.↗

A scaling law for energy transfer by inelastic electron-molecule collisions in mixtures

The equation governing the electron energy distribution in the presence of a spatially uniform electric field in a weakly ionized gas was reformulated into an integral equation for the logarithmic slope of the distribution function. For gas mixtures in which the dominant electron energy loss mechanism is by vibrational excitation of the molecules, this equation is suitable for approximate analysis and exact numerical solution by iteration. Superelastic collisions are easily included in this formulation, and do not seriously effect the convergence of the numerical scheme. The approximate analytical results are only qualitatively correct, but suggest appropriate parameters which correlate the exact numerical results very well. The distribution function as well as certain gross properties such as net energy transfer into vibration, mean energy, and drift velocity depend primarily on a single nondimensional parameter involving only E/N and the cross sections.

Bienkowski, G. K.↗

Application of the AMI C sub l sub max prediction method to a number of airfoils

A method for calculating the flow about airfoils up to and beyond the stall is described. It is an iterative procedure between potential flow and boundary layer solutions. The separated region is modeled in the potential flow analysis using free vortex sheets which require an inner iteration to establish their shapes. The free vortex sheet length is an important parameter in the potential flow calculation. Results so far indicate a possible correlation between wake length and airfoil thickness/chord ratio. Calculated and experimental results are compared for a series of airfoils.

Dvorak, F. A.↗

Analytical study of STOL Aircraft in ground effect. Part 1: Nonplanar, nonlinear wing/jet lifting surface method

A nonlinear, nonplanar three dimensional jet flap analysis, applicable to the ground effect problem, is presented. Lifting surface methodology is developed for a wing with arbitrary planform operating in an inviscid and incompressible fluid. The classical, infintely thin jet flap model is employed to simulate power induced effects. An iterative solution procedure is applied within the analysis to successively approximate the jet shape until a converged solution is obtained which closely satisfies jet and wing boundary conditions. Solution characteristics of the method are discussed and example results are presented for unpowered, basic powered and complex powered configurations. Comparisons between predictions of the present method and experimental measurements indicate that the improvement of the jet with the ground plane is important in the analyses of powered lift systems operating in ground proximity. Further development of the method is suggested in the areas of improved solution convergence, more realistic modeling of jet impingement and calculation efficiency enhancements.

Shollenberger, C. A.↗

On the structure of parallelism in a highly concurrent PDE solver

A parallel multigrid algorithm for solving elliptic partial differential equations is developed and evaluated. A V-cycle multigrid method is altered to increase the degree of parallelism. A numerical analysis of the resulting concurrent-iteration multigrid algorithm is performed; its architectural implications are considered; highly parallel systems without shared memory are examined (including mesh-connected arrays, mesh-shuffle-connected systems, permutation networks, and direct VLSI embeddings); and the results of numerical experiments are presented in tables and graphs.

Gannon, D.↗

Computer-Aided Design Of Turbine Blades And Vanes

Quasi-three-dimensional method for determining aerothermodynamic configuration of turbine uses computer-interactive analysis and design and computer-interactive graphics. Design procedure executed rapidly so designer easily repeats it to arrive at best performance, size, structural integrity, and engine life. Sequence of events in aerothermodynamic analysis and design starts with engine-balance equations and ends with boundary-layer analysis and viscous-flow calculations. Analysis-and-design procedure interactive and iterative throughout.

Hsu, Wayne Q.↗