Engineering Papers⌕ Search

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.

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

At least 1,693 records · Page 94
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Multidisciplinary Optimization of an Electric Quadrotor Urban Air Mobility Aircraft

Urban Air Mobility (UAM) vehicles have the potential to augment urban transportation systems, allowing passengers to skip the traffic below for a fee. This emerging market is opening up the design space for a new class of Urban Air Mobility (UAM) vehicles which could be powered by electric propulsion systems to be economical and environmentally friendly. However, development of these UAM concepts presents several additional challenges in the design process. First, these concept designs require including new disciplinary models for subsystems including the electric motors, cables, batteries and thermal management systems. Second, correctly designing and evaluating these various subsystems requires tight coupling between the discipline models to capture interactions. This paper presents the continued development of a multidisciplinary design optimization environment to aid in the development of these vehicle concepts. The multidisciplinary environment fully couples the various subsystem models allowing for the full vehicle to be designed and optimized simultaneously. In this research, the developed modeling approach is demonstrated in the analysis of a small, all-electric quadrotor UAM concept. Results from these studies show that numerous disciplines can be tightly coupled and producing improved overall vehicle designs.

Multidisciplinary Optimization↗
NASA NTRSDate not supplied

Quality Control Methods of Tower Data at Kennedy Space Center’s Launch Complex 39-B and the USAF Cape Canaveral 500 Foot Tower

The National Aeronautics and Space Administration (NASA) has long used meteorological data from weather towers located at Kennedy Space Center (KSC) and the United States Air Force’s Eastern Range(ER) in support of their various launch vehicles and climate studies. Some of the most valuable data is gathered at the Launch Complex 39-B (LC39-B) and the 500 foot tower located approximately two miles west of LC39-B and three miles north of the Vehicle Assembly Building (VAB). The data from LC39-B gives NASA Space Launch System (SLS) engineers valuable insight into the weather the vehicle can expect at the pad, while the 500 foot tower has a longer period of record and has been used by the Shuttle and SLS programs. However, numerous data quality (DQ) control issues have arisen when utilizing this data for analysis. In addition, instrumentation tends to deteriorate faster than normal due to the corrosive nature of the high salt content in the air. NASA’s Natural Environments branch has developed several QC databases of these towers, but these studies were for fixed periods of records. EV44 identified a need for a continually updated and QC’ed database of tower data not only to provide along term QC’ed database for vehicle and climate analyses, but also to provide the capability to investigate recent weather events such as downbursts and other high wind events. This study follows the methods developed in prior Natural Environment QC tower databases, which includes variable specific QC thresholds and checks used to generate QC flags. However, this study developed new techniques such as the development and design of the graphical user interface (GUI) for manual verification, and the creation of the final monthly QC files.

Quality Control↗
NASA NTRSDate not supplied

Coastal Land Change Due to Tectonic Processes and Implications for Relative Sea-Level Rise in the Samoan Islands

Of the major coastal land change mechanisms responsible for relative sea-level change, tectonic subsidence is generally quoted as ranging from < mm/yr to 1 cm/yr. However, we documented coseismic and ongoing post-earthquake surface displacements from continuous GPS and tide gauge/altimetry data that indicated rapid subsidence on two of the major Samoan Islands of 12 - 20 cm during and following the 8.1 2009 Tonga-Samoa earthquake. Earlier results and our modeling of GRACE-derived gravimetric data provided a preliminary forecast of future relative sea-level rise through rapid land subsidence [Han et al., 2019]. Of course these numerical forecasts of time-dependent deformation are only as good as our input observations and our assumed rheological models. As part of our current NASA Earth Surface and Interior study, we are obtaining a wider range of data to constrain and test alternate models of ongoing postseismic deformation across American Samoa and Upolu, Samoa: (1) times series of altimetry plus tide gauge data processed to complement the cGPS data available to provide high-temporal resolution, point measurements of uplift/subsidence, (2) InSAR derived observations of surface deformation across the highly vegetated Samoan Islands, (3) evaluating and using NASA satellite lidar data (ICESat-I & ICESat-II, GEDI) for fusion with multi-source topographic data sets and for estimating topographic change on the decadal time scale. We are evaluating and using these new observations to better understand and separate out local, island-wide, and multi-island subsidence patterns and to evaluate the high impact of rising sea-level in a tectonically active region.

