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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,567 records · Page 87

All sources NASA NTRS DOE OSTI

NASA NTRS2012

Absolute Radiometric Calibration of Narrow-Swath Imaging Sensors with Reference to Non-Coincident Wide-Swath Sensors

An inter-calibration method is developed to provide absolute radiometric calibration of narrow-swath imaging sensors with reference to non-coincident wide-swath sensors. The method predicts at-sensor radiance using non-coincident imagery from the reference sensor and knowledge of spectral reflectance of the test site. The imagery of the reference sensor is restricted to acquisitions that provide similar view and solar illumination geometry to reduce uncertainties due to directional reflectance effects. Spectral reflectance of the test site is found with a simple iterative radiative transfer method using radiance values of a well-understood wide-swath sensor and spectral shape information based on historical ground-based measurements. At-sensor radiance is calculated for the narrow-swath sensor using this spectral reflectance and atmospheric parameters that are also based on historical in situ measurements. Results of the inter-calibration method show agreement on the 2 5 percent level in most spectral regions with the vicarious calibration technique relying on coincident ground-based measurements referred to as the reflectance-based approach. While the variability of the inter-calibration method based on non-coincident image pairs is significantly larger, results are consistent with techniques relying on in situ measurements. The method is also insensitive to spectral differences between the sensors by transferring to surface spectral reflectance prior to prediction of at-sensor radiance. The utility of this inter-calibration method is made clear by its flexibility to utilize image pairings with acquisition dates differing in excess of 30 days allowing frequent absolute calibration comparisons between wide- and narrow-swath sensors.

McCorkel, Joel↗

NASA NTRS1997

Iterative edge- and wavelet-based image registration of AVHRR and GOES satellite imagery

Most automatic registration methods are either correlation-based, feature-based, or a combination of both. Examples of features which can be utilized for automatic image registration are edges, regions, corners, or wavelet-extracted features. In this paper, we describe two proposed approaches, based on edge or edge-like features, which are very appropriate to highlight regions of interest such as coastlines. The two iterative methods utilize the Normalized Cross-Correlation of edge and wavelet features and are applied to such problems as image-to-map registration, landmarking, and channel-to-channel co-registration, utilizing test data, AVHRR data, as well as GOES image data.

LeMoigne, Jacqueline↗

DOE OSTI2013

Assessment of outdoor radiofrequency electromagnetic field exposure through hotspot localization using kriging-based sequential sampling

In this study, a novel methodology is proposed to create heat maps that accurately pinpoint the outdoor locations with elevated exposure to radiofrequency electromagnetic fields (RF-EMF) in an extensive urban region (or, hotspots), and that would allow local authorities and epidemiologists to efficiently assess the locations and spectral composition of these hotspots, while at the same time developing a global picture of the exposure in the area. Moreover, no prior knowledge about the presence of radiofrequency radiation sources (e.g., base station parameters) is required. After building a surrogate model from the available data using kriging, the proposed method makes use of an iterative sampling strategy that selects new measurement locations at spots which are deemed to contain the most valuable information—inside hotspots or in search of them—based on the prediction uncertainty of the model. The method was tested and validated in an urban subarea of Ghent, Belgium with a size of approximately 1 km{sup 2}. In total, 600 input and 50 validation measurements were performed using a broadband probe. Five hotspots were discovered and assessed, with maximum total electric-field strengths ranging from 1.3 to 3.1 V/m, satisfying the reference levels issued by the International Commission on Non-Ionizing Radiation Protection for exposure of the general public to RF-EMF. Spectrum analyzer measurements in these hotspots revealed five radiofrequency signals with a relevant contribution to the exposure. The radiofrequency radiation emitted by 900 MHz Global System for Mobile Communications (GSM) base stations was always dominant, with contributions ranging from 45% to 100%. Finally, validation of the subsequent surrogate models shows high prediction accuracy, with the final model featuring an average relative error of less than 2 dB (factor 1.26 in electric-field strength), a correlation coefficient of 0.7, and a specificity of 0.96. -- Highlights: • We present an iterative measurement and modeling method for outdoor RF-EMF exposure. • Hotspots are rapidly identified, and accurately characterized. • An accurate graphical representation, or heat map, is created, using kriging. • Random validation shows good correlation (0.7) and low relative errors (2 dB)

61 RADIATION PROTECTION AND DOSIMETRY↗

