Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Iterative methods”

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 541 records · Page 30
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

BigNeuron: a resource to benchmark and predict performance of algorithms for automated tracing of neurons in light microscopy datasets

BigNeuron is an open community bench-testing platform with the goal of setting open standards for accurate and fast automatic neuron tracing. We gathered a diverse set of image volumes across several species that is representative of the data obtained in many neuroscience laboratories interested in neuron tracing. Here, we report generated gold standard manual annotations for a subset of the available imaging datasets and quantified tracing quality for 35 automatic tracing algorithms. The goal of generating such a hand-curated diverse dataset is to advance the development of tracing algorithms and enable generalizable benchmarking. Together with image quality features, we pooled the data in an interactive web application that enables users and developers to perform principal component analysis, t-distributed stochastic neighbor embedding, correlation and clustering, visualization of imaging and tracing data, and benchmarking of automatic tracing algorithms in user-defined data subsets. The image quality metrics explain most of the variance in the data, followed by neuromorphological features related to neuron size. Furthermore, we observed that diverse algorithms can provide complementary information to obtain accurate results and developed a method to iteratively combine methods and generate consensus reconstructions. The consensus trees obtained provide estimates of the neuron structure ground truth that typically outperform single algorithms in noisy datasets. However, specific algorithms may outperform the consensus tree strategy in specific imaging conditions. Finally, to aid users in predicting the most accurate automatic tracing results without manual annotations for comparison, we used support vector machine regression to predict reconstruction quality given an image volume and a set of automatic tracings.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Development of a continuous synthesis process for carbamazepine using validated in-line Raman spectroscopy and kinetic modelling for disturbance simulation

Mitigation of failure modes in the continuous synthesis (CS) of a drug substance (DS) has the potential to widen the adoption of continuous manufacturing (CM) technologies by the pharmaceutical industry. Here, this work demonstrates the development of a robust continuous process for the synthesis of carbamazepine (CBZ), an essential medicine as per the World Health Organization (WHO), facilitated by kinetic modelling and monitored by in-line Raman spectroscopy. Accurate kinetic modelling and the use of validated process analytical technology (PAT) models for quantitative measurement were found to play an important role in developing CS of drug substances. Kinetic data for the formation of CBZ from iminostilbene (ISB) were collected by batch reaction sampling and high-performance liquid chromatography (HPLC) analysis. A non-linear solver and iterative method was applied to determine two sets of Arrhenius parameters simultaneously for the reaction system by minimizing the standard error of the model fit. The start-up and dynamic equilibrium stages for the CS of CBZ using a continuous stirred tank reactor (CSTR) were modelled based on the batch kinetic data and employed to optimize conversion and simulate process disturbances. An in-line Raman spectroscopy method was successfully developed, validated, and integrated to determine the concentrations of CBZ and ISB within the operating range for the CS. The CS kinetic model was evaluated experimentally from startup to dynamic equilibrium over 10 residence times with monitoring by HPLC and in-line Raman spectroscopy. The developed kinetic model in tandem with in-line Raman spectroscopy successfully predicted disturbances due to changes in process variables and can serve as a useful tool in the future design of advanced process control strategies for the continuous synthesis of CBZ.

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

Stress-induced Effects Caused by 3D IC TSV Packaging in Advanced Semiconductor Device Performance

Potential challenges with managing mechanical stress and the consequent effects on device performance for advanced 3D through-silicon-via (TSV) based technologies are outlined. The paper addresses the growing need in a simulation-based design verification flow capable to analyze a design of 3D IC stacks and to determine across-die out-of-spec variations in device electrical characteristics caused by the layout and through-silicon-via (TSV)/package-induced mechanical stress. The limited characterization/measurement capabilities for 3D IC stacks and a strict ''good die'' requirement make this type of analysis critical for the achievement of an acceptable level of functional and parametric yield and reliability. The paper focuses on the development of a design-for-manufacturability (DFM) type of methodology for managing mechanical stresses during a sequence of designs of 3D TSV-based dies, stacks and packages. A set of physics-based compact models for a multi-scale simulation to assess the mechanical stress across the device layers in silicon chips stacked and packaged with the 3D TSV technology is proposed. A calibration technique based on fitting to measured stress components and electrical characteristics of the test-chip devices is presented. A strategy for generation of a simulation feeding data and respective materials characterization approach are proposed, with the goal to generate a database for multi-scale material parameters of wafer-level and package-level structures. For model validation, high-resolution strain measurements in Si channels of the test-chip devices are needed. At the nanoscale, the transmission electron microscopy (TEM) is the only technique available for sub-10 nm strain measurements so far.

