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,027 records · Page 57
All sourcesNASA NTRSDOE OSTI
NASA NTRS1976

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.↗
NASA NTRS1994

Aspects of unstructured grids and finite-volume solvers for the Euler and Navier-Stokes equations

Basic algorithms for unstructured mesh generation and fluid flow calculation are discussed. In particular the following are addressed: preliminaries of graphs and meshes; duality and data structures; basic graph operations important in CFD (Computational Fluid Dynamics); triangulation methods, including Varonoi diagrams and Delaunay triangulation; maximum principle analysis; finite volume schemes for scalar conservation law equations; finite volume schemes for the Euler and Navier-Stokes equations; and convergence acceleration for steady state calculations.

Barth, T. J.↗
NASA NTRS1996

Quasi-Optimal Schwarz Methods for the Conforming Spectral Element Discretization

Fast methods are proposed for solving the system K(sub N)x = b resulting from the discretization of self-adjoint elliptic equations in three dimensional domains by the spectral element method. The domain is decomposed into hexahedral elements, and in each of these elements the discretization space is formed by polynomials of degree N in each variable. Gauss-Lobatto-Legendre (GLL) quadrature rules replace the integrals in the Galerkin formulation. This system is solved by the preconditioned conjugate gradients method. The conforming finite element space on the GLL mesh consisting of piecewise Q(sub 1) elements produces a stiffness matrix K(sub h) that is spectrally equivalent to the spectral element stiffness matrix K(sub N). The action of the inverse of K(sub h) is expensive for large problems, and is therefore replaced by a Schwarz preconditioner B(sub h) of this finite element stiffness matrix. The preconditioned operator then becomes B(sub h)(exp -l)K(sub N). The technical difficulties stem from the nonregularity of the mesh. Tools to estimate the convergence of a large class of new iterative substructuring and overlapping Schwarz preconditioners are developed. This technique also provides a new analysis for an iterative substructuring method proposed by Pavarino and Widlund for the spectral element discretization.

Casarin, Mario↗
NASA NTRS1999

Realistic NLTE Radiative Transfer for Modeling Stellar Winds

This NASA grant supported the development of codes to solve the non-LTE multi-level spherical radiative transfer problem in the presence of velocity fields. Much of this work was done in collaboration with Graham Harper (CASA, University of Colorado). These codes were developed for application to the cool, low-velocity winds of evolved late-type stars. Particular emphasis was placed on modeling the wind of lambda Velorum (K4 lb), the brightest K supergiant in the sky, based on extensive observations of the ultraviolet spectrum with the HST/GHRS from GO program 5307. Several solution techniques were examined, including the Eddington factor Approach described in detail by Bennett & Harper (1997). An Eddington factor variant of Harper's S-MULTI code (Harper 1994) for stationary atmospheres was developed and implemented, although full convergence was not realized. The ratio of wind terminal velocity to turbulent velocity is large (approx. 0.3-0.5) in these cool star winds so this assumption of stationarity provides reasonable starting models. Final models, incorporating specified wind laws, were converged using the comoving CRD S-MULTI code. Details of the solution procedure were published by Bennett & Harper (1997). Our analysis of the wind of lambda Vel, based on wind absorption superimposed on chromospheric emission lines in the ultraviolet, can be found in Carpenter et al. (1999). In this paper, we compare observed wind absorption features to an exact CRD calculation in the comoving frame, and also to a much quicker, but approximate, method using the SEI (Sobolev with Exact Integration) code of Lamers, Cerruti-Sola, & Perinotto (1987). Carpenter et al. (1999) provide detailed comparisons of the exact CRD and approximate SEI results and discuss when SEI is adequate to use for computing wind line profiles. Unfortunately, the observational material is insufficient to unambiguously determine the wind acceleration law for lambda Vel. Relatively few unblended Fe II lines of optical depth sensitive to the wind acceleration region are present in the existing HST/GHRS data set. Most of the Fe II lines are either too optically thick (resulting in a board, black wind absorption profile) or too optically thin (resulting in no wind absorption feature present). Also, most of the ultraviolet spectra obtained from HST GO-5307 was at medium resolution (R approx. 40,000, corresponding to a velocity resolution of 7.5 km/s). This is simply inadequate to resolve the turbulence in the outer wind; a key parameter in theoretical wind models. We can now say that an unambiguous determination of the wind velocity law in lambda Vel will require complete coverage of the ultraviolet spectrum at high dispersion (R approx. 10(exp 5), or 3 km/s). This is now feasible usin, the STIS echelle spectrograph on-board HST.

