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 1,801 records · Page 100

All sources NASA NTRS DOE OSTI

NASA NTRS1987

NAS-current status and future plans

The Numerical Aerodynamic Simulation (NAS) has met its first major milestone, the NAS Processing System Network (NPSN) Initial Operating Configuration (IOC). The program has met its goal of providing a national supercomputer facility capable of greatly enhancing the Nation's research and development efforts. Furthermore, the program is fulfilling its pathfinder role by defining and implementing a paradigm for supercomputing system environments. The IOC is only the begining and the NAS Program will aggressively continue to develop and implement emerging supercomputer, communications, storage, and software technologies to strengthen computations as a critical element in supporting the Nation's leadership role in aeronautics.

Bailey, F. R.↗

NASA NTRS1994

Ice nucleation in the upper troposphere: Sensitivity to aerosol number density, temperature, and cooling rate

We have investigated the processes that control ice crystal nucleation in the upper troposphere using a numerical model. Nucleation of ice resulting from cooling was simulated for a range of aerosol number densities, initial temperatures, and cooling rates. In contrast to observations of stratus clouds, we find that the number of ice crystals that nucleate in cirrus is relatively insensitive to the number of aerosols present. The ice crystal size distribution at the end of the nucleation process is unaffected by the assumed initial aerosol number density. Essentially, nucleation continues until enough ice crystals are present such that their deposition growth rapidly depletes the vapor and shuts off any further nucleation. However, the number of ice crystals nucleated increases rapidly with decreasing initial temperature and increasing cooling rate. This temperature dependence alone could explain the large ice crystal number density observed in very cold tropical cirrus.

Jensen, E. J.↗

NASA NTRS2004

Technology Development and Demonstration Concepts for the Space Elevator

During the 1990s several discoveries and advances in the development of carbon nano-tube (CNT) materials indicated that material strengths many times greater than common high-strength composite materials might be possible. Progress in the development of this material led to renewed interest in the space elevator concept for construction of a tether structure from the surface of the Earth through a geostationary orbit (GEO) and thus creating a new approach to Earth-to-orbit transportation infrastructures. To investigate this possibility the author, in 1999, managed for NASA a space elevator work:hop at the Marshall Space Flight Center to explore the potential feasibility of space elevators in the 21 century, and to identify the critical technologies and demonstration missions needed to make development of space elevators feasible. Since that time, a NASA Institute for Advanced Concepts (NIAC) funded study of the Space Elevator proposed a concept for a simpler first space elevator system using more near-term technologies. This paper will review some of the latest ideas for space elevator development, the critical technologies required, and some of the ideas proposed for demonstrating the feasibility for full-scale development of an Earth to GEO space elevator. Critical technologies include CNT composite materials, wireless power transmission, orbital object avoidance, and large-scale tether deployment and control systems. Numerous paths for technology demonstrations have been proposed utilizing ground experiments, air structures. LEO missions, the space shuttle, the international Space Station, GEO demonstration missions, demonstrations at the lunar L1 or L2 points, and other locations. In conclusion, this paper finds that the most critical technologies for an Earth to GEO space elevator include CNT composite materials development and object avoidance technologies; that lack of successful development of these technologies need not preclude continued development of space elevator systems in general; and that the critical technologies required for the Earth to GEO space elevator are not required for similar systems at the Moon, Mars, Europa, or for orbital tether systems at GEO, Luna, and other locations.

Smitherman, David V., Jr.↗

NASA NTRS2015

Expanding the Operational Use of Total Lightning Ahead of GOES-R

