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 577 records · Page 32
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Low-order preconditioning of the Stokes equations

A well-known strategy for building effective preconditioners for higher-order discretizations of some PDEs, such as Poisson's equation, is to leverage effective preconditioners for their low-order analogs. In this work, we show that high-quality preconditioners can also be derived for the Taylor–Hood discretization of the Stokes equations in much the same manner. In particular, we investigate the use of geometric multigrid based on the Q 1 iso Q 2 /Q 1 discretization of the Stokes operator as a preconditioner for the Q 2 /Q 1 discretization of the Stokes system. We utilize local Fourier analysis to optimize the damping parameters for Vanka and Braess–Sarazin relaxation schemes and to achieve robust convergence. Furthermore, these results are then verified and compared against the measured multigrid performance. While geometric multigrid can be applied directly to the Q 2 /Q 1 system, our ultimate motivation is to apply algebraic multigrid within solvers for Q 2 /Q 1 systems via the Q 1 iso Q 2 /Q 1 discretization, which will be considered in a companion paper.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Assessing the Impact of a Novel TBC Material on Heat Transfer in a Spark Ignition Engine through 3D CFD-FEA Co-Simulation Routine

Thermal barrier coatings (TBCs) have been of interest since the 1970s for application in internal combustion (IC) engines. Thin TBCs exhibit a temperature swing phenomenon wherein wall temperatures dynamically respond to the transient working-gas temperature throughout the engine cycle, thus reducing the temperature difference driving the heat transfer. Determining these varying wall temperatures is necessary to evaluate and study the effect of coatings on wall heat transfer. This study focuses on developing a 3D computational fluid dynamics (CFD)-finite element analysis (FEA) coupled simulation, or co-simulation, routine to determine the wall temperatures of a piston coated with a thin TBC layer subject to spark ignition combustion heat flux. A CONVERGE 3D-CFD model was used to simulate the combustion process in a single-cylinder, light-duty experimental spark ignition (SI) engine. Transient piston heat transfer analysis was conducted using ABAQUS, a FEA package, under the simulated combustion heat flux load. The effect of the temperature swing phenomenon due to this TBC layer was observed in a CFD simulation by implementing the FEA results as the piston thermal boundary conditions. The boundary conditions were passed between the CFD and FEA tools until a quasi-steady state solution was achieved. Furthermore, a reduction in wall heat transfer was observed due to a reduced temperature difference between the wall and the working gas.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2022

Shock compression behavior of stainless steel 316L octet-truss lattice structures

Lattice structures offer desirable mechanical properties for applications of energy absorption and impact mitigation but limited research has been carried out on their shock compression behavior. In this work, the shock compression behavior of stainless steel 316L (SS316L) octet-truss lattice structures was investigated through experimental techniques and numerical simulations. Plate impact experiments with high-speed imaging were conducted at impact velocities of 270 – 390 m/s on lattice specimens with 5x5x10 unit cell geometries additively manufactured (AM) using direct metal laser sintering. High-speed imaging together with digital image correlation was used to extract full-field measurements and define a two-wave structure consisting of an elastic wave and planar compaction (shock) wave which propagated along the impact direction. A linear shock velocity versus particle velocity relation was found to approximate the measurements with a unit slope and a linear fit constant equal to the crushing speed. Furthermore, the shock velocity versus particle velocity relation, full-field measurements, and elastic limit together with the Eulerian form of the Rankine-Hugoniot jump conditions were used to find relations for the stress and internal energy behind the shock. Stress behind the shock increased with relative density and particle velocity, and specific internal energy converged to a single curve similar to that of bulk AM SS316L. Explicit finite element analysis using the Johnson-Cook constitutive model demonstrated similar shock behavior observed in experiments and a linear shock velocity versus particle velocity relation and corresponding Hugoniot calculations were found to be in agreement with experimental results. Numerical simulations confirmed negligible effects of exterior versus interior measurements and further validated the application of one-dimensional shock theory.

42 ENGINEERING↗
DOE OSTI2020

Dynamic compressed sensing for real-time tomographic reconstruction

