SEARCH · Engineering Papers

Results for “numerical continuation”

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 793 records · Page 44

All sources NASA NTRS DOE OSTI

NASA NTRS1981

Flow and far field acoustic amplification properties of heated and unheated jets

The interaction of an acoustic pulse with the experimentally determined mean flow field of a spreading jet is simulated numerically. The simulation is obtained through solving the Euler equations linearized about the spreading jet. The model reveals a small, sustained oscillation long after the original pulse has passed. This remnant is considered a continual shedding of vortices from the nozzle lip, together with the generation of acoustic ripples. IT is shown that the jet also acts as an amplifier of sound. This amplification is traced to the jet's stability characteristics. It is demonstrated that some of the observed differences in the spectra of heated and unheated jets can be attributed to differences in the stability characteristics of the jets.

Maestrello, L.↗

NASA NTRS1986

Study of a homotopy continuation method for early orbit determination with the Tracking and Data Relay Satellite System (TDRSS)

A recent mathematical technique for solving systems of equations is applied in a very general way to the orbit determination problem. The study of this technique, the homotopy continuation method, was motivated by the possible need to perform early orbit determination with the Tracking and Data Relay Satellite System (TDRSS), using range and Doppler tracking alone. Basically, a set of six tracking observations is continuously transformed from a set with known solution to the given set of observations with unknown solutions, and the corresponding orbit state vector is followed from the a priori estimate to the solutions. A numerical algorithm for following the state vector is developed and described in detail. Numerical examples using both real and simulated TDRSS tracking are given. A prototype early orbit determination algorithm for possible use in TDRSS orbit operations was extensively tested, and the results are described. Preliminary studies of two extensions of the method are discussed: generalization to a least-squares formulation and generalization to an exhaustive global method.

Smith, R. L.↗

NASA NTRS1989

Fourier transform methods in local gravity modeling

New algorithms were derived for computing terrain corrections, all components of the attraction of the topography at the topographic surface and the gradients of these attractions. These algoriithms utilize fast Fourier transforms, but, in contrast to methods currently in use, all divergences of the integrals are removed during the analysis. Sequential methods employing a smooth intermediate reference surface were developed to avoid the very large transforms necessary when making computations at high resolution over a wide area. A new method for the numerical solution of Molodensky's problem was developed to mitigate the convergence difficulties that occur at short wavelengths with methods based on a Taylor series expansion. A trial field on a level surface is continued analytically to the topographic surface, and compared with that predicted from gravity observations. The difference is used to compute a correction to the trial field and the process iterated. Special techniques are employed to speed convergence and prevent oscillations. Three different spectral methods for fitting a point-mass set to a gravity field given on a regular grid at constant elevation are described. Two of the methods differ in the way that the spectrum of the point-mass set, which extends to infinite wave number, is matched to that of the gravity field which is band-limited. The third method is essentially a space-domain technique in which Fourier methods are used to solve a set of simultaneous equations.

Harrison, J. C.↗

NASA NTRS1989

Deformation analysis of boron/aluminum specimens by moire interferometry

Whole-field surface deformations were measured for two slotted tension specimens from multiply laminates, one with 0 deg fiber orientation in the surface ply and the other with 45 deg orientation. Macromechanical and micromechanical details were revealed using high-sensitivity moire interferometry. Although global deformations of all plies were essentially equal, numerous random or anomalous features were observed. Local deformations of adjacent 0 deg and 45 deg plies were very different, both near the slot and remote from it, requiring large interlaminar shear strains for continuity. Shear strains were concentrated in the aluminum matrix. For 45 deg plies, a major portion of the deformation was by shear; large plastic slip of matrix occurred at random locations in 45 deg plies, wherein groups of fibers slipped relative to other groups. Shear strains in the interior, between adjacent fibers, were larger than the measured surface strains.

Post, Daniel↗

NASA NTRS1990

Turbulence in astrophysics