Bennett, Philip D.↗
NASA NTRS2005

Quantifying Exact Motions Along Lineaments on Europa

Evaluating the precise motions along lineaments on the surface of Jupiter's icy moon, Europa, is a valuable tool for interpreting the development and history of lineaments of various morphologies. Such morphologies include strike-slip faults, dilational bands, ridges, and convergence zones. However, the exact mode of origin and kinematic behavior of these various lineaments are not obvious based on morphology alone. In fact, the apparent motions implied by displaced crosscut features can provide misleading indications of true motions along lineaments. Identifying the precise motions (combinations of sliding and opening/closing) is critical to the accurate characterization and interpretation of each of these lineament types. Lineaments of interest (i.e., those having displaced relatively older features in some manner) are identified on Galileo spacecraft images and measurements are made of the total offset, the separation, and relative orientations of crosscut features with respect to the lineament of interest. Specifically, by using these measured quantities and a series of trigonometric equations, the precise motions (i.e., dilation, convergence, strike-slip, or a combination of strike-slip and dilation or convergence) can be determined. These measurements are, however, limited by the resolution of the available images. This study focuses on motion analysis techniques for Europan lineaments and the precise characterization of fault-orthogonal and/or strike-slip motion along lineaments of varying morphologies. We highlight potential pitfalls of cursory analyses of motion indicators. For example, lineaments with obvious lateral offsets have typically been identified simply as strike-slip faults. This assumption may actually be incorrect, as fault-orthogonal motions may contribute to apparent lateral displacements (offsets or separations). Also, variability in the amount of fault motion along the trace length should theoretically be identifiable using the outlined technique. Strike-slip faults on Europa have conventionally been presumed to have a constant slip magnitude along their lengths; however, we posit that variable distributions of slip may be common, as is typical along terrestrial strike-slip faults.

Vetter, J. C.↗
NASA NTRS2016

Closing the Contrast Gap Between Testbed and Model Prediction with WFIRST-CGI Shaped Pupil Coronagraph

JPL has recently passed an important milestone in its technology development for a proposed NASA WFIRST mission: demonstration of a better than 1x10-8 contrast over broad bandwidth (10%) on both shaped pupil coronagraph (SPC) and hybrid Lyot coronagraph (HLC) testbeds with WFIRST telescope optics. Challenges remain, however, in the technology readiness for the proposed mission. One such is the discrepancy between the achieved contrast on testbeds and their corresponding model predictions. A series of testbed diagnosis and modeling activity were planned and (some) carried out on SPC testbed in order to close the gap. A very useful tool we developed was to analyze the measured testbed Jacobian and compare with its model version that was used to control image plane speckle pattern. The difference between these two is the error in the control Jacobian. When the error, which includes both amplitude and phase, was inserted into a separate prediction model with nominal Jacobian, the model prediction became closely matching the SPC testbed behavior in both contrast floor and contrast convergence speed. This offers a new perspective and is a step closer toward model validation for high contrast coronagraph. Further Jacobian analysis and modeling provided clues to the possible sources for the mismatch: the DM misregistration and the underrepresented testbed optics wavefront error (WFE) and the deformable mirror (DM) setting for correcting this WFE. These analyses suggested that high contrast coronagraph has a tight tolerance in the accuracy of its control Jacobian. Modifications to both testbed control model as well as prediction model are being implemented and future works discussed.

diffraction modeling↗
NASA NTRS2020

Passive Roll Stabilization of the Near Earth Asteroid Scout Solar Sail Mission

The Near Earth Asteroid (NEA) Scout is a small satellite, solar sail mission set to launch on Artemis I. Analysis of the NEA Scout solar radiation pressure model determined that, for certain solar incidence angles, there exists at least one locally stable equilibrium point about the “roll” axis, or the axis normal to the solar sail plane. This analysis is extended to three other geometrically similar solar sail models, and a Lyapunov stability analysis is conducted demonstrating that the stable equilibria are locally stable for the majority of roll angles in a wide range of solar incidence angles. Under certain assumptions, it is possible that only two axis control is required for non-spinning solar sails. If the roll axis is controlled, an active rate damping controller can be designed to asymptotically stabilize the system. Simulation results are presented that follow the stability analysis, where in the homogenous case, the system enters a limit cycle and, in the actively damped case, the system asymptotically converges to the equilibrium point.