NASA's Short‐term Prediction Research and Transition Center (SPoRT) has been transitioning real‐time total lightning observations from ground‐based lightning mapping arrays since 2003. This initial effort was with the local Weather Forecast Offices (WFO) that could use the North Alabama Lightning Mapping Array (NALMA). These early collaborations established a strong interest in the use of total lightning for WFO operations. In particular the focus started with warning decision support, but has since expanded to include impact‐based decision support and lightning safety. SPoRT has used its experience to establish connections with new lightning mapping arrays as they become available. The GOES‐R / JPSS Visiting Scientist Program has enabled SPoRT to conduct visits to new partners and expand the number of operational users with access to total lightning observations. In early 2014, SPoRT conducted the most recent visiting scientist trips to meet with forecast offices that will used the Colorado, Houston, and Langmuir Lab (New Mexico) lightning mapping arrays. In addition, SPoRT met with the corresponding Center Weather Service Units (CWSUs) to expand collaborations with the aviation community. These visits were an opportunity to learn about the forecast needs of each office visited as well as to provide on‐site training for the use of total lightning, setting the stage for a real‐time assessment during May‐July 2014. With five lightning mapping arrays covering multiple geographic locations, the 2014 assessment has demonstrated numerous uses of total lightning in varying situations. Several highlights include a much broader use of total lightning for impact‐based decision support ranging from airport weather warnings, supporting fire crews, and protecting large outdoor events. The inclusion of the CWSUs has broadened the operational scope of total lightning, demonstrating how these data can support air traffic management, particularly in the Terminal Radar Approach Control Facilities (TRACON) region around an airport. These collaborations continue to demonstrate, from the operational perspective, the utility of total lightning and the importance of continued training and preparation in advance of the Geostationary Lightning Mapper.

GLM↗

NASA NTRS2017

Near Real-Time Probabilistic Damage Diagnosis Using Surrogate Modeling and High Performance Computing

This work investigates novel approaches to probabilistic damage diagnosis that utilize surrogate modeling and high performance computing (HPC) to achieve substantial computational speedup. Motivated by Digital Twin, a structural health management (SHM) paradigm that integrates vehicle-specific characteristics with continual in-situ damage diagnosis and prognosis, the methods studied herein yield near real-time damage assessments that could enable monitoring of a vehicle's health while it is operating (i.e. online SHM). High-fidelity modeling and uncertainty quantification (UQ), both critical to Digital Twin, are incorporated using finite element method simulations and Bayesian inference, respectively. The crux of the proposed Bayesian diagnosis methods, however, is the reformulation of the numerical sampling algorithms (e.g. Markov chain Monte Carlo) used to generate the resulting probabilistic damage estimates. To this end, three distinct methods are demonstrated for rapid sampling that utilize surrogate modeling and exploit various degrees of parallelism for leveraging HPC. The accuracy and computational efficiency of the methods are compared on the problem of strain-based crack identification in thin plates. While each approach has inherent problem-specific strengths and weaknesses, all approaches are shown to provide accurate probabilistic damage diagnoses and several orders of magnitude computational speedup relative to a baseline Bayesian diagnosis implementation.

Warner, James E.↗

NASA NTRS1983

A numerical model of acoustic choking. I - Shock free solutions

The phenomenon of acoustic choking in near sonic flows is investigated using the one-dimensional gasdynamic equations for an ideal gas. It is pointed out that this approach eliminates the need to make the classical small disturbance assumption. A finite difference scheme is elaborated to approximate these equations. Boundary conditions that correctly model the physics of the acoustics problems are chosen. The use of a fast Fourier transform routine allows the results to be compared with those obtained by using the conventional harmonic assumptions. Only solutions in which shocks have not fully developed are reported. Continuous solutions are sought, and these are compared with the Fubini solution and with certain results of Myers and Callegari (1977, 1978, 1979).

Walkington, N. J.↗

DOE OSTI2026

Harnessing the Second-Order Metal−Insulator Transition for Neuromorphic Computing

Vanadium oxides are widely studied phase change materials for brain-inspired computing architectures. Systems like VO 2 and V 2 O 3 exhibit first-order metal−insulator transitions (MITs) with hysteresis and percolative switching, increasing stochasticity and device variability. Here, we focus on the less explored Magnéli phase V 4 O 7 , which undergoes a continuous, non-hysteretic, second-order MIT. This surprisingly enables highly reproducible volatile resistive switching in spiking-neuron-type devices. We synthesize V 4 O 7 films, characterize their structural and transport properties, and demonstrate voltage and current-driven threshold switching with electrothermal feedback. In a Pearson–Anson oscillator, V 4 O 7 devices produce stable, tunable spiking across 20–200 kHz, with consistent operation among multiple devices. We introduce a numerical analog leaky-integrate-and-fire (aLIF) model that captures waveform shapes and their dependence on load resistance, temperature, and voltage. Furthermore, these findings suggest that second-order MIT materials like V 4 O 7 are promising for deterministic, scalable spiking neuron arrays for neuromorphic computing.

