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 829 records · Page 46
All sourcesNASA NTRSDOE OSTI
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 OSTI2019

Domain decomposition and multiscale mortar mixed finite element methods for linear elasticity with weak stress symmetry

Two non-overlapping domain decomposition methods are presented for the mixed finite element formulation of linear elasticity with weakly enforced stress symmetry. The methods utilize either displacement or normal stress Lagrange multiplier to impose interface continuity of normal stress or displacement, respectively. By eliminating the interior subdomain variables, the global problem is reduced to an interface problem, which is then solved by an iterative procedure. The condition number of the resulting algebraic interface problem is analyzed for both methods. A multiscale mortar mixed finite element method for the problem of interest on non-matching multiblock grids is also studied. It uses a coarse scale mortar finite element space on the non-matching interfaces to approximate the trace of the displacement and impose weakly the continuity of normal stress. A priori error analysis is performed. It is shown that, with appropriate choice of the mortar space, optimal convergence on the fine scale is obtained for the stress, displacement, and rotation, as well as some superconvergence for the displacement. Computational results are presented in confirmation of the theory of all proposed methods.

Mathematics↗
DOE OSTI2014

Electronic state spectroscopy of diiodomethane (CH{sub 2}I{sub 2}): Experimental and computational studies in the 30 000–95 000 cm{sup −1} region

The electronic absorption spectrum of diiodomethane in the 30 000–95 000 cm{sup −1} region is investigated using synchrotron radiation; the spectrum in the 50 000–66 500 cm{sup −1} region is reported for the first time. The absorption bands in the 30 000–50 000 cm{sup −1} region are attributed to valence transitions, while the vacuum ultraviolet (VUV) spectrum (50 000–95 000 cm{sup −1}) is dominated by several Rydberg series converging to the first four ionization potentials of CH{sub 2}I{sub 2} at 9.46, 9.76, 10.21, and 10.56 eV corresponding to the removal of an electron from the outermost 3b{sub 2}, 2b{sub 1}, 1a{sub 2}, and 4a{sub 1} non-bonding orbitals, respectively. Rydberg series of ns, np, and nd type converging to each of the four ionization potentials are assigned based on a quantum defect analysis. Time dependent density functional theory calculations of excited states support the analysis and help in interpretation of the Rydberg and valence nature of observed transitions. Density functional theory calculations of the neutral and ionic ground state geometries and vibrational frequencies are used to assign the observed vibronic structure. Vibronic features accompanying the Rydberg series are mainly due to excitation of the C-I symmetric stretch (ν{sub 3}) and CH{sub 2} wag (ν{sub 8}) modes, with smaller contributions from the C-H symmetric stretch (ν{sub 1}). UV absorption bands are assigned to low lying valence states 1{sup 1}B{sub 2}, 1{sup 1}B{sub 1}, 2{sup 1}A{sub 1}, 3{sup 1}A{sub 1}, 2{sup 1}B{sub 1}, and 2{sup 1}B{sub 2} and the unusually high underlying intensity in parts of the VUV spectrum is attributed to valence states with high oscillator strength. This is the first report of a comprehensive Rydberg series and vibronic analysis of the VUV absorption spectrum of CH{sub 2}I{sub 2} in the 50 000–85 000 cm{sup −1} region. The VUV absorption spectrum of CD{sub 2}I{sub 2} which serves to verify and consolidate spectral assignments is also reported here for the first time.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2018

Dynamic compaction of nickel powder examined by x-ray phase contrast imaging

Understanding the shock response of porous, granular materials is important for many scientific applications. In this work, a propagation-based x-ray phase contrast imaging technique was used to examine the shock compaction response of Ni powder (particle size: 30 and 45 µm), encapsulated in PMMA cylinders, in situ and in real time. The propagating shock wave in the PMMA cylinder and the deformation of the Ni powder column were recorded, but the compaction wave in the Ni powder could not be observed due to insufficient penetration of x-ray photons through the sample. The overall shape of the deformed Ni column downstream of the PMMA shock looked qualitatively similar for both particle sizes: it deformed inward and then outward similar to a converging-diverging nozzle, followed by a ‘mushroom’-like structure at the impact face. A preliminary analysis of our data using shock-polar method provided insight into the observed shock and flow deflection angles at the PMMA/Ni interface, and permitted evaluation of the Ni powder EOS used in this work.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2019

A Homotopy-Based Method for Robust Computation of Controlling Unstable Equilibrium Points

This paper proposes a homotopy-based method for efficiently computing and precisely determining controlling unstable equilibrium points (CUEPs) in transient stability analysis and screening using direct methods. The proposed method is intended to overcome potential numerical convergence problems associated with computing CUEPs. These convergence problems are mostly due to irregularities in the region of convergence of CUEPs which make it difficult to determine suitable initial points. The most commonly used initial points in computing CUEPs are minimum gradient points (MGPs) on the stability boundary which in turn rely on the determination of exit points (EPs-points at which sustained-fault trajectories exit system stability boundaries). However, determination of MGPs and EPs are computationally involved and sophisticated numerical methods have been introduced such as stability-boundary-following procedure and shadowing method to correct the solution trajectory. The proposed method does not depend on MGPs and also does not require accurate EPs. The proposed method maps the solution from EPs to the CUEPs. The proposed method is applied on several known systems including the NE 39 bus system and the reduced WECC system and the results are provided. The results are compared with the traditional methods such as Boundary of stability region based CUEP method (BCU method).

42 ENGINEERING↗
DOE OSTI2016

Computing and Applying Atomic Regulons to Understand Gene Expression and Regulation

Understanding gene function and regulation is essential for the interpretation, prediction, and ultimate design of cell responses to changes in the environment. A multitude of technologies, abstractions, and interpretive frameworks have emerged to answer the challenges presented by genome function and regulatory network inference. Here, we propose a new approach for producing biologically meaningful clusters of coexpressed genes, called Atomic Regulons (ARs), based on expression data, gene context, and functional relationships. We demonstrate this new approach by computing ARs for Escherichia coli , which we compare with the coexpressed gene clusters predicted by two prevalent existing methods: hierarchical clustering and k-means clustering. We test the consistency of ARs predicted by all methods against expected interactions predicted by the Context Likelihood of Relatedness (CLR) mutual information based method, finding that the ARs produced by our approach show better agreement with CLR interactions. We then apply our method to compute ARs for four other genomes: Shewanella oneidensis , Pseudomonas aeruginosa , Thermus thermophilus , and Staphylococcus aureus . We compare the AR clusters from all genomes to study the similarity of coexpression among a phylogenetically diverse set of species, identifying subsystems that show remarkable similarity over wide phylogenetic distances. We also study the sensitivity of our method for computing ARs to the expression data used in the computation, showing that our new approach requires less data than competing approaches to converge to a near final configuration of ARs. We go on to use our sensitivity analysis to identify the specific experiments that lead most rapidly to the final set of ARs for E. coli . As a result, this analysis produces insights into improving the design of gene expression experiments.

59 BASIC BIOLOGICAL SCIENCES↗
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.↗