36 MATERIALS SCIENCE↗
DOE OSTI2012

Plasma screening effects on the energies of hydrogen atom

A more general exponential cosine screened Coulomb potential is used for the first time to investigate the screening effects on the hydrogen atom in plasmas. This potential is examined for four different cases that correspond to four different type potentials when the different parameters are used in the potential within the framework of the well-known asymptotic iteration method. By solving the corresponding the radial Schroedinger equation with the screened and exponential cosine screened Coulomb potentials and comparing the obtained energy eigenvalues with the results of other studies, the applicability of the method to this kind of plasma physics problem is shown. The energy values of more general exponential cosine screened Coulomb potential are presented for various parameters in the potential. One of the advantages of the present potential is that it exhibits stronger screening effect than that of the exponential cosine screened Coulomb potential and it is also reduced to screened Coulomb and exponential cosine screened Coulomb as well as Coulomb potentials for special values of parameters. The parameters in the potential would be useful to model screening effects which cause an increase or decrease in the energy values of hydrogen atom in both Debye and quantum plasmas and in this manner this potential would be useful for the investigations of the atomic structure and collisions in plasmas.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Plasma screening effects on the energies of hydrogen atom under the influence of velocity-dependent potential

In order to examine the plasma screening and velocity-dependent potential effects on the hydrogen atom, the Schrödinger equation including a more general exponential cosine screened Coulomb and velocity-dependent potential is solved numerically in the framework asymptotic iteration method. The more general exponential cosine screened Coulomb potential is used to model Debye and quantum plasma for the specific values of the parameters in its structure. However, in order to examine effects of velocity-dependent potential on energy values of hydrogen atom in Debye and quantum plasma, the isotropic form factor of velocity-dependent potential is given as harmonic oscillator type, ρ(r)=ρ{sub o}r{sup 2}. Then, the energies of s and p states are calculated numerically without any approximation. In order to investigate thoroughly plasma screening effects and contribution of velocity-dependent potential on energy values of hydrogen atom, the corresponding calculations are carried out by using different values of parameters of more general exponential cosine screened Coulomb potential and isotropic dependence, results of which are discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

The hydrogen atom in plasmas with an external electric field

We numerically solve the Schrödinger equation, using a more general exponential cosine screened Coulomb (MGECSC) potential with an electric field, in order to investigate the screening and weak external electric field effects on the hydrogen atom in plasmas. The MGECSC potential is examined for four different cases, corresponding to different screening parameters of the potential and the external electric field. The influences of the different screening parameters and the weak external electric field on the energy eigenvalues are determined by solving the corresponding equations using the asymptotic iteration method (AIM). It is found that the corresponding energy values shift when a weak external electric field is applied to the hydrogen atom in a plasma. This study shows that a more general exponential cosine screened Coulomb potential allows the influence of an applied, weak, external electric field on the hydrogen atom to be investigated in detail, for both Debye and quantum plasmas simultaneously. This suggests that such a potential would be useful in modeling similar effects in other applications of plasma physics, and that AIM is an appropriate method for solving the Schrödinger equation, the solution of which becomes more complex due to the use of the MGECSC potential with an applied external electric field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2014

Analytical solutions of a generalized non-central potential in N-dimensions

