Search NASA⌕ Search

SEARCH · Search NASA

Results for “Local Method”

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 307 records · Page 17

Reliability enhancement of Navier-Stokes codes through convergence enhancement

Reduction of total computing time required by an iterative algorithm for solving Navier-Stokes equations is an important aspect of making the existing and future analysis codes more cost effective. Several attempts have been made to accelerate the convergence of an explicit Runge-Kutta time-stepping algorithm. These acceleration methods are based on local time stepping, implicit residual smoothing, enthalpy damping, and multigrid techniques. Also, an extrapolation procedure based on the power method and the Minimal Residual Method (MRM) were applied to the Jameson's multigrid algorithm. The MRM uses same values of optimal weights for the corrections to every equation in a system and has not been shown to accelerate the scheme without multigriding. Our Distributed Minimal Residual (DMR) method based on our General Nonlinear Minimal Residual (GNLMR) method allows each component of the solution vector in a system of equations to have its own convergence speed. The DMR method was found capable of reducing the computation time by 10-75 percent depending on the test case and grid used. Recently, we have developed and tested a new method termed Sensitivity Based DMR or SBMR method that is easier to implement in different codes and is even more robust and computationally efficient than our DMR method.

Choi, K.-Y.↗

Underlying Fundamentals of Kalman Filtering for River Network Modeling

The grand challenge of producing hydrometeorological estimates every time and everywhere has motivated the fusion of sparse observations with dense numerical models, with a particular interest on discharge in river modeling. Ensemble methods are largely preferred as they enable the estimation of error properties, but at the expense of computational load and generally with underestimations. These imperfect stochastic estimates motivate the use of correction methods, that is, error localization and inflation, although the physical justifications for their optimality are limited. The purpose of this study is to use one of the simplest forms of data assimilation when applied to river modeling and reveal the underlying mechanisms impacting its performance. Our framework based on assimilating daily averaged in situ discharge measurements to correct daily averaged runoff was tested over a 4-yr case study of two rivers in Texas. Results show that under optimal conditions of inflation and localization, discharge simulations are consistently improved such that the mean values of Nash–Sutcliffe efficiency are enhanced from211.32 to 0.55 at observed gauges and from212.24 to21.10 at validation gauges. Yet, parameters controlling the inflation and the localization have a large impact on the performance. Further investigations of these sensitivities showed that optimal inflation occurs when compensating exactly for discrepancies in the magnitude of errors while optimal localization matches the distance traveled during one assimilation window. These results may be applicable to more advanced data assimilation methods as well as for larger applications motivated by upcoming river-observing satellite missions, such as NASA’s Surface Water and Ocean Topography mission.

Streamflow↗

An Amplitude-Based Estimation Method for International Space Station (ISS) Leak Detection and Localization Using Acoustic Sensor Networks

The development of a robust and efficient leak detection and localization system within a space station environment presents a unique challenge. A plausible approach includes the implementation of an acoustic sensor network system that can successfully detect the presence of a leak and determine the location of the leak source. Traditional acoustic detection and localization schemes rely on the phase and amplitude information collected by the sensor array system. Furthermore, the acoustic source signals are assumed to be airborne and far-field. Likewise, there are similar applications in sonar. In solids, there are specialized methods for locating events that are used in geology and in acoustic emission testing that involve sensor arrays and depend on a discernable phase front to the received signal. These methods are ineffective if applied to a sensor detection system within the space station environment. In the case of acoustic signal location, there are significant baffling and structural impediments to the sound path and the source could be in the near-field of a sensor in this particular setting.

Tian, Jialin↗

Sound propagation in a duct of periodic wall structure

A boundary condition, which accounts for the coupling in the sections behind the duct boundary, is given for the sound-absorbing duct with a periodic structure of the wall lining and using regular partition walls. The soundfield in the duct is suitably described by the method of differences. For locally active walls this renders an explicit approximate solution for the propagation constant. Coupling may be accounted for by the method of differences in a clear manner. Numerical results agree with measurements and yield information which has technical applications.