Electron tomography has achieved higher resolution and quality at reduced doses with recent advances in compressed sensing. Compressed sensing (CS) exploits the inherent sparse signal structure to efficiently reconstruct three-dimensional (3D) volumes at the nanoscale from undersampled measurements. However, the process bottlenecks 3D reconstruction with computation times that run from hours to days. Here we demonstrate a framework for dynamic compressed sensing that produces a 3D specimen structure that updates in real-time as new specimen projections are collected. Researchers can begin interpreting 3D specimens as data is collected to facilitate high-throughput and interactive analysis. Using scanning transmission electron microscopy (STEM), we show that dynamic compressed sensing accelerates the convergence speed by ~3-fold while also reducing its error by 27% for a Au/SrTiO 3 nanoparticle specimen. Before a tomography experiment is completed, the 3D tomogram has interpretable structure within ~33% of completion and fine details are visible as early as ~66%. Upon completion of an experiment, a high-fidelity 3D visualization is produced without further delay. Additionally, we report reconstruction parameters that tune data fidelity can be manipulated throughout the computation without re-running the entire process.

47 OTHER INSTRUMENTATION↗
DOE OSTI2020

An Investigation of Poorly Studied Open Cluster NGC 4337 Using Multicolor Photometric and Gaia DR2 Astrometric Data

We present a comprehensive analysis (photometric and kinematical) of the poorly studied open cluster NGC 4337 using 2MASS, WISE, APASS, and Gaia DR2 databases. By determining the membership probabilities of stars, we identified the 624 most probable members with membership probability higher than 50% by using proper motion and parallax data taken from Gaia DR2. The mean proper motion of the cluster is obtained as μ{sub x}=−8.83±0.01 and μ {sub y} = 1.49 ± 0.006 mas yr{sup −1}. We find the normal interstellar extinction toward the cluster region. The radial distribution of members provides a cluster radius of 7.′75 (5.63 pc). The estimated age of 1600 ± 180 Myr indicates that NGC 4337 is an old open cluster with a bunch of red giant stars. The overall mass function slope for main-sequence stars is found as 1.46±0.18 within the mass range 0.75−2.0 M{sub ⊙}, which is in fair agreement with Salpeter’s value (x = 1.35) within uncertainty. The present study demonstrates that NGC 4337 is a dynamically relaxed open cluster. Using the Galactic potential model, Galactic orbits are obtained for NGC 4337. We found that this object follows a circular path around the Galactic center. Under the kinematical analysis, we compute the apex coordinates (A, D) by using two methods: (i) the classical convergent point method and (ii) the AD-diagram method. The obtained coordinates are (A {sub conv}, D {sub conv}) = (96.°27 ± 0.°10, 13.°14 ± 0.°27) and (A {sub ◦}, D {sub ◦}) = (100.°282 ± 0.°10, 9.°577 ± 0.°323) respectively. We also computed the Velocity Ellipsoid Parameters, matrix elements (μ {sub ij}), direction cosines (l {sub j}, m {sub j}, n {sub j}), and the Galactic longitude of the vertex (l {sub 2}).

79 ASTRONOMY AND ASTROPHYSICS↗
NASA NTRS1977

Viking satellite orbit determination

During the summer of 1976, the two Viking spacecraft, each consisting of an orbiter-lander combination, were inserted into orbit about Mars. The paper describes the experiences of the Viking Satellite Orbit Determination Team in determining Mars centered ephemerides of the orbiters and positions of the landers from the two-way Doppler and range data, and synthesizes the different phases of the navigation plan which involves pre-flight modeling and error analysis for all Viking navigation functions from launch through landing. The problem of initial orbit convergence is solved by using DPODP's (Double Precision Orbit Determination Program) square-root batch data filter, and gravity models for both Viking I and Viking II were produced from the combination of short arc estimates. The importance of synchronous Viking orbits (with the rotational period of Mars) is stressed and future extended missions of the spacecraft are outlined.

Hildebrand, C. E.↗
NASA NTRS1978

Rain-rate estimates for an attenuating radar

An iterative method has been used to produce rain-rate estimates for an attenuating frequency radar. An infinite number of higher-order estimates are shown to converge in the limit to the Hitschfeld-Bordan solution under certain conditions. An error analysis was performed by a model accounting for the randomness of the radar return power, the k-Z, Z-R relations, and offsets in the radar calibration constant. Since the behavior of the estimates strongly depends on system errors, the choice of the best estimate requires a knowledge of the variance and range of offsets in meteorological data and the calibration constant. As the errors increase, the use of the lower orders avoids significant overpredictions. In order to obtain reliable rain-rate predictions in the presence of realistic errors, an antenna pointing angles away from the horizontal at frequencies in the lower end of the X-band may be used. Such antenna configurations insure low attenuation.

