Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical approximation & analysis”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
NASA NTRS1992

Optimal matrix approximants in structural identification

Problems of model correlation and system identification are central in the design, analysis, and control of large space structures. Of the numerous methods that have been proposed, many are based on finding minimal adjustments to a model matrix sufficient to introduce some desirable quality into that matrix. In this work, several of these methods are reviewed, placed in a modern framework, and linked to other previously known ideas in computational linear algebra and optimization. This new framework provides a point of departure for a number of new methods which are introduced here. Significant among these is a method for stiffness matrix adjustment which preserves the sparsity pattern of an original matrix, requires comparatively modest computational resources, and allows robust handling of noisy modal data. Numerical examples are included to illustrate the methods presented herein.

Beattie, C. A.↗
DOE OSTI2026

Computational capacity in hydrodynamic real-time hybrid simulation applied to simulate the dynamic response of floating offshore wind turbines

Real-time hybrid simulation (RTHS) mitigates similitude distortions in model-scale tests of floating offshore wind turbines (FOWTs) by coupling physical experiments with numerical models in real time. The coupling requires faster-than-real-time numerical computations to satisfy temporal similitude with the physical experiment, presenting a bottleneck for using more complex numerical models in RTHS. This paper presents a hydrodynamic-RTHS (hydro-RTHS) framework for FOWTs that simulates the hydrodynamics physically and the aerodynamics numerically with sensor feedback from the physical testing. The framework adapts the three-loop hardware architecture to leverage greater computational resources and mitigate strict temporal requirements, enabling more computationally demanding numerical analyses in hydro-RTHS. The three-loop hardware architecture integrates multiple machines, each dedicated to either numerical analysis or RTHS controls, with a rate-transition algorithm to synchronize the tasks executed across the different machine processors. Virtual and physical tests verified and validated the hydro-RTHS framework, respectively. The ”virtual” tests, which approximates the physical domain numerically, verified the RTHS framework with respect to a numerical full-scale complete FOWT model simulated in the open-source software, OpenFAST. The virtual tests were able to maintain comparable control signals while enabling greater computational resources for the numerical calculations. Real-world physical tests demonstrated that the hydro-RTHS framework computes aerodynamic forces similar to the complete OpenFAST model, validating the hydro-RTHS framework using the three-loop hardware architecture. Findings show that the hydro-RTHS framework with the three-loop hardware architecture is computationally efficient, with reserve capacity to simulate more complex problems due to the customized software, hardware, and rate-transition algorithm.

17 WIND ENERGY↗
NASA NTRS1994

Analysis of the early spectra and light curve of SN 1987A

Numerical modeling of supernova spectra, light curves, and hydrodynamics requires physical inputs, numerical techniques, approximations, and assumptions which must be thoroughly understood in order to study the details of supernova explosions. Here, we discuss some of these in the context of the early evolution of supernova 1987A. Gray radiation-hydrodynamics is used to calculate the bolometric light curve and the hydrodynamic evolution of the supernova. Synthetic spectra are then obtained for the resulting density and velocity structure. The spectrum calculations are performed using a special-relativistic treatment of the radiative transfer equation in the comoving frame, line blanketing by about 10(exp 5) spectral lines, and departures from local thermodynamic equilibrium (LTE) for H I, He I, Mg II, and Ca II. We find that we are able to simultaneously fit the early light curve and spectra reasonably well, using a progenitor model from Arnett (1991a), without fine-tuning the free parameters. Temperature structures and radiative equilibrium, non-LTE effects, homologous expansion, and mean opacities are discussed.

Hauschildt, Peter H.↗
DOE OSTI2020

Application of a Null-Space Monte Carlo Flow Model Set to the Composite Analysis Base Case Fate and Transport Modeling