Ivan Rodrigues Bertaska↗
DOE OSTI2020

MULTI-DIMENSIONAL MODELING OF THE CFR ENGINE FOR THE INVESTIGATION OF SI NATURAL GAS COMBUSTION AND CONTROLLED END-GAS AUTOIGNITION

Engine knock and misfire are barriers to pathways leading to high-efficiency Spark-Ignited (SI) Natural Gas engines. The general tendency to knock is highly dependent on engine operating conditions and the fuel reactivity. The problem is further complicated by low emission limits and the wide range of chemical reactivity in pipeline quality natural gas. Depending on the region and the source of the natural gas, its reactivity, described by its methane number (analogous to the octane number for liquid SI fuels) can span from 65 - 95. In order to realize diesel-like efficiencies, SI natural gas engines must be designed to operate at high BMEP near knock limits over a wide range of fuel reactivity. This requires a deep understanding regarding the combustion-engine interactions pertaining to flame propagation and end-gas autoignition (EGAI). However, EGAI, if controlled, provides an opportunity to increase SI natural gas engine efficiency by increasing combustion rate and the total burned fuel, mitigating the effects of the slow flame speeds of natural gas fuels which generally reduce BMEP and increase unburned hydrocarbon emissions. For this reason, in order to study EGAI phenomenon, the present work highlights multi-dimensional computational fluid dynamics (CFD) models of the Cooperative Fuel Research (CFR) engine. The CFR engine models are used to investigate fuel-engine interactions that lead to EGAI with natural gas, including effects of fuel reactivity, engine operating parameters, and exhaust gas recirculation (EGR). A Three-Pressure Analysis, performed with GT-Power, was used to estimate initial and boundary conditions for the three-dimensional CFD model. CONVERGE CFD v2.4 was used for the three-dimensional CFD modeling where the level set G-Equation model and SAGE detailed chemical kinetics solver were used. An assessment of the different modeling approaches is also provided to evaluate their limitations, advantages and disadvantages, and for which situations they are most applicable. Model validation was performed with experimental data taken with a CFR engine over varying compression ratio, CA50, EGR fraction, and IMEP and shows good agreement in Peak Cylinder Pressure (PCP), PCP crank angle, and the location of the 10%, 50%, and 90% mass fraction burned (CA10, CA50, and CA90, respectively). The models can predict the onset crank angle and pressure rise rate for light, medium, and heavy EGAI under a variety of fuel reactivities and engine operating conditions.

Bestel, Diego↗
DOE OSTI2013

Fission matrix-based Monte Carlo criticality analysis of fuel storage pools

Standard Monte Carlo transport procedures experience difficulties in solving criticality problems in fuel storage pools. Because of the strong neutron absorption between fuel assemblies, source convergence can be very slow, leading to incorrect estimates of the eigenvalue and the eigenfunction. This study examines an alternative fission matrix-based Monte Carlo transport method that takes advantage of the geometry of a storage pool to overcome this difficulty. The method uses Monte Carlo transport to build (essentially) a fission matrix, which is then used to calculate the criticality and the critical flux. This method was tested using a test code on a simple problem containing 8 assemblies in a square pool. The standard Monte Carlo method gave the expected eigenfunction in 5 cases out of 10, while the fission matrix method gave the expected eigenfunction in all 10 cases. In addition, the fission matrix method provides an estimate of the error in the eigenvalue and the eigenfunction, and it allows the user to control this error by running an adequate number of cycles. Because of these advantages, the fission matrix method yields a higher confidence in the results than standard Monte Carlo. We also discuss potential improvements of the method, including the potential for variance reduction techniques. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2012

PROBING PRIMORDIAL NON-GAUSSIANITY WITH WEAK-LENSING MINKOWSKI FUNCTIONALS