Meneghini, R.↗
NASA NTRS1980

Interlaminar stress singularities at a straight free edge in composite laminates

A quasi three dimensional finite element analysis was used to analyze the edge stress problem in four-ply, composite laminates. Convergence studies were made to explore the existence of stress singularities near the free edge. The existence of stress singularities at the intersection of the interface and the free edge is confirmed.

Raju, I. S.↗
NASA NTRS1980

Geomagnetic field modeling by optimal recursive filtering

Data sets selected for mini-batches and the software modifications required for processing these sets are described. Initial analysis was performed on minibatch field model recovery. Studies are being performed to examine the convergence of the solutions and the maximum expansion order the data will support in the constant and secular terms.

Source record↗
NASA NTRS1981

On the warming of the polar troposphere during the stratospheric warming of January 1977

In January 1977, a warming of the polar troposphere occurred following a major stratospheric warming. The temperature increase in the polar troposphere is found to be larger below 400 mb than in the upper troposphere. Analysis of the heat budget of the polar troposphere shows that the eddy heat flux convergence is the major factor responsible for the tropospheric warming.

Hansen, A. R.↗
NASA NTRS1983

Aerothermal modeling program, phase 1

Aerothermal submodels used in analytical combustor models are analyzed. The models described include turbulence and scalar transport, gaseous full combustion, spray evaporation/combustion, soot formation and oxidation, and radiation. The computational scheme is discussed in relation to boundary conditions and convergence criteria. Also presented is the data base for benchmark quality test cases and an analysis of simple flows.

Srinivasan, R.↗
NASA NTRS1983

Case study of Pampa, Texas, multicell storms

Analysis of the Pampa, TX, multicell storms showed prestorm environmental conditions of a strong horizontal convergence of moisture, mixing ratio 12 g/kg at the 850 mb level in the Texas Panhandle areas and strong winds with shear veering with height. Rapid-scan satellite imagery showed that the clouds penetrated above the tropopause and cloud-top temperatures were at least 4 to 9 C colder than the temperature of the tropopause. This formation and collapsing of high-mass-density overshooting cloud tops above the tropopause is characteristic of tornadoes in the middle portion of the United States.

Hung, R. J.↗
NASA NTRS1987

A quadrilateral shell element using a mixed formulation

A simple quadrilateral shell element consisting of five nodes, four corner nodes and a central node, is developed for linear elastic analysis of thin as well as moderately thick shells. Based on a modified Hellinger-Reissner principle, finite element equations are derived from the assumed displacement and strain fields. By carefully choosing appropriate strain terms, all kinematic deformation modes are suppressed. Although the present element is similar to a displacement-based degenerated shell, no locking is experienced when it is applied to thin shell problems. Five examples are given to illustrate the analysis capability of the shell element. Numerical results indicate that the element shows fast mesh convergence and gives excellent stress predictions.

Saleeb, A. F.↗
NASA NTRS1987

Exact image theory for the problem of dielectric/magnetic slab

Exact image method, recently introduced for the exact solution of electromagnetic field problems involving homogeneous half spaces and microstrip-like geometries, is developed for the problem of homogeneous slab of dielectric and/or magnetic material in free space. Expressions for image sources, creating the exact reflected and transmitted fields, are given and their numerical evaluation is demonstrated. Nonradiating modes, guided by the slab and responsible for the loss of convergence of the image functions, are considered and extracted. The theory allows, for example, an analysis of finite ground planes in microstrip antenna structures.

Lindell, I. V.↗
NASA NTRS1989

Mesh refinement in finite element analysis by minimization of the stiffness matrix trace