Kurze, U.↗

The theoretical accuracy of Runge-Kutta time discretizations for the initial boundary value problem: A careful study of the boundary error

The conventional method of imposing time dependent boundary conditions for Runge-Kutta (RK) time advancement reduces the formal accuracy of the space-time method to first order locally, and second order globally, independently of the spatial operator. This counter intuitive result is analyzed in this paper. Two methods of eliminating this problem are proposed for the linear constant coefficient case: (1) impose the exact boundary condition only at the end of the complete RK cycle, (2) impose consistent intermediate boundary conditions derived from the physical boundary condition and its derivatives. The first method, while retaining the RK accuracy in all cases, results in a scheme with much reduced CFL condition, rendering the RK scheme less attractive. The second method retains the same allowable time step as the periodic problem. However it is a general remedy only for the linear case. For non-linear hyperbolic equations the second method is effective only for for RK schemes of third order accuracy or less. Numerical studies are presented to verify the efficacy of each approach.

Carpenter, Mark H.↗

Rigorous numerical treatment of the no-slip condition in a vorticity formulation

It is shown that a rigorous mathematical treatment of the boundary conditions for rigid no-slip surfaces requires, in general, that two passes be made in solving Poisson's equation. A Fourier solution procedure was the key to this realization. It is also shown that a corresponding method exists for finite and, for the same local approximation, the usual one step method gives the same result for the geometry considered here. The more fundamental two step method gives a new insight into the physical significance of boundary vorticity.

Gazdag, J.↗

Mode localization in multi-span beams

The influence of numerous effects on mode localization in multi-span beams is investigated. Finite-element methods are used to study localization as a function of: Timoshenko beam effects; beam end conditions; span length, mass, and stiffness imperfection; viscous damping; axial force; transverse support and rotational coupling stiffness; and modeling resolution. Three configurations are studied, commencing with two different two-span models, and culminating in a ten-span configuration resembling lattice-type large space structures. Results indicate that, in addition to the ratio of imperfection to coupling stiffness being an important localization parameter, transverse support stiffness and Timoshenko beam effects greatly affect the tendency of a structure to exhibit localized modes.

Lust, S. D.↗

The Effect of Fiber Strength Stochastics and Local Fiber Volume Fraction on Multiscale Progressive Failure of Composites

Continuous fiber unidirectional polymer matrix composites (PMCs) can exhibit significant local variations in fiber volume fraction as a result of processing conditions that can lead to further local differences in material properties and failure behavior. In this work, the coupled effects of both local variations in fiber volume fraction and the empirically-based statistical distribution of fiber strengths on the predicted longitudinal modulus and local tensile strength of a unidirectional AS4 carbon fiber/ Hercules 3502 epoxy composite were investigated using the special purpose NASA Micromechanics Analysis Code with Generalized Method of Cells (MAC/GMC); local effective composite properties were obtained by homogenizing the material behavior over repeating units cells (RUCs). The predicted effective longitudinal modulus was relatively insensitive to small (~8%) variations in local fiber volume fraction. The composite tensile strength, however, was highly dependent on the local distribution in fiber strengths. The RUC-averaged constitutive response can be used to characterize lower length scale material behavior within a multiscale analysis framework that couples the NASA code FEAMAC and the ABAQUS finite element solver. Such an approach can be effectively used to analyze the progressive failure of PMC structures whose failure initiates at the RUC level. Consideration of the effect of local variations in constituent properties and morphologies on progressive failure of PMCs is a central aspect of the application of Integrated Computational Materials Engineering (ICME) principles for composite materials.

multiscale analysis↗

A multilevel control system for the large space telescope

