SEARCH · Engineering Papers

Results for “convergence rate”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

SourceResearch connectionFrom yearThrough year

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

1,310 records · Page 73

All sources NASA NTRS DOE OSTI

NASA NTRS1997

Dynamics Questions Associated with the AERCam Sprint Free-Flyer

The International Space Station will require the development of small robotic vehicles for such tasks as external inspection, monitoring of extravehicular activities (EVA's) and station build-up, and providing additional lighting at EVA worksites. The Autonomous EVA Robotic Camera (AERCam) family of free-flyers is currently being developed at NASA Johnson Space Center to perform these functions; the first member of the family is the AERCam Sprint vehicle. Many interesting dynamical questions are associated with the Sprint free-flyer. For instance, the reaction of a vehicle which is nearly spherically symmetric (such as Sprint) to a stuck-on thruster is significantly more complicated than that obtained for an idealized, perfectly spherical, spacecraft model. In particular, the real spacecraft will experience a form of forced nutation, with convergence towards either its major or minor principal axis, depending on both the applied torque and the mass properties of the vehicle. Furthermore, the body-fixed jet force vector may have a significant component along this principal axis, so giving rise to a considerable net linear acceleration of the spacecraft. The large velocity that can result is very important, as it may lead to collision with the nearby Orbiter, and is completely overlooked in the idealized analysis. This report will firstly briefly describe the stuck-on thruster dynamics of the real vehicle, and outline how the small products of inertia of the spacecraft determine the time constants of the motion. Secondly, the dynamical effects of a failed-off jet on the Sprint free-flyer will be described in more detail, and compared with the stuck-on thruster case. This will help to show whether the two malfunctions should be dealt with differently in flight. Finally, the stuck-on thruster detection software (known as the uncommanded motion algorithm) that is proposed to be flown on the Sprint vehicle will be analyzed, and all possible perturbation sources that may tend to give rise to false stuck-on thruster alarms quantified. It will be shown that false alarms can be triggered by discrete events such as a failed-off thruster, a glancing collision with Orbiter structure, or thruster saturation. They can also be triggered by a combination of errors introduced by constantly present sources such as thrust level errors, thruster misalignment, inertial cross-coupling, angular rate sensor noise, and structural vibrations induced by thruster firings. Based on this fact, the original plan to have the spacecraft enter safe mode whenever the uncommanded motion algorithm is triggered does not appear to be advisable. A better approach seems to be to provide the Sprint pilot with a warning whenever uncommanded motion is detected, and then allow the crew to determine whether safing the vehicle is appropriate.

Williams, Trevor↗

NASA NTRS2010

Investigation of Long-Term Impacts of Urbanization when Considering Global Warming for a Coastal Tropical Region

The overachieving goal of this project is to gain a better understanding of the climate impacts caused by the combined effects of land cover and land use (LCLU) changes and increasing global concentrations of green house gases (GHG) in tropical coastal areas, regions where global, regional and local climate phenomena converge, taking as the test case the densely populated northeast region of the Caribbean island of Puerto Rico. The research uses an integrated approach of high-resolution remote sensing information linked to a high resolution Regional Atmospheric Modeling System (RAMS), which was employed to perform ensembles of climate simulations (combining 2-LCLU and 2-GHG concentration scenarios). Reconstructed agricultural maps are used to define past LCLU, and combined with reconstructed sea surface temperatures (SST) for the same period form the PAST climate scenario (1951-1956); while the PRESENT scenario (2000-2004) was additionally supported by high resolution remote sensing data (10-m-res). The climate reconstruction approach is validated with available observed climate data from surface weather stations for both periods of time simulated. The selection of the past and present climate scenarios considers large-scale biases (i.e. ENSO/NAO) as reflected in the region of interest. Direct and cross comparison of the results is allowing quantifying single, combined, and competitive effects. Results indicate that global GHG have dominant effects on minimum temperatures (following regional tendencies), while urban sprawl dominates maximum temperatures. To further investigate impacts of land use the Bowen Ratio and the thermal response number (TRN) are analyzed. The Bowen ratio indicates that forestation of past agricultural high areas have an overwhelmingly mitigation effect on increasing temperatures observed in different LCLU scenarios, but when abandoned agricultural lands are located in plains, the resulting shrub/grass lands produce higher surface temperatures. The TRN (J/m^2/degC) is a surface property defined as the ratio of the surface net radiation to the rate of change in surface temperature, expresses how those fluxes are reacting to radiant energy inputs. Natural vegetated surfaces have a greater TRN than urban and barren surfaces because the net radiation processed by them is mostly used for latent heat and thermal storage heat rather than sensible heat (heating the air). Significant changes in TRN were observed in the metropolitan area of San Juan for the two analyzed periods reflecting a reduction of this variable in the present from the past consistent with increasing in thermal mass, or intense urbanization.

Gonalez, Jorge E.↗

NASA NTRS2013

Tuning and Robustness Analysis for the Orion Absolute Navigation System

The Orion Multi-Purpose Crew Vehicle (MPCV) is currently under development as NASA's next-generation spacecraft for exploration missions beyond Low Earth Orbit. The MPCV is set to perform an orbital test flight, termed Exploration Flight Test 1 (EFT-1), some time in late 2014. The navigation system for the Orion spacecraft is being designed in a Multi-Organizational Design Environment (MODE) team including contractor and NASA personnel. The system uses an Extended Kalman Filter to process measurements and determine the state. The design of the navigation system has undergone several iterations and modifications since its inception, and continues as a work-in-progress. This paper seeks to show the efforts made to-date in tuning the filter for the EFT-1 mission and instilling appropriate robustness into the system to meet the requirements of manned space ight. Filter performance is affected by many factors: data rates, sensor measurement errors, tuning, and others. This paper focuses mainly on the error characterization and tuning portion. Traditional efforts at tuning a navigation filter have centered around the observation/measurement noise and Gaussian process noise of the Extended Kalman Filter. While the Orion MODE team must certainly address those factors, the team is also looking at residual edit thresholds and measurement underweighting as tuning tools. Tuning analysis is presented with open loop Monte-Carlo simulation results showing statistical errors bounded by the 3-sigma filter uncertainty covariance. The Orion filter design uses 24 Exponentially Correlated Random Variable (ECRV) parameters to estimate the accel/gyro misalignment and nonorthogonality. By design, the time constant and noise terms of these ECRV parameters were set to manufacturer specifications and not used as tuning parameters. They are included in the filter as a more analytically correct method of modeling uncertainties than ad-hoc tuning of the process noise. Tuning is explored for the powered-flight ascent phase, where measurements are scarce and unmodelled vehicle accelerations dominate. On orbit, there are important trade-off cases between process and measurement noise. On entry, there are considerations about trading performance accuracy for robustness. Process Noise is divided into powered flight and coasting ight and can be adjusted for each phase and mode of the Orion EFT-1 mission. Measurement noise is used for the integrated velocity measurements during pad alignment. It is also used for Global Positioning System (GPS) pseudorange and delta- range measurements during the rest of the flight. The robustness effort has been focused on maintaining filter convergence and performance in the presence of unmodeled error sources. These include unmodeled forces on the vehicle and uncorrected errors on the sensor measurements. Orion uses a single-frequency, non-keyed GPS receiver, so the effects due to signal distortion in Earth's ionosphere and troposphere are present in the raw measurements. Results are presented showing the efforts to compensate for these errors as well as characterize the residual effect for measurement noise tuning. Another robustness tool in use is tuning the residual edit thresholds. The trade-off between noise tuning and edit thresholds is explored in the context of robustness to errors in dynamics models and sensor measurements. Measurement underweighting is also presented as a method of additional robustness when processing highly accurate measurements in the presence of large filter uncertainties.

Holt, Greg N.↗

DOE OSTI2023

Cost-constrained adaptive simulations of transient spray combustion in a gas turbine combustor

Predictive high-fidelity simulations of turbulent spray combustion must capture the combined effects of complex chemistry, multiphase evaporating flow and spray-flame interactions to achieve physical accuracy. Finite-rate chemistry (FRC) combined with a realistic chemical mechanism is a combustion model well-suited for this purpose, but has a high computational cost due to the large number and stiffness of transported chemical species. In contrast, flamelet-based models achieve lower cost by transporting a small number of quantities of reduced stiffness, but assumptions regarding local flame topology, boundary conditions and inter-phase coupling limit their physical accuracy. Recently, the Pareto-efficient combustion (PEC) framework was developed to dynamically assign combustion models based on local cost and accuracy metrics in gas-phase reacting flows. In this work, we extend this PEC framework to spray combustion through the rigorous analysis of the multiphase coupling terms in the governing equations. The derivation shows that spray evaporation causes errors in the prediction of species mass fractions for flamelet-based models due to the sensitivity of the local thermo-chemical state to changes in composition caused by fuel vaporization across combustion regimes present in practical spray combustion devices. Sub-model assignment is formulated as a multiple-choice knapsack problem, where computational cost is directly controlled through the fraction of the domain assigned to the FRC sub-model. The extended PEC formulation is applied to the simulation of a realistic rich-quench-lean gas turbine combustor at steady-state conditions, as well as transient operation resulting in lean blow-out (LBO). Analysis of transient simulations during LBO demonstrates the extended PEC formulation’s capacity to dynamically adapt to changing conditions within the combustor. Transient combustor dynamics are shown to approach convergence with limited increases in computational cost, while retaining substantial computational cost reduction compared to monolithic FRC simulations. Through PEC simulations with increasing fractions of the domain assigned to FRC, monolithic flamelet simulations are shown to over-predict flame stability during LBO. Finally, the extended PEC formulation is thus shown to overcome deficiencies of monolithic models by controlling modeling error for multiphase combustion modeling.

33 ADVANCED PROPULSION SYSTEMS↗

DOE OSTI2025

2D CFD of lean premixed hydrogen–air flame quenching under locomotive engine conditions

Hydrogen (H 2 ) is a promising fuel for reducing emissions in heavy-duty internal combustion engines (ICEs), but its low quenching distance increases the risk of flame propagation into narrow crevice regions, such as piston-liner gaps. This work uses detailed CFD simulations with finite-rate chemistry to investigate premixed H 2 –air flame quenching in a two-dimensional (2D) region consistent with the piston-liner gap of a diesel ICE. Model accuracy was assessed by comparison with experimentally measured quenching distances in an annular stepwise diverging tube (ASDT). A parametric study was conducted to assess the influence of crevice width (0.05–1.18 mm), pressure (50–150 bar), unburned gas temperature (431–573 K), and equivalence ratio ( Φ= 0.3–0.6). Results show that the critical Péclet number for flame survival is no greater than 3.25, consistent with prior literature, even under ultra-lean and high-pressure conditions (Φ ≤ 0.3, > 50 bar). Additionally, flames with Péclet numbers exceeding 6.15 exhibited front wrinkling, suggesting the onset of velocity-driven instabilities and enhanced flame robustness. These findings help define thresholds for flame quenching in confined geometries and support the safe design of H 2 fueled ICEs.

08 HYDROGEN↗

NASA NTRS2023

Broader Impacts and Outreach Videos on Ocean Physics and Biogeochemistry

This is just a preliminary NF1676 in advance of filming for a project....A convergent team of experts in biology, chemistry, physics, engineering, mathematics, and computational modeling examine how dynamic and coupled phytoplankton-pathogen-particle-predator linkages coalesce to explain the observed high spatial variability in the efficiency of the export of particulate organic carbon (POC) to the deep ocean. They elucidate and quantify the linkages between viruses and ballast minerals to increase understanding of carbon cycling in the oceanic biological carbon pump and the impact of viruses within it. By providing knowledge than can be used to improve the parameterization of carbon export in Earth system models, the project will help reduce uncertainty in regional marine biogeochemical projections, potentially improving marine ecosystem and fisheries management on timescales from seasons to decades. The project includes activities that provide teaching resources and hands-on training to educators within a ‘Tools of Science’ program that provides a simple and succinct way to communicate the process of scientific research to students in a way that is useful to teachers. The project couples laboratory-based experiments on model host-virus-grazer systems with extensive field based observational and manipulative studies on natural populations of diatoms and coccolithophores, the two phytoplankton groups that account for most of the estimated particulate organic matter flux to the deep ocean. Experiments and measurements integrate diagnostic biological and chemical controls on infection and particle coagulation theory with microscale physics and grazing to quantify links to each hypothesized export mechanism under field-relevant turbulent conditions. Cutting-edge engineering and analytical tools are used to diagnose and track infection dynamics while characterizing and quantifying particle aggregation and disaggregation, mineral dissolution, sinking dynamics, grazing rates, and fecal pellet production at unprecedented resolution and under well-defined, microscale physical regimes. Field campaigns elucidate the relative efficiency of hypothesized mechanisms in stimulating POC export in natural blooms, while providing bulk and size-resolved estimates of POC flux.

Daniel B Whitt↗

DOE OSTI2025

Large‐Scale 2D Perovskite Nanocrystals Photodetector Array via Ultrasonic Spray Synthesis

Abstract 2D perovskite (PVSK) single crystals have received significant attention due to their unique optical and optoelectronic properties. However, current synthesis methods face limitations, particularly in large‐area fabrication, which remain critical barriers to practical applications. In this study, the synthesis of red/green/purple‐blue‐colored 2D PVSK nanocrystals over a large area (4‐inch wafer) and the fabrication of high‐performance photodetector arrays are presented via a facile yet efficient spray‐coating approach with a liquid‐bridge transport effect. The photodetector array achieves 100% working yield, high photo‐responsivity (1.5 × 10 6 A W −1 ) and specific‐detectivity (1.1 × 10 16 Jones) with competitive photomapping characteristics. An intelligent vision system for automatic shape recognition is further demonstrated with a recognition rate exceeding 90%. This study provides significant advances in the scalable synthesis of nanoscale 2D PVSK crystals, their integration into large‐area optoelectronic devices, and their potential use in artificial‐intelligence systems.

Lee, Yoon Ho [Davidson School of Chemical Engineer↗

NASA NTRS2014

Mean Flow Augmented Acoustics in Rocket Systems

Oscillatory motion in solid rocket motors and liquid engines has long been a subject of concern. Many rockets display violent fluctuations in pressure, velocity, and temperature originating from the complex interactions between the combustion process and gas dynamics. The customary approach to modeling acoustic waves inside a rocket chamber is to apply the classical inhomogeneous wave equation to the combustion gas. The assumption of a linear, non-dissipative wave in a quiescent fluid remains valid while the acoustic amplitudes are small and local gas velocities stay below Mach 0.2. The converging section of a rocket nozzle, where gradients in pressure, density, and velocity become large, is a notable region where this approach is not applicable. The expulsion of unsteady energy through the nozzle of a rocket is identified as the predominate source of acoustic damping for most rocket systems. An accurate model of the acoustic behavior within this region where acoustic modes are influenced by the presence of a steady mean flow is required for reliable stability predictions. Recently, an approach to address nozzle damping with mean flow effects was implemented by French [1]. This new approach extends the work originated by Sigman and Zinn [2] by solving the acoustic velocity potential equation (AVPE) formulated by perturbing the Euler equations [3]. The acoustic velocity potential (psi) describing the acoustic wave motion in the presence of an inhomogeneous steady high-speed flow is defined by, (del squared)(psi) − (lambda/c)(exp 2)(psi) − M(dot)[M(dot)(del)(del(psi))] − 2(lambda(M/c) + (M(dot)del(M))(dot)del(psi)−2(lambda)(psi)[M(dot)del(1/c)]=0 (1) with M as the Mach vector, c as the speed of sound, and lambda as the complex eigenvalue. French apply the finite volume method to solve the steady flow field within the combustion chamber and nozzle with inviscid walls. The complex eigenvalues and eigenvector are determined with the use of the ARPACK eigensolver. The present study employs the COMSOL Multphysics framework to solve the coupled eigenvalue problem using the finite element approach. The study requires one way coupling of the CFD High Mach Number Flow (HMNF) and mathematics module. The HMNF module evaluated the gas flow inside of a solid rocket motor using St. Robert's law modeling solid propellant burn rate, slip boundary conditions, and the supersonic outflow condition. Results from the HMNF model are used by the coefficient form of the mathematics module to determine the eigenvalues of the AVPE. The mathematics model is truncated at the nozzle sonic line, where a zero flux boundary condition is self-satisfying. The remaining boundaries are modeled with a zero flux boundary condition, assuming zero acoustic absorption on all surfaces. Pertinent results from these analyses are the complex valued eigenvalue and eigenvectors. Comparisons are made to the French results to evaluate the modeling approach. A comparison of the French results with that of the present analysis is displayed in figures 1 and 2, respectively. The graphic shows the first tangential eigenvector's real (a) and imaginary (b) values.

Fischbach, Sean R.↗

NASA NTRS1997

Effects of Buoyancy on Laminar and Turbulent Premixed V-Flame

Turbulent combustion occurs naturally in almost all combustion systems and involves complex dynamic coupling of chemical and fluid mechanical processes. It is considered as one of the most challenging combustion research problems today. Though buoyancy has little effect on power generating systems operating under high pressures (e.g., IC engines and turbines), flames in atmospheric burners and the operation of small to medium furnaces and boilers are profoundly affected by buoyancy. Changes in burner orientation impacts on their blow-off, flash-back and extinction limits, and their range of operation, burning rate, heat transfer, and emissions. Theoretically, buoyancy is often neglected in turbulent combustion models. Yet the modeling results are routinely compared with experiments of open laboratory flames that are obviously affected by buoyancy. This inconsistency is an obstacle to reconciling experiments and theories. Consequently, a fundamental understanding of the coupling between turbulent flames and buoyancy is significant to both turbulent combustion science and applications. The overall effect of buoyancy relates to the dynamic interaction between the flame and its surrounding, i.e., the so-called elliptical problem. The overall flame shape, its flowfield, stability, and mean and local burning rates are dictated by both upstream and downstream boundary conditions. In steady propagating premixed flames, buoyancy affects the products region downstream of the flame zone. These effects are manifested upstream through the mean and fluctuating pressure fields to influence flame stretch and flame wrinkling. Intuitively, the effects buoyancy should diminish with increasing flow momentum. This is the justification for excluding buoyancy in turbulent combustion models that treats high Reynolds number flows. The objectives of our experimental research program is to elucidate flame-buoyancy coupling processes in laminar and turbulent premixed flames, and to characterize microgravity (micro g) premixed flames. The results are used to derive appropriate scaling parameters for guiding the development of theoretical models to include the effects of buoyancy. Knowledge gain from the analysis will also contribute to further understanding of the elliptical nature of premixed flames. Our current emphasis is to examine the momentum limit above which the effects of buoyancy would become insignificant. This is accomplished by comparing the flowfields and the mean properties of normal gravity flames (+g), and reversed gravity flames (-g, up-side-down flames) at different flow velocities and turbulence intensities. Microgravity (micro g) flames experiments provide the key reference data to reconcile the differences between flames in +g and -g. As flame configuration has significant impact on premixed flames characteristics we have studied axi-symmetric conical flames and plane-symmetric rod-stabilized v-flames. The two configurations produce distinct features that dictates how the flames couple with buoyancy. In a conical flame, the hot products plume completely envelopes the flame cone and shields the flame from direct interaction with the ambient air. The plume originates at the burner rim and generates a divergent flowfield. In comparison, the products region of v-flames forms between the twin flame sheets and it is convergent towards the center-plane. Interaction with ambient air is limited to the two end regions of the stabilized rod and beyond the flame sheets.

Cheng, Robert K.↗

NASA NTRS2013

Issues and Strategies in Solving Multidisciplinary Optimization Problems

Optimization research at NASA Glenn Research Center has addressed the design of structures, aircraft and airbreathing propulsion engines. The accumulated multidisciplinary design activity is collected under a testbed entitled COMETBOARDS. Several issues were encountered during the solution of the problems. Four issues and the strategies adapted for their resolution are discussed. This is followed by a discussion on analytical methods that is limited to structural design application. An optimization process can lead to an inefficient local solution. This deficiency was encountered during design of an engine component. The limitation was overcome through an augmentation of animation into optimization. Optimum solutions obtained were infeasible for aircraft and airbreathing propulsion engine problems. Alleviation of this deficiency required a cascading of multiple algorithms. Profile optimization of a beam produced an irregular shape. Engineering intuition restored the regular shape for the beam. The solution obtained for a cylindrical shell by a subproblem strategy converged to a design that can be difficult to manufacture. Resolution of this issue remains a challenge. The issues and resolutions are illustrated through a set of problems: Design of an engine component, Synthesis of a subsonic aircraft, Operation optimization of a supersonic engine, Design of a wave-rotor-topping device, Profile optimization of a cantilever beam, and Design of a cylindrical shell. This chapter provides a cursory account of the issues. Cited references provide detailed discussion on the topics. Design of a structure can also be generated by traditional method and the stochastic design concept. Merits and limitations of the three methods (traditional method, optimization method and stochastic concept) are illustrated. In the traditional method, the constraints are manipulated to obtain the design and weight is back calculated. In design optimization, the weight of a structure becomes the merit function with constraints imposed on failure modes and an optimization algorithm is used to generate the solution. Stochastic design concept accounts for uncertainties in loads, material properties, and other parameters and solution is obtained by solving a design optimization problem for a specified reliability. Acceptable solutions can be produced by all the three methods. The variation in the weight calculated by the methods was found to be modest. Some variation was noticed in designs calculated by the methods. The variation may be attributed to structural indeterminacy. It is prudent to develop design by all three methods prior to its fabrication. The traditional design method can be improved when the simplified sensitivities of the behavior constraint is used. Such sensitivity can reduce design calculations and may have a potential to unify the traditional and optimization methods. Weight versus reliability traced out an inverted-S-shaped graph. The center of the graph corresponded to mean valued design. A heavy design with weight approaching infinity could be produced for a near-zero rate of failure. Weight can be reduced to a small value for a most failure-prone design. Probabilistic modeling of load and material properties remained a challenge.

Patnaik, Surya↗

DOE OSTI2023

Constraints on simulated past Arctic amplification and lapse rate feedback from observations

The Arctic has warmed more rapidly than the global mean during the past few decades. The lapse rate feedback (LRF) has been identified as being a large contributor to the Arctic amplification (AA) of climate change. This particular feedback arises from the vertically non-uniform warming of the troposphere, which in the Arctic emerges as strong near-surface and muted free-tropospheric warming. Stable stratification and meridional energy transport are two characteristic processes that are evoked as causes for this vertical warming structure. Our aim is to constrain these governing processes by making use of detailed observations in combination with the large climate model ensemble of the sixth Coupled Model Intercomparison Project (CMIP6). We build on the result that CMIP6 models show a large spread in AA and Arctic LRF, which are positively correlated for the historical period of 1951–2014. Thus, we present process-oriented constraints by linking characteristics of the current climate to historical climate simulations. In particular, we compare a large consortium of present-day observations to co-located model data from subsets that show a weak and strong simulated AA and Arctic LRF in the past. Our analyses suggest that the vertical temperature structure of the Arctic boundary layer is more realistically depicted in climate models with weak (w) AA and Arctic LRF (CMIP6/w) in the past. In particular, CMIP6/w models show stronger inversions in the present climate for boreal autumn and winter and over sea ice, which is more consistent with the observations. These results are based on observations from the year-long Multidisciplinary Drifting Observatory for the Study of Arctic Climate (MOSAiC) expedition in the central Arctic, long-term measurements at the Utqiagvik site in Alaska, USA, and dropsonde temperature profiling from aircraft campaigns in the Fram Strait. In addition, the atmospheric energy transport from lower latitudes that can further mediate the warming structure in the free troposphere is more realistically represented by CMIP6/w models. In particular, CMIP6/w models systemically simulate a weaker Arctic atmospheric energy transport convergence in the present climate for boreal autumn and winter, which is more consistent with fifth generation reanalysis of the European Centre for Medium-Range Weather Forecasts (ERA5). We further show a positive relationship between the magnitude of the present-day transport convergence and the strength of past AA. With respect to the Arctic LRF, we find links between the changes in transport pathways that drive vertical warming structures and local differences in the LRF. This highlights the mediating influence of advection on the Arctic LRF and motivates deeper studies to explicitly link spatial patterns of Arctic feedbacks to changes in the large-scale circulation.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2022

Robust Carbon Dioxide Plume Imaging Using Joint Tomographic Inversion of Seismic Onset Time and Distributed Pressure and Temperature Measurements (Final Report)

We develop and demonstrate rapid and cost-effective methodologies for spatiotemporal tracking of CO2 plumes during geologic sequestration using joint inversion of seismic data and distributed pressure and temperature measurements. Key elements of our methodology are: (a) a computationally efficient approach to pressure and temperature propagation, (b) analysis of time lapse seismic data using a novel ‘seismic onset time’ approach to detect fluid front propagation, and (c) data assimilation and uncertainty assessment via joint inversion of pressure, temperature and time lapse seismic data, and (d) validating the numerical tomographic inversion using a CO2 injection demonstration projects, specifically data collected from the from the Petra Nova Parish Holdings CCUS project in the West Ranch Field, Texas and the Chester-16 reef CO2 injection site in Northern Michigan which is part of the DOE Midwestern Carbon Sequestration Project. The research team is led by Texas A&M University and includes Battelle as a subcontractor with support from Shell, Anadarko, Chevron and JX Nippon. A carbon dioxide (CO2) water-alternating-gas (WAG) pilot was conducted to gain insights into tertiary oil recovery potential via CO2 flood in the West Ranch Field as part of the Petra Nova project, the world’s largest post-combustion CO2 capture and utilization initiative. With a fluvial formation geology and large contrasts in permeability, this is a challenging and novel application of CO2 enhanced oil recovery (EOR). We build a predictive dynamic model of the subsurface that incorporates the multiphase and compositional data acquired during the pilot operation. The calibrated model is used for the carbon dioxide plume imaging. The study began with an initialization of the pilot sector model extracted from a calibrated full-field model. The pilot model calibration follows a two-step hierarchical workflow. First, we performed a large-scale update of the permeability distribution by integrating available bottomhole pressure and multiphase production data. In the second step, local permeability field is fine-tuned using a streamline-based method to match CO2 breakthrough times at the producers. The predictive capability of the calibrated model was verified through two blind validation tests: (1) the model showed good agreement with saturation logs acquired at two observation wells; and (2) the model reproduced the CO2 recovery as a fraction of the injected CO2. The use of seismic onset times has shown great promise for integrating near-continuous seismic surveys for updating geologic models. In this study, we analyze the impact of seismic survey frequency on the onset time approach aiming to extend the application of onset time to infrequent seismic surveys. In addition, we quantitatively examine the nonlinearity of the onset time method and compare it to the commonly used amplitude inversion method. We carry out a sensitivity analysis of seismic survey frequency based on the complete seismic survey data (over 175 surveys) of steam injection in a heavy oil reservoir (Peace River Unit) in Canada. Our results show that an adequate onset time map can be obtained from the infrequent seismic surveys by interpolation between seismic surveys as long as there is no change in the dominant underlying physics between the successive surveys. The study also shows that nonlinearity of the onset time method can be -smaller than that of the amplitude inversion method by several orders of magnitude. Application to the Brugge benchmark case shows that the onset time method obtains comparable permeability update as the traditional seismic amplitude inversion method with faster computation and improved convergence characteristics. We extend the streamline-based data integration approach to incorporate distributed temperature sensor (DTS) data using the concept of thermal tracer travel time. Then, a hierarchical workflow composed of evolutionary and streamline methods is employed to jointly history match the DTS and pressure data. Finally, CO2 saturation and streamline maps are used to visualize the CO2 plume movement during the sequestration process. The hierarchical workflow is applied to a carbon sequestration project in a carbonate reef reservoir within the Northern Niagaran Pinnacle Reef Trend in Michigan, USA. The monitoring data set consists of distributed temperature sensing (DTS) data acquired at the injection well and a monitoring well, flowing bottom-hole pressure data at the injection well, and time-lapse pressure measurements at several locations along the monitoring well. The history matching results indicate that the CO2 movement is mostly restricted to the intended zones of injection which is consistent with an independent warm-back analysis of the temperature data. In addition to employing simulation models and inverse methods for CO2 plume imaging, we also initialized a data-driven technology for detecting inter-well connectivity based on production and pressure data. Our machine-learning framework is built on the statistical recurrent unit (SRU) model and interprets well-based injection/production data into inter-well connectivity without relying on a geologic model. We test it on synthetic and field-scale CO2 EOR projects utilizing the water-alternating-gas (WAG) process. The validation of the proposed data-driven inter-well connectivity assessment is performed using synthetic data from simulation models where inter-well connectivity can be easily measured using the streamline-based flux allocation. The SRU model is shown to offer excellent prediction performance on the synthetic case. Despite significant measurement noise and frequent well shut-ins imposed in the field-scale case, the SRU model offers good prediction accuracy, the overall relative error of the phase production rates at most producers ranges from 10% to 30%. It is shown that the dominant connections identified by the data-driven method and streamline method are in close agreement. Texas A&M University, the lead organization in the project, was primarily responsible for the development of tomographic approaches for CO2 plume mapping in conjunction with distributed pressure, temperature and seismic onset time data. Battelle, as a subcontractor, was primarily responsible for the development of analytical and empirical methods for analyzing transient injection rate and pressure data from point/line sources such as injection and monitoring wells. An additional area of emphasis for Battelle was the use of machine learning for such tasks as inferring reservoir connectivity information from injection-production data, and identifying variable importance for machine learning-based proxy models developed from full-physics simulations. The two organizations also collaborated on the application of the tomographic inversion methodology for a field data set.

02 PETROLEUM↗

NASA NTRSDate not supplied

NASA Earth Systems Digital Twins (ESDT)

"Similarly to artificial intelligence, which is now revolutionizing many aspects of our daily lives, Earth system digital twin technologies have the potential to revolutionize the way Earth Science research will be conducted in the future, and how results and knowledge from this research will provide information to support decision making and yield impactful societal benefits. An Earth System Digital Twin or ESDT is a dynamic and interactive information system that first provides a digital replica of the past and current states of the Earth or Earth system as accurately and timely as possible; second, allows for computing forecasts of future states under nominal assumptions and based on the current replica; and third, offers the capability to investigate many hypothetical scenarios under varying impact assumptions. In other words, an ESDT provides the integrated What-Now, What-Next, and What-If pictures of the Earth or Earth system, by continuously ingesting newly observed data and by leveraging multiple interconnected models, machine learning as well advanced computing and visualization capabilities. Digital twins have been developed in engineering since 2002, but the interest in digital twins for the Earth domain is more recent and stems from the convergence of several developments: - The huge amount of diverse data that has now been collected continuously for more than 50 years, and that is becoming more and more difficult to access, understand, and utilize. - At the same time, because of climate change and its impacts the information produced by all of this data is becoming of interest to many new non-traditional users for analyzing and predicting various phenomena. - Because of advances in computational and visualization capabilities and the parallel unprecedented development of machine learning (ML), extracting relevant information from these large amounts of data and running complex models faster has become possible. As a result, it is becoming necessary and possible to build intuitive and interactive frameworks that will enable users with various skill levels and/or organizational hierarchy levels to easily access large amounts of targeted information along with the relevant tools and models (Earth system and human activity models), to support them in analyzing and visualizing this information, to help them understand interactions among models, to visualize the potential outcomes of various impacts, and to support decision or policy making. The full power of digital twins is that, through an integrated representation and standardized tools and software technologies, the same digital replica can address the needs of multiple users at various resolutions (spatial and temporal) and for various applications (science, economic, policy, etc.) – “from farmer to scientist”. With all these interests at stake, the challenges of building optimal digital twins are many and complex. The first challenge is to determine if a Digital Twin should be global or local, and multi-domain or thematic. For example, some domains such as Climate or Weather will require a global Digital Twin or Digital Twin capabilities while science areas such as Biodiversity might be more local. We can also envision that multiple thematic ESDTs, e.g., Air Quality, Wildfires, Hydrology could be federated or provide input to other ESDTs, either on a regional level or to a more global ESDT. Overall, we can imagine a future “web” of Digital Twins co-existing in a hierarchy or in a network, and capable of being connected or federated depending on the needs. This last point brings up the very important challenge of interoperability, including standards and protocols that will need to be built into these systems from the beginning. Each individual digital twin would have full flexibility in internal construction but would need standards-based interfaces (input and output) or hooks to make it compatible with others. Another challenge when building digital twins will be to decide how to organize each digital replica. Based on the applications targeted by the DT under implementation, various amounts and types of raw data, Analysis Ready Data (ARD) and information will need to be incorporated. Depending on the required latencies and needs of the users, various solutions can be considered, including Data Cubes, Data Lakes, pointers, or computing information on demand. We envision that each ESDT will choose a solution adapted to its specific objectives. Another important challenge is the type(s) of visualization that will be used, as well as the level of interactivity and refresh rate that will be required. Again, this will depend on the objectives of the ESDT, but also on the various users’ needs. In most cases, several types of visualizations and human interfaces will need to be offered depending on the projected users of that system. In parallel to the challenges highlighted above, there are also many tools and technologies that will need to be developed or improved for all types of digital twins. Among those are improved machine learning technologies, for example providing explainability, but also ML techniques for causality and providing a better integration of physics models. Additionally, reliable uncertainty quantification methods will be needed for all ESDT components, from validating data fusion and assimilation to assessing the accuracy of ML models and weighing the values of decisions supported by those systems. This presentation introduces the ESDT concept, presents several ESDT use cases, and a proposed ESDT architecture framework, as well as various technologies being developed by the Advanced Information Systems Technology (AIST) Program."

Earth Science Remote Sensing; Information Systems↗

NASA NTRS2020

AERoBOND Project Summary

Under NASA’s Convergent Aeronautics Solutions (CAS) project, the Adhesive-Free Bonding of Complex Composites (AERoBOND) project investigated off-stoichiometric epoxy polymers for fast, reliable assembly of epoxy matrix composite structures. The project goal was to demonstrate feasibility of the AERoBOND joining method by demonstrating mechanical properties greater than 80% of conventional co-cured materials while reducing structure weight by 1%. The project consisted of three convergent research areas: material and process development, systems analysis, and material and process modeling. Material and process development was the largest component of AERoBOND with approximately 6 FTE and 1WYE of support to formulate and characterize new resins, prepare carbon fiber prepregs, fabricate laminates, measure mechanical properties, analyze failure results, and select material and process improvements. The systems analysis activity estimated the potential reduction in part count and aircraft weight by comparing models of composite wing boxes with no fasteners (co-cured structure), fasteners in major joints (co-cured stringers), and fasteners in all joints. The materials and process modeling activity included a molecular model of the AERoBOND materials system to predict mechanical properties of resins with offset stoichiometry and a process model to predict the effect of resin formulation and processing conditions on the extent of mixing and degree of cure in a finished joint. As the number of airline passenger trips doubles in the next 20 years (IATA/Tourism Economics Air Passenger Forecasts, April 2019), the increased demand for new commercial aircraft is now the single greatest technical challenge to the airframe manufacturing industry. To meet efficiency requirements, new aircraft must be fabricated primarily from high performance structural composites, but manufacturing processes are inherently slow with the largest bottleneck attributed to assembly and installation of fasteners (NASA/TM–2019-220428). Manufactures of commercial transport aircraft are compelled to install more than 100,000 redundant fasteners into bonded joints to prevent failures due to unpredictable weak bonds. In structural adhesive bonds, the interface between adherend and adhesive is nearly two-dimensional making it susceptible to minute quantities of contamination, which can cause weak bonds. Currently, bond strength assessment is only possible through destructive testing (i.e., breaking the joint). For these reasons, regulatory organizations such as the Federal Aviation Administration (FAA) often require redundant load paths in secondary-bonded, primary-structures to alleviate concerns with bond performance. The AERoBOND process enables reflow of matrix resin during assembly to eliminate the material discontinuity at the interface, thereby eliminating the dependence of mechanical performance on interfacial adhesion. The AERoBOND joint is equivalent to the interlaminar region obtained during a co-cure process, so joint performance depends on the cohesive properties of the matrix resin. Conventional co-cured structures, although too costly and complex for large-scale manufacturing, are trusted by manufacturers and regulators, and are certified for flight with few or no redundant fasteners.Systems analysis performed on a composite wing model at the scale of a single-aisle commercial transport aircraft indicated that >20,000 redundant fasteners per wing could be eliminated by implementing the AERoBOND joining method. A total weight reduction of 15% was predicted in a wing box by eliminating fasteners and thinning components that must no longer support localized fastener loads and accommodate fastener dimensions. Interlaminar shear fracture toughness measured by the end-notched flexure test was greater than 1 kJ/m2 (nearly 140% of the co-cured benchmark property), which is greatly in excess of the project goals for mechanical properties. Testing was planned to measure interlaminar tensile fracture toughness as well as interlaminar tensile and shear strengths using the same AERoBOND configuration, but was delayed due to closure of LaRC facilities during the COVID-19 pandemic. The AERoBOND process model is partially validated and available for experimental use. It allows the user to input AERoBOND process parameters such as material composition, laminate configuration, and cure cycle to predict the final cure state of the AERoBOND joint. A preliminary, multi-scale material model was developed to predict AERoBOND joint mechanical properties (stiffness and strength) based on the cure state of the joint provided by the process model. The timing for transition of this technology within NASA is excellent as NASA initiates new enduring projects to address composites manufacturing rate challenges. AERoBOND technology is well suited to AAVP/AATT objectives for rapid manufacturing of a composite wing. A minimal effort (1 FTE/$15k procurement/0 WYE) is proposed in FY21 to continue a minor mechanical testing effort and maintain a SAA with ASX composites to develop commercial quality prepreg material. An RFI with the composites industry is suggested to quantify the technology gap between the current TRL and the TRL needed for transition to industry. A moderate effort [3-4 FTE/$150k/1 WYE (~$115k)] is proposed in FY22 for the “high rate composites manufacturing” project currently in planning. The partnership with ASX Composites will be expanded to produce material for sub-element/element-scale “panel-off” activities. Industry partnerships with airframe manufacturers is an expected component to explore damage tolerance and environmental stability. Further development of multi-scale modeling tools (process model, meso-scale model, and molecular model) is planned to enhance and deliver tools for rapid manufacturing infusion.

Frank Louis Palmieri↗