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,495 records · Page 83

Application of dynamical systems theory to nonlinear aircraft dynamics

Dynamical systems theory has been used to study nonlinear aircraft dynamics. A six degree of freedom model that neglects gravity has been analyzed. The aerodynamic model, supplied by NASA, is for a generic swept wing fighter and includes nonlinearities as functions of the angle of attack. A continuation method was used to calculate the steady states of the aircraft, and bifurcations of these steady states, as functions of the control deflections. Bifurcations were used to predict jump phenomena and the onset of periodic motion for roll coupling instabilities and high angle of attack maneuvers. The predictions were verified with numerical simulations.

Culick, Fred E. C.↗

Active Control of the Forced and Transient Response of a Finite Beam

When studying structural vibrations resulting from a concentrated source, many structures may be modelled as a finite beam excited by a point source. The theoretical limit on cancelling the resulting beam vibrations by utilizing another point source as an active controller is explored. Three different types of excitation are considered, harmonic, random, and transient. In each case, a cost function is defined and minimized for numerous parameter variations. For the case of harmonic excitation, the cost function is obtained by integrating the mean squared displacement over a region of the beam in which control is desired. A controller is then found to minimize this cost function in the control interval. The control interval and controller location are continuously varied for several frequencies of excitation. The results show that control over the entire beam length is possible only when the excitation frequency is near a resonant frequency of the beam, but control over a subregion may be obtained even between resonant frequencies at the cost of increasing the vibration outside of the control region. For random excitation, the cost function is realized by integrating the expected value of the displacement squared over the interval of the beam in which control is desired. This is shown to yield the identical cost function as obtained by integrating the cost function for harmonic excitation over all excitation frequencies. As a result, it is always possible to reduce the cost function for random excitation whether controlling the entire beam or just a subregion, without ever increasing the vibration outside the region in which control is desired. The last type of excitation considered is a single, transient pulse. A cost function representative of the beam vibration is obtained by integrating the transient displacement squared over a region of the beam and over all time. The form of the controller is chosen a priori as either one or two delayed pulses. Delays constrain the controller to be causal. The best possible control is then examined while varying the region of control and the controller location. It is found that control is always possible using either one or two control pulses. The two pulse controller gives better performance than a single pulse controller, but finding the optimal delay time for the additional controllers increases as the square of the number of control pulses.

Post, John Theodore↗

Computation of turbine flowfields with a Navier-Stokes code

A new technique has been developed for the solution of the incompressible Navier-Stokes equations. The numerical technique, derived from a pressure substitution method (PSM), overcomes many of the deficiencies of the pressure crrection method. This technique allows for the direct solution of the actual pressure in the form of a Poisson equation which is derived from the pressure weighted substitution of the full momentum equations into the continuity equation. In two-dimensions a turbine flowfield, including heat transfer, has been computed with this method and the prediction of the cascade performance is presented. The extension of the pressure correction method for the solution of three-dimensional flows is also presented for laminar flow in an S-shaped duct and turbulent flow in the end-wall region of a turbine cascade.

Hobson, G. V.↗

A pressure based method for the solution of viscous incompressible turbomachinery flows

A new technique was developed for the solution of the incompressible Navier-Stokes equations. The numerical technique, derived from a pressure substitution method (PSM), overcomes many of the deficiencies of the pressure correction method. This technique allows for the direct solution of the actual pressure in the form of a Poisson equation which is derived from the pressure weighted substitution of the full momentum equations into the continuity equation. Two dimensional internal flows are computed with this method. The prediction of cascade performance is presented. The extention of the pressure correction method for the solution of three dimensional flows is also presented.

Hobson, Garth Victor↗

The Use of a Block Diagram Simulation Language for Rapid Model Prototyping

The research performed this summer was a continuation of work performed during the 1995 NASA/ASEE Summer Fellowship. The focus of the work was to expand previously generated predictive models for liquid oxygen (LOX) loading into the external fuel tank of the shuttle. The models which were developed using a block diagram simulation language known as VisSim, were evaluated on numerous shuttle flights and found to well in most cases. Once the models were refined and validated, the predictive methods were integrated into the existing Rockwell software propulsion advisory tool (PAT). Although time was not sufficient to completely integrate the models developed into PAT, the ability to predict flows and pressures in the orbiter section and graphically display the results was accomplished.

Whitlow, Johnathan E.↗

NASA's Current Earth Science Program

