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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
NASA NTRS2002

Trajectories for High Specific Impulse High Specific Power Deep Space Exploration

Flight times and deliverable masses for electric and fusion propulsion systems are difficult to approximate. Numerical integration is required for these continuous thrust systems. Many scientists are not equipped with the tools and expertise to conduct interplanetary and interstellar trajectory analysis for their concepts. Several charts plotting the results of well-known trajectory simulation codes were developed and are contained in this paper. These charts illustrate the dependence of time of flight and payload ratio on jet power, initial mass, specific impulse and specific power. These charts are intended to be a tool by which people in the propulsion community can explore the possibilities of their propulsion system concepts. Trajectories were simulated using the tools VARITOP and IPOST. VARITOP is a well known trajectory optimization code that involves numerical integration based on calculus of variations. IPOST has several methods of trajectory simulation; the one used in this paper is Cowell's method for full integration of the equations of motion. An analytical method derived in the companion paper was also evaluated. The accuracy of this method is discussed in the paper.

Polsgrove, Tara↗
NASA NTRSDate not supplied

Evidence of High-Angle Ejecta from Impacts Iinto Various Granular Materials

The dominant mechanism for resurfacing airless bodies is the excavation and ejection of material during the formation of impact craters. Ejecta that is deposited proximal to the rim of craters formed in the gravity-regime (where gravitational forces are greater than material strength) form a relatively smooth and continuous deposit. Aside from far-field secondary craters, craters observed on continuous ejecta deposits are mostly assumed to be representative of the primary impactor flux after the time of deposition. Sustained primary bombardment should produce a homogeneous crater population across this surface, thus providing an opportunity to determine crater size-frequency distributions and surface chronologies. Deviations from this ideal have been observed as different crater-size distributions appearing between the continuous deposits and their superimposed melt ponds. Numerous explanations for this discrepancy have been proposed, one of which is the occurrence of self-secondary craters, which would be created by ejecta launched at sufficiently high-angles to re-impact the fresh continuous ejecta deposits. Here, we present results from cratering experiments to explore the influence of density, porosity, and friction on ejecta kinematics and report on the observation of such high-angle ejecta in cohesionless granular targets.

Impact crater↗
NASA NTRS1995

Annual Research Briefs: 1995

This report contains the 1995 annual progress reports of the Research Fellows and students of the Center for Turbulence Research (CTR). In 1995 CTR continued its concentration on the development and application of large-eddy simulation to complex flows, development of novel modeling concepts for engineering computations in the Reynolds averaged framework, and turbulent combustion. In large-eddy simulation, a number of numerical and experimental issues have surfaced which are being addressed. The first group of reports in this volume are on large-eddy simulation. A key finding in this area was the revelation of possibly significant numerical errors that may overwhelm the effects of the subgrid-scale model. We also commissioned a new experiment to support the LES validation studies. The remaining articles in this report are concerned with Reynolds averaged modeling, studies of turbulence physics and flow generated sound, combustion, and simulation techniques. Fundamental studies of turbulent combustion using direct numerical simulations which started at CTR will continue to be emphasized. These studies and their counterparts carried out during the summer programs have had a noticeable impact on combustion research world wide.

Source record↗
DOE OSTI2025

Artificial Intelligence and Machine Learning Applications in Modern Power Systems

Machine learning (ML) and artificial intelligence (AI) algorithms offer valuable tools for the analysis and interpretation of large datasets. These tools have the capability to uncover insights that may not be readily apparent within these datasets. In recent years, the integration of ML and AI has become increasingly prevalent in various applications within the power system domain. One of the earliest instances of machine learning in power systems can be traced back to demand forecasting, where artificial neural networks were employed for short-term load forecasting. In contemporary power systems, an abundance of high-resolution geospatial and temporal data is generated at various time intervals, ranging from sub-seconds (Phasor Measurement Units or PMUs) to seconds (Supervisory Control and Data Acquisition or SCADA), minutes (Process Information or PI), and extending to days, months, and years. These datasets contain valuable information concerning system reliability and performance. This information holds the potential to offer critical insights into system operations, as well as solutions for predicting and mitigating contingencies to prevent cascading outages. Despite the immense power of machine learning tools, system operators, planners, and utilities often exhibit hesitancy in fully embracing AI-enabled system operations and planning. This cautious approach persists, even as numerous diverse applications of machine learning continue to emerge in the realm of power systems. In this chapter, our focus will delve deep into ML and AI applications tailored for power systems. These applications aim to furnish system operators with enhanced situational awareness and augment their decision-making capabilities, especially during challenging operating conditions. Specific areas of interest encompass root cause analyses of electricity market datasets and the strategic selection of representative samples from vast power system databases for training ML/AI models. Finally, the chapter will conclude with a short discussion on the future of ML/AI in power systems and possible directions that the industry is moving towards.