Some of the astrophysical scenarios in which turbulence plays an important role are discussed in view of the comparative advantages of currently available models of turbulence phenomena; attention is given to a specific model that has undergone continuous refinement since 1985. The desideratum in this inquiry is a turbulence model which incorporates the best features of an a priori deterministic model, as well as a redundant set of results from full numerical simulations for a wide variety of turbulent flows; there should also be a simplification of the former, and an enlargement of the complexities of the latter.

Canuto, V. M.↗

NASA NTRS1993

A Spectral Algorithm for Envelope Reduction of Sparse Matrices

The problem of reordering a sparse symmetric matrix to reduce its envelope size is considered. A new spectral algorithm for computing an envelope-reducing reordering is obtained by associating a Laplacian matrix with the given matrix and then sorting the components of a specified eigenvector of the Laplacian. This Laplacian eigenvector solves a continuous relaxation of a discrete problem related to envelope minimization called the minimum 2-sum problem. The permutation vector computed by the spectral algorithm is a closest permutation vector to the specified Laplacian eigenvector. Numerical results show that the new reordering algorithm usually computes smaller envelope sizes than those obtained from the current standard algorithms such as Gibbs-Poole-Stockmeyer (GPS) or SPARSPAK reverse Cuthill-McKee (RCM), in some cases reducing the envelope by more than a factor of two.

Barnard, Stephen T.↗

NASA NTRS1998

Algorithms and Libraries

This exploratory study initiated our inquiry into algorithms and applications that would benefit by latency tolerant approach to algorithm building, including the construction of new algorithms where appropriate. In a multithreaded execution, when a processor reaches a point where remote memory access is necessary, the request is sent out on the network and a context--switch occurs to a new thread of computation. This effectively masks a long and unpredictable latency due to remote loads, thereby providing tolerance to remote access latency. We began to develop standards to profile various algorithm and application parameters, such as the degree of parallelism, granularity, precision, instruction set mix, interprocessor communication, latency etc. These tools will continue to develop and evolve as the Information Power Grid environment matures. To provide a richer context for this research, the project also focused on issues of fault-tolerance and computation migration of numerical algorithms and software. During the initial phase we tried to increase our understanding of the bottlenecks in single processor performance. Our work began by developing an approach for the automatic generation and optimization of numerical software for processors with deep memory hierarchies and pipelined functional units. Based on the results we achieved in this study we are planning to study other architectures of interest, including development of cost models, and developing code generators appropriate to these architectures.

Dongarra, Jack↗

NASA NTRS2001

Modeling and Analysis of Wrinkled Membranes: An Overview

Thin-film membranes are basic elements of a variety of space inflatable/deployable structures. Wrinkling degrades the performance and reliability of these membrane structures, and hence has been a topic of continued interest. Wrinkling analysis of membranes for general geometry and arbitrary boundary conditions is quite challenging. The objective of this presentation is two-fold. Firstly, the existing models of wrinkled membranes and related numerical solution methods are reviewed. The important issues to be discussed are the capability of a membrane model to characterize taut, wrinkled and slack states of membranes in a consistent and physically reasonable manner; the ability of a wrinkling analysis method to predict the formation and growth of wrinkled regions, and to determine out-of-plane deformation and wrinkled waves; the convergence of a numerical solution method for wrinkling analysis; and the compatibility of a wrinkling analysis with general-purpose finite element codes. According to this review, several opening issues in modeling and analysis of wrinkled membranes that are to be addressed in future research are summarized, The second objective of this presentation is to discuss a newly developed membrane model of two viable parameters (2-VP model) and associated parametric finite element method (PFEM) for wrinkling analysis are introduced. The innovations and advantages of the proposed membrane model and PFEM-based wrinkling analysis are: (1) Via a unified stress-strain relation; the 2-VP model treat the taut, wrinkled, and slack states of membranes consistently; (2) The PFEM-based wrinkling analysis has guaranteed convergence; (3) The 2-VP model along with PFEM is capable of predicting membrane out-of-plane deformations; and (4) The PFEM can be integrated into any existing finite element code. Preliminary numerical examples are also included in this presentation to demonstrate the 2-VP model and PFEM-based wrinkling analysis approach.