Jeanne Sauber-Rosenberg↗
DOE OSTI2021

Entanglement transitions from restricted Boltzmann machines

The search for novel entangled phases of matter has lead to the recent discovery of a new class of “en- tanglement transitions,” exemplified by random tensor networks and monitored quantum circuits. Most known examples can be understood as some classical ordering transitions in an underlying statistical mechanics model, where entanglement maps onto the free-energy cost of inserting a domain wall. In this paper we study the possibility of entanglement transitions driven by physics beyond such statistical mechanics mappings. Motivated by recent applications of neural-network-inspired variational Ansätze, we investigate under what conditions on the variational parameters these Ansätze can capture an entanglement transition. We study the entanglement scaling of short-range restricted Boltzmann machine (RBM) quantum states with random phases. For uncorrelated random phases, we analytically demonstrate the absence of an entanglement transition and reveal subtle finite-size effects in finite-size numerical simulations. Introducing phases with correlations decaying as 1/r α in real space, we observe three regions with a different scaling of entanglement entropy depending on the exponent α. We study the nature of the transition between these regions, finding numerical evidence for critical behavior. Furthermore, our work establishes the presence of long-range correlated phases in RBM-based wave functions as a required ingredient for entanglement transitions.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

A face-upwinded spectral element method

Here we present a new high-order accurate discretisation on unstructured meshes of quadrilateral elements. Our Face Upwinded Spectral Element (FUSE) method uses the same node distribution as a high-order continuous Galerkin (CG) method, but with a particular choice of node locations within each element and an upwinded stencil on the face nodes. This results in a number of benefits, including fewer degrees of freedom and straight-forward integration with CG. We present the derivation of the scheme and the analysis of its properties, in particular showing stability using von Neumann analysis. We show numerical evidence for its accuracy and efficiency on multiple classes of problems including convection-dominated flows, Poisson's equation, and the incompressible Navier-Stokes equations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Towards a classification of bifurcations in Vlasov equations

Here, we propose a classification of bifurcations of Vlasov equations, based on the strength of the resonance between the unstable mode and the continuous spectrum on the imaginary axis. We then identify and characterize a new type of generic bifurcation where this resonance is weak, but the unstable mode couples with a stable mode and a Casimir invariant of the system to form a size-3 Jordan block. We derive a three-dimensional reduced noncanonical Hamiltonian system describing this bifurcation. Comparison of the reduced dynamics with direct numerical simulations on a test case gives excellent agreement. We finally discuss the relevance of this bifurcation to specific physical situations.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS1979

Two-dimensional zonal mean flow model

The interactions of dynamics, photochemistry and radiation in the stratosphere can be described by the continuity, momentum, and energy equations. A more comprehensive model for stratospheric transport theory was developed for the purpose of aiding predictions of changes in the stratospheric ozone content as a consequence of natural and anthropogenic processes. This model is time dependent and the dependent variables are zonal means of the relevant meteorological quantities which are functions of latitude and height. The detailed formulation of a numerical model both in physics and mathematics is given. A set of fundamental dynamical equations, and the numerical method used in the integration are described.

Chen, H. C.↗
NASA NTRS1990

A current disruption mechanism in the neutral sheet - A possible trigger for substorm expansions