power system applications, machine learning (ML), ↗
DOE OSTI2021

A duality-based coupling of Cosserat crystal plasticity and phase field theories for modeling grain refinement

High-rate deformation processes of metals entail intense grain refinement and special attention needs to be paid to capture the evolution of microstructure. In this article, a new formulation for coupling Cosserat crystal plasticity and phase field is developed. A common approach is to penalize kinematic incompatibility between lattice orientation and displacement-based elastic rotation. However, this can lead to significant solution sensitivity to the penalty parameter, resulting in low accuracy and convergence rates. To address these issues, a duality-based formulation is developed which directly imposes the rotational kinematic compatibility. A weak inf-sup-based skew-symmetric stress projection is introduced to suppress instabilities present in the dual formulation. An additional least squares stabilization is introduced to suppress the spurious lattice rotation with a suitable parameter range derived analytically and validated numerically. The required high-order continuity is attained by the reproducing kernel approximation. It is observed that equal order displacement-rotation-phase field approximations are stable, which allows efficient employment of the same set of shape functions for all independent variables. The proposed formulation is shown to yield superior accuracy and convergence with marginal parameter sensitivity compared to the penalty-based approach and successfully captures the dominant rotational recrystallization mechanism including block dislocation structures and grain boundary migration.

42 ENGINEERING↗
DOE OSTI2024

D2NO: Efficient handling of heterogeneous input function spaces with distributed deep neural operators

Neural operators have been applied in various scientific fields, such as solving parametric partial differential equations, dynamical systems with control, and inverse problems. However, challenges arise when dealing with input functions that exhibit heterogeneous properties, requiring multiple sensors to handle functions with minimal regularity. To address this issue, discretization-invariant neural operators have been used, allowing the sampling of diverse input functions with different sensor locations. However, existing frameworks still require an equal number of sensors for all functions. We propose a novel distributed approach to further relax the discretization requirements and solve the heterogeneous dataset challenges. Our method involves partitioning the input function space and processing individual input functions using independent and separate neural networks. A centralized neural network is used to handle shared information across all output functions. This distributed methodology reduces the number of gradient descent back-propagation steps, improving efficiency while maintaining accuracy. Here, we demonstrate that the corresponding neural network is a universal approximator of continuous nonlinear operators and present three numerical examples to validate its performance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Assessment of chemo-mechanical impacts of CO 2 sequestration on the caprock formation in Farnsworth oil field, Texas

This study evaluates the chemo-mechanical influence of injected CO 2 on the Morrow B sandstone reservoir and the upper Morrow shale caprock utilizing data from the inverted 5-spot pattern centered on Well 13-10A within the Farnsworth unit (FWU). This study also seeks to evaluate the integrity of the caprock and the long-term CO 2 storage capability of the FWU. The inverted 5-spot pattern was extracted from the field-scale model and tuned with the available field observed data before the modeling work. Two coupled numerical simulation models were utilized to continue the study. First, a coupled hydro-geochemical model was constructed to simulate the dissolution and precipitation of formation minerals by modeling three intra-aqueous and six mineral reactions. In addition, a coupled hydro-geomechanical model was constructed and employed to study the effects of stress changes on the caprock’s porosity, permeability, and ground displacement. The Mohr–Coulomb circle and failure envelope were used to determine caprock failure. In this work, the CO 2 -WAG injection is followed by the historical field-observed strategy. During the forecasting period, a Water Alternating Gas (WAG) injection ratio of 1:3 was utilized with a baseline bottom-hole pressure constraint of 5500 psi for 20 years. A post-injection period of 1000 years was simulated to monitor the CO 2 plume and its effects on the CO 2 storage reservoir and caprock integrity. The simulation results indicated that the impacts of the geochemical reactions on the porosity of the caprock were insignificant as it experienced a decrease of about 0.0003% at the end of the 1000-year post-injection monitoring. On the other hand, the maximum stress-induced porosity change was about a 1.4% increase, resulting in about 4% in permeability change. It was estimated that about 3.3% of the sequestered CO 2 in the formation interacted with the caprock. Despite these petrophysical property alterations and CO 2 interactions in the caprock, the caprock still maintained its elastic properties and was determined to be far from its failure.

04 OIL SHALES AND TAR SANDS↗
DOE OSTI2021

Chemical and structural characterization of EUV photoresists as a function of depth by standing-wave x-ray photoelectron spectroscopy