We present that N-dimensional non-relativistic wave equation for the generalized non-central potential with arbitrary angular momentum is analytically solvable in the hyperspherical coordinates. Asymptotic iteration method as a different approach is applied to obtain N-dimensional energy eigenvalues and the corresponding eigenfunctions. In hyperspherical coordinates, the wave function solutions are obtained in terms of hypergeometric functions and Jacobi polynomials. The bound states of quantum systems under consideration for some special cases, such as Hartmann and Makarov potentials, have been discussed in N-dimensions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2015

Confinement effects of magnetic field on two-dimensional hydrogen atom in plasmas

In this study, for the first time, the Schrödinger equation with more general exponential cosine screened Coulomb (MGECSC) potential is solved numerically in the presence and in the absence of an external magnetic field within two-dimensional formalism using the asymptotic iteration method. The MGECSC potential includes four different potential forms when considering different sets of the parameters in the potential. The plasma screening effects in the weak and strong magnetic field regimes as well as the confinement effects of magnetic field on the two-dimensional hydrogen atom in Debye and quantum plasmas are investigated by solving the corresponding equations. It is found that applying a uniform magnetic field on the hydrogen atom embedded in a plasma leads to change in the profile of the total interaction potential. Thus, confinement effects of magnetic field on hydrogen atom embedded in Debye and quantum plasmas modeled by a MGECSC potential lead to shift bound state energies. This effect would be important to isolate the plasma from the external environment in the experimental applications of plasma physics.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2015

Effects of laser radiation field on energies of hydrogen atom in plasmas

In this study, for the first time, the Schrödinger equation with more general exponential cosine screened Coulomb (MGECSC) potential is solved numerically in the presence of laser radiation field within the Ehlotzky approximation using the asymptotic iteration method. The MGECSC potential includes four different potential forms in consideration of different sets of the parameters in the potential. By applying laser field, the total interaction potential of hydrogen atom embedded in plasmas converts to double well-type potential. The plasma screening effects under the influence of laser field as well as confinement effects of laser field on hydrogen atom in Debye and quantum plasmas are investigated by solving the Schrödinger equation with the laser-dressed MGECSC potential. It is resulted that since applying a monochromatic laser field on hydrogen atom embedded in a Debye and quantum plasma causes to shift in the profile of the total interaction potential, the confinement effects of laser field on hydrogen atom in plasmas modeled by the MGECSC potential change localizations of energy states.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2019

Electrostatic equilibria of non-neutral plasmas confined in a Penning trap with axially varying magnetic field

A procedure for computing the electrostatic equilibria of non-neutral plasmas in a Penning trap with a nonuniform magnetic field by solving Poisson's equation using an iterative method is described. Plasma equilibria in a model Penning trap with high and low field regions are computed. The plasma is assumed to follow the Boltzmann density distribution along magnetic field lines. Correspondence with previous investigations examining similar configurations analytically and using particle-in-cell simulations is found. The relationship between the plasma density in low and high field regions is examined for various plasma temperatures, densities, magnetic mirror ratios, and plasma and electrode radii. An analytical description of the radial density profile in the high field region is developed and compared with the computed equilibria. A concept is described for cooling a positron plasma with laser-cooled ions trapped axially within a high magnetic field region, while antiprotons are trapped axially separated from the laser-cooled ions within a low field region, and the positron plasma extends to both regions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2011

TEMPERATURE STRUCTURE OF PROTOPLANETARY DISKS UNDERGOING LAYERED ACCRETION