NASA's Earth science program is a scientific endeavor whose goal is to provide long-term understanding of the Earth as an integrated system of land, water, air and life. A highly developed scientific knowledge of the Earth system is necessary to understand how the environment affects humanity, and how humanity may be affecting the environment. The remote sensing technologies used to gather the global environmental data used in such research also have numerous practical applications. Current applications of remote sensing data demonstrate their practical benefits in areas such as the monitoring of crop conditions and yields, natural disasters and forest fires; hazardous waste clean up; and tracking of vector-borne diseases. The long-term availability of environmental data is essential for the continuity of important research and applications efforts. NASA's Earth observation program has undergone many changes in the recent past.

Charles, Leslie Bermann↗

Calculation of the Full Scattering Amplitude without Partial Wave Decomposition: Inclusion of Exchange - 2

The development of a practical method of accurately calculating the full scattering amplitude, without making a partial wave decomposition is continued. The method is developed in the context of electron-hydrogen scattering, and here exchange is dealt with by considering e-H scattering in the static exchange approximation. The Schroedinger equation in this approximation can be simplified to a set of coupled integro-differential equations. The equations are solved numerically for the full scattering wave function. The scattering amplitude can most accurately be calculated from an integral expression for the amplitude; that integral can be formally simplified, and then evaluated using the numerically determined wave function. The results are essentially identical to converged partial wave results.

Shertzer, Janine↗

Comparisons of Crosswind Velocity Profile Estimates Used in Fast-Time Wake Vortex Prediction Models

Five methods for estimating crosswind profiles used in fast-time wake vortex prediction models are compared in this study. Previous investigations have shown that temporal and spatial variations in the crosswind vertical profile have a large impact on the transport and time evolution of the trailing vortex pair. The most important crosswind parameters are the magnitude of the crosswind and the gradient in the crosswind shear. It is known that pulsed and continuous wave lidar measurements can provide good estimates of the wind profile in the vicinity of airports. In this study comparisons are made between estimates of the crosswind profiles from a priori information on the trajectory of the vortex pair as well as crosswind profiles derived from different sensors and a regional numerical weather prediction model.

Pruis, Mathew J.↗

Collaborative Platforms Aid Emergency Decision Making

Terra. Aqua. Cloudsat. Landsat. NASA runs and partners in many missions dedicated to monitoring the Earth, and the tools used in these missions continuously return data on everything from shifts in temperature to cloud formation to pollution levels over highways. The data are of great scientific value, but they also provide information that can play a critical role in decision making during times of crisis. Real-time developments in weather, wind, ocean currents, and numerous other conditions can have a significant impact on the way disasters, both natural and human-caused, unfold. "NASA has long recognized the need to make its data from real-time sources compatible and accessible for the purposes of decision making," says Michael Goodman, who was Disasters Program manager at NASA Headquarters from 2009-2012. "There are practical applications of NASA Earth science data, and we d like to accelerate the use of those applications." One of the main obstacles standing in the way of eminently practical data is the fact that the data from different missions are collected, formatted, and stored in different ways. Combining data sets in a way that makes them useful for decision makers has proven to be a difficult task. And while the need for a collaborative platform is widely recognized, very few have successfully made it work. Dave Jones, founder and CEO of StormCenter Communications Inc., which consults with decision makers to prepare for emergencies, says that "when I talk to public authorities, they say, If I had a nickel for every time someone told me they had a common operating platform, I d be rich. But one thing we ve seen over the years is that no one has been able to give end users the ability to ingest NASA data sets and merge them with their own."

Source record↗

Compatibility Between Exploration EVA System and Exploration Spacecraft

Over the life of the Extravehicular Mobility Unit (EMU), numerous products detailing “how to build Extravehicular Activity (EVA) System hardware”, “how to interface with EVA System hardware” and “how to design hardware EVA will access and utilize” were generated to provide interoperability between a suited crewmember and the specified vehicle’s EVA task. Since the inception of these products, some have continued to receive updates due to the necessity of an on-going program while others remained unchanged for years and have led to discrepancies between the current accepted values and those considered outdated. For EVA-suited crewmember tasks beyond Low Earth Orbit (LEO), new vehicles need a single consolidated location for the best practices and lessons learned from the EVA Community. This paper outlines what common EVA compatibility design requirements are expected of an Exploration spacecraft that has interactions between the vehicle and an EVA suited crewmember for and an approach for standardizing EVA compatibility across various vehicles at various destinations. The approach of standardization allows for flexibility by tailoring the applicability to meet the EVA tasks required for that vehicle’s operation beyond Low Earth Orbit. This paper will also describe the broad difference between microgravity and partial gravity EVA compatibility and how those requirements were identified and will be informed.