V4O7↗

DOE OSTI2023

On fast simulation of dynamical system with neural vector enhanced numerical solver

The large-scale simulation of dynamical systems is critical in numerous scientific and engineering disciplines. However, traditional numerical solvers are limited by the choice of step sizes when estimating integration, resulting in a trade-off between accuracy and computational efficiency. To address this challenge, we introduce a deep learning-based corrector called Neural Vector (NeurVec), which can compensate for integration errors and enable larger time step sizes in simulations. Our extensive experiments on a variety of complex dynamical system benchmarks demonstrate that NeurVec exhibits remarkable generalization capability on a continuous phase space, even when trained using limited and discrete data. NeurVec significantly accelerates traditional solvers, achieving speeds tens to hundreds of times faster while maintaining high levels of accuracy and stability. Moreover, NeurVec’s simple-yet-effective design, combined with its ease of implementation, has the potential to establish a new paradigm for fast-solving differential equations based on deep learning.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2021

An experimental database of cell performance for vanadium redox flow battery

The continual growth in energy demand has resulted in the deployment of renewable energy generators to reduce the impact of fossil fuel dependence. However, these generators often suffer from intermittency and require energy storage when there is over-generation and the subsequent release of this stored energy at high demand. One promising energy storage technology which can provide a solution to improve energy management and grid stability, is the redox flow battery. Among the numerous flow battery systems, vanadium redox flow battery is the most iconic solution to large scale energy storage, giving a more efficient link between energy production, especially from renewables, and energy demand. The aim of the current database is to characterize the performance of the cell design and to provide training/validation data for physical model or data-driven model. The database includes hundreds of experimental cell performance data of vanadium redox flow battery with various current densities for multiple charge-discharge cycles. All the cell parameters, chemical parameters, material parameters, operation parameters and thermodynamic parameters of the cell system are listed. Coulomb, voltaic and energy efficiencies are also provided. The database will be helpful for researchers in the field of redox flow batteries.

Gao, Peiyuan↗

NASA NTRS1986

A comparison of the regional slope characteristics of Venus and earth - Implications for geologic processes on Venus

The range of 3 degree by 3 degree regional slopes of the Earth and Venus is similar (approximately 0.0-2.4 degrees), although the surface distribution of these values differs significantly. On earth, cratonic and abyssal plains form extensive regions of 0.0 degree slope. Within these regions a variety of features (mid-ocean ridges, volcanic island chains, subduction zones, and floded mountains) have regional slope characteristics influenced by seafloor spreading and plate recycling, as well as an active weathering regime. The plains provinces of Venus are much more rugged than earth's plains and are marked by numerous closely spaced circular and linear features (0.1-0.2 degree regional slope) concentrated into broad linear zones of global extent. Although Venus highlands are bounded by narrow zones of relatively steep slope, the margins of Aphrodite Terra and Beta Regio are not as steep as earth's continental margins and appear to be best developed parallel to the trends of major chasmata within these regions. Ishtar Terra's margins are significantly steeper and more continuous than other highland margins and are comparable to passive margins on earth. The Venus highlands do not contain appreciable smooth, flat interior regions, implying that highland topography is not significantly modified by erosion or deposition.

Sharpton, V. L.↗

NASA NTRS1989

NASA/MSFC FY88 Global Scale Atmospheric Processes Research Program Review

Interest in environmental issues and the magnitude of the environmental changes continues. One way to gain more understanding of the atmosphere is to make measurements on a global scale from space. The Earth Observation System is a series of new sensors to measure globally atmospheric parameters. Analysis of satellite data by developing algorithms to interpret the radiance information improves the understanding and also defines requirements for these sensors. One measure of knowledge of the atmosphere lies in the ability to predict its behavior. Use of numerical and experimental models provides a better understanding of these processes. These efforts are described in the context of satellite data analysis and fundamental studies of atmospheric dynamics which examine selected processes important to the global circulation.

Wilson, Greg S.↗