NASA NTRS1977

Calculation of transonic potential flowfields about complex, three-dimensional configurations

Methods for extending iterative, finite-difference calculations of transonic potential flowfields to complex three-dimensional configurations are discussed. One particularly attractive approach is to use relatively simple conformal mappings in combination with shearing transformations to generate computational domains that are nearly-conformally mapped from the physical space in one family of coordinate surfaces, and which map the complex boundaries to grid surfaces. The application of such a method to a general wing-body combination or to a multi-bladed fan is discussed. A transformation to map the wing-fuselage or fan-hub combination to a convenient computational domain is proposed. The transformation is useful in its own right for treating the two-dimensional problems of flow past a profile in a wind tunnel or through a cascade. Some results of preliminary calculations are presented.

Caughey, D. A.↗

NASA NTRS2015

A Principled Approach to the Specification of System Architectures for Space Missions

Modern space systems are increasing in complexity and scale at an unprecedented pace. Consequently, innovative methods, processes, and tools are needed to cope with the increasing complexity of architecting these systems. A key systems challenge in practice is the ability to scale processes, methods, and tools used to architect complex space systems. Traditionally, the process for specifying space system architectures has largely relied on capturing the system architecture in informal descriptions that are often embedded within loosely coupled design documents and domain expertise. Such informal descriptions often lead to misunderstandings between design teams, ambiguous specifications, difficulty in maintaining consistency as the architecture evolves throughout the system development life cycle, and costly design iterations. Therefore, traditional methods are becoming increasingly inefficient to cope with ever-increasing system complexity. We apply the principles of component-based design and platform-based design to the development of the system architecture for a practical space system to demonstrate feasibility of our approach using SysML. Our results show that we are able to apply a systematic design method to manage system complexity, thus enabling effective data management, semantic coherence and traceability across different levels of abstraction in the design chain. Just as important, our approach enables interoperability among heterogeneous tools in a concurrent engineering model based design environment.

platform-based design↗

NASA NTRS1989

Domain decomposition methods for systems of conservation laws: Spectral collocation approximations

Hyperbolic systems of conversation laws are considered which are discretized in space by spectral collocation methods and advanced in time by finite difference schemes. At any time-level a domain deposition method based on an iteration by subdomain procedure was introduced yielding at each step a sequence of independent subproblems (one for each subdomain) that can be solved simultaneously. The method is set for a general nonlinear problem in several space variables. The convergence analysis, however, is carried out only for a linear one-dimensional system with continuous solutions. A precise form of the error reduction factor at each iteration is derived. Although the method is applied here to the case of spectral collocation approximation only, the idea is fairly general and can be used in a different context as well. For instance, its application to space discretization by finite differences is straight forward.

Quarteroni, Alfio↗

NASA NTRS2004

Aerodynamic Optimization of Rocket Control Surface Geometry Using Cartesian Methods and CAD Geometry

Aerodynamic design is an iterative process involving geometry manipulation and complex computational analysis subject to physical constraints and aerodynamic objectives. A design cycle consists of first establishing the performance of a baseline design, which is usually created with low-fidelity engineering tools, and then progressively optimizing the design to maximize its performance. Optimization techniques have evolved from relying exclusively on designer intuition and insight in traditional trial and error methods, to sophisticated local and global search methods. Recent attempts at automating the search through a large design space with formal optimization methods include both database driven and direct evaluation schemes. Databases are being used in conjunction with surrogate and neural network models as a basis on which to run optimization algorithms. Optimization algorithms are also being driven by the direct evaluation of objectives and constraints using high-fidelity simulations. Surrogate methods use data points obtained from simulations, and possibly gradients evaluated at the data points, to create mathematical approximations of a database. Neural network models work in a similar fashion, using a number of high-fidelity database calculations as training iterations to create a database model. Optimal designs are obtained by coupling an optimization algorithm to the database model. Evaluation of the current best design then gives either a new local optima and/or increases the fidelity of the approximation model for the next iteration. Surrogate methods have also been developed that iterate on the selection of data points to decrease the uncertainty of the approximation model prior to searching for an optimal design. The database approximation models for each of these cases, however, become computationally expensive with increase in dimensionality. Thus the method of using optimization algorithms to search a database model becomes problematic as the number of design variables is increased.

Nelson, Andrea↗

DOE OSTI2015