Compatibility↗

Numerical simulation of interacting, three-dimensional vortex filaments

Unsteady three-dimensional fluid flows which are characterized by low viscosity and the presence of distinct regions of high vorticity embedded in an otherwise irrotational flow are considered, taking into account the case in which an intermingling of vortex filaments appears. Incompressible flows of the considered type are simulated with the aid of an approach in which the vorticity distribution is modeled in terms of continuous closed filaments. It is attempted to track these filaments in a Lagrangian reference frame.

Leonard, A.↗

Theoretical properties of Omega-loops in the convective zone of the Sun. 3: Extended updrafts

It was pointed out in an earlier paper that the continuing emergence of Omega-loops at localized sites on the surface of the Sun indicates a continuing updraft at those sites. The updraft evidently extends all the way from the base of the convective zone to somewhat near (approximately 10(exp 9)cm) the surface. We pointed out that such updrafts enhance the convective heat transport to the surface, accounting for a major part of the increased solar brightness or irradiance during times of solar activity. The problem is to work out, as much as possible, the dynamical nature of the extended updrafts, initiated as the wakes of successive rising Omega-loops and driven therafter by the convective forces. The question is, does the updraft take on a long-lived columnar form of its own, or does it never devlop beyond a sequence of rising wakes, resembling beads on a string? The dynamics of a columnar updraft is complicated by both the large Reynolds number and the strong stratication of the atmosphere, and by a total lack of direct observational information. Extended slender updrafts are not a spontaneous occurrence in numerical simulations of thermal convection in a stratified atmosphere, although slender concentrated downdrafts commonly occur. This paper examines several aspects of a columnar updrft in a convective atmosphere under various idealized circumstances to investigate to what extent that state can be maintained against the diminishing vorticity and expansion in the updraft. It appears that the successive passage of Omega-loops from the bottom to the top of the convective zone is an essential feature of the continuing existence of the updraft.

Parker, E. N.↗

Growth of Jupiter: Formation in disks of gas and solids and evolution to the present epoch

The formation of Jupiter is modeled via core-nucleated accretion, and the planet's evolution is simulated up to the present epoch. Throughout the phases when the planet acquires most of its heavy-element content, the calculation of solids' accretion accounts for interactions with an evolving disk of planetesimals. The phase of growth from an embryo of a few hundred kilometers in radius until the time when the accretion of gas overtakes solids' accretion was presented by D'Angelo et al., and the same numerical methods are applied here. Those calculations followed the formation for about 4 × 10 5 years, until the epoch when the heavy-element and hydrogen/helium masses were $M_Z ≈ 7.3$ and $M_{XY} ≈ 0.15$ Earth's masses ($M_⊕$), respectively, and $\dot {M} _{XY} ≈ $$\dot {M}_Z$. In this work, the calculation is continued through the phase when $M_{XY}$ grows to equal $M_Z$ , at which age, about 2.4 × 10 6 years, the total mass of the planet is $M_p ≈ 20 M_⊕$. About 9 × 10 5 years later, $M_p$ is approximately $60 M_⊕$ and $M_Z ≈ 16 M_⊕$, three-quarters of which are delivered by planetesimals larger than 10 km in radius. Around this epoch, the contraction of the envelope dictates gas accretion rates a few times 10 –3 $M_⊕$ per year, initiating the regime of disk-limited accretion, whereby the planet can accrete all the gas provided by the disk, and its evolution is therefore tied to disk's evolution. Growth is continued by constructing simplified models of protosolar accretion disks that evolve through viscous diffusion, winds, and accretion on the planet. Jupiter's formation ends after ≈ 3.4–4.2 Myr, depending on the applied disk viscosity parameter, when nebula gas disperses. The young Jupiter is 4.5–5.5 times as voluminous as it is presently and thousands of times as luminous, ~10 –5 $L_⊙$. The heavy-element mass is ≈ 20 $M_⊕$. The evolution proceeds through the cooling and contraction phase, in isolation except for solar irradiation. After 4570 Myr, the age of the solar system, radius and luminosity of the planet are within 10$\%$ of current values, accounting also for uncertainties in the power absorbed from the Sun. During formation, and soon thereafter, the planet exhibits features, e.g., luminosity and effective temperature, that may probe aspects of the latter stages of formation, if observable. These possibly distinctive features, however, seem to disappear within a few tens of Myr.

