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 415 records · Page 23
All sourcesNASA NTRSDOE OSTI
NASA NTRS2007

Response of the Earth's Magnetosphere to Changes in the Solar Wind

The solar wind couples to the magnetosphere via dynamic pressure and electric field. Pressure establishes the size and shape of the system, while the electric field transfers energy, mass, and momentum to the magnetosphere. When the interplanetary magnetic field (IMF) is antiparallel to the dayside magnetic field, magnetic reconnection connects the IMF to the dipole field. Solar wind transport of the newly opened field lines to the nightside creates an internal convection system. These open field lines must ultimately be closed by reconnection on the nightside. For many decades, it was thought that a magnetospheric substorm was the process for accomplishing this and that all magnetic activity was a consequence of substorms. It is now recognized that there are a variety of modes of response of the magnetosphere to the solar wind. In this paper, we briefly describe these modes and the conditions under which they occur. They include substorms, pseudo-breakups, poleward boundary intensifications (PBI), steady magnetospheric convection (SMC), sawtooth injection events, magnetic storms, high-intensity long-duration continuous AE activities (HILDCAAs), and storm-time activations. There are numerous explanations for these different phenomena, some of which do not involve magnetic reconnection. However, we speculate that it is possible to interpret each mode in terms of differences in the way magnetic reconnection occurs on the nightside.

McPherron, Robert L.↗
NASA NTRS2011

Shuttle Propulsion Overview - The Design Challenges

The major elements of the Space Shuttle Main Propulsion System include two reusable solid rocket motors integrated into recoverable solid rocket boosters, an expendable external fuel and oxidizer tank, and three reusable Space Shuttle Main Engines. Both the solid rocket motors and space shuttle main engines ignite prior to liftoff, with the solid rocket boosters separating about two minutes into flight. The external tank separates, about eight and a half minutes into the flight, after main engine shutdown and is safely expended in the ocean. The SSME's, integrated into the Space Shuttle Orbiter aft structure, are reused after post landing inspections. The configuration is called a stage and a half as all the propulsion elements are active during the boost phase, with only the SSME s continuing operation to achieve orbital velocity. Design and performance challenges were numerous, beginning with development work in the 1970's. The solid rocket motors were large, and this technology had never been used for human space flight. The SSME s were both reusable and very high performance staged combustion cycle engines, also unique to the Space Shuttle. The multi body side mount configuration was unique and posed numerous integration and interface challenges across the elements. Operation of the system was complex and time consuming. This paper describes the design challenges and key areas where the design evolved during the program.

Owen, James W.↗
DOE OSTI2025

Reducing NOx Emissions in Ammonia Combustors

Ammonia continues to attract growing interest as a carbonneutral replacement fuel, motivating numerous research efforts toward understanding fundamental ammonia combustion characteristics. A major challenge for the use of ammonia is the development of combustor technologies for mitigating potentially high NOx emissions from the fuel-bound nitrogen chemical pathways to acceptable levels. Our work focuses on a staged RQL combustor architecture for minimizing the NOx emission levels through burning fuel-rich in the primary stage to formcombustion products containing significant levels of hydrogen in addition to nitrogen and water with minimal NOx formation. The subsequent quench and burnout stages of the combustor must then quickly burn residual hydrogen with flame-temperatures moderated by nitrogen and water forming in the first stage. Chemical Reactor Network modeling was used to understand and identify optimal stoichiometry and residence times in each stage for minimizing NOx emissions and to quantify pressure and temperature effects. Reducing the overall NOx emissions requires relatively long residence times in the primary stage to achieve near equilibrium NO levels due to kinetically controlling processes. For conditions relevant to gas turbines (e.g., 30 atm), our work indicates that NOx emissions below 20 ppm are theoretically achievable in a staged RQL combustor architecture. However, these emission predictions significantly depend on the accuracies of currently available chemical kinetic mechanisms which have not been extensively validated under elevated pressure and temperature conditions relevant to gas turbines.

10 SYNTHETIC FUELS↗
DOE OSTI2022

Removing numerical pathologies in a turbulence parameterization through convergence testing