A linear analysis is performed to investigate the kinetic cross-field streaming instability in the earth's magnetotail neutral sheet region. Numerical solution of the dispersion equation shows that the instability can occur under conditions expected for the neutral sheet just prior to the onset of substorm expansion. The excited waves are obliquely propagating whistlers with a mixed polarization in the lower hybrid frequency range. The ensuing turbulence of this instability can lead to a local reduction of the cross-tail current causing it to continue through the ionosphere to form a substorm current wedge. A substorm expansion onset scenario is proposed based on this instability in which the relative drift between ions and electrons is primarily due to unmagnetized ions undergoing current sheet acceleration in the presence of a cross-tail electric field. The required electric field strength is within the range of electric field values detected in the neutral sheet region during substorm intervals. The skew in local time of substorm onset location and the three conditions under which substorm onset is observed can be understood on the basis of the proposed scenario.

Lui, A. T. Y.↗
NASA NTRS1997

Extinguishment of a Diffusion Flame Over a PMMA Cylinder by Depressurization in Low-Gravity

The behavior of flames in low-speed flows in low-gravity is relevant to spacecraft fire safety. Previous work has shown that flames in the presence of low-speed forced flows in low-gravity may be more flammable than in a forced flow of the same magnitude in normal gravity. Additionally, fire suppression plans for the International Space Station include the use of venting (depressurization) as an emergency option for extinguishing fires. This procedure would induce flows in the affected compartment that could temporarily intensify the fire, as was observed in flammability tests of solids conducted on board Skylab. Despite a general Understanding, current knowledge of the combined effects of reduced pressure and forced flow on a burning solid in low-gravity is inadequate for the design of a venting extinguishment system. Previous studies in low-g have examined flammability limits for thermally thin solids. However, there are differences when burning thick materials because the interior solid-phase temperature continuously changes, which affects the percentage of gas-phase heat feedback to the solid-phase. Changes in the heat feedback to the solid-phase can affect the flammability characteristics of the material. In the current work, the extinction of a diffusion flame burning over PMMA (polymethyl methacrylate) cylinders during depressurization with a low-speed cross flow was examined experimentally and via numerical simulations.

Goldmeer, Jeffrey S.↗
NASA NTRS2002

On the Numerical Convergence to Steady State of Hypersonic Flows Over Bodies with Concavities

Two recent numerical studies of hypersonic flows over bodies with concavities revealed problems with convergence to a steady state with an oft-used application of local-time-stepping. Both simulated flows showed a time-like, periodic shedding of vortices in a subsonic domain bounded by supersonic external flow although the simulations, using local-time-stepping, were not time accurate. Simple modifications to the numerical algorithm were implemented to enable implicit, first-order accurate in time simulations. Subsequent time-accurate simulations of the two test problems converged to a steady state. The baseline algorithm and modifications for temporal accuracy are described. The requirement for sub-iterations to achieve convergence is demonstrated. Failure to achieve convergence without time accuracy is conjectured to arise from temporal errors being continuously refocused into a subsonic domain.

Gnoffo, Peter A.↗
NASA NTRS2015

Landsat-8 Sensor Characterization and Calibration

Landsat-8 was launched on 11 February 2013 with two new Earth Imaging sensors to provide a continued data record with the previous Landsats. For Landsat-8, pushbroom technology was adopted, and the reflective bands and thermal bands were split into two instruments. The Operational Land Imager (OLI) is the reflective band sensor and the Thermal Infrared Sensor (TIRS), the thermal. In addition to these fundamental changes, bands were added, spectral bandpasses were refined, dynamic range and data quantization were improved, and numerous other enhancements were implemented. As in previous Landsat missions, the National Aeronautics and Space Administration (NASA) and United States Geological Survey (USGS) cooperated in the development, launch and operation of the Landsat- 8 mission. One key aspect of this cooperation was in the characterization and calibration of the instruments and their data. This Special Issue documents the efforts of the joint USGS and NASA calibration team and affiliates to characterize the new sensors and their data for the benefit of the scientific and application users of the Landsat archive. A key scientific use of Landsat data is to assess changes in the land-use and land cover of the Earth's surface over the now 43-year record. In order to perform these analyses and avoid confusing sensor changes with Earth surface changes, a solid understanding of the sensors' performance, consistent geolocation and radiometry are essential. Particularly with the significant changes in the Landsat-8 sensors relative to previous Landsat missions, this characterization becomes all the more important.

