SEARCH · Engineering Papers

Results for “solve time”

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 469 records · Page 26

All sources NASA NTRS DOE OSTI

NASA NTRS1979

Laminar and turbulent flows over spherically blunted cone and hyperboloid with massive surface blowing

Numerical solutions are presented for the flow over a spherically blunted cone and hyperboloid with massive surface blowing. Time-dependent viscous shock-layer equations are used to describe the flow field. The boundary conditions on the body surface include a prescribed blowing-rate distribution. The governing equations are solved by a time-asymptotic finite-difference method. Results presented here are only for a perfect gas-type flow at zero angle of attack. Both laminar and turbulent flow solutions are obtained. It is found that the effect of the surface blowing on the laminar flow field is to smooth out the curvature discontinuity at the sphere-cone juncture point, which results in a positive pressure gradient over the body. The shock slope increases on the downstream portion of the body as the surface blowing rate is increased. The turbulent flow with surface blowing is found to redevelop a boundary-layer-like region near the surface. The effects of this boundary-layer region on the flow field and heating rates are discussed.

Kumar, A.↗

NASA NTRS1983

A fully vectorized numerical solution of the incompressible Navier-Stokes equations

A vectorizable algorithm is presented for the implicit finite difference solution of the incompressible Navier-Stokes equations in general curvilinear coordinates. The unsteady Reynolds averaged Navier-Stokes equations solved are in two dimension and non-conservative primitive variable form. A two-layer algebraic eddy viscosity turbulence model is used to incorporate the effects of turbulence. Two momentum equations and a Poisson pressure equation, which is obtained by taking the divergence of the momentum equations and satisfying the continuity equation, are solved simultaneously at each time step. An elliptic grid generation approach is used to generate a boundary conforming coordinate system about an airfoil. The governing equations are expressed in terms of the curvilinear coordinates and are solved on a uniform rectangular computational domain. A checkerboard SOR, which can effectively utilize the computer architectural concept of vector processing, is used for iterative solution of the governing equations.

Patel, N.↗

NASA NTRS1994

Advances in parameter estimation techniques applied to flexible structures

In this work, various parameter estimation techniques are investigated in the context of structural system identification utilizing distributed parameter models and 'measured' time-domain data. Distributed parameter models are formulated using the PDEMOD software developed by Taylor. Enhancements made to PDEMOD for this work include the following: (1) a Wittrick-Williams based root solving algorithm; (2) a time simulation capability; and (3) various parameter estimation algorithms. The parameter estimations schemes will be contrasted using the NASA Mini-Mast as the focus structure.

Maben, Egbert↗

DOE OSTI2024

Local enhancement in transient absorption spectroscopy by gating the resonance in the time domain

The concept of resonant perfect absorption, enabled by the combined action of pulse propagation and an auxiliary gate pulse, was recently proposed and demonstrated in a group of two-level systems [Y. He , ]. Here we exploit this method in a more realistic scenario by solving the coupled time-dependent Schrödinger equation and the Maxwell wave equation in helium. Through emptying the population of the 1 s 2 p excited state after its excitation, we explore the evolution of the spectral profile with time delay and propagation distance and link the observations to the controlled interference between the original field and the gated new field. We find that resonant absorption for higher-lying states can also be strongly enhanced, in spite of the congestion of multiple resonances and the presence of complex laser-induced couplings. Our results show that interferometric control of absorption using intense laser fields can be applied selectively in both the temporal and spatial domains. Published by the American Physical Society 2024

He, Yu (ORCID:0000000152120176)↗

NASA NTRS1994

Direct numerical simulation of sheared turbulent flow