Iterative stability analysis of spatial domain decomposition based on block Jacobi algorithm for the diamond-difference scheme

Here, we study convergence of the integral transport matrix method (ITMM) based on a parallel block Jacobi (PBJ) iterative strategy for solving particle transport problems. The ITMM is a spatial domain decomposition method proposed for massively parallel computations. AFourier analysis of the PBJ-based iterations applied to SN diamond-difference equations in 1D slab and 2D Cartesian geometries is performed. It is carried out for infinite-medium problems with homogeneous material properties. To analyze the performance of the ITMM with the PBJ algorithm and evaluate its potential in scalability we consider a limiting case of one spatial cell per subdomain. The analysis shows that in such limit the spectral radius of the iteration method is one without regard to values of the scattering ratio and optical thickness of the spatial cells. This implies lack of convergence in infinite medium. Numerical results of finite-medium problems are presented. They demonstrate effects of finite size of spatial domain on the performance of the iteration algorithm as well as its asymptotic behavior when the extent of the spatial domain increases. Finally, these numerical experiments also show that for finite domains iterative convergence to a finite criterion is achievable in a multiple of the sum of number of cells in each dimension.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

DOE OSTI2018

A dissipation-based method for improving the accuracy of computational fluid dynamics simulations of high level non-Newtonian wastes

Here this paper investigates a new viscosity model for laminar, transitional, and turbulent pipe flow simulations of a non-Newtonian slurry using the Star-CCM+ software. Herschel-Bulkley (HB) rheology equation was used in its original and modified forms to model the fluid viscosity. The realizable k-ε model was used for modeling of turbulence. It was observed that simulations using the HB model predicted inaccurate profiles of U (mean velocity) and U + (mean velocity in wall units). Therefore, we developed alternative values of viscosity obtained directly from non-linearity of stress–strain variation (α) in each iteration. The developed method, “Alpha method”, applied changes in places where dissipation rate of turbulent kinetic energy occurred and was designed based on local or global dissipation rate values. Inverse values of α were used to test effects of more viscosity reduction for all flow regimes. In prediction of the mean velocities, significant improvements were obtained through global adjustment of viscosity in inverse form during simulation of the laminar and turbulent flow cases. For the transitional flow, obtaining improvements in prediction of U and U + from the same model version faced challenges. While the global adjustments in direct form were found significantly effective in improving the prediction of U, they failed to predict the U + accurately. For this purpose, Local viscosity adjustments were found significantly effective.

42 ENGINEERING↗

DOE OSTI2014

Density cumulant functional theory from a unitary transformation: N-representability, three-particle correlation effects, and application to O{sub 4}{sup +}

A new approach to density cumulant functional theory is developed that derives density cumulant N-representability conditions from an approximate Fock space unitary transformation. We present explicit equations for the third- and fourth-order two-particle cumulant N-representability, as well as the second-order contributions that depend on the connected three-particle density cumulant. These conditions are used to formulate the ODC-13 method and the non-iterative (λ{sub 3}) correction that employ an incomplete description of the fourth-order two-particle cumulant N-representability and the second-order three-particle correlation effects, respectively. We perform an analysis of the ODC-13 N-representability description for the dissociation of H{sub 2} and apply the ODC-13 method and the (λ{sub 3}) correction to diatomic molecules with multiple bond character and the symmetry-breaking tetraoxygen cation (O{sub 4}{sup +}). For the O{sub 4}{sup +} molecule, the vibrational frequencies of the ODC-13(λ{sub 3}) method do not exhibit spatial symmetry breaking and are in a good agreement with the recent infrared photodissociation experiment. We report the O{sub 4}{sup +} equilibrium structure, harmonic frequencies, and dissociation energy computed using ODC-13(λ{sub 3}) with a diffuse, core-correlated aug-cc-pCVTZ basis set.

74 ATOMIC AND MOLECULAR PHYSICS↗

NASA NTRSDate not supplied

History and Status of the CDISC Aerodynamic Design Method

This paper will review the development and application of the CDISC aerodynamic design method that began in the early 1980s and continues in active use today. The method uses an iterative, knowledge-based approach that has been coupled with numerous flow solvers utilizing a variety of grid types and levels of flow physics. It has been applied to configurations varying in complexity from 2-D airfoils to full 3-D aircraft at flow conditions ranging from low-speed, high-lift to supersonic cruise. The knowledge-based approach provides a very rapid design process, with designs typically completed in approximately 1-2x the time required for an initial converged analysis of the baseline. The evolution of the method as it tracked the expansive growth of grid and flow solver technology, along with some major applications over 5 decades, are given, followed by a description of the philosophy and process used in design. Finally, an example of its use in a preliminary design environment is given.