Radiometry↗
DOE OSTI2020

Wash Penalty Factors: Insight into the Effects of Water Wash Down of Building Materials for Radionuclide/Nuclear (Rn) Decontamination - 20153

The release of radiological material to the urban environment can have a devastating socio-economic impact. In an urban release scenario, areas around the dispersal epicenter and any areas contaminated via atmospheric deposition may require extensive decontamination. Most recently, our group reported the decontamination efficacy of pressurized water washing of common porous building materials, and compared the penetration depth of the contaminations before and after exposure to pressurized tap water. Here, we continue evaluation of the decontamination efficacy of wash solutions based on potassium (K{sup +}) salt for cesium (Cs), strontium (Sr), and europium (Eu) contaminations on common building surfaces. The surfaces included: the cement used in concrete, the aggregate used in concrete, intact concrete, roadway asphalt, asphalt shingles, and latex paint. We deposited solutions of {sup 137}Cs, {sup 85}Sr, and {sup 152}Eu onto numerous common building materials and exposed these materials to a static bath or low-pressure flow of tap water, 0.1 M potassium chloride (KCl), and 0.5 M KCl. The decontamination efficacy and the depth profile for residual contamination were measured to determine the conditions under which applying a wash solution has benefit compared to physically removing the surface material. To aid in these analyses, we introduced a 'wash penalty factor' to quantify the degree to which contamination migrated into the bulk material by the action of the wash method. On asphalt, 70-80% of the radionuclides were found to be within 0.02 mm of the surface. Cs{sup +} bonded to negative charge sites on the surface of the mineral aggregate (used in concrete) so decontamination was improved by ion exchange with K{sup +}. Sr{sup 2+} is soluble and removed by water, while not bonding to aggregate or demonstrating solubility in the liquid asphalt. Eu{sup 3+} is precipitated and not removed in gentle wash but displays mobility within the pore space that suggests its presence as a colloid. Concrete is more porous than asphalt, and 80% of the radionuclides were within 0.2 mm of the surface for {sup 137}Cs and {sup 152}Eu and 50- 80% for {sup 85}Sr. On concrete, Cs{sup +} bonds similarly to asphalt. Sr{sup 2+} bonds to the cement or precipitates from the alkaline pore water. Eu{sup 3+} sorbs to the cement and then precipitates. In brick, the results for Cs{sup +} and Sr{sup 2+} suggested that diffusion of ions from deeper in the coupon was important. Eu{sup 3+} is removed poorly due to precipitation. We will explain and report example wash penalty factors that vary from 0 to > 1 for the various methods and the implications for deciding when to implement wash down methods. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2026

Advancing the HERO WEC Through Integrated Modeling, Testing, and Field Deployments: Preprint

The Hydraulic and Electric Reverse Osmosis Wave Energy Converter (HERO WEC) was developed by the National Laboratory of the Rockies as a modular platform for advancing wave-powered desalination technologies. Designed to operate in either a hydraulic or electric configuration, the system enables direct comparison of mechanical-to-water and electrical-to-water conversion pathways using a single hardware architecture. Deployments at the Jennette's Pier test site in 2022 and 2024 demonstrated freshwater production between 60 and 300 gallons per day, including successful operation in wave heights as low as 0.25 m. A structured evaluation approach combining numerical modeling of hydrodynamic and PTO response, controlled laboratory testing, and in-ocean field deployments has been used to characterize and refine system performance. Insights from these efforts are now informing the second-generation HERO WEC (V2), which incorporates improved drivetrain components, refined geometry, enhanced control systems, and design updates aimed at increasing robustness and long-duration survivability. The HERO WEC platform continues to serve as an open-access research asset supporting broader marine energy and desalination development.