Most finite element packages provide means to generate meshes automatically. However, the user is usually confronted with the problem of not knowing whether the mesh generated is appropriate for the problem at hand. Since the accuracy of the finite element results is mesh dependent, mesh selection forms a very important step in the analysis. Indeed, in accurate analyses, meshes need to be refined or rezoned until the solution converges to a value so that the error is below a predetermined tolerance. A-posteriori methods use error indicators, developed by using the theory of interpolation and approximation theory, for mesh refinements. Some use other criterions, such as strain energy density variation and stress contours for example, to obtain near optimal meshes. Although these methods are adaptive, they are expensive. Alternatively, a priori methods, until now available, use geometrical parameters, for example, element aspect ratio. Therefore, they are not adaptive by nature. An adaptive a-priori method is developed. The criterion is that the minimization of the trace of the stiffness matrix with respect to the nodal coordinates, leads to a minimization of the potential energy, and as a consequence provide a good starting mesh. In a few examples the method is shown to provide the optimal mesh. The method is also shown to be relatively simple and amenable to development of computer algorithms. When the procedure is used in conjunction with a-posteriori methods of grid refinement, it is shown that fewer refinement iterations and fewer degrees of freedom are required for convergence as opposed to when the procedure is not used. The mesh obtained is shown to have uniform distribution of stiffness among the nodes and elements which, as a consequence, leads to uniform error distribution. Thus the mesh obtained meets the optimality criterion of uniform error distribution.

Kittur, Madan G.↗
NASA NTRS1994

Neural networks: Application to medical imaging

The research mission is the development of computer assisted diagnostic (CAD) methods for improved diagnosis of medical images including digital x-ray sensors and tomographic imaging modalities. The CAD algorithms include advanced methods for adaptive nonlinear filters for image noise suppression, hybrid wavelet methods for feature segmentation and enhancement, and high convergence neural networks for feature detection and VLSI implementation of neural networks for real time analysis. Other missions include (1) implementation of CAD methods on hospital based picture archiving computer systems (PACS) and information networks for central and remote diagnosis and (2) collaboration with defense and medical industry, NASA, and federal laboratories in the area of dual use technology conversion from defense or aerospace to medicine.

Clarke, Laurence P.↗
NASA NTRS2001

Seasonal Variation in Precipitation Patterns to the Global Ocean: An Analysis of the GPCP Version 2 Data Set

An analysis of temporal and spatial variation of oceanic precipitation was conducted on the GPCP version two data set. While the precipitation patterns observed are generally similar to the previous climatologies, new features and greater detail of global precipitation were revealed from out analysis of the GPCP data set. High precipitation waw observed in the inter-tropical convergence zone, the South Pacific convergence zone, and the storm tracks in the North Pacific and Atlantic Oceans. Low precipitation was observe in the Polar regions and in the subtropics of the East Pacific, East Atlantic, and the Southeast and Northwest Indian Ocean. The spatial coverage of these high and low precipitation regions changed thruough the year. A strong seasonal cycle or precipitation was observed for the Northern and the Southern Hemispheres and for each ocean basin. Global precipitation also varied significantly with both latitude and longitude, with a latitudinal maximum at 56 degrees South, 39 degrees South, 4 degrees South, 6 degrees North, 39 degrees North, and 56 degress North, and a longitudinal maxiumum over each ocean. The seasonal varying precipitation patterns are a foundation for evaluating the effect of wet deposition on ocean circulation, flux of chemical species, and its effect on marine ecosystems.

Miller, Richard↗
NASA NTRS2018

Hollow Cathode Thermal Modelling and Self-Consistent Plasma Solution: Work Function Evaluation for a LaB6 Cathode.

Self-heating hollow cathodes are central components in modern electric thrusters. The plasma discharge inside these devices heats the internal components, thus maintaining the temperatures required for electron emission. Precise knowledge of the physical phenomena governing hollow cathode operation is key to predict their lifetime, specifically, their thermionic emission characteristics. A simulation platform has been built to couple plasma and thermal models of the self-heating hollow cathode to produce a self-consistent solution. A self-consistent solution has been found for a LaB6 hollow cathode operating at 25A and 13 sccm where the work function is assumed to be spatially uniform along the emitter with a value which is allowed to vary as the coupled model iterates to a self-consistent solution. The emitter temperature from the converged solution does not agree with experimental temperature measurements, however. The results of a sensitivity analysis suggest that none of the tolerances in the measurement are responsible for the discrepancy. We hypothesize that either the work function needs to be a function of position along the emitting surfaces or the heat fluxes have been overestimated in the plasma solver.

Meiron, Daniel I.↗