CDISC↗

NASA NTRSDate not supplied

History and Status of the CDISC Aerodynamic Design Method

CDISC↗

DOE OSTI2013

Excitation energies with cost-reduced variant of the active-space EOMCCSDT method: the EOMCCSDt-(3) over-bar approach

In this paper we discuss the performance of the several simplified variants of equation-of-motion coupled cluster method (EOMCC) with iterative inclusion of singles, doubles and active-space triples (EOMCCSDt). In particular, we explore simplified EOMCCSDt approaches which enable one to generate the triply excited amplitudes in on-the-fly manner. The original EOMCCSDt formulation has already demonstrated a great success in encapsulating the most important excited-state correlation effects due to triples. In analogy to the original EOMCCSDT formulation, the proposed approach can by-pass the typical bottlenecks associated with the need for storing triply-excited amplitudes. In this paper, we illustrate the performance of several approximate EOMCCSDt methods, named EOMCCSDt-3 and EOMCCSdt-3x, on typical benchmark systems including C2, N2, and the ozone molecules. These new methods yield excitation energies close to the EOMCCSDt ones. The extrapolation of excitation energies for basis sets ranging from cc-pVDZ to cc-pV6Z for N2 and C2 shows very good convergence to the experimental results for states dominated by single excitations.

Coupled cluster↗

DOE OSTI2011

Exact exchange potential evaluated from occupied Kohn-Sham and Hartree-Fock solutions

The reported algorithm determines the exact exchange potential v{sub x} in an iterative way using energy shifts (ESs) and orbital shifts (OSs) obtained with finite-difference formulas from the solutions (occupied orbitals and their energies) of the Hartree-Fock-like equation and the Kohn-Sham-like equation, the former used for the initial approximation to v{sub x} and the latter for increments of ES and OS due to subsequent changes of v{sub x}. Thus, the need for solution of the differential equations for OSs, used by Kuemmel and Perdew [Phys. Rev. Lett. 90, 043004 (2003)], is bypassed. The iterated exchange potential, expressed in terms of ESs and OSs, is improved by modifying ESs at odd iteration steps and OSs at even steps. The modification formulas are related to the optimized-effective-potential equation (satisfied at convergence) written as the condition of vanishing density shift (DS). They are obtained, respectively, by enforcing its satisfaction through corrections to approximate OSs and by determining the optimal ESs that minimize the DS norm. The proposed method, successfully tested for several closed-(sub)shell atoms, from Be to Kr, within the density functional theory exchange-only approximation, proves highly efficient. The calculations using the pseudospectral method for representing orbitals give iterative sequences of approximate exchange potentials (starting with the Krieger-Li-Iafrate approximation) that rapidly approach the exact v{sub x} so that, for Ne, Ar, and Zn, the corresponding DS norm becomes less than 10{sup -6} after 13, 13, and 9 iteration steps for a given electron density. In self-consistent density calculations, orbital energies of 10{sup -4} hartree accuracy are obtained for these atoms after, respectively, 9, 12, and 12 density iteration steps, each involving just two steps of v{sub x} iteration, while the accuracy limit of 10{sup -6} to 10{sup -7} hartree is reached after 20 density iterations.

74 ATOMIC AND MOLECULAR PHYSICS↗

DOE OSTI2015

Investigation and quantification of nonlinearity using surrogate data in a glow discharge plasma

Detection of nonlinearity has been carried out in periodic and aperiodic floating potential fluctuations of DC glow discharge plasma by generating surrogate data using iterative amplitude adjusted Fourier transform method. We introduce “delay vector variance” analysis (DVV) for the first time, which allows reliable detection of nonlinearity and provides some easy to interpret diagram conveying information about the nature of the experimental floating potential fluctuations (FPF). The method of false nearest neighbourhood is deployed on the FPF's to find a good embedding so as to be acquainted with the precise knowledge of m, which is desirable for carrying out DVV analysis. The emergence of nonlinearity with increase in discharge voltage has been ensured by taking into consideration the total energy present in different band of frequencies excited due to nonlinear processes. Rejection of null hypothesis has been verified by performing the rank test method that confirms the presence of nonlinearity quantitatively.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2026