The Plateau to River Groundwater Model (P2R Model) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by CH2M HILL Plateau Remediation Company at the Hanford Site in Washington State. Figure 1-1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model is utilized in the Composite Analysis (CA) for the Hanford Site as the computational engine for computing F&T predictions as described in CP-60406, Hanford Site Composite Analysis Technical Approach Description: Groundwater. The model simulates contaminants of concern within the saturated zone of the uppermost aquifer beneath the Central Plateau and downgradient to the Columbia River. CP-57037, Model Package Report for the Plateau to River Model Version 8.3, documents the current version of the P2R Model including a description of the conceptual site model, model development and calibration, and limitations to the model application. The overall objective of the saturated zone modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within and downgradient of the Central Plateau of the Hanford Site. Specifically, the purpose of this environmental calculation is to describe the application of the hydraulic property fields and recharge parameters documented in ECF-HANFORD-20-0027, Null Space Monte Carlo Evaluation of the Plateau to River Model, to the CA flow and fate and transport simulation results to quantify the uncertainty in the simulated results due to input parameter selection. Use of numerical groundwater models is always accompanied with uncertainty in the results produced by a model because models are approximations of reality. Thus, by definition, models lack the detail to fully represent observed behavior. Use of numerical techniques, such as a NSMC analysis, can help in identifying and quantifying the potential uncertainties associated with a numerical model such as the P2R Model. The result of NSCM analysis is a set of F&T simulations that provide an estimate of the range of possible outcomes that are used to quantify the uncertainty in simulated concentrations produced using the base case simulations. The simulated concentrations from all simulations will support calculation of the uncertainty of the total dose calculated in a separate calculation.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1993

Program Computes Flows Of Fluids And Heat

SINDA'85/FLUINT incorporates lumped-parameter-network and one-dimensional-flow mathematical models. System enables analysis of mutual influences of thermal and flow phenomena. Offers two finite-difference numerical solution techniques: forward-difference explicit approximation and Crank-Nicholson approximation. Enables simulation of nonuniform heating and facilitates mathematical modeling of thin-walled heat exchangers. Ability to model nonequilibrium behavior within two-phase volumes included. Recent changes in program improve modeling of real evaporator pumps and other capillary-assist evaporators. Written in FORTRAN 77.

Cullimore, Brent↗
NASA NTRS1993

Influence of test configuration on the combustion characteristics of polymers as ignition sources

The experimental evaluation of polymers as ignition sources for metals was accomplished at the NASA White Sands Test Facility (WSTF) using a standard promoted combustion test. These tests involve the transient burning of materials in high-pressure oxygen environments. They have provided data from which design decisions can be made; data include video recordings of ignition and non-ignition for specific combinations of metals and polymers. Other tests provide the measured compositions of combustion products for polymers at select burn times and an empirical basis for estimating burn rates. With the current test configuration, the detailed analysis of test results requires modeling a three-dimensional, transient convection process involving fluid motion, thermal conduction and convection, the diffusion of chemical species, and the erosion of sample surface. At the high pressure extremes, it even requires the analysis of turbulent, transient convection where the physics of the problem are not well known and the computation requirements are not practical at this time. An alternative test configuration that can be analyzed with a relatively-simple convection model was developed during the summer period. The principal change constitutes replacing a large-diameter polymer disk at the end of the metal test rod with coaxial polymer cylinders that have a diameter nearer to that of the metal rod. The experimental objective is to assess the importance of test geometries on the promotion of metal ignition by testing with different lengths of the polymer and, with an extended effort, to analyze the surface combustion in the redesigned promoted combustion tests through analytical modeling of the process. The analysis shall use the results of cone-calorimeter tests of the polymer material to model primary chemical reactions and, with proper design of the promoted combustion test, modeling of the convection process could be conveniently limited to a quasi-steady boundary layer analysis where the economical solution of parabolic equations is involved. The products for the summer period are: (1) a conceptual-level redesign of the test apparatus, and (2) the development and use of an approximate integral boundary layer analysis to demonstrate the influence of geometry changes prior to testing. A computer code STAN5, an accurate numerical boundary layer model whose earlier versions were developed for the NASA Lewis Research Center by the Fellow, also was installed and validated on the WSTF and New Mexico State University computer systems as a starting point in the development of a more detailed fluid mechanics and combustion model.

Julien, Howard L.↗
NASA NTRS1998

A Hybrid SBR/MoM Technique for Analysis of Scattering from Small Protrusions on a Large Conducting Body