The summer assignment to study sheared turbulent flow was divided into three phases which were: (1) literature survey, (2) computational familiarization, and (3) pilot computational studies. The governing equations of fluid dynamics or Navier-Stokes equations describe the velocity, pressure, and density as functions of position and time. In principle, when combined with conservation equations for mass, energy, and thermodynamic state of the fluid a determinate system could be obtained. In practice the Navier-Stokes equations have not been solved due to the nonlinear nature and complexity of these equations. Consequently, the importance of experiments in gaining insight for understanding the physics of the problem has been an ongoing process. Reasonable computer simulations of the problem have occured as the computational speed and storage of computers has evolved. The importance of the microstructure of the turbulence dictates the need for high resolution grids in extracting solutions which contain the physical mechanisms which are essential to a successful simulation. The recognized breakthrough occurred as a result of the pioneering work of Orzag and Patterson in which the Navier-Stokes equations were solved numerically utilizing a time saving toggling technique between physical and wave space, known as a spectral method. An equally analytically unsolvable problem, containing the same quasi-chaotic nature as turbulence, is known as the three body problem which was studied computationally as a first step this summer. This study was followed by computations of a two dimensional (2D) free shear layer.

Harris, Vascar G.↗

DOE OSTI2024

The Effective Fragment Molecular Orbital Method: Achieving High Scalability and Accuracy for Large Systems

The effective fragment molecular orbital (EFMO) method has been developed to predict the total energy of a very large molecular system accurately (with respect to the underlying quantum mechanical method) and efficiently by taking advantage of the locality of strong chemical interactions and employing a two-level hierarchical parallelism. The accuracy of the EFMO method is partly attributed to the accurate and robust intermolecular interaction prediction between distant fragments, in particular, the many-body polarization and dispersion effects, which require the generation of static and dynamic polarizability tensors by solving the coupled perturbed Hartree–Fock (CPHF) and time-dependent HF (TDHF) equations, respectively. Solving the CPHF and TDHF equations is the main EFMO computational bottleneck due to the inefficient (serial) and I/O-intensive implementation of the CPHF and TDHF solvers. In this work, the efficiency and scalability of the EFMO method are significantly improved with a new CPU memory-based implementation for solving the CPHF and TDHF equations that are parallelized by either message passing interface (MPI) or hybrid MPI/OpenMP. Here, the accuracy of the EFMO method is demonstrated for both covalently bonded systems and noncovalently bound molecular clusters by systematically examining the effects of basis sets and a key distance-related cutoff parameter, R cut . R cut determines whether a fragment pair (dimer) is treated by the chosen ab initio method or calculated using the effective fragment potential (EFP) method (separated dimers). Decreasing the value of Rcut increases the number of separated (EFP) dimers, thereby decreasing the computational effort. It is demonstrated that excellent accuracy (<1 kcal/mol error per fragment) can be achieved when using a sufficiently large basis set with diffuse functions coupled with a small R cut value. With the new parallel implementation, the total EFMO wall time is substantially reduced, especially with a high number of MPI ranks. Given a sufficient workload, nearly ideal strong scaling is achieved for the CPHF and TDHF parts of the calculation. For the first time, EFMO calculations with the inclusion of long-range polarization and dispersion interactions on a hydrated mesoporous silica nanoparticle with explicit water solvent molecules (more than 15k atoms) are achieved on a massively parallel supercomputer using nearly 1000 physical nodes. In addition, EFMO calculations on the carbinolamine formation step of an amine-catalyzed aldol reaction at the nanoscale with explicit solvent effects are presented.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS2003

"Virtual Cockpit Window" for a Windowless Aerospacecraft

A software system processes navigational and sensory information in real time to generate a three-dimensional-appearing image of the external environment for viewing by crewmembers of a windowless aerospacecraft. The design of the particular aerospacecraft (the X-38) is such that the addition of a real transparent cockpit window to the airframe would have resulted in unacceptably large increases in weight and cost. When exerting manual control, an aircrew needs to see terrain, obstructions, and other features around the aircraft in order to land safely. The X-38 is capable of automated landing, but even when this capability is utilized, the crew still needs to view the external environment: From the very beginning of the United States space program, crews have expressed profound dislike for windowless vehicles. The wellbeing of an aircrew is considerably promoted by a three-dimensional view of terrain and obstructions. The present software system was developed to satisfy the need for such a view. In conjunction with a computer and display equipment that weigh less than would a real transparent window, this software system thus provides a virtual cockpit window. The key problem in the development of this software system was to create a realistic three-dimensional perspective view that is updated in real time. The problem was solved by building upon a pre-existing commercial program LandForm C3 that combines the speed of flight-simulator software with the power of geographic-information-system software to generate real-time, three-dimensional-appearing displays of terrain and other features of flight environments. In the development of the present software, the pre-existing program was modified to enable it to utilize real-time information on the position and attitude of the aerospacecraft to generate a view of the external world as it would appear to a person looking out through a window in the aerospacecraft. The development included innovations in realistic horizon-limit modeling, three-dimensional stereographic display, and interfaces for utilization of data from inertial-navigation devices, Global Positioning System receivers, and laser rangefinders.