NASA NTRS1979

Numerical simulation of turbulent flows with a three-dimensional vortex-in-cell method

A three-dimensional vortex-in-cell method has been developed for the evaluation of local flow fields due to a family of vortex filaments which employs the principles and architecture of a code developed for magnetic field evaluation in plasma simulations. The computational effort in the new method, as compared to 'vortex pushing' by direct Biot-Savarat interaction, increases directly, rather than quadratically, with the number of vortex elements. The method is well suited for studying large number of vortex filaments or rings and can be used for simulating continuous vorticity.

Buneman, O.↗

DOE OSTI2026

MFC 5.0: An exascale many-physics flow solver

Many problems of interest in engineering, medicine, and the fundamental sciences rely on high-fidelity flow simulation, making performant computational fluid dynamics solvers a mainstay of the open-source software community. Previous work MFC 3.0 was made a published, documented, and open-source solver via Bryngelson et al. Comp. Phys. Comm. (2021) with numerous physical features, numerical methods, and scalable infrastructure. MFC 5.0 is a significant update to MFC 3.0, featuring a broad set of well-established and novel physical models and numerical methods, as well as the introduction of GPU and APU (or superchip) acceleration. Here, we exhibit state-of-the-art performance and ideal scaling on the first two exascale supercomputers, OLCF Frontier and LLNL El Capitan. Combined with MFC’s single-accelerator performance, MFC achieves exascale computation in practice, and achieved the largest-to-date public CFD simulation at 200 trillion grid points as a 2025 ACM Gordon Bell Prize finalist. New physical features include the immersed boundary method, N-fluid phase change, Euler–Euler and Euler–Lagrange sub-grid bubble models, fluid-structure interaction, hypo- and hyper-elastic materials, chemically reacting flow, two-material surface tension, magnetohydrodynamics (MHD), and more. Numerical techniques now represent the current state-of-the-art, including general relaxation characteristic boundary conditions, WENO variants, Strang splitting for stiff sub-grid flow features, and low Mach number treatments. Weak scaling to tens of thousands of GPUs on OLCF Summit and Frontier and LLNL El Capitan achieves efficiencies within 5% of ideal to over 90% of their respective system sizes. Strong scaling results for a 16-times increase in device count show parallel efficiencies over 90% on OLCF Frontier. MFC’s software stack has undergone further improvements, including continuous integration, which ensures code resilience and correctness through over 300 regression tests; metaprogramming, which reduces code length while maintaining performance portability; and code generation for computing chemical reactions

Computational fluid dynamics↗

DOE OSTI2025

A new tubular ironless linear induction motor using non-interrupted winding method

The ironless linear induction motor device can be used for any application which requires a large motion span or rotation restrictedly by the ironless feature due to the strong magnetic field environment. One of the possible applications is ITER diagnostics shutters. The tokamak reactor systems have a severe radiation issue during the operation. The traditional Kapton-insulated wires have significantly lower radiation resistance than an expected large radiation dose at ITER Motion Stark Effect (MSE) shutter location. Furthermore, the traditional single-sided linear induction motor winds the wire to make the N-turn pitch coils independently and then connect them in series to make a one-phase winding. Therefore, it is a big challenge to make so many radiation-protection connectors in such a limited space as inside the vacuum vessel of fusion reactors. The approved mineral-insulated kinds of cables need a large bend radius for the end turns, so the traditional single-sided linear motor with two-layer three-phase windings need to be upgraded to the tubular linear motors. The new proposed tubular ironless linear induction motor can not only overcome the large bend radius issue with the mineral-insulated cables but also allow the motor coils wound continuously around the tubular-shaped stator, forward and backward N times to make a N-turn winding for each phase, in order to resolve the radiation protection issue with the winding connection. The 3D transient Maxwell electromagnetic model with a tubular three-phase linear induction motor has been analyzed. Furthermore, this paper will present the detailed design and winding method of ironless linear induction motor actuators, and results of the 3D transient numerical simulation of the actuator performance.

Ironless↗

DOE OSTI2023

Finite-difference time-domain methods