Yang, B.↗

NASA NTRS2003

VLA Imaging of Protoplanetary Environments

Our program uses the techniques of millimeter interferometry to make high resolution observations of dust continuum emission to study the structure of protoplanetary disks around nearby stars. We continue to use the Very Large Array (VLA) of the National Radio Astronomy Observatories (NRAO) to make produce high angular resolution images, and we are improving our capabilities to interpret these observations with detailed numerical models. We have also made the first millimeter interferometer observations of protoplanetary disks in the southern hemisphere, using the recently upgraded Australia Telescope Compact Array (ATCA), taking advantage of the techniques we have developed in our work on northern sources.

Wilner, David J.↗

NASA NTRS2008

LiGa(OTf)(sub 4) as an Electrolyte Salt for Li-Ion Cells

Lithium tetrakis(trifluoromethane sulfo - nato)gallate [abbreviated "LiGa(OTf)4" (wherein "OTf" signifies trifluoro - methanesulfonate)] has been found to be promising as an electrolyte salt for incorporation into both liquid and polymer electrolytes in both rechargeable and non-rechargeable lithium-ion electrochemical cells. This and other ingredients have been investigated in continuing research oriented toward im proving the performances of rechargeable lithium-ion electrochemical cells, especially at low temperatures. This research at earlier stages, and the underlying physical and chemical principles, were reported in numerous previous NASA Tech Briefs articles. As described in more detail in those articles, lithiumion cells most commonly contain nonaqueous electrolyte solutions consisting of lithium hexafluorophosphate (LiPF6) dissolved in mixtures of cyclic and linear alkyl carbonates, including ethylene carbonate (EC), propylene carbonate (PC), dimethyl carbonate (DMC), diethyl carbonate (DEC), and ethyl methyl carbonate (EMC). Although such LiPF6-based electrolyte solutions are generally highly ionically conductive and electrochemically stable, as needed for good cell performance, there is interest in identifying alternate lithium electrolyte salts that, relative to LiPF6, are more resilient at high temperature and are less expensive. Experiments have been performed on LiGa(OTf)4 as well as on several other candidate lithium salts in pursuit of this interest. As part of these experiments, LiGa(OTf)4 was synthesized by the reaction of Ga(OTf)3 with an equimolar portion of LiOTf in a solvent consisting of anhydrous acetonitrile. Evaporation of the solvent yielded LiGa(OTf)4 as a colorless crystalline solid. The LiGa(OTf)4 and the other salts were incorporated into solutions with PC and DMC. The resulting electrolyte solutions exhibited reasonably high ionic conductivities over a relatively wide temperature range down to 40 C (see figure). In cyclic voltammetry measurements, LiGa(OTf)4 and the other salts exhibited acceptably high electrochemical stability over the relatively wide potential window of 0 to 5 V versus Li+/Li. 13C nuclear-magneticresonance measurements yielded results that suggested that in comparison with the other candidate salts, LiGa(OTf)4 exhibits less ion pairing. Planned further development will include optimization of the salt and solvent contents of such electrolyte solutions and incorporation of LiGa(OTf)4 into gel and solid-state polymer electrolytes. Of the salts, LiGa(OTf)4 is expected to be especially desirable for incorporation into lithium polymer electrolytes, wherein decreased ion pairing is advantageous and the large delocalized anions can exert a plasticizing effect.

Reddy, V. Prakash↗

NASA NTRS2009

The Capillary Flow Experiments Aboard the International Space Station: Increments 9-15