Abernathy, Michael F.↗

NASA NTRS1974

Water-management models in Florida from ERTS-1 data

A prototype multiparameter data acquisition network, installed and operated by the U.S. Geological Survey is a viable approach for obtaining near real-time data needed to solve hydrologic problems confronting nearly 2.5 million residents of south Florida. Selected water quantity and quality data obtained from ground stations are transmitted for relay via ERTS-1 to NASA receiving stations in virtual real time. This data relay system has been very reliable and, by coupling the ground information with ERTS imagery, a modeling technique is available for water resource management in south Florida. An ecological model has been designed for the Shark River Slough in Everglades National Park.

Higer, A. L.↗

NASA NTRS1971

Time dependent shear stress and temperature distribution over an insulated flat plate moving at hypersonic speed.

The laminar two-dimensional flow over a stepwise accelerated flat plate moving with hypersonic speed at zero angle of attack is analysed. The governing equations in the self-similar form are linearized and solved numerically for small times. The solutions obtained are the deviations of the velocity and the temperature profiles from those of steady state. The presented results may be used to find the first order boundary layer induced pressure on the plate.

Rodkiewicz, C. M.↗

DOE OSTI2023

Frequency Domain Reflectometry (FDR) Simulation Techniques for Digital Twin Representation of an Electrical Cable

Simulation of cable system response to frequency domain reflectometry (FDR) tests can be instrumental to understanding these tests and the nature and influence of various cable anomalies on test signatures. Reflectometry simulations are based upon a finite element representation of cable conductors and insulation to produce an S-parameter at each evaluated frequency. The aggregate collection of cable model S-parameters can simulate the influence on a test signal injected into a physical cable. Such an approach was undertaken in this work to produce a digital twin simulation of a low-voltage electrical cable. The electrical cable digital twin examined the influence of test simulation parameters and the relative influence of cable anomalies, including thermal aging, water or moisture exposure, water or moisture ingress, and other anomalies. The digital twin in this work included modeling of the conductors, insulation, jacket, and surrounding environment (air, water, etc.). The digital twin could be expanded to include cable bends, junctions and splices, branch or T systems, and termination impedances of motors or instruments. Observations and conclusions of this work include: 1. Fully 3D digital twin simulation of an electrical cable using an FDR approach is possible. However, there are tradeoffs between simulation fidelity and solution time, which must be balanced to ensure the simulation solves in an adequate amount of time (e.g., less than 20 minutes). Simulation parameters to balance include frequency bandwidth, number of frequencies, mesh density, connection impedance, and permittivity tolerance. 2. The digital twin simulation can explain FDR sensitivity to various cable anomalies, including entry and exit from an oven or water bath. 3. The digital twin simulation FDR response attenuates with distance along the cable and is further affected by the frequency bandwidth, which is similar to that observed with physical measurements. 4. The resolution of the digital twin FDR peaks increased with increasing bandwidth and with increasing number of frequencies, again similar to physical measurements. 5. The presence of multiple anomalies in the digital twin does not substantially attenuate the FDR response to anomalies located beyond the first encountered anomaly and impedance mismatch. 6. Spectral variation of the permittivity did not have a significant effect on the FDR response compared to a fixed nominal value. 7. Extension of the digital twin to 1000 ft still allowed for detection of distal anomalies near the far end of the electrical cable from the instrument connection point. 8. The ARENA test bed facilitates efficient NDE evaluations of well understood cable anomalies with various NDE methods without risking actual plant damage.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NASA NTRS1983

A finite element method for diffusion dominated unsteady viscous flows

A general conforming finite element scheme for computing viscous flows is presented which is of second-order accuracy in space and time. Viscous terms are treated implicitly and advection terms are treated explicitly in the time marching segment of the algorithm. A method for solving the algebraic equations at each time step is given. The method is demonstrated on two test problems, one of them being a plane vortex flow for which asymptotic methods are used to obtain suitable numerical boundary conditions at each time step.