We study the cosmological information contained in the Minkowski functionals (MFs) of weak gravitational lensing convergence maps. We show that the MFs provide strong constraints on the local-type primordial non-Gaussianity parameter f {sub NL}. We run a set of cosmological N-body simulations and perform ray-tracing simulations of weak lensing to generate 100 independent convergence maps of a 25 deg{sup 2} field of view for f {sub NL} = -100, 0 and 100. We perform a Fisher analysis to study the degeneracy among other cosmological parameters such as the dark energy equation of state parameter w and the fluctuation amplitude {sigma}{sub 8}. We use fully nonlinear covariance matrices evaluated from 1000 ray-tracing simulations. For upcoming wide-field observations such as those from the Subaru Hyper Suprime-Cam survey with a proposed survey area of 1500 deg{sup 2}, the primordial non-Gaussianity can be constrained with a level of f {sub NL} {approx} 80 and w {approx} 0.036 by weak-lensing MFs. If simply scaled by the effective survey area, a 20,000 deg{sup 2} lensing survey using the Large Synoptic Survey Telescope will yield constraints of f {sub NL} {approx} 25 and w {approx} 0.013. We show that these constraints can be further improved by a tomographic method using source galaxies in multiple redshift bins.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2015

KINEMATICS IN PARTIALLY IONIZED MOLECULAR CLOUDS: IMPLICATIONS FOR THE TRANSITION TO COHERENCE

A previous paper by Bailey and Basu shows analysis of density and mass-to-flux ratio maps for simulations with either an ionization profile which takes into account photoionization (step-like profile) or a cosmic ray only ionization profile. We extend this study to analyze the effect of these ionization profiles on velocity structures, kinematics, and synthetic spectra. Clump regions are found to occur at the convergence of two flows with a low velocity region and velocity direction transition occurring at the junction. Models with evident substructure show that core formation occurs on the periphery of these velocity valleys. Analysis of synthetic spectra reveals the presence of large non-thermal components within low-density gas, especially for models with the step-like ionization profile. All cores show small, sub-thermal relative motions compared to background gas. Large deviations within this analysis are due to the line of sight intersecting low- and high-density regions across the velocity switch transition. Positive deviations correspond to a foreground core moving away from the observer while negative deviations correspond to a background core moving toward the observer. Comparison of velocities resulting from different ionization profiles suggest that high ionization fractions yield supersonic velocities, up to two times the sound speed, while regions with low ionization fractions tend to be subsonic or mildly supersonic. This suggests that the transition to coherence within cores could be a transition between high and low ionization fractions within the gas.

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS2018

Analyzing the Relationship between Tracking and Covariance for Satellite Collision Avoidance

This presentation will display statistical analysis performed with over 150,000 Conjunction Data Messages (CDMs) received for the Earth Observing System (EOS) Aqua, Aura and Terra satellites within the three year period of March 2015 through February 2018. The analysis performed demonstrates the relationship between the covariance of a secondary object and the amount of tracks it receives, as a function of Time to Closest Approach (TCA). The examination indicates that higher tracking frequencies result in lower uncertainties inclusively from the time the data is screened up until TCA. Further investigation displays the impacts of how additional observations of secondary objects lower their position uncertainties and what period of time prior to TCA convergence of the covariance can be expected depending on how well the objects are tracked. Data points from several recent EOS High Interest Events (HIEs) are used to exemplify this analysis.

Collision Avoidance↗
NASA NTRS1994

Multiblock analysis for Shuttle Orbiter reentry heating from Mach 24 to Mach 12

A multiblock, laminar heating analysis for the shuttle orbiter at three trajectory points ranging from Mach 24.3 to Mach 12.86 on reentry is described. The analysis is performed using the Langley Aerothermodynamic Upwind Relaxation Algorithm with a seven species chemical nonequilibrium model. A finite-catalytic-wall model appropriate for shuttle tiles at a radiative equilibrium wall temperature is applied. Computed heating levels are generally in good agreement with the flight data, although a few rather large discrepancies remain unexplained. The multiblock relaxation strategy partitions the flowfield into manageable blocks requiring a fraction of the computational resources (time and memory) required by a full domain approach. In fact, the computational cost for a solution at even a single trajectory point would be prohibitively expensive at the given resolution without the multiblock approach. Converged blocks are reassembled to enable a fully coupled converged solution over the entire vehicle, starting from a nearly converged initial condition.

LANGLEY AEROTHERMODYNAMIC UPWI↗
NASA NTRS1984

Analysis and testing for rotordynamic coefficients of turbulent annular seals with different, directionally-homogeneous surface-roughness treatment for rotor and stator elements