Discretized numerical models of the atmosphere are usually intended to faithfully represent an underlying set of continuous equations, but this necessary condition is violated sometimes by subtle pathologies that have crept into the discretized equations. Such pathologies can introduce undesirable artifacts, such as sawtooth noise, into the model solutions. The presence of these pathologies can be detected by numerical convergence testing. This study employs convergence testing to verify the discretization of the Cloud Layers Unified By Binormals (CLUBB) model of clouds and turbulence. That convergence testing identifies two aspects of CLUBB's equation set that contribute to undesirable noise in the solutions. First, numerical limiters (i.e. clipping) used by CLUBB introduce discontinuities or slope discontinuities in model fields. Second, this noise can be amplified by an advective term in CLUBB's background diffusion. Smoothing the limiters and removing the advective component of the background diffusion reduces the noise and restores the expected first-order convergence in CLUBB's solutions. These model reformulations improve the results at coarser, near-operational grid spacing and time step in cumulus cloud and dry turbulence tests. In addition, convergence testing is proved to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

Mathematics and Computing, Geosciences↗
DOE OSTI2023

Removing Numerical Pathologies in a Turbulence Parameterization Through Convergence Testing

Abstract Discretized numerical models of the atmosphere are usually intended to faithfully represent an underlying set of continuous equations, but this necessary condition is violated sometimes by subtle pathologies that have crept into the discretized equations. Such pathologies can introduce undesirable artifacts, such as sawtooth noise, into the model solutions. The presence of these pathologies can be detected by numerical convergence testing. This study employs convergence testing to verify the discretization of the Cloud Layers Unified By Binormals (CLUBB) model of clouds and turbulence. That convergence testing identifies two aspects of CLUBB's equation set that contribute to undesirable noise in the solutions. First, numerical limiters (i.e., clipping) used by CLUBB introduce discontinuities or slope discontinuities in model fields. Second, nonlinear artificial diffusion employed for improving numerical stability can introduce unintended small‐scale features into the solution of the model equations. Smoothing the limiters and using linear artificial diffusion reduces the noise and restores the expected first‐order convergence in CLUBB's solutions. These model reformulations enhance our confidence in the trustworthiness of solutions from CLUBB by eliminating the unphysical oscillations in high‐resolution simulations. The improvements in the results at coarser, near‐operational grid spacing and timestep are also seen in cumulus cloud and dry turbulence tests. In addition, convergence testing is proven to be a valuable tool for detecting pathologies, including unintended discontinuities and grid dependence, in the model equation set.

58 GEOSCIENCES↗
NASA NTRS1996

Numerical Simulation of High-Speed Turbulent Reacting Flows

The purpose of this research is to continue our efforts in advancing the state of knowledge in large eddy simulation (LES) methods for computational analysis of high-speed reacting turbulent flows. We have just completed the first year of Phase 3 of this research.

Jaberi, F. A.↗
DOE OSTI2021

Oxygenated Fuel Blending Effects in a Mixing-Controlled Compression-Ignition Engine Equipped with Ducted Fuel Injection [Slides]

Continued creation of harmful emissions such as NO x and soot from compression-ignition engines utilizing mixing-controlled combustion systems (i.e., diesel engines) remains a problem and is the subject of on-going research. The inherently high efficiency, relatively low cost, and numerous other desirable attributes of such engines, coupled with a widely supported infrastructure, motivates their continued advancement. Recently, a scientifically distinct and mechanically simple technology called ducted fuel injection (DFI) has shown a robust ability to allow such engines to operate with simultaneously low engine-out soot and NO x emissions when it is employed with simulated exhaust-gas recirculation. To better understand the property ranges of sustainable, oxygenated-fuel blending stocks that will most improve engine performance, two oxygenated blendstocks were separately blended with a commercial diesel base fuel and tested within a heavy-duty diesel optical engine equipped with a four-duct DFI configuration. Conventional and crank-angle-resolved optical diagnostics were used to elucidate the effects of fuel ignition quality, oxygenate molecular structure, and overall oxygen content on engine performance.

33 ADVANCED PROPULSION SYSTEMS↗
NASA NTRS1978

Turbulence models from the point of view of an industrial user

Since there appears to be little prospect of the development of methods for the solution of the full, time dependent Navier-Stokes equations in the near future, the need continues for turbulence models to approximate the Reynolds stress terms that appear in the time averaged Navier-Stokes equations. Even if methods were available for solving the full equations, this would not necessarily be the optimum choice in all cases. As the cost of numerical computations decreases, the trend toward the use of more complex methods is likely to continue, but there will always be a need for a range of methods, depending on the accuracy and detail required from the calculation.

Birch, S. F.↗
NASA NTRS1981

The potential impact of scatterometry on oceanography - A wave forecasting case

A series of observing system simulation experiments have been performed in order to assess the potential impact of marine surface wind data on numerical weather prediction. In addition to conventional data, the experiments simulated the time-continuous assimilation of remotely sensed marine surface wind or temperature sounding data. The wind data were fabricated directly for model grid points intercepted by a Seasat-1 scatterometer swath and were assimilated into the lowest active level (945 mb) of the model using a localized successive correction method. It is shown that Seasat wind data can greatly improve numerical weather forecasts due to better definition of specific features. The case of the QE II storm is examined.