The success in the miniaturization of the electronic device constituents depends mostly on the photolithographic techniques. Recently, to achieve patterning at the sub-10-nm node, extreme ultraviolet (EUV) lithography has been introduced into high volume production. Continued scaling of EUV via increased numerical aperture to achieve nodes at 3-nm and below requires the development of fundamentally new patterning materials and new characterization methods. Current EUV-resist film thicknesses are in the 20- to 40-nm range, and further thickness reduction is required for the next generation. Therefore, interfaces become exceedingly important, and the properties of the resist film would be dominated by top and bottom interfacial effects. X-ray photoelectron spectroscopy (XPS) combined with standing-wave excitation (SW-XPS), a fairly new method in the EUV lithography field, previously had been largely applied in multilayers and superlattices for characterizing the composition and electronic structure of buried layers and interfaces as a function of depth. We applied the SW-XPS method to organic/inorganic photoresists to provide depth-selective information on their structural and chemical conditions of as a function of temperature, EUV exposure, different underlayers, and other fundamental parameters. As a first attempt, we perform an SW-XPS feasibility study on self-assembled monolayer (SAM) films after exposure to an electron beam. By SW-XPS, we determined that the interface between the Al 2 O 3 underlayer and the SAMs is smooth, with a mean roughness of about 0.2 nm. Moreover, we determined that the SAM chains are, on average, tilted by 1/430 deg off the sample normal. The SW-XPS results also suggest that the SAM is not a perfectly aligned and uniform monolayer, with some areas having thickness higher than a single monolayer. We demonstrated that SW-XPS can provide useful information on ultrathin materials with high potential for being used as a characterization method of organic/inorganic photoresists.

36 MATERIALS SCIENCE↗
NASA NTRS1972

6300 A quantum efficiency of the recombination mechanism in the night-time F layer.

Simultaneous airglow and electron content measurements made at Hawaii are used to infer the number of 6300 and 6364 A quanta produced per electron lost in the nighttime F layer of the ionosphere. The equation of continuity of electrons is then solved numerically to obtain the electron density profile, and the amount of quenching is estimated. This leads to the number of excitations of O(super-1 D) per O2(+) recombination (epsilon). We find, for an exospheric temperature of 1100 K, epsilon is equal to 1.1 plus or minus 0.6, in good agreement with Zipf's laboratory measurement at 300 K.

Brown, W. E.↗
NASA NTRS1975

Liquid jet pumped by rising gas bubbles

A two-phase mathematical model is proposed for calculating the induced turbulent vertical liquid flow. Bubbles provide a large buoyancy force and the associated drag on the liquid moves the liquid upward. The liquid pumped upward consists of the bubble wakes and the liquid brought into the jet region by turbulent entrainment. The expansion of the gas bubbles as they rise through the liquid is taken into account. The continuity and momentum equations are solved numerically for an axisymmetric air jet submerged in water. Water pumping rates are obtained as a function of air flow rate and depth of submergence. Comparisons are made with limited experimental information in the literature.

Hussain, N. A.↗
NASA NTRS1975

Liquid jet pumped by rising gas bubbles

From observations of a stream of gas bubbles rising through a liquid, a two-phase mathematical model is proposed for calculating the induced turbulent vertical liquid flow. The bubbles provide a large buoyancy force and the associated drag on the liquid moves the liquid upward. The liquid pumped upward consists of the bubble wakes and the liquid brought into the jet region by turbulent entrainment. The expansion of the gas bubbles as they arise through the liquid is taken into account. The continuity and momentum equations are solved numerically for an axisymmetric air jet submerged in water. Water pumping rates are obtained as a function of air flow rate and depth of submergence. Comparisons are made with limited experimental information in the literature.

Hussain, N. A.↗
NASA NTRS1975

A computer program for fitting smooth surfaces to three-dimensional aircraft configurations

A computer program developed to fit smooth surfaces to the component parts of three-dimensional aircraft configurations was described. The resulting equation definition of an aircraft numerical model is useful in obtaining continuous two-dimensional cross section plots in arbitrarily defined planes, local tangents, enriched surface plots and other pertinent geometric information; the geometry organization used as input to the program has become known as the Harris Wave Drag Geometry.

Craidon, C. B.↗
NASA NTRS1976

A limiting analysis for edge effects in angle-ply laminates

A zeroth order solution for edge effects in angle ply composite laminates using perturbation techniques and a limiting free body approach was developed. The general method of solution for laminates is developed and then applied to the special case of a graphite/epoxy laminate. Interlaminar stress distributions are obtained as a function of the laminate thickness to width ratio h/b and compared to existing numerical results. The solution predicts stable, continuous stress distributions, determines finite maximum tensile interlaminar normal stress for two laminates, and provides mathematical evidence for singular interlaminar shear stresses.