Gunzburger, M. D.↗

NASA NTRS1990

A tool for modeling concurrent real-time computation

Real-time computation is a significant area of research in general, and in AI in particular. The complexity of practical real-time problems demands use of knowledge-based problem solving techniques while satisfying real-time performance constraints. Since the demands of a complex real-time problem cannot be predicted (owing to the dynamic nature of the environment) powerful dynamic resource control techniques are needed to monitor and control the performance. A real-time computation model for a real-time tool, an implementation of the QP-Net simulator on a Symbolics machine, and an implementation on a Butterfly multiprocessor machine are briefly described.

Sharma, D. D.↗

NASA NTRS1996

Satisfiability Test with Synchronous Simulated Annealing on the Fujitsu AP1000 Massively-Parallel Multiprocessor

Solving the hard Satisfiability Problem is time consuming even for modest-sized problem instances. Solving the Random L-SAT Problem is especially difficult due to the ratio of clauses to variables. This report presents a parallel synchronous simulated annealing method for solving the Random L-SAT Problem on a large-scale distributed-memory multiprocessor. In particular, we use a parallel synchronous simulated annealing procedure, called Generalized Speculative Computation, which guarantees the same decision sequence as sequential simulated annealing. To demonstrate the performance of the parallel method, we have selected problem instances varying in size from 100-variables/425-clauses to 5000-variables/21,250-clauses. Experimental results on the AP1000 multiprocessor indicate that our approach can satisfy 99.9 percent of the clauses while giving almost a 70-fold speedup on 500 processors.

Sohn, Andrew↗

DOE OSTI2024

Parallel computing for power system climate resiliency: Solving a large-scale stochastic capacity expansion problem with mpi-sppy

Here we propose a nodal stochastic generation and transmission expansion planning model that incorporates the output from high-resolution global climate models through load and generation availability scenarios. We implement our model in Pyomo and perform computational studies on a realistically-sized test case of the California electric grid in a high performance computing environment. We propose model reformulations and algorithm tuning to efficiently solve this large problem using a variant of the Progressive Hedging Algorithm. We utilize the parallelization capabilities and overall versatility of mpi-sppy, exploiting its hub-and-spoke architecture to concurrently obtain inner and outer bounds on an optimal expansion plan. Initial results show that instances with 360 representative days on a system with over 8,000 buses can be solved to within 5% of optimality in under 4 h of wall clock time, a first step towards solving a large-scale power system expansion planning problem across a wide range of climate-informed operational scenarios.

24 POWER TRANSMISSION AND DISTRIBUTION↗

DOE OSTI2020

Parallel-in-Time Multigrid Methods for Hyperbolic Problems, with a Focus on the Shallow Water Equations

The coming massive parallelism of exascale computing presents a pressing challenge for the many DOE simulations of time-dependent partial differential equations, which typically use traditional sequential time stepping methods. Since this traditional approach is inherently serial, it presents a sequential bottleneck when moving to exascale computing, because future performance gains will come through greater concurrency, not faster clock speeds. Thus, the goal of this work is to research parallelism in time, i.e., methods that compute multiple time values simultaneously, not sequentially. The focus will be on hyperbolic problems of programmatic interest to DOE, with the goal of enabling scalable simulations of time-dependent hyperbolic problems on future architectures. The difficulty lies in the fact that hyperbolic problems are well-known to be difficult for parallel-in-time methods, with the most common parallel-in-time method, Parareal, diverging in many cases. Here, the chosen methodology for scalably solving these hyperbolic space-time equations is multigrid-reductionin-time (MGRIT), because multigrid (when it works) is a powerful, optimal, and scalable solver for discretized PDEs. To further research in this area, this project will explore a model hyperbolic problem (the shallow water equations) in the context of recent advances (e.g., by Wingate and Haut) in constructing improved coarse time-grid time-propagators by using the slow asymptotic structure of the equations. In particular, we will investigate if such coarse time-propagators offer significant advantages in an MGRIT setting, and explore any broader insights gained into parallel-in-time for hyperbolic problems.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1982