This report provides a summary of the experimental, analytical, and numerical results of the Capillary Flow Experiment (CFE) performed aboard the International Space Station (ISS). The experiments were conducted in space beginning with Increment 9 through Increment 16, beginning August 2004 and ending December 2007. Both primary and extra science experiments were conducted during 19 operations performed by 7 astronauts including: M. Fincke, W. McArthur, J. Williams, S. Williams, M. Lopez-Alegria, C. Anderson, and P. Whitson. CFE consists of 6 approximately 1 to 2 kg handheld experiment units designed to investigate a selection of capillary phenomena of fundamental and applied importance, such as large length scale contact line dynamics (CFE-Contact Line), critical wetting in discontinuous structures (CFE-Vane Gap), and capillary flows and passive phase separations in complex containers (CFE-Interior Corner Flow). Highly quantitative video from the simply performed flight experiments provide data helpful in benchmarking numerical methods, confirming theoretical models, and guiding new model development. In an extensive executive summary, a brief history of the experiment is reviewed before introducing the science investigated. A selection of experimental results and comparisons with both analytic and numerical predictions is given. The subsequent chapters provide additional details of the experimental and analytical methods developed and employed. These include current presentations of the state of the data reduction which we anticipate will continue throughout the year and culminate in several more publications. An extensive appendix is used to provide support material such as an experiment history, dissemination items to date (CFE publication, etc.), detailed design drawings, and crew procedures. Despite the simple nature of the experiments and procedures, many of the experimental results may be practically employed to enhance the design of spacecraft engineering systems involving capillary interface dynamics.

Jenson, Ryan M.↗

NASA NTRS2011

Collocation and Galerkin Time-Stepping Methods

We study the numerical solutions of ordinary differential equations by one-step methods where the solution at tn is known and that at t(sub n+1) is to be calculated. The approaches employed are collocation, continuous Galerkin (CG) and discontinuous Galerkin (DG). Relations among these three approaches are established. A quadrature formula using s evaluation points is employed for the Galerkin formulations. We show that with such a quadrature, the CG method is identical to the collocation method using quadrature points as collocation points. Furthermore, if the quadrature formula is the right Radau one (including t(sub n+1)), then the DG and CG methods also become identical, and they reduce to the Radau IIA collocation method. In addition, we present a generalization of DG that yields a method identical to CG and collocation with arbitrary collocation points. Thus, the collocation, CG, and generalized DG methods are equivalent, and the latter two methods can be formulated using the differential instead of integral equation. Finally, all schemes discussed can be cast as s-stage implicit Runge-Kutta methods.

Huynh, H. T.↗

NASA NTRS2011

Inversion of Magnetic Measurements of the CHAMP Satellite Over the Pannonian Basin

The Pannonian Basin is a deep intra-continental basin that formed as part of the Alpine orogeny. In order to study the nature of the crustal basement we used the long-wavelength magnetic anomalies acquired by the CHAMP satellite. The anomalies were distributed in a spherical shell, some 107,927 data recorded between January 1 and December 31 of 2008. They covered the Pannonian Basin and its vicinity. These anomaly data were interpolated into a spherical grid of 0.5 x 0.5, at the elevation of 324 km by the Gaussian weight function. The vertical gradient of these total magnetic anomalies was also computed and mapped to the surface of a sphere at 324 km elevation. The former spherical anomaly data at 425 km altitude were downward continued to 324 km. To interpret these data at the elevation of 324 km we used an inversion method. A polygonal prism forward model was used for the inversion. The minimum problem was solved numerically by the Simplex and Simulated annealing methods; a L2 norm in the case of Gaussian distribution parameters and a L1 norm was used in the case of Laplace distribution parameters. We INTERPRET THAT the magnetic anomaly WAS produced by several sources and the effect of the sable magnetization of the exsolution of hemo-ilmenite minerals in the upper crustal metamorphic rocks.

Kis, K. I.↗

NASA NTRS1983

A comparison of numerical results of Arctic Sea ice modeling with satellite images