We calculate the temperature structures of protoplanetary disks (PPDs) around T Tauri stars heated by both incident starlight and viscous dissipation. We present a new algorithm for calculating the temperatures in disks in hydrostatic and radiative equilibrium, based on Rybicki's method for iteratively calculating the vertical temperature structure within an annulus. At each iteration, the method solves for the temperature at all locations simultaneously, and converges rapidly even at high (>>10{sup 4}) optical depth. The method retains the full frequency dependence of the radiation field. We use this algorithm to study for the first time disks evolving via the magnetorotational instability. Because PPD midplanes are weakly ionized, this instability operates preferentially in their surface layers, and disks will undergo layered accretion. We find that the midplane temperatures T{sub mid} are strongly affected by the column density {Sigma}{sub a} of the active layers, even for fixed mass accretion rate M-dot . Models assuming uniform accretion predict midplane temperatures in the terrestrial planet forming region several x 10{sup 2} K higher than our layered accretion models do. For M-dot < 10{sup -7} M{sub sun} yr{sup -1} and the column densities {Sigma}{sub a} < 10 g cm{sup -2} associated with layered accretion, disk temperatures are indistinguishable from those of a passively heated disk. We find emergent spectra are insensitive to {Sigma}{sub a}, making it difficult to observationally identify disks undergoing layered versus uniform accretion.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2012

POLARIZED PARTIAL FREQUENCY REDISTRIBUTION IN SUBORDINATE LINES. II. SOLUTION OF THE TRANSFER EQUATION WITH RAYLEIGH SCATTERING

It is quite common in line formation theory to treat scattering in subordinate lines under the assumption of complete frequency redistribution (CRD). The partial frequency redistribution (PRD) in subordinate lines cannot always be approximated by CRD, especially when the polarization state of the line radiation is taken into account. Here we investigate the PRD effects in subordinate lines including scattering polarization. The line formation is described by a polarized non-LTE line transfer equation based on a two-level atom model. We use the recently derived subordinate line redistribution matrix. We devise polarized approximate lambda iteration methods to solve the concerned transfer problem. The linear polarization profiles of subordinate lines formed in non-magnetic (Rayleigh) scattering atmospheres are discussed. We consider one-dimensional isothermal planar model atmospheres. We show that in the polarized line transfer calculations of subordinate lines, PRD plays as important of a role as it does in the case of resonance lines. We also study the effect of collisions on linear polarization profiles of subordinate lines.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2013

EXPLICIT-IMPLICIT SCHEME FOR RELATIVISTIC RADIATION HYDRODYNAMICS

We propose an explicit-implicit scheme for numerically solving special relativistic radiation hydrodynamic equations, which ensures a conservation of total energy and momentum (matter and radiation). In our scheme, zeroth and first moment equations of the radiation transfer equation are numerically solved without employing a flux-limited diffusion approximation. For an hyperbolic term, of which the timescale is the light crossing time when the flow velocity is comparable to the speed of light, is explicitly solved using an approximate Riemann solver. Source terms describing an exchange of energy and momentum between the matter and the radiation via the gas-radiation interaction are implicitly integrated using an iteration method. The implicit scheme allows us to relax the Courant-Friedrichs-Lewy condition in optically thick media, where heating/cooling and scattering timescales could be much shorter than the dynamical timescale. We show that our numerical code can pass test problems of one- and two-dimensional radiation energy transport, and one-dimensional radiation hydrodynamics. Our newly developed scheme could be useful for a number of relativistic astrophysical problems. We also discuss how to extend our explicit-implicit scheme to the relativistic radiation magnetohydrodynamics.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Recent development of Monte Carlo shell model and its application to no-core calculations

One of the major challenges in nuclear theory is to reproduce and to predict nuclear structure from ab initio calculations with realistic nuclear forces. As the current limitation of direct diagonalization of Hamiltonian matrices by Lanczos iteration method is around the order of matrix dimensionality 10 10 in shell-model calculations, it is difficult to access heavier nuclei beyond the p shell with sufficiently large basis spaces. It is possible to overcome this difficulty by utilizing efficient approximate methods to reproduce full ab initio solutions with good precision and quantified uncertainties. Following the major success of the Monte Carlo shell model (MCSM) with an assumed inert core in the sd- and pf-shell regions and also by recent developments in the MCSM algorithm, the no-core MCSM is expected to be one of the most powerful tools to meet these conditions. We have performed benchmark calculations in the p-shell region. Results of energies are compared with those in the full configuration interaction and no-core full configuration methods. These are found to be consistent with each other within quoted uncertainties when they could be quantified. We also compare and discuss the radial density of the helium-4 ground state extracted from the MCSM and FCI many-body wave functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2011