Hsu, P. W.↗
NASA NTRS1980

The compatible balancing approach to initialization, and four-dimensional data assimilation

A unified theoretical approach to both the four-dimensional assimilation of asynoptic data and the initialization problem is attempted. This approach relies on the derivation of certain relationships between geopotential tendencies and tendencies of the horizontal velocity field in primitive-equation models of atmospheric flow. The approach is worked out and analyzed in detail for some simple barotropic models. Certain independent results of numerical experiments for the time-continuous assimilation of real asynoptic meteorological data into a complex, baroclinic weather prediction model are discussed in the context of the present approach. Tentative inferences are drawn for practical assimilation procedures.

Ghil, M.↗
NASA NTRS1982

Seasonal variations of the high-latitude F region for strong convection

A plasma convection model is combined with an ionospheric-atmospheric composition model in order to study the seasonal variations of the high-latitude F region for geomagnetic conditions leading to strong convection. In a model calculation, a field tube of plasma is followed as it moves along a convection trajectory through a moving neutral atmosphere. Altitude profiles of the ion densities are obtained by solving the appropriate continuity, momentum, and energy equations including numerous high-latitude processes. It is found that the high-latitude ionosphere exhibits a significant UT variation both during the winter and summer. In general, the electron density at high-latitudes is lower in winter than in summer. In both summer and winter, the major region of low electron density is associated with the 'main' or mid-latitude' trough.

Sojka, J. J.↗
NASA NTRS1984

Transformation of two and three-dimensional regions by elliptic systems

The research during this period continued to expand the class of numerical algorithms that can be accurately and efficiently implemented on overlapping grids. Whereas previous calculations have been used to solve elliptic equations and to find the steady-state solution of parabolic equations, the present work is aimed towards developing time-accurate solution techniques for parabolic and hyperbolic equations. The primary difficulty here is in the correct treatment of the interior boundary nodes that must be updated at each iteration. The implementation of explicit methods is straightforward. However, the common practice of lagging these values when using an implicit methods leads to inconsistencies in the difference equation. One way to avoid this problem is to alternately calculate with an implicit and an explicit method on each subgrid. With this procedure, the explicit method generates boundary values at the next time level which are then used by the implicit step. It can be shown that when a backward implicit method is combined with a forward explicit method, the composite method is second order accurate and unconditionally stable for linear problems. A second area in which progress can be reported is in the distribution of grid points on curves and surfaces.

Mastin, E. W.↗
NASA NTRS1991

Charge collected by diffusion from an ion track under mixed boundary conditions

Charge-carrier diffusion from an ion track in a silicon substrate at least a few hundred microns thick is analyzed. The substrate upper surface is treated as reflective except for a small section, intended to represent a reverse-biased junction, which is treated as a sink. Total charge collected by the sink is calculated by assuming transport to be governed by an ambipolar diffusion equation with temporally constant and spatially uniform carrier lifetime and diffusion coefficient. Present results apply to a normally incident track but could easily be generalized to arbitrary track direction. The collected charge is found to depend on track length and on the electrostatic capacitance, rather than the area, of the sink. Theoretical predictions are compared to the results of a numerical simulation called the Poisson and Continuity Equation Solver (PISCES) for three cases and are found to agree within a factor of two in the worst case.

Edmonds, Larry D.↗
NASA NTRS1993

Effects of a parallel electric field and the geomagnetic field in the topside ionosphere on auroral and photoelectron energy distributions

The consequences of electric field acceleration and an inhomogeneous magnetic field on auroral electron energy distributions in the topside ionosphere are investigated. The one-dimensional, steady state electron transport equation includes elastic and inelastic collisions, an inhomogeneous magnetic field, and a field-aligned electric field. The case of a self-consistent polarization electric field is considered first. The self-consistent field is derived by solving the continuity equation for all ions of importance, including diffusion of O(+) and H(+), and the electron and ion energy equations to derive the electron and ion temperatures. The system of coupled electron transport, continuity, and energy equations is solved numerically. Recognizing observations of parallel electric fields of larger magnitude than the baseline case of the polarization field, the effect of two model fields on the electron distribution function is investigated. In one case the field is increased from the polarization field magnitude at 300 km to a maximum at the upper boundary of 800 km, and in another case a uniform field is added to the polarization field. Substantial perturbations of the low energy portion of the electron flux are produced: an upward directed electric field accelerates the downward directed flux of low-energy secondary electrons and decelerates the upward directed component. Above about 400 km the inhomogeneous magnetic field produces anisotropies in the angular distribution of the electron flux. The effects of the perturbed energy distributions on auroral spectral emission features are noted.

Min, Q.-L.↗