For analysis of large-scale electromagnetic scattering problems, high-frequency asymptotic methods are fast but approximate, whereas low-frequency numerical methods are accurate but slow. Neither can produce an efficient and accurate solution to scattering by large bodies containing small structures. A promising approach is to combine the best features of both types of methods to produce a hybrid technique that is sufficiently fast, reasonably accurate, and applicable to a class of unsolvable problems such as the scatterers mentioned above. There are two extremes for this type of hybridization. One is simply to superimpose solutions from asymptotic and numerical methods. While this approach is most widely used in practical applications, it neglects the interactions between the two solutions, which can be significant in many problems. The other extreme is to combine an asymptotic and a numerical method in an exact manner. In this approach, the effect of a large body is included by incorporating its diffraction into the Green's function in the integral equation for the small structures, which accounts for all interactions. While this approach is accurate, it is difficult to be implemented in a general-purpose computer code because of its complex nature. A more practical approach is to develop a technique that can include all significant interactions and neglect all trivial interactions. The resulting hybrid technique can produce sufficient accuracy and can be implemented in a general-purpose computer code. In this paper, we develop a technique that combines the shooting and- bouncing-ray (SBR) method and the method of moments (MoM) to solve for the scattering by large conducting bodies with small structures mounted on their surfaces.

Jin, J. M.↗
NASA NTRS2023

A Finite Element Approach for Simplified 2D Nonlinear Dynamic Contact/Impact Analysis

In this paper, a simplified numerical approach for finite element dynamic analysis of an inelastic solid structure subjected to solid object impact is presented. The approach approximates the impacting solid as the selected multiple nodes, for which mass of the impactor is distributed. The node-to-segment contact formulation with the penalty constraint technique incorporated is employed to impose contact conditions between the nodes and the surface of the receiver structure. The node-to-segment algorithm is integrated into Newton–Raphson time integration scheme and the Lagrange multiplier technique is applied to enforce the identical displacements for the selected nodes throughout the analysis process. The approach is verified using two-dimensional plane strain models considering elastic-perfectly-plastic material behavior. The results obtained using the proposed approach are in a good agreement with those simulated using a commercial finite element code, ABAQUS dynamic/implicit, in terms of displacements and stress distribution fields. The proposed approach is shown to be computationally superior to general finite element method-based contact/impact analysis without significantly sacrificing the accuracy.

Seungwook Seok↗
DOE OSTI2023

Exponential acceleration of macroscopic quantum tunneling in a Floquet Ising model

The exponential suppression of macroscopic quantum tunneling (MQT) in the number of elements to be reconfigured is an essential element of broken symmetry phases. This suppression is also a core bottleneck in quantum algorithms, such as traversing an energy landscape in optimization, and adiabatic state preparation more generally. In this work, we demonstrate exponential acceleration of MQT through Floquet engineering with the application of a uniform, high frequency transverse drive field. Using the ferromagnetic phase of the transverse field Ising model in one and two dimensions as a prototypical example, we identify three phenomenological regimes as a function of drive strength. For weak drives, the system exhibits exponentially decaying tunneling rates but robust magnetic order; in the crossover regime at intermediate drive strength, we find polynomial decay of tunnelling alongside vanishing magnetic order; and at very strong drive strengths both the Rabi frequency and time-averaged magnetic order are approximately constant with increasing system size. We support these claims with extensive full wavefunction and tensor network numerical simulations, and theoretical analysis. An experimental test of these results presents a technologically important and novel scientific question accessible on NISQ-era quantum computers.

Grattan, George↗
NASA NTRS1991

The asymptotic spectra of banded Toeplitz and quasi-Toeplitz matrices

Toeplitz matrices occur in many mathematical, as well as, scientific and engineering investigations. This paper considers the spectra of banded Toeplitz and quasi-Toeplitz matrices with emphasis on non-normal matrices of arbitrarily large order and relatively small bandwidth. These are the type of matrices that appear in the investigation of stability and convergence of difference approximations to partial differential equations. Quasi-Toeplitz matrices are the result of non-Dirichlet boundary conditions for the difference approximations. The eigenvalue problem for a banded Toeplitz or quasi-Toeplitz matrix of large order is, in general, analytically intractable and (for non-normal matrices) numerically unreliable. An asymptotic (matrix order approaches infinity) approach partitions the eigenvalue analysis of a quasi-Toeplitz matrix into two parts, namely the analysis for the boundary condition independent spectrum and the analysis for the boundary condition dependent spectrum. The boundary condition independent spectrum is the same as the pure Toeplitz matrix spectrum. Algorithms for computing both parts of the spectrum are presented. Examples are used to demonstrate the utility of the algorithms, to present some interesting spectra, and to point out some of the numerical difficulties encountered when conventional matrix eigenvalue routines are employed for non-normal matrices of large order. The analysis for the Toeplitz spectrum also leads to a diagonal similarity transformation that improves conventional numerical eigenvalue computations. Finally, the algorithm for the asymptotic spectrum is extended to the Toeplitz generalized eigenvalue problem which occurs, for example, in the stability of Pade type difference approximations to differential equations.