79 ASTRONOMY AND ASTROPHYSICS↗

Regional seismic velocity changes following the 2019 M w 7.1 Ridgecrest, California earthquake from autocorrelations and P / S converted waves

SUMMARY We examine regional transient changes of seismic velocities generated by the 2019 Mw 7.1 Ridgecrest earthquake in California, using autocorrelations of moving time windows in continuous waveforms recorded at regional stations. We focus on traveltime differences in a prominent phase generated by an interface around 2 km depth, associated with transmitted Pp waves and converted Ps waves from the ongoing microseismicity. Synthetic tests demonstrate the feasibility of the method for monitoring seismic velocity changes. Taking advantage of the numerous aftershocks in the early period following the main shock, we obtain a temporal resolution of velocity changes up to 20 min in the early post-main-shock period. The results reveal regional coseismic velocity drops in the top 1–3 km with an average value of ∼2 per cent over distances up to 100 km from the Ridgecrest event. These average velocity drops are likely dominated by larger changes in the shallow materials and are followed by rapid recoveries on timescales of days. Around the north end of the Ridgecrest rupture and the nearby Coso geothermal region, the observed coseismic velocity drops are up to ∼8 per cent. The method allows monitoring temporal changes of seismic velocities with high temporal resolution, fast computation and precise spatial mapping of changes. The results suggest that significant temporal changes of seismic velocities of shallow materials are commonly generated on a regional scale by large events.

58 GEOSCIENCES↗

Hybrid quantum simulations with qubits and qumodes on trapped-ion platforms

We explore the feasibility of gate-based hybrid quantum computing using both discrete (qubit) and continuous (qumode) variables on trapped-ion platforms. Trapped-ion systems have demonstrated record one- and two-qubit gate fidelities and long qubit coherence times, while qumodes, which can be represented by the collective vibrational modes of the ion chain, have remained relatively unex- plored for their use in computing. Using numerical simulations, we show that high-fidelity hybrid gates and measurement operations can be achieved for existing trapped-ion quantum platforms. As an exemplary application, we consider quantum simulations of the Jaynes-Cummings-Hubbard model, which is given by a one-dimensional chain of interacting spin and boson degrees of free- dom. Using classical simulations, we study its real-time evolution and develop a suitable variational quantum algorithm for ground state preparation. Furthermore, our results motivate further studies of hybrid quantum computing in this context, which may lead to direct applications in condensed matter and fundamental particle and nuclear physics.

Lower-dimensional field theories↗

Nonperturbative phase diagram of two-dimensional N = ( 2 , 2 ) super-Yang-Mills theory

We consider two-dimensional N = ( 2 , 2 ) Yang-Mills theory with gauge group SU ( N ) in Euclidean signature compactified on a torus with thermal fermion boundary conditions imposed on one cycle. We perform nonperturbative lattice analyses of this theory for large 12 ≤ N ≤ 20 . Although no holographic dual of this theory is yet known, we conduct numerical investigations to check for features similar to the two-dimensional N = ( 8 , 8 ) Yang–Mills theory, which has a well-defined gravity dual. We perform lattice field theory calculations to determine the phase diagram, observing a spatial deconfinement transition similar to the maximally supersymmetric case. However, the transition does not continue to low temperature, implying the absence of a topology-changing transition between black hole geometries in any holographic dual for this four-supercharge theory. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Dynamic Earth Energy Storage: Terawatt-year, Grid-scale Energy Storage Using Planet Earth as a Thermal Battery (GeoTES): Phase I Project (Final Report)