Electric quadrupole transitions of the Bohr Hamiltonian with the Morse potential

Eigenfunctions of the collective Bohr Hamiltonian with the Morse potential have been obtained by using the asymptotic iteration method (AIM) for both {gamma}-unstable and rotational structures. B(E2) transition rates have been calculated and compared to experimental data. Overall good agreement is obtained for transitions within the ground-state band, while some interband transitions appear to be systematically underpredicted in {gamma}-unstable nuclei and overpredicted in rotational nuclei.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Lagrangian displacement field estimators in cosmology

The late-time nonlinear Lagrangian displacement field is highly correlated with the initial field, so reconstructing it could enable us to extract primordial cosmological information. Our previous work carefully studied the displacement field reconstructed from the late-time density field using the iterative method proposed and found that it does not fully converge to the true, underlying displacement field (e.g., ~8 % offset at k ~ 0.2 h Mpc –1 at z = 0.6 ). We also constructed the Lagrangian perturbation theory model for the reconstructed field, but the model could not explain the discrepancy between the true and the reconstructed fields in the previous work. Here, the main sources of the discrepancy were speculated to be a numerical artifact in the displacement estimator due to the discreteness of the sample. In this paper, we develop two new estimators of the displacement fields to reduce such a numerical discreteness effect, the normalized momentum estimator and the rescaled resumed estimator. We show that the discrepancy reported is not due to the numerical artifacts. We conclude that the method cannot fully reconstruct the shape of the nonlinear displacement field at the redshift we studied, while it is still an efficient baryon acoustic oscillation reconstruction method. In parallel, by properly accounting for the UV-sensitive term in a reconstruction procedure with an effective field theory approach, we improve the theoretical model for the reconstructed displacement field, by almost 5 times, from ~15% to the level of a few percent at k ~ 0.2 h Mpc –1 at the redshift z = 0.6.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2018

Indigo: A Domain-Specific Language for Fast, Portable Image Reconstruction

© 2018 IEEE. Linear operators used in iterative methods like conjugate gradient have typically been implemented either as ''matrix-driven'' subroutines backed by explicit sparse or dense matrices, or as ''matrix-free'' subroutines that implement specific linear operations directly (e.g. FFTs). The matrix-driven approach is generally more portable because it can target widely-Available BLAS libraries, but it can be inefficient in terms of time and space complexity. In contrast, the matrix-free approach is more performant because it leverages structure in operations, but it requires each operator be re-implemented on each new platform. To increase performance and portability, we propose a hybrid approach that represents linear operators as expression trees. Leaf nodes in the tree are either matrix-free or matrix-driven operators, and interior nodes represent mathematical compositions (sums, products, transposes) or structural compositions (stacks, block diagonals, etc.) of the leaf operators. This representation enables expert-guided reordering and fusion transformations that can improve performance or reduce memory pressure. We implement our approach in a domain-specific language called Indigo. We assess Indigo on image reconstruction problems arising in four application areas: magnetic resonance imaging, ptychography, magnetic particle imaging, and fluorescent microscopy. We give performance results from vendor BLAS libraries, and we introduce specializations to Sparse BLAS routines that achieve near-Roofline performance on multi-core, many-core, and GPU systems.

Driscoll, Michael↗
DOE OSTI2015

Software-defined Radio Based Wireless Tomography: Experimental Demonstration and Verification

This letter presents an experimental demonstration of software-defined-radio-based wireless tomography using computer-hosted radio devices called Universal Software Radio Peripheral (USRP). This experimental brief follows our vision and previous theoretical study of wireless tomography that combines wireless communication and RF tomography to provide a novel approach to remote sensing. Automatic data acquisition is performed inside an RF anechoic chamber. Semidefinite relaxation is used for phase retrieval, and the Born iterative method is utilized for imaging the target. Experimental results are presented, validating our vision of wireless tomography.

Bonior, Jason D↗