Beam, Richard M.↗
NASA NTRS1991

Mission and system optimization of nuclear electric propulsion vehicles for lunar and Mars missions

The detailed mission and system optimization of low thrust electric propulsion missions is a complex, iterative process involving interaction between orbital mechanics and system performance. Through the use of appropriate approximations, initial system optimization and analysis can be performed for a range of missions. The intent of these calculations is to provide system and mission designers with simple methods to assess system design without requiring access or detailed knowledge of numerical calculus of variations optimizations codes and methods. Approximations for the mission/system optimization of Earth orbital transfer and Mars mission have been derived. Analyses include the variation of thruster efficiency with specific impulse. Optimum specific impulse, payload fraction, and power/payload ratios are calculated. The accuracy of these methods is tested and found to be reasonable for initial scoping studies. Results of optimization for Space Exploration Initiative lunar cargo and Mars missions are presented for a range of power system and thruster options.

Gilland, James H.↗
DOE OSTI2024

A Provably Accurate Randomized Sampling Algorithm for Logistic Regression

In statistics and machine learning, logistic regression is a widely-used supervised learning technique primarily employed for binary classification tasks. When the number of observations greatly exceeds the number of predictor variables, we present a simple, randomized sampling-based algorithm for logistic regression problem that guarantees high-quality approximations to both the estimated probabilities and the overall discrepancy of the model. Our analysis builds upon two simple structural conditions that boil down to randomized matrix multiplication, a fundamental and well-understood primitive of randomized numerical linear algebra. We analyze the properties of estimated probabilities of logistic regression when leverage scores are used to sample observations, and prove that accurate approximations can be achieved with a sample whose size is much smaller than the total number of observations. To further validate our theoretical findings, we conduct comprehensive empirical evaluations. Overall, our work sheds light on the potential of using randomized sampling approaches to efficiently approximate the estimated probabilities in logistic regression, offering a practical and computationally efficient solution for large-scale datasets.

Chowdhury, Agniva↗
NASA NTRS1994

Approximate analysis for repeated eigenvalue problems with applications to controls-structure integrated design

A method for eigenvalue and eigenvector approximate analysis for the case of repeated eigenvalues with distinct first derivatives is presented. The approximate analysis method developed involves a reparameterization of the multivariable structural eigenvalue problem in terms of a single positive-valued parameter. The resulting equations yield first-order approximations to changes in the eigenvalues and the eigenvectors associated with the repeated eigenvalue problem. This work also presents a numerical technique that facilitates the definition of an eigenvector derivative for the case of repeated eigenvalues with repeated eigenvalue derivatives (of all orders). Examples are given which demonstrate the application of such equations for sensitivity and approximate analysis. Emphasis is placed on the application of sensitivity analysis to large-scale structural and controls-structures optimization problems.

Kenny, Sean P.↗
DOE OSTI2022

Droplet formation simulation using mixed finite elements

Droplet formation happens in finite time due to the surface tension force. The linear stability analysis is useful to estimate the size of a droplet but fails to approximate the shape of the droplet. This is due to a highly nonlinear flow description near the point where the first pinch-off happens. A one-dimensional axisymmetric mathematical model was first developed by Eggers and Dupont [“Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] using asymptotic analysis. This asymptotic approach to the Navier–Stokes equations leads to a universal scaling explaining the self-similar nature of the solution. Numerical models for the one-dimensional model were developed using the finite difference [Eggers and Dupont, “Drop formation in a one-dimensional approximation of the Navier–Stokes equation,” J. Fluid Mech. 262, 205–221 (1994)] and finite element method [Ambravaneswaran et al., “Drop formation from a capillary tube: Comparison of one-dimensional and two-dimensional analyses and occurrence of satellite drops,” Phys. Fluids 14, 2606–2621 (2002)]. The focus of this study is to provide a robust computational model for one-dimensional axisymmetric droplet formation using the Portable, Extensible Toolkit for Scientific Computation. The code is verified using the Method of Manufactured Solutions and validated using previous experimental studies done by Zhang and Basaran [“An experimental study of dynamics of drop formation,” Phys. Fluids 7, 1184–1203 (1995)]. The present model is used for simulating pendant drops of water, glycerol, and paraffin wax, with an aspiration of extending the application to simulate more complex pinch-off phenomena.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1982