Frequency and Time Domain Modeling of Acoustic Liner Boundary Conditions

As part of a research program directed at the acoustics of advanced subsonic propulsion systems undertaken at NASA Langley, Duke University was funded to develop a boundary condition model for bulk-reacting nacelle liners. The overall objective of the Langley program was to understand and predict noise from advanced subsonic transport engines and to develop related noise control technology. The overall technical areas included: fan and propeller source noise, acoustics of ducts and duct liners, interior noise, subjective acoustics, and systems noise prediction. The Duke effort was directed toward duct liner acoustics for the development of analytical methods to characterize liner behavior in both frequency domain and time domain. A review of duct acoustics and liner technology can be found in Reference [1]. At that time, NASA Langley was investigating the propulsion concept of an advanced ducted fan, with a large diameter housed inside a relatively short duct. Fan diameters in excess of ten feet were proposed. The lengths of both the inlet and exhaust portions of the duct were to be short, probably less than half the fan diameter. The nacelle itself would be relatively thin-walled for reasons of aerodynamic efficiency. The blade-passage frequency was expected to be less than I kHz, and very likely in the 200 to 300 Hz range. Because of the design constraints of a short duct, a thin nacelle, and long acoustic wavelengths, the application of effective liner technology would be especially challenging. One of the needs of the NASA Langley program was the capability to accurately and efficiently predict the behavior of the acoustic liner. The traditional point impedance method was not an adequate model for proposed liner designs. The method was too restrictive to represent bulk reacting liners and to allow for the characterization of many possible innovative liner concepts. In the research effort at Duke, an alternative method, initially developed to handle bulk reacting layers as described in Reference [2], was extended to apply to a broad range of liner types. This method included the effect of local gradients along the liner surface, and was particularly appropriate for situations with flow over the liner and grazing incidence acoustic fields. In order to utilize time domain computational methods to solve for the propfan acoustic field, corresponding liner boundary conditions were developed for time domain solutions rather than frequency domain solutions.

Bliss, Donald B.↗

NASA NTRS1993

Real-time trajectory optimization on parallel processors

A parallel algorithm has been developed for rapidly solving trajectory optimization problems. The goal of the work has been to develop an algorithm that is suitable to do real-time, on-line optimal guidance through repeated solution of a trajectory optimization problem. The algorithm has been developed on an INTEL iPSC/860 message passing parallel processor. It uses a zero-order-hold discretization of a continuous-time problem and solves the resulting nonlinear programming problem using a custom-designed augmented Lagrangian nonlinear programming algorithm. The algorithm achieves parallelism of function, derivative, and search direction calculations through the principle of domain decomposition applied along the time axis. It has been encoded and tested on 3 example problems, the Goddard problem, the acceleration-limited, planar minimum-time to the origin problem, and a National Aerospace Plane minimum-fuel ascent guidance problem. Execution times as fast as 118 sec of wall clock time have been achieved for a 128-stage Goddard problem solved on 32 processors. A 32-stage minimum-time problem has been solved in 151 sec on 32 processors. A 32-stage National Aerospace Plane problem required 2 hours when solved on 32 processors. A speed-up factor of 7.2 has been achieved by using 32-nodes instead of 1-node to solve a 64-stage Goddard problem.

Psiaki, Mark L.↗

NASA NTRS1985

Three dimensional CSCM method for the compressible Navier-Stokes equations with application to a multi-nozzle exhaust flowfield

The effective design of three-dimensional flow systems such as a clustered multinozzle shrouded exhaust system requires detailed analysis of the internal flowfield. It is pointed out that new numerical techniques on the most modern vector computers are making these analyses practicable. In this paper, the upwind differenced implicit conservative supracharacteristics method (CSCM) for three-dimensional viscous flow is presented. The new method is an implicit 'method of planes' symmetric Gauss-Seidel relaxation scheme. The data is conveniently stored on successive planes along the streamwise coordinate, and the system of equations is solved twice in each successive plane of the streamwise coordinate. The developed new three-dimensional compressible Navier-Stokes algorithm combines the best features of the storage and computationally efficient space marching procedures with the generality of time dependent techniques to solve flows with elliptic and streamwise separated flows.

Bardina, J.↗