Cane, M. A.↗
NASA NTRS1993

Calculation methods for rotor wake formation in forward flight

The primary focus of this research centered on improving the vorticity-capturing properties of Eulerian flow solvers in convection-dominated flows. The work was motivated by the fact that excessive numerical diffusion observed in current Euler solvers on grids of feasible size is a major obstacle to accurately computing the flow about helicopter rotors, especially in the forward flight case. The work was originally supported by NASA Grant NAG2-421 for two years. The research then continued for four more years under an ONR contract. The research goal of reduction of numerical diffusion in Euler solvers remained the same, although the target application was now the prediction of internal vortical secondary flows. This was part of a larger ONR-sponsored research program at MIT aimed at reducing the noise of turbopumps. The research work culminated in the PhD thesis titled 'A Coupled Eulerian/Lagrangian Method for the Solution of Three-Dimensional Vortical Flows'.

Drela, Mark↗
DOE OSTI2024

Anti-symmetric and positivity preserving formulation of a spectral method for Vlasov-Poisson equations

We analyze the anti-symmetric properties of a spectral discretization for the one-dimensional Vlasov-Poisson equations. The discretization is based on a spectral expansion in velocity with the symmetrically weighted Hermite basis functions, central finite differencing in space, and an implicit Runge Kutta integrator in time. The proposed discretization preserves the anti-symmetric structure of the advection operator in the Vlasov equation, resulting in a stable numerical method. We apply such discretization to two formulations: the canonical Vlasov-Poisson equations and their continuously transformed square-root representation. The latter preserves the positivity of the particle distribution function. We derive analytically the conservation properties of both formulations, including particle number, momentum, and energy, which are verified numerically on the following benchmark problems: manufactured solution, linear and nonlinear Landau damping, two-stream instability, bump-on-tail instability, and ion-acoustic wave.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

EERC Final Topical Report – Findings on Subtask 3.2 – Produced Water Management Through Geologic Homogenization, Conditions, and Reuse