A multilevel scheme was proposed for control of Large Space Telescope (LST) modeled by a three-axis-six-order nonlinear equation. Local controllers were used on the subsystem level to stabilize motions corresponding to the three axes. Global controllers were applied to reduce (and sometimes nullify) the interactions among the subsystems. A multilevel optimization method was developed whereby local quadratic optimizations were performed on the subsystem level, and global control was again used to reduce (nullify) the effect of interactions. The multilevel stabilization and optimization methods are presented as general tools for design and then used in the design of the LST Control System. The methods are entirely computerized, so that they can accommodate higher order LST models with both conceptual and numerical advantages over standard straightforward design techniques.

D D Siljak↗

Efficient reanalysis of structures by a direct modification method

A procedure for the local stiffness modifications of large structures is described. It enables structural modifications without an a priori definition of the changes in the original structure and without loss of efficiency due to multiple loading conditions. The solution procedure, implemented in NASTRAN, involved the decomposed stiffness matrix and the displacement vectors of the original structure. It solves the modified structure exactly, irrespective of the magnitude of the stiffness changes. In order to investigate the efficiency of the present procedure and to test its applicability within a design environment, several real and large structures were solved. The results of the efficiency studies indicate that the break-even point of the procedure varies between 8% and 60% stiffness modifications, depending upon the structure's characteristics and the options employed.

Raibstein, A. I.↗

Grid generation for time dependent problems: Criteria and methods

The problem of generating local mesh refinements when solving time dependent partial differential equations was examined. The problem of creating an appropriate grid, given a mesh function h defined over the spatial domain is discussed. A data structure which permits efficient use of the resulting grid is described. A good choice for h is an estimate of the local truncation error, and several ways to estimate it are discussed. The efficiency and implementation problems of these error estimates were compared.

Berger, M.↗

Viscoplastic response of structures for intense local heating

A thermoviscoplastic finite element method employing the Bodner-Partom constitutive model is used to investigate the response of simplified thermal-structural models to intense local heating. With rapid rises of temperature, the nickel alloy structures display initially higher yield stresses due to strain rate effects. As temperatures approach elevated values, yield stress and stiffness degrade rapidly and pronounced plastic deformation occurs.

Thornton, Earl A.↗

Free and forced response of nearly periodic multi-span beams and multi-bay trusses

The influence of several effects on mode localization in free-free multispan beams is investigated. Finite-element methods are used to study localization as a function of Timoshenko beam effects, imperfection strength, coupling strength, and imperfection type. The multispan beam study is conducted with random realizations of a forty-span free-free beam. In addition, free and forced response of free-free trusses and frames is investigated. Results indicate that, for the trusses and frames considered, the dynamic behavior can be approximated by continuum modeling approaches for only the lowest few global modes. Also, for reasonable values of modal damping, localized modes play a subtle role in the overall structural response.

Lust, S. D.↗

Non-LTE line-blanketed model atmospheres of hot stars. 1: Hybrid complete linearization/accelerated lambda iteration method

A new munerical method for computing non-Local Thermodynamic Equilibrium (non-LTE) model stellar atmospheres is presented. The method, called the hybird complete linearization/accelerated lambda iretation (CL/ALI) method, combines advantages of both its constituents. Its rate of convergence is virtually as high as for the standard CL method, while the computer time per iteration is almost as low as for the standard ALI method. The method is formulated as the standard complete lineariation, the only difference being that the radiation intensity at selected frequency points is not explicity linearized; instead, it is treated by means of the ALI approach. The scheme offers a wide spectrum of options, ranging from the full CL to the full ALI method. We deonstrate that the method works optimally if the majority of frequency points are treated in the ALI mode, while the radiation intensity at a few (typically two to 30) frequency points is explicity linearized. We show how this method can be applied to calculate metal line-blanketed non-LTE model atmospheres, by using the idea of 'superlevels' and 'superlines' introduced originally by Anderson (1989). We calculate several illustrative models taking into accont several tens of thosands of lines of Fe III to Fe IV and show that the hybrid CL/ALI method provides a robust method for calculating non-LTE line-blanketed model atmospheres for a wide range of stellar parameters. The results for individual stellar types will be presented in subsequent papers in this series.