A combined analytical-computational method has been developed to calculate the transient pressure field and rotordynamic coefficients for high-pressure annular seals. The solution procedure applies to constant-clearance or the convergent tapered geometries of multistage centrifugal pumps which may have different surface roughness conditions on the stator and rotor seal elements. In experimental calculations with the method, the turbulent equations of Hirs (1973) are modified slightly to account for different surface roughness conditions, and a perturbation analysis is employed to develop zeroth and first order perturbation equations. Zeroth equations are also used to define both the leakage and circumferential flow due to shear stresses around the stator and rotor surfaces. The solution to analytical equations for four different surface roughnesses confirm the predicted net damping for the seals. A round-holed stator pattern yielded the highest net damping and lowest leakage of all the seals tested.

Childs, D. W.↗
NASA NTRS1975

Simultaneous calculation of aircraft design loads and structural member sizes

A design process which accounts for the interaction between aerodynamic loads and changes in member sizes during sizing of aircraft structures is described. A simultaneous iteration procedure is used wherein both design loads and member sizes are updated during each cycle yielding converged, compatible loads and member sizes. A description is also given of a system of programs which incorporates this process using lifting surface theory to calculate aerodynamic pressure distributions, using a finite-element method for structural analysis, and using a fully stressed design technique to size structural members. This system is tailored to perform the entire process with computational efficiency in a single computer run so that it can be used effectively during preliminary design. Selected results, considering maneuver, taxi, and fatigue design conditions, are presented to illustrate convergence characteristics of this iterative procedure.

Giles, G. L.↗
NASA NTRS2014

Summary and Statistical Analysis of the First AIAA Sonic Boom Prediction Workshop

A summary is provided for the First AIAA Sonic Boom Workshop held 11 January 2014 in conjunction with AIAA SciTech 2014. Near-field pressure signatures extracted from computational fluid dynamics solutions are gathered from nineteen participants representing three countries for the two required cases, an axisymmetric body and simple delta wing body. Structured multiblock, unstructured mixed-element, unstructured tetrahedral, overset, and Cartesian cut-cell methods are used by the participants. Participants provided signatures computed on participant generated and solution adapted grids. Signatures are also provided for a series of uniformly refined workshop provided grids. These submissions are propagated to the ground and loudness measures are computed. This allows the grid convergence of a loudness measure and a validation metric (dfference norm between computed and wind tunnel measured near-field signatures) to be studied for the first time. Statistical analysis is also presented for these measures. An optional configuration includes fuselage, wing, tail, flow-through nacelles, and blade sting. This full configuration exhibits more variation in eleven submissions than the sixty submissions provided for each required case. Recommendations are provided for potential improvements to the analysis methods and a possible subsequent workshop.

Park, Michael A.↗
DOE OSTI2017

Analysis of Monte Carlo accelerated iterative methods for sparse linear systems: Analysis of Monte Carlo accelerated iterative methods for sparse linear systems

Here, we consider hybrid deterministic-stochastic iterative algorithms for the solution of large, sparse linear systems. Starting from a convergent splitting of the coefficient matrix, we analyze various types of Monte Carlo acceleration schemes applied to the original preconditioned Richardson (stationary) iteration. We expect that these methods will have considerable potential for resiliency to faults when implemented on massively parallel machines. We also establish sufficient conditions for the convergence of the hybrid schemes, and we investigate different types of preconditioners including sparse approximate inverses. Numerical experiments on linear systems arising from the discretization of partial differential equations are presented.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Constraining accuracy of the pairwise velocities in N -body simulations using scale-free models

We present a continuation of an analysis that aims to quantify resolution of N-body simulations by exploiting large (up to N = 4096 3 ) simulations of scale-free cosmologies run using abacus. Here, we focus on radial pairwise velocities of the matter field, both by direct estimation and through the cumulative two-point correlation function (using the pair conservation equation). We find that convergence at the 1 per cent level of the mean relative pairwise velocity can be demonstrated over a range of scales, evolving from a few times the grid spacing at early times to slightly below this scale at late times. We show the analysis of two different box sizes as well as from averaging results from the smaller boxes, and compare the power of the two aforementioned estimators in constraining accuracy at each scale. Down to scales of the order of the smoothing parameter, convergence is obtained at ∼5 per cent precision, and shows a behaviour indicating asymptotic stable clustering. In conclusion, we also infer for LCDM simulations conservative estimates on the evolution of the lower cut-off to resolution (at 1 and 5 per cent precision) as a function of redshift.

79 ASTRONOMY AND ASTROPHYSICS↗