The Energy & Environmental Research Center (EERC) was awarded a contract by the North Dakota Industrial Commission (NDIC) Oil and Gas Research Program (NDIC No. G-051-101) to conduct a study on the recycling of water used in oil and gas operations, also known as produced water, from oil- and gas-producing regions of North Dakota as directed by Section 19 of North Dakota House Bill 1014. This final report provides a compilation of results of the study, which include regulatory, scientific, technological, and feasibility methods and considerations associated with North Dakota produced water management. The report provides a synopsis of this project’s previously submitted produced water assessment report entitled “Produced Water Management and Recycling Options in North Dakota” (Energy & Environmental Research Center, 2020), with updated values provided as appropriate. The report provides the results from the investigation of a novel produced water management strategy, referred to as geologic homogenization, conditioning, and reuse (GHCR), which aims to use a subsurface geologic formation as a natural medium for managing produced water recycling and reuse. Water management is a significant technical and economic challenge for sustainable oil and gas production, and water volumes are intrinsically linked to oil production volumes. North Dakota oil production rose to over 1.5 million barrels (MMbbl)/day in 2019, and despite a downturn in oil price in early 2020, North Dakota oil production has recovered to 1.1 MMbbl/day as of August 2021. Bakken petroleum system development between 2008 and 2020 has resulted in a nearly fourfold increase in produced water volumes to 642 MMbbl/yr in 2020 after peaking at 740 MMbbl/yr in 2019 and a fivefold increase in saltwater disposal (SWD) volumes to 565 MMbbl/yr in 2020 after a peak of 682 MMbbl/yr in 2019. Produced water and SWD volumes are forecasted to double by 2030. SWD is the primary method of produced water management used in North Dakota, with approximately 95% of the SWD volume occurring through subsurface injection into sandstones of the Dakota Group (Dakota). Localized pressurization of the Dakota resulting from SWD and projected increases in produced water volumes could impact the economics of North Dakota oil production. As a result, there is an emerging need to pursue alternative produced water management approaches, including recycling and reuse. While produced water recycling is not yet widespread, commercial operators are making strides in overcoming the technical challenges of using high salinity produced water in completion operations (Marathon Oil, 2020). As water management continues to be a key focal point in companies’ environmental, social, and governance (ESG) initiatives, focus on water management, including recycling, will likely continue to increase. Laboratory column testing, field sample collection, geologic modeling and numerical simulation, and techno-economic analysis all indicate that GHCR could feasibly be implemented as a potential water management option. Laboratory column testing and field sample collection indicate that the Inyan Kara sandstone and native formation fluid are capable of homogenizing with the Bakken produced water to a point where the fluid composition appears to stabilize. Extracting that stabilized fluid could be considered homogeneous and capable of providing individual batches of hydraulic fracturing fluid. Numerical simulation results indicate that extraction of fluids from the Inyan Kara in a GHCR implementation scenario is capable of reducing formation pressure, which would help ease localized pressurization of the Inyan Kara and extend the available capacity for nearby existing SWD wells. Economic analysis indicates that there are scenarios where GHCR implementation can be a competitive or even lower-cost option than a conventional water management approach. Site-specific conditions will dictate the economic potential of GHCR, but potentially attractive sites for GHCR implementation will be those that are located above a pressurized zone of the Inyan Kara, need six or more Bakken infill wells, and face high costs for conventional SWD and/or freshwater. Based on the regulatory review, drilling into the Inyan Kara for SWD and to harness as a source water for industrial use have precedent in North Dakota, and a workable regulatory solution for GHCR seems likely. However, restrictions in the state regarding surface storage and transport of produced fluids may limit some activities, which will affect how GHCR could ultimately be implemented. In summary, this study reveals pursuing GHCR can be a viable approach to water management in North Dakota. The GHCR concept addresses some of the challenges that hinder the more traditional approaches to recycling in the industry. Furthermore, an assessment of the current landscape of water management within the state reveals the ongoing trend of increasing volumes of produced water and SWD. Projections reveal that the volumes of produced water that need to be managed are expected to double over the next decade (Energy & Environmental Research Center, 2020). With the continued development of the Bakken and continuing driving factors related to ESG initiatives, implementing a practice such as GHCR is a feasible approach to adding recycling of produced water to industry within the state. This subtask was cofunded through the EERC–U.S. Department of Energy Joint Program on Research and Development for Fossil Energy-Related Resources Cooperative Agreement No. DE-FE0024233. Nonfederal funding was provided by the North Dakota Industrial Commission Oil and Gas Research Program. References Energy & Environmental Research Center, 2020, Produced water management and recycling options in North Dakota: Final Report for North Dakota Legislative Management Energy Development and Transmission Committee and North Dakota Industrial Commission. Marathon Oil, 2020, Sustainability report: https://cdn.sanity.io/files/ghcnw9z2/website/ 91744eb6ef8fbe59505a911c6b8d2e8dd9a537fa.pdf?dl (accessed November 2021).

02 PETROLEUM↗
NASA NTRS2001

Nonlinear Thermomechanical Response of Composite Panels with Continuous and Terminated Stiffeners

A two-phase approach and a computational procedure are used for predicting the variability of the response of stiffened composite panels associated with variations in the geometric and material parameters of the structures. In the first phase, hierarchical sensitivity analysis is used to identify the major parameters that have the most effect on the response quantities of interest. In the second phase, the major parameters are taken to be fuzzy parameters, and a fuzzy set analysis is used to determine the range of variation of the response, associated with preselected variations in the major parameters. Numerical results are presented showing the variability of the response of panels with both continuous and terminated stiffeners associated with variations in the micro mechanical and geometric parameters. Both flat and curved panels are considered.

Noor, Ahmed K.↗
DOE OSTI2024

Light-ray wave functions and integrability

Using integrability, we construct (to leading order in perturbation theory) the explicit form of twist-three light-ray operators in planar $\mathcal{N}$ = 4 SYM. This construction allows us to directly compute analytically continued CFT data at complex spin. We derive analytically the “magic” decoupling zeroes previously observed numerically. Using the Baxter equation, we also show that certain Regge trajectories merge together into a single unifying Riemann surface. Perhaps more surprisingly, we find that this unification of Regge trajectories is not unique. If we organize twist-three operators differently into what we call “cousin trajectories” we find infinitely more possible continuations. We speculate about which of these remarkable features of twist-three operators might generalize to other operators, other regimes and other theories.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1993

A global analytical representation of the magnetic field produced by the region 2 Birkeland currents and the partial ring current