Higher-Order Panel Method for Aerodynamic Flow Analysis

PANAIR uses high-order panel method to predict inviscid subsonic or supersonic flows about arbitrary configuration. Panel method solves linear partial differential equation numerically by approximating configuration surface with panels on which unknown "singularity strengths" are defined. PANAIR includes advanced software technology as well as advanced aerodynamic technology.

Erickson, L.↗
DOE OSTI2020

Null Space Monte Carlo Evaluation of the Plateau to River Model

The Plateau to River Groundwater Model (P2R Model) is a groundwater flow and contaminant fate and transport (F&T) simulation model used to support remedial activities conducted by CH2M HILL Plateau Remediation Company at the Hanford Site in Washington State. Figure 1-1 illustrates the P2R Model extents, discretization, and boundary conditions. The P2R Model provides a computational framework to simulate the F&T of contaminants in groundwater associated with the 200-PO-1, 200-UP-1, 200-BP-5, and 200-ZP-1 Groundwater Operable Units (OUs) in the Hanford Site Central Plateau. In addition, the model includes adjacent areas and facilities (e.g., the State Approved Land Disposal Site). Intended and anticipated uses of the model include calculating water levels, hydraulic gradients, and groundwater flows throughout the model domain (encompassing the 200 West and 200 East Areas) for use in subsequent F&T calculations for contaminants of concern and developing scale-appropriate, telescopic-mesh refinement models for detailed evaluation of areas within the model domain where required. The overall objective of the modeling effort is to provide a basis for making informed remedial action decisions based on descriptions of current and expected future contaminant concentrations in groundwater at decision points within the OU boundaries. The objective for the model development phase is to create a common modeling platform that can be used for investigations of the Central Plateau groundwater OUs and areas downgradient toward the Columbia River. The P2R Model calibration to historical data observed at the Hanford Site is documented in CP-57037, Model Package Report for the Plateau to River Model Version 8.3. The purpose of this environmental calculation is to describe a null space Monte Carlo (NSMC) evaluation was conducted with the historic calibration of the P2R model. Use of numerical groundwater models is always accompanied with uncertainty in the results produced by a model because models are approximations of reality. Thus, by definition, lack the detail to fully represent observed behavior. Use of numerical techniques, such as a NSMC analysis, can help in identifying and quantifying the potential uncertainties associated with a numerical model such as the P2R Model. Use of the NSMC approach results in 100 groundwater flow models that are variants of the calibrated P2R Model. These variant models can be used to evaluate uncertainty in model predictions made by the calibrated P2R Model for other analyses. A secondary purpose of the environmental calculation is to establish these variant models for use with other applications.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
NASA NTRS1981

Numerical simulation of flows in curved diffusers with cross-sectional transitioning using a three-dimensional viscous analysis

A three dimensional analysis for fully viscous, subsonic, compressible flow is evaluated. An approximate form of the Navier Stokes equations is solved by an implicit spatial marching technique. Calculations were made for flow in a circular S duct and in the F 16 inlet duct. The computed total pressure contours and secondary flow velocity vectors are presented. Qualitative comparisons with experiment are shown for both ducts. The analysis is used to show how the cross section transitioning in the F 16 inlet suppresses the development of a secondary flow vortex.

Towne, C. E.↗
NASA NTRS1985

Numerical Estimation of the Curvature of Biological Surfaces

Many biological systems may profitably be studied as surface phenomena. A model consisting of isotropic growth of a curved surface from a flat sheet is assumed. With such a model, the Gaussian curvature of the final surface determines whether growth rate of the surface is subharmonic or superharmonic. These properties correspond to notions of convexity and concavity, and thus to local excess growth and local deficiency of growth. In biological models where the major factors controlling surface growth are intrinsic to the surface, researchers thus gained from geometrical study information on the differential growth undergone by the surface. These ideas were applied to an analysis of the folding of the cerebral cortex, a geometrically rather complex surface growth. A numerical surface curvature technique based on an approximation to the Dupin indicatrix of the surface was developed. A metric for comparing curvature estimates is introduced, and considerable numerical testing indicated the reliability of this technique.

Todd, P. H.↗