Observations made from 1972 to 1976 with the Electrically Scanning Microwave Radiometer on board the Nimbus 5 Satellite provide sequential synoptic information of the Arctic sea ice cover. This 4 year data set was used to construct a fairly continuous series of 3 day average 19-GHz passive microwave images which has become a valuable source of polar information, yielding many anticipated and unanticipated discoveries of the sea ice canopy observed in its entirety through the clouds and during the polar night. Interpretation of the passive microwave satellite data set was performed by comparing selected sequential images with correspondin microwave profiles and images acquired by the NASA CV-990 airborne laboratory and with various in situ microwave and physical observations.

Ling, C. H.↗

NASA NTRSDate not supplied

Turbulence, Transition, and Numerical Method Technologies

The close collaboration in the validation experimentation effort was excellent, the complementary research efforts related tofundamental numerical model development, and the activities related to effective HPC utilization on near term architectures that are coming down the pipeline were exactly what NASA should be doing. Even though the TQR panel cannot say that the success criteria were fully met, we do feel that the research funded by TTT/RCA was successful and useful in determining the next suite of CFD validation tests and the problematic flow conditions that continue to be difficult to predict using our current suite of computational tools. Given the TQR Panel’s review criteria: 1. The deliverables meet the technical requirements; 2.The appropriate technical approaches were followed in producing the deliverables; and 3. All technical caveats and concerns have been identified and the deliverables: a.Development of more accurate physics-based methods (e.g., higher moment closure); b. Large eddy simulation (LES); c. Advanced numerical methods; d. Transition prediction and modeling; e.Validation experiments; f. Multidisciplinary analysis and design (high fidelity). TQR panel feels that the technical processes were all completed and that the review criteria and all the deliverables were met. Thus, the Milestone TACP-2016-001 has been completed.

Mujeeb R Malik↗

NASA NTRS1981

Vertical-structure effects on planetary microwave brightness temperature measurements - Applications to the lunar regolith

An attempt to constrain the effects of vertical variations in dielectric properties on lunar microwave observations is presented. A numerical approach for deriving the reflectivity and microwave weighting function of a vertically varying half-space is used, assuming variance in the dielectric properties with depth only, and negligible magnetic effects. The cases of continuous and stratified models of vertical structures are discussed, and a concentration of emitted energy in upper layers is found. The total emitted energy oscillates, varying with the thickness of the upper soil layer, but averaging out interference effects due to random variations in the substrate depth. Consideration is also given to the vertical structure effects on the lunation-mean disk-center brightness temperature, its variations, and the regolith electrical loss, and predicted reflectivity effects by feasible models of the lunar regolith dielectric profile.

Keihm, S. J.↗

NASA NTRS1983

A spectral numerical method for the Navier-Stokes equations with applications to Taylor-Couette flow

A new spectral method for solving the incompressible Navier-Stokes equations in a plane channel and between concentric cylinders is presented. The method uses spectral expansions which inherently satisfy the boundary conditions and the continuity equation and yield banded matrices which are efficiently solved at each time step. In addition, the number of dependent variables is reduced, resulting in a reduction in computer memory requirements. Several test problems have been computed for the channel flow and for flow between concentric cylinders, including Taylor-Couette flow with axisymmetric Taylor vortices and wavy vortices. In all cases, agreement with available experimental and theoretical results is very good.

Moser, R. D.↗

NASA NTRS1986

Modelling short channel mosfets for use in VLSI

In an investigation of metal oxide semiconductor field effect transistor (MOFSET) devices, a one-dimensional mathematical model of device dynamics was prepared, from which an accurate and computationally efficient drain current expression could be derived for subsequent parameter extraction. While a critical review revealed weaknesses in existing 1-D models (Pao-Sah, Pierret-Shields, Brews, and Van de Wiele), this new model in contrast was found to allow all the charge distributions to be continuous, to retain the inversion layer structure, and to include the contribution of current from the pinched-off part of the device. The model allows the source and drain to operate in different regimes. Numerical algorithms used for the evaluation of surface potentials in the various models are presented.

Klafter, Alex↗