A quantitative model is developed of the magnetic field produced by the electric current system of region 2 Birkeland currents, closed via the partial ring current. The distribution of j-perpendicular is computed from a given axially asymmetric spatial distribution of hot isotropic magnetospheric plasma over an infinitely thin L shell in an axisymmetric purely dipolar geomagnetic field, while the field-aligned current density is found from the continuity of the net electric current. The magnetic field distribution is derived by a Biot-Savart integral over the electric current system. An assumed cosine dependence of the plasma pressure on local time makes it possible to reduce the problem of analytical representation of the B field to two dimensions. The obtained numerical fits for the partial ring current/region 2 Birkeland current magnetic field are relatively simple, continuous, and valid throughout the whole extraterrestrial space from ionospheric heights up to tens of earth radii.

Tsyganenko, N. A.↗
NASA NTRS1995

An Experimental and Theoretical Study of Radiative Extinction of Diffusion Flames

The objective of this research was to experimentally and theoretically investigate the radiation-induced extinction of gaseous diffusion flames in microgravity. The microgravity conditions were required because radiation-induced extinction is generally not possible in 1-g but is highly likely in microgravity. In 1-g, the flame-generated particulates (e.g. soot) and gaseous combustion products that are responsible for flame radiation, are swept away from the high temperature reaction zone by the buoyancy-induced flow and a steady state is developed. In microgravity, however, the absence of buoyancy-induced flow which transports the fuel and the oxidizer to the combustion zone and removes the hot combustion products from it enhances the flame radiation due to: (1) transient build-up of the combustion products in the flame zone which increases the gas radiation, and (2) longer residence time makes conditions appropriate for substantial amounts of soot to form which is usually responsible for most of the radiative heat loss. Numerical calculations conducted during the course of this work show that even non-radiative flames continue to become "weaker" (diminished burning rate per unit flame area) due to reduced rates of convective and diffusive transport. Thus, it was anticipated that radiative heat loss may eventually extinguish the already "weak" microgravity diffusion flame. While this hypothesis appears convincing and our numerical calculations support it, experiments for a long enough microgravity time could not be conducted during the course of this research to provide an experimental proof. Space shuttle experiments on candle flames show that in an infinite ambient atmosphere, the hemispherical candle flame in microgravity will burn indefinitely. It was hoped that radiative extinction can be experimentally shown by the aerodynamically stabilized gaseous diffusion flames where the fuel supply rate was externally controlled. While substantial progress toward this goal was made during this project, identifying the experimental conditions for which radiative extinction occurs for various fuels requires further study. Details concerning this research which are discussed in published articles are included in the appendices.

Atreya, Arvind↗
DOE OSTI2026

OC6 Project Phase IV: Validation of Numerical Models for Novel Floating Offshore Wind Support Structures

The fourth phase of the Offshore Code Comparison, Collaboration, Continued, with Correlation and unCertainty (OC6) aims to assess the ability of the tools to model a novel floating wind design. The design that was selected for this project is the TetraSpar designed by Stiesdal Offshore Technologies. The platform is a unique design that features a spar like stability with a very low center of gravity, but does not require deep-water ports. This combination is achieved with a system made of two separate steel structures, a hull at the water’s surface that supports the wind turbine tower, and a keel that is suspended below the main body. This lower keel has a very high density and drives the total platform’s low center of gravity. Scaled model tests of the TetraSpar design were performed in 2018-2019 by the University of Maine. The tests were performed with 1:43 Froude scaling. The data collected in the model tests provide a detailed description of the response of the platform in various controlled and relatively well understood conditions. Comparison with these data provides a good opportunity for validation of numerical models. The OC6 Phase IV work will focus on the specific parameters used in these scaled model tests. The load cases (LC) considered in this project are as follows: LC 1.X - Equilibrium and Static Offset Tests LC 2.X - Platform Free Decays LC 3.X - Turbine Thrust (Wind Only) LC 4.X - Regular and Irregular Waves (Wave only) LC 5.X - Combined Wind and Wave Loading

17 WIND ENERGY↗
NASA NTRS1978

Large eddy simulation of incompressible turbulent channel flow

The three-dimensional, time-dependent primitive equations of motion were numerically integrated for the case of turbulent channel flow. A partially implicit numerical method was developed. An important feature of this scheme is that the equation of continuity is solved directly. The residual field motions were simulated through an eddy viscosity model, while the large-scale field was obtained directly from the solution of the governing equations. An important portion of the initial velocity field was obtained from the solution of the linearized Navier-Stokes equations. The pseudospectral method was used for numerical differentiation in the horizontal directions, and second-order finite-difference schemes were used in the direction normal to the walls. The large eddy simulation technique is capable of reproducing some of the important features of wall-bounded turbulent flows. The resolvable portions of the root-mean square wall pressure fluctuations, pressure velocity-gradient correlations, and velocity pressure-gradient correlations are documented.

Moin, P.↗