The finite-difference time-domain (FDTD) method is a widespread numerical tool for full-wave analysis of electromagnetic fields in complex media and for detailed geometries. Applications of the FDTD method cover a range of time and spatial scales, extending from subatomic to galactic lengths and from classical to quantum physics. Technology areas that benefit from the FDTD method include biomedicine — bioimaging, biophotonics, bioelectronics and biosensors; geophysics — remote sensing, communications, space weather hazards and geolocation; metamaterials — sub-wavelength focusing lenses, electromagnetic cloaks and continuously scanning leaky-wave antennas; optics — diffractive optical elements, photonic bandgap structures, photonic crystal waveguides and ring-resonator devices; plasmonics — plasmonic waveguides and antennas; and quantum applications — quantum devices and quantum radar. This Primer summarizes the main features of the FDTD method, along with key extensions that enable accurate solutions to be obtained for different research questions. Additionally, hardware considerations are discussed, plus examples of how to extract magnitude and phase data, Brillouin diagrams and scattering parameters from the output of an FDTD model. Furthermore, the Primer ends with a discussion of ongoing challenges and opportunities to further enhance the FDTD method for current and future applications.

42 ENGINEERING↗

NASA NTRS1990

Numerical simulation of a catastrophe model for coronal mass ejections

Van Tend and Kuperus (1978) proposed a simple catastrophe model for magnetically driving coronal mass ejections, prominence eruptions, and two-ribbon flares. Their model, which is based on simple circuit concepts, suggests that a stable configuration containing a current filament will lose equilibrium when the filament current exceeds a critical value. Here, a two-dimensional numerical simulation is used to test how the Van Tend-Kuperus model works in an ideal MHD fluid. The simulation exhibits the expected loss of mechanical equilibrium near the predicted critical value, but the current filament moves only a short distance upward before coming to rest at a new equilibrium. However, this new equilibrium contains a current sheet which is resistively unstable to magnetic reconnection, and if magnetic reconnection occurs rapidly, the filament can continue to move upward at Alfvenic speeds.

Forbes, T. G.↗

NASA NTRS1991

Robustness in the presence of mixed parametric uncertainty and unmodeled dynamics

Continuing the development of the structured singular value approach to robust control design, the authors investigate the problem of computing the structured singular value (mu) in the case of mixed real parametric and complex uncertainty. The problem is shown to be equivalent to a smooth constrained finite-dimensional optimization problem. In view of the fact that the functional to be maximized may have several local extrema, an upper bound on mu whose computation is numerically tractable is established; this leads to a sufficient condition of robust stability and performance. A historical perspective on the development of the mu theory is included.

Fan, Michael K. H.↗

NASA NTRS2000

Application of Three-Component PIV to a Hovering Rotor Wake

The key to accurate predictions of rotorcraft aerodynamics, acoustics, and dynamics lies in the accurate representation of the rotor wake. The vortical wake computed by rotorcraft CFD analyses typically suffer from numerical dissipation before the first blade passage. With some a priori knowledge of the wake trajectory, grid points can be concentrated along the trajectory to minimize the dissipation. Comprehensive rotorcraft analyses based on lifting-line theory rely on classical vortex models and/or semi-empirical information about the tip vortex structure. Until the location, size, and strength of the trailed tip vortex can be measured over a range of wake ages, the analyses will continue to be adjusted on a trial and error basis in order to correctly predict blade airloads, acoustics, dynamics, and performance. Using the laser light sheet technique, tip vortex location can be acquired in a straightforward manner. Measuring wake velocities and vortex core size, however, has been difficult and tedious using point-measurement techniques such as laser velocimetry. Recently, the Particle Image Velocimetry (PIV) technique has proven to be an efficient method for acquiring velocity measurements over relatively large areas and volumes of a rotor wake. The work reported to date, however, has been restricted to 2-component velocity measurements of the rotor wake. Three-component velocity measurements of a hovering rotor wake were acquired at NASA Ames Research Center in May 1999. This experiment represents a major step toward understanding the detailed structure of a rotor wake. This paper will focus primarily on the experimental technique used in acquiring this data. The accuracy and limitations of the current technique will also be discussed. Representative velocity field measurements will be included.

Yamauchi, Gloria K.↗

Refine your search to explore more results.