16 TIDAL AND WAVE POWER↗
NASA NTRS1982

Geometrically derived difference formulae for the numerical integration of trajectory problems

An initial value problem for the autonomous system of ordinary differential equations dy/dt = f(y), where y is a vector, is considered. In a number of practical applications the interest lies in obtaining the curve traced by the solution y. These applications include the computation of trajectories in mechanical problems. The term 'trajectory problem' is employed to refer to these cases. Lambert and McLeod (1979) have introduced a method involving local rotation of the axes in the y-plane for the two-dimensional case. The present investigation continues the study of difference schemes specifically derived for trajectory problems. A simple geometrical way of constructing such methods is presented, and the local accuracy of the schemes is investigated. A circularly exact, fixed-step predictor-corrector algorithm is defined, and a variable-step version of a circularly exact algorithm is presented.

Mcleod, R. J. Y.↗
DOE OSTI2020

An energy-conserving and asymptotic-preserving charged-particle orbit implicit time integrator for arbitrary electromagnetic fields

Here we present a new implicit asymptotic preserving time integration scheme for charged-particle orbit computation in arbitrary electromagnetic fields. The scheme is built on the Crank-Nicolson integrator and continues to recover full-orbit motion in the small time-step limit, but also recovers all the first-order guiding center drifts as well as the correct gyroradius when stepping over the gyration time-scale. In contrast to previous efforts in this direction, the new scheme also features exact energy conservation. In the derivation of the scheme, we find that a new numerical time-scale is introduced. This scale is analyzed and the resulting restrictions on time-step are derived. Based on this analysis, we develop an adaptive time-stepping strategy the respects these constraints while stepping over the gyration scale when physically justified. It is shown through numerical tests on single-particle motion that the scheme's energy conservation property results in tremendous improvements in accuracy, and that the scheme is able to transition smoothly between magnetized and unmagnetized regimes as a result of the adaptive time-stepping.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1972

Subsonic, transonic, and supersonic nozzle flow by the inverse technique.

The inverse technique is used to obtain a mathematically and physically consistent solution of the flowfield in a nozzle from the mass generation surface through the supersonic region. The inverse method employs an assumed centerline function which is of the Cauchy type in that the values and the derivatives of the function are known. Since the Cauchy boundary conditions can give rise to numerical instabilities, the governing gasdynamic equations for rotational steady flow were transformed into a form which puts the geometry into a rectangular shape, and which spaces the network of interior points more finely in regions of the greatest gradients of the dependent variables. For arbitrarily specified centerline data, the solution of the governing flow equations may not exist, and if it does it may not depend continuously on the data.

Norton, D. J.↗
NASA NTRS1974

Plastic buckling

The present article is concerned mainly with the post-bifurcation and imperfection-sensitivity aspects of plastic buckling. A simple two-degree-of-freedom model is used to introduce post-bifurcation behavior and a second model illustrates features of the behavior of continuous solids and structures. Hill's bifurcation criterion for a class of three-dimensional solids is applied to the Donnell-Mushtari-Vlasov (DMV) theory of plates and shells. A general treatment of the initial post-bifurcation behavior of plates and shells is given within the context of the DMV theory. This is illustrated by problems involving columns and circular plates under radial compression. Numerical results are given for a column under axial compression, a circular plate under radial compression, and spherical and cylindrical shells.

Hutchinson, J. W.↗
NASA NTRS1974

Interpretation of air pollution data as measured by an airborne remote sensor

The investigation described is a continuation of the work reported by Smith et al. (1974) in which a single source was studied. In the current study, multiple sources of known location are considered. The study is concerned with the strength of each source and the resulting pollution concentration field. The characteristics of the remotely sensed data are discussed along with the parameter estimation procedure, the estimation of pollution parameters, and a numerical example.

Smith, G. L.↗