Hubeny, I.↗

Accuracy/Computation Performance of a New Trilateration Scheme for GPS-Style Localization

We recently introduced a new geometric trilateration (GT) method for GPS-style positioning. Preliminary singlepoint analysis using simplistic error assumptions indicates that the new scheme delivers almost indistinguishable localization accuracy as the traditional Newton-Raphson (NR) approach. Also, the same computation procedure can be used to perform high-accuracy relative positioning between a reference vehicle and an arbitrary number of target vehicles. This scheme has the potential to enable a) new mission concepts in collaborative science, b) in-situ navigation services for human Mars missions, and c) lower cost and faster acquisition of GPS signals for consumer-grade GPS products. The new GT scheme differs from the NR scheme as follows: 1. The new scheme is derived from Pythagoras Theorem, whereas the NR method is based on the principle of linear regression. 2. The NR method uses the absolute locations (xi, yi, zi)’s of the GPS satellites as input to each step of the localization computation. The GT method uses the Directional Cosines Ui’s from Earth’s center to the GPS satellite Si. 3. Both the NR method and the GT method iterate to converge to a localized solution. In each iteration step, multiple matrix operations are performed. The NR method constructs a different matrix in each iterative step, thus requires performing a new set of matrix operations in each step. The GT scheme uses the same matrix in each iteration, thus requiring computing the matrix operations only once for all subsequent iterations. In this paper, we perform an in-depth comparison between the GT scheme and the NR method in terms of a) GPS localization accuracy in the GPS operation environment, b) its sensitivity with respect to systematic errors and random errors, and c) computation load required to converge to a localization solution.

Lee, Charles↗

Simplified method for measuring mean opacities and effective absorption coefficients.

Use of wide bandwidths and a number of path lengths to obtain mean opacities and effective absorption coefficients. This method is applicable if the temperatures and pressures at which these qualities are desired can be duplicated experimentally. Additional restrictions are that the gas in the experiments must be essentially isothermal, homogeneous, nonscattering, and in local thermodynamic equilibrium. The method is most advantageous in shock tube experiments and experiments utilizing a pulsed arc of cylindrical symmetry viewed axially.

Patch, R. W.↗

Local wall heat flux/temperature meter for convective flow and method of utilizing same

According to one embodiment of the invention, a method includes providing a conduit having a fluid flowing therethrough, disposing a plurality of temperature measurement devices inside a wall of the conduit, positioning at least some of the temperature measurement devices proximate an inside surface of the wall of the conduit, positioning at least some of the temperature measurement devices at different radial positions at the same circumferential location within the wall, measuring a plurality of temperatures of the wall with respective ones of the temperature measurement devices to obtain a three-dimensional temperature topology of the wall, determining the temperature dependent thermal conductivity of the conduit, and determining a multi-dimensional thermal characteristic of the inside surface of the wall of the conduit based on extrapolation of the three-dimensional temperature topology and the temperature dependent thermal conductivities.

Boyd, Ronald D.↗

Ballistic-Electron-Emission Microscopy Techniques for Nanometer-scale Characterization of Interfaces

Semiconductor interface properties are among the most important phenomena in materials science and technology. The study of metal/semiconductor Schottky barrier interfaces has been the primary focus of a large research and development community for decades. Throughout the long history of interface investigation, the study of interface defect electronic properties have been seriously hindered by the fundamental experimental difficulty of probing subsurface structures. A new method, Ballistic-Electron-Emission Microscopy (BEEM), has been developed which not only enables spectroscopic probing of subsurface interface properties, but also, provides nanometer-resolution imaging capabilities. BEEM employs Scanning Tunneling Microscopy (STM) and a unique spatially localized ballistic electron spectroscopy method...

Bell, L. D.↗