A Scalable Reduced‐Order Model for the Steady Navier–Stokes Equations

Scaling up new scientific technologies from laboratory to industry often involves demonstrating performance on a larger scale. Computer simulations can accelerate design and predictions in the deployment process, though traditional numerical methods are computationally intractable even for intermediate pilot plant scales. Recently, the component reduced order modeling method has been developed to tackle this challenge by combining projection reduced order modeling and discontinuous Galerkin domain decomposition. However, while many scientific or engineering applications involve nonlinear physics, this method has only been demonstrated for various linear systems. In this work, the component reduced order modeling method is extended to steady Navier–Stokes flow, with application to general nonlinear physics in view. The large‐scale, global domain is decomposed into a combination of small‐scale unit component. Linear subspaces for flow velocity and pressure are identified via proper orthogonal decomposition over sample snapshots collected from each small‐scale unit component. Velocity bases are augmented with a pressure supremizer to satisfy the inf–sup condition for stable pressure prediction. Two different nonlinear reduced order modeling methods are employed and compared for efficient evaluation of nonlinear advection: A third‐order tensor projection operator and the empirical quadrature procedure. The proposed method is demonstrated on the flow over arrays of five different unit objects, achieving a 23‐fold speedup with less than 4% relative error in domains up to 256 times larger than the unit components. Furthermore, a numerical experiment with the pressure supremizer strongly indicates the need for a supremizer for stable pressure prediction. A comparison between the tensorial approach and the empirical quadrature procedure revealed a slight advantage of the empirical quadrature procedure. The framework is compared with an alternating Schwarz‐based reduced‐order approach, demonstrating improved efficiency and robustness for the DG‐based global solver while retaining flexibility for sub‐scale iterative solvers. The method is further extended to a coupled advection–diffusion and Navier–Stokes system, illustrating its applicability to multi‐physics problems and its potential for more general, inter‐coupled nonlinear systems.

42 ENGINEERING↗

DOE OSTI2018

Distributed Link Removal Using Local Estimation of Network Topology

In this paper, we consider the problem of network structure manipulation in the absence of information on the global network topology. In particular, the problem of removing some of links is investigated in order to slow or stop the spread of disease in a network while preserving its connectivity. Existing methods solve this combinatorial problem in a centralized manner and they require the global information of network structure. In this paper, we propose a distributed design algorithm to compute a suboptimal solution to this problem efficiently by mimicking gradient based method, namely by iteratively removing one or multiple links at a time from the network. Specifically, using matrix perturbation analysis we formulate an optimization problem involving the eigenvector associated with the largest eigenvalue of the adjacency matrix and whose solution is equal to a suboptimal solution to the original problem. This strategy also enables us to overcome the combinatorial issue of the problem. Distributed algorithms to estimate the eigenvector and to verify network’s connectivity are then proposed which facilitate us to solve the new optimization problem. In addition, topological insights into the proposed algorithm and optimality of its solution are also discussed. Lastly, the proposed distributed design method is demonstrated and evaluated via several numerical examples.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2015

On the Use of Finite Difference Matrix-Vector Products in Newton-Krylov Solvers for Implicit Climate Dynamics with Spectral Elements

Efficient solution of global climate models requires effectively handling disparate length and time scales. Implicit solution approaches allow time integration of the physical system with a time step dictated by accuracy of the processes of interest rather than by stability governed by the fastest of the time scales present. Implicit approaches, however, require the solution of nonlinear systems within each time step. Usually, a Newton s method is applied for these systems. Each iteration of the Newton s method, in turn, requires the solution of a linear model of the nonlinear system. This model employs the Jacobian of the problem-defining nonlinear residual, but this Jacobian can be costly to form. If a Krylov linear solver is used for the solution of the linear system, the action of the Jacobian matrix on a given vector is required. In the case of spectral element methods, the Jacobian is not calculated but only implemented through matrix-vector products. The matrix-vector multiply can also be approximated by a finite-difference which may show a loss of accuracy in the overall nonlinear solver. In this paper, we review the advantages and disadvantages of finite-difference approximations of these matrix-vector products for climate dynamics within the spectral-element based shallow-water dynamical-core of the Community Atmosphere Model (CAM).

climate modeling↗