Grid-scale energy storage has been identified by the U.S. Department of Energy’s (DOE) Energy Storage Grand Challenge as a necessary technology to support the continued build-out of intermittent renewable energy resources required to attain a carbon-free energy future. To meet this goal, the 2018 Department of Energy Research and Innovation Act mandated the creation of a comprehensive program to accelerate the development and commercialization of next-generation energy storage technologies. One of numerous energy storage technology options is the storage of excess energy as heated geothermal brine in suitable geologic formations. This concept, known as reservoir thermal energy storage (RTES), geologic thermal energy storage (GeoTES), aquifer thermal energy storage (ATES), etc., relies on the storage of thermal energy in geologic formations for recovery and use in large-scale direct use geothermal (e.g., district heating, industrial processes, etc.) and electrical power generation applications. This thermal energy is derived from excess or waste heat from any high-temperature heat source, such as concentrated solar or from conventional thermal/nuclear generation. As such, RTES can potentially play a significant role in meeting the energy storage shortfall in the coming decades. RTES can provide energy arbitrage through both the storage and production of thermal energy stored in geologic formations for direct use applications and can serve as a source of hot fluids that can be used to generate electricity to support peak demand ramping, thus easing stress on transmission and distribution. This energy storage option has geographic benefits in that energy can be stored locally or regionally depending on the various needs/loads. RTES can also be located across an enormous geographic area, without the need for a traditional hydrothermal resource but where thermal gradients and hydrogeology allow economic exploitation of subsurface heat. The work conducted for this project includes (1) a review of lessons learned from past high-temperature RTES international projects; (2) geochemical experimental investigation and numerical simulations of potential domestic sedimentary reservoirs and (3) development of a thermo-hydrological-mechanical (THM) numerical simulation tool for optimizing formation properties and design parameters to maximize thermal energy storage performance.

15 GEOTHERMAL ENERGY↗

Finding ways to reduce nuclear waste: searching for the unknown one step at a time

In my home country, Venezuela, research has been stagnant. Due to the political turmoil and the crisis, many educated people have left the country in search for a better life. This has caused a deficit in any technological and scientific advances, making Venezuela one of the first South American countries to have its rate of publications decline by 29% in 2013. Currently, Venezuela lacks the infrastructure and the means to keep up with the research progress as compared to other countries in South America, such as Brazil. Since coming to the United States (US), and currently working for a national laboratory, the active research environment endorses a wide range of careers and engineering programs that allow researchers to thrive at any given field. Researchers have access to funds and tools to succeed in developing materials for the future. There are 17 national laboratories in the US, and all of these have a different research focus/objective. As examples, Los Alamos National Laboratory and Sandia National Laboratory focus is on national homeland security, weapon science, radiation effects, among others. Argonne National Laboratory focuses on nuclear energy, energy storage, high performance computing, etc. At Idaho National Laboratory (INL) the research focuses on innovating nuclear energy and clean energy resources, critical infrastructure materials, along with fuel cycle solutions to manage, dispose and find ways to recycle current and future radiological waste. Compared to other national laboratories, INL focuses slightly more on applied processes and how nuclear energy can be innovated to next reactor design and technologies. The research being conducted at INL made me apply for a Seaborg distinguished postdoctoral position. For the position itself, the researcher must submit a proposal related to actinide chemistry on a research field area. In this position, 50% of my time will be focused on my own proposal. The proposal that I am working on is focused on the innovation of nuclear energy and fuel cycle recycling, which is why I was mainly interested on working at this national laboratory. To give a bit more context of what my proposal is about, a little bit of background is necessary: After the nuclear fuel (UO2) is used in a reactor, the fuel matrix is then characterized by various fission products (FP). Among these FP (including rare earth elements, alkali/alkaline earths, and actinides), many can potentially be recovered through nuclear reprocessing technologies. In pyroprocessing, the used nuclear fuel undergoes electrochemical dissolution into a molten chloride salt mixture in an electrorefiner. Initially, uranium is reduced onto an inert cathode by applied potentials. However, numerous remaining FPs accumulate in the melt and pose challenges for recovery by an inert electrode, particularly the rare earth elements (e.g., Nd, Gd, Pr, Sm) due to their multivalent oxidation states and tendencies toward side reactions, leading to their dissolution in the electrolyte. These recovery challenges result in inefficiencies and necessitate the continual discarding of the molten chloride salt, thereby generating additional waste. Furthermore, the presence of rare earth elements and other fission products in the molten salt electrolyte alters its physical and chemical properties, affecting both uranium recovery efficiency and the longevity of the molten chloride salt. To improve the recovery efficiency of the FP, specifically rare earth elements, I am investigating the fundamental interactions between rare earth elements in the molten chloride salt and their metallic form. The kinetic pathways and the chemical reactions of these elements will give insights on how the recovery efficiency can be improved. The interactions and speciation of these elements are being studied by spectro-electrochemistry at high temperature environments in quartz and other ceramic materials (e.g., alumina crucibles). Some of the challenges I am facing specifically relates the reactivity of some of these elements with different glass and crucible materials. Although my research focuses on fundamental science, it will benefit the applied process by generating new scientific knowledge and closing the gap for an efficient recycling of the waste: one step at a time.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