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.

At least 937 records · Page 52

All sources NASA NTRS DOE OSTI

NASA NTRS1999

Climatic Variability of Precipitation from the Seasonal Cycle to ENSO Using GPCP's Merged Data Product and SSM/I-Based Microwave Estimates

Satellite estimates and gauge observations of precipitation are useful in understanding the water cycle, analyzing climatic variability, and validating climate models. The Global Precipitation Climatology Project (GPCP) released a community merged precipitation data set for the period July 1987 through the present, and has recently extended that data set back to 1986. One objective of this study is to use GPCP estimates to describe and quantify the seasonal variation of precipitation, with emphasis on the Asian summer monsoon. Another focus is the 1997-98 El Nino Southern Oscillation (ENSO) and associated extreme precipitation events. The summer monsoon tends to be drier than normal in El Nino ears. This was not observed for 1997 or 1998, while for 1997 the NCEP model produced the largest summer rain rates over India in years. This inconsistency will be examined. The average annual global precipitation rate is 2.7 mm day as estimated by GPCP, which is similar to values computed from long-term climatologies. From 30 deg N to 30 deg S the average precipitation rate is 2.7 mm day over land with a maximum in the annual cycle occurring in February-March, when the Amazon basin receives abundant rainfall. The average precipitation rate is 3.1 mm day over the tropical oceans, with a peak earlier in the season (November-December), corresponding with the transition from a strong Pacific Intertropical Convergence Zone (ITCZ) from June to November to a strong South Pacific Convergence Zone (SPCZ) from December to March. The seasonal evolution of C, C, the Asian summer monsoon stands out with rains in excess of 15 mm day off the coast of Burma in June. The GPROF pentad data also captures the onset of the tropical Pacific rainfall patterns associated with the 1997-98 ENSO. From February to October 1997 at least four rain-producing systems traveled from West to East in the equatorial corridor. A rapid transition from El Nino to La Nina conditions occurred in May-June 1998. GPCP and GPROF were used to construct precipitation-based ENSO indices to monitor El Ninos (EL) and La Ninas and (LI).

Curtis, Scott↗

DOE OSTI2023

Performance Characteristics of a Rotating Detonation Combustor Exiting Into a Pressurized Plenum to Simulate Gas Turbine Inlet

Abstract The present study aims to experimentally characterize the performance of a rotating detonation combustion (RDC) system integrated with a pressurized downstream plenum to simulate the high-pressure inlet conditions of power-generating gas turbines. A thorough understanding of the operational behavior including wave mode behavior, static pressure profile along the combustor length, and dynamic features of pressure fluctuations is crucial for successful integration of RDC with the turbine. In this study, two RDC configurations are investigated, RDC with a constant area annulus and RDC with a converging nozzle. In both cases, the RDC flow exited into a plenum chamber kept at pressures varying from 155 kPa to 330 kPa. RDC was operated on methane and oxygen-enriched air to represent reactants used in land-based power generation. Experiments were conducted for the two RDCs configurations operated at three reactant mass flow rates (0.23, 0.32, and 0.46 kg/s). The RDC performance is characterized by time-averaged static pressure measurements, and wave velocity determined by ionization probes. In addition, dynamic pressure measurements were recorded both inside and near the exit of RDC channel to investigate wave interactions between RDC and downstream plenum. Results show that the RDC with the converging nozzle achieved superior performance while minimizing detrimental interactions with the reflected shock and/or acoustic waves from the downstream plenum.

Engineering↗

NASA NTRS1999

Efficient Multi-Stage Time Marching for Viscous Flows via Local Preconditioning

A new method has been developed to accelerate the convergence of explicit time-marching, laminar, Navier-Stokes codes through the combination of local preconditioning and multi-stage time marching optimization. Local preconditioning is a technique to modify the time-dependent equations so that all information moves or decays at nearly the same rate, thus relieving the stiffness for a system of equations. Multi-stage time marching can be optimized by modifying its coefficients to account for the presence of viscous terms, allowing larger time steps. We show it is possible to optimize the time marching scheme for a wide range of cell Reynolds numbers for the scalar advection-diffusion equation, and local preconditioning allows this optimization to be applied to the Navier-Stokes equations. Convergence acceleration of the new method is demonstrated through numerical experiments with circular advection and laminar boundary-layer flow over a flat plate.

Kleb, William L.↗

NASA NTRS1999

Efficient Multi-Stage Time Marching for Viscous Flows via Local Preconditioning

Kleb, William L.↗

DOE OSTI2024

Modeling plasticity-mediated void growth at the single crystal scale: A physics-informed machine learning approach

Modeling the evolution of voids during plastic flow as well as their effects on plastic dissipation is critical for both component manufacturing and lifetime estimation purposes. To this end, we propose a rate-dependent constitutive model to homogenize the effects of semi-randomly distributed voids on single crystal plasticity whilst capturing void interaction and plastic anisotropy. Here, this present work focuses on the case of face centered cubic crystals to introduce an anisotropic gauge function applicable within the crystal plasticity formalism. The approach combines analytical methods to describe the micromechanics of the system in combination with symbolic regression to capture analytically intractable mechanisms from data. The hybrid framework uses a physics-informed genetic programming-based symbolic regression algorithm to solve a multiform optimization problem simultaneously producing a new gauge function and a new strain rate equation. This is also a multi-objective optimization problem with many competing objectives. A new search and selection step is introduced to the genetic algorithm that promotes convergence toward a global solution that better satisfies all the objectives. Overall, the symbolic equations produced leverage data-driven methods to achieve greater accuracy than comparable alternatives on an analytically intractable problem while maintaining model transparency.

36 MATERIALS SCIENCE↗

DOE OSTI2020

Multiresolution GPC-Structured Control of a Single-Loop Cold-Flow Chemical Looping Testbed

Chemical looping is a near-zero emission process for generating power from coal. It is based on a multi-phase gas-solid flow and has extremely challenging nonlinear, multi-scale dynamics with jumps, producing large dynamic model uncertainty, which renders traditional robust control techniques, such as linear parameter varying H ∞ design, largely inapplicable. This process complexity is addressed in the present work through the temporal and the spatiotemporal multiresolution modeling along with the corresponding model-based control laws. Namely, the nonlinear autoregressive with exogenous input model structure, nonlinear in the wavelet basis, but linear in parameters, is used to identify the dominant temporal chemical looping process dynamics. The control inputs and the wavelet model parameters are calculated by optimizing a quadratic cost function using a gradient descent method. The respective identification and tracking error convergence of the proposed self-tuning identification and control schemes, the latter using the unconstrained generalized predictive control structure, is separately ascertained through the Lyapunov stability theorem. The rate constraint on the control signal in the temporal control law is then imposed and the control topology is augmented by an additional control loop with self-tuning deadbeat controller which uses the spatiotemporal wavelet riser dynamics representation. The novelty of this work is three-fold: (1) developing the self-tuning controller design methodology that consists in embedding the real-time tunable temporal highly nonlinear, but linearly parametrizable, multiresolution system representations into the classical rate-constrained generalized predictive quadratic optimal control structure, (2) augmenting the temporal multiresolution loop by a more complex spatiotemporal multiresolution self-tuning deadbeat control loop, and (3) demonstrating the effectiveness of the proposed methodology in producing fast recursive real-time algorithms for controlling highly uncertain nonlinear multiscale processes. The latter is shown through the data from the implemented temporal and augmented spatiotemporal solutions of a difficult chemical looping cold flow tracking control problem.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

NASA NTRS1993

Global positioning system reoccupation of early triangulation sites - Tectonic deformation of the Southern Coast Ranges

We study tectonic deformation in the Southern Coast Range, California. We use triangulation and astronomic azimuth data collected since 1875, trilateration since 1970, and global positioning system data collected from 1986 to 1987. Two modeling techniques have been used. An elastic block-fault model is applied to study the tectonic motion of the San Andreas Fault and the San Gregorio-Hosgri Fault. Station velocities are modeled to study regional deformations. Results show that the regional deformation is predominantly controlled by deep strike-slip motion along the San Andreas Fault, at a rate of 33 +/- 2 mm/yr. Deep slip along the San Gregorio-Hosgri Fault is about 0-4 mm/yr, assuming a locked suit to a depth of 20 km. Convergence normal to the San Andreas Fault in the Southern Coast Ranges is not greater than 0.02 microrad/yr.

Shen, Zheng-Kang↗

NASA NTRS1998

Transient Simulation of Pressure Oscillations in the Fuel Feedline of the Fastrac Engine Thrust Chamber

During mainstage testing of the 60,000 lbf thrust Fastrac thrust chamber at MSFC's Test Stand 116 (TS 116), sustained, large amplitude oscillations near 530 Hz were observed in the pressure data. These oscillations were detected both in the RP-1 feedline, downstream of the cavitating venturi, and in the combustion chamber. The driver of the instability is believed to be feedline excitation driven by either periodic cavity collapse at the exit of the cavitating venturi or combustion instability. In covitating venturi, static pressure drops as the flow passes through a constriction resembling a converging-diverging nozzle until the vapor pressure is reached. At the venturi throat, the flow is essentially choked, which is why these devices are typically used for mass flow rate control and disturbance isolation. Typically, a total pressure drop of 15% or more across the venturi is required for cavitation. For much larger pressure differentials, unstable cavities can form and subsequently collapse downstream of the throat. Although the disturbances generated by cavitating venturis is generally considered to be broad-band, this type of phenomena could generate periodic behavior capable of exciting the feedline. An excitation brought about by combustion instability would result from the coupling of a combustion chamber acoustic mode and a feedline resonance frequency. This type of coupling is referred to as "buzz" and is not uncommon for engines in this thrust range.

Bullard, Brad↗

NASA NTRS2005

Modelling the Neutral Atmosphere and Plasma Environment of Saturn

The first year of this contract has resulted in two publications with the P.I. and co-I Jurac as lead authors and two publications where these team members are co-authors. These papers discuss modeling work undertaken in preparation for Cassini; the goal was to summarize our current best knowledge of the ion and neutrals sources and distributions. One of the major goals of this project is to improve models of the plasma and neutral environment near Saturn. The paper "A self-consistent model of plasma and neutrals at Saturn: Neutral cloud morphology" [Jurac and Richardson, 20051 presents results on the neutral clouds near Saturn using a model which for the first times treats the ions and neutrals self-consistently. We also for the first time include a directly sputtered H source. The Voyager and HST observations are used as model constraints. The neutral source is adjusted to give a good match to the HST observations of OH. For this initial run the ion parameters from Richardson et al. are used; charge exchange with ions is a major neutral loss process. The neutral profile derived from the model is then used in a model of plasma transport and chemistry (with the plasma diffusion rate the only free parameter). This model gives new values of the ion composition which are then fed back into the neutral model. This iteration continues until the values converge.

Richardson, John D.↗

NASA NTRSDate not supplied

Modeling of Cavity Residence Time via Modeled Lagrangian Particle Tracking

One of the critical components in a scramjet is the cavity flame holder. In evaluating the sizing of this component, the residence time of the fuel in the cavity is a common parameter used. However, to obtain the residence time via CFD is an expensive undertaking. One method is to perform an LES simulation of the cavity and use Lagrangian particle tracking to obtain both the mean and probability density function of the time that the particles remain in the cavity. While this method affords a significant amount of information, it is very costly in CPU run time. A simpler method involves performing a time accurate RANS simulation of the cavity, and when converged, overlay a traceable scalar in the cavity region, and then record the decay of this scalar’s flux at a downstream location. The time scale of the decay rate can then be estimated and a residence time inferred. This method is less computationally expensive than the LES approach, but still requires a time-accurate solution while only providing a mean cavity residence time evaluation. In this paper a method of determining the mean cavity residence time as well as the Probability Density Function (PDF) of the residence time is proposed and demonstrated. This method is based on Lagrangian particle tracking modeled by the Generalized Langevin Equation and is implemented as a post-processing step for a steady RANS solution. Computational cost is on the order of minutes and the results are in agreement with values obtained by the LES and scalar decay methods.

Langevin↗

NASA NTRS2017

Improvement of Automated POST Case Success Rate Using Support Vector Machines

During early conceptual design of complex systems, concept down selection can have a large impact upon program life-cycle cost. Therefore, any concepts selected during early design will inherently commit program costs and affect the overall probability of program success. For this reason it is important to consider as large a design space as possible in order to better inform the down selection process. For conceptual design of launch vehicles, trajectory analysis and optimization often presents the largest obstacle to evaluating large trade spaces. This is due to the sensitivity of the trajectory discipline to changes in all other aspects of the vehicle design. Small deltas in the performance of other subsystems can result in relatively large fluctuations in the ascent trajectory because the solution space is non-linear and multi-modal [1]. In order to help capture large design spaces for new launch vehicles, the authors have performed previous work seeking to automate the execution of the industry standard tool, Program to Optimize Simulated Trajectories (POST). This work initially focused on implementation of analyst heuristics to enable closure of cases in an automated fashion, with the goal of applying the concepts of design of experiments (DOE) and surrogate modeling to enable near instantaneous throughput of vehicle cases [2]. Additional work was then completed to improve the DOE process by utilizing a graph theory based approach to connect similar design points [3]. The conclusion of the previous work illustrated the utility of the graph theory approach for completing a DOE through POST. However, this approach was still dependent upon the use of random repetitions to generate seed points for the graph. As noted in [3], only 8% of these random repetitions resulted in converged trajectories. This ultimately affects the ability of the random reps method to confidently approach the global optima for a given vehicle case in a reasonable amount of time. With only an 8% pass rate, tens or hundreds of thousands of reps may be needed to be confident that the best repetition is at least close to the global optima. However, typical design study time constraints require that fewer repetitions be attempted, sometimes resulting in seed points that have only a handful of successful completions. If a small number of successful repetitions are used to generate a seed point, the graph method may inherit some inaccuracies as it chains DOE cases from the non-global-optimal seed points. This creates inherent noise in the graph data, which can limit the accuracy of the resulting surrogate models. For this reason, the goal of this work is to improve the seed point generation method and ultimately the accuracy of the resulting POST surrogate model. The work focuses on increasing the case pass rate for seed point generation.

Zwack, Matthew R.↗

DOE OSTI2022

Adaptive elasticity policies for staging-based in situ visualization

In situ processing aims to alleviate the growing gap between computation and I/O capabilities by performing data processing close to the data source. In situ processing is widely used to process data generated by multiple data sources, including observation data from edge devices or scientific observational facilities and the simulation data generated by scientific computation on a high-performance computing (HPC) platform. For a scientific workflow that is run on an HPC platform and composed of a simulation program and an in situ data analytics or visualization (abbreviated as ana/vis) task, there is an implicit assumption that the computing resources assigned to the workflow keep static during the workflow execution. However, with the converging trend between the HPC and cloud computing platform, running the in situ ana/vis task in an elastic way is promising to decrease its overhead and improve its resource utilization rate. Resource elasticity represents the ability to change resource configurations such as the number of computing nodes/processes during workflow execution. An elastic job may dynamically adjust resource configurations; it may use a few resources at the beginning and more resources toward the end of the job when interesting data appear. However, it is hard to predict a priori how many computing nodes/processes need to be added/removed during the workflow execution to adapt to changing workflow needs. How to efficiently guide elasticity operations, such as growing or shrinking the number of processes used for in situ analysis during workflow execution, is an open-ended research question. In this article, we present adaptive elasticity policies that adopt workflow runtime information collected during workflow execution to predict how to trigger the addition/removal of processes in order to minimize in situ processing overhead. Taking in situ visualization tasks as an example, we integrate the presented elasticity policies into a staging-based elastic workflow and evaluate its efficiency in multiple elasticity scenarios. Compared with the situation without elasticity or with a static elasticity policy that uses a fixed number of processes for each rescaling operation, the adaptive elasticity policy can save overhead in finding a proper resource configuration and improve resource utilization efficiency. Furthermore, one experiment illustrates that the adaptive elasticity policy saves 41% of core-hours compared with the situation without the resource elasticity.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2020

Methodology for inferring reactor core power distribution from an optical fiber based gamma thermometer array

Here, a 4-step data analytic methodology has been devised for the purpose of inferring the distribution of power in a reactor core, based on the response of an array of optical fiber based gamma thermometers (OFBGTs). This data analytic methodology is crucial for the development of a system of OFBGTs for the purpose of calibrating local power range monitors in boiling water reactors. Such a system would be an improvement to the present calibration system in boiling water reactors, in terms of safety, efficiency, and permanence. The first step of this methodology is to establish an energy balance method. In this method, one uses MCNP to determine response functions, which allow one to convert from gamma thermometer response to power. The gamma thermometers and the reactor core are segmented, such that each gamma thermometer segment provides an estimate of power for each reactor core segment. The estimates of power for the reactor core segments (hereafter referred to as fuel assembly segments) are calculated based on the measured response of the gamma thermometers. The second step of the methodology is to employ a weighting scheme to combine the estimates of the power of the various fuel assembly segments, based on the response of the various OFGBT segments and the incremental dose rates of the various fuel assembly segments to the OFBGT segments. The third step is to iteratively calculate the estimates of the power of the fuel assembly segments, until a convergence criterion is met, which indicates that the calculation has converged. The fourth and final step of the methodology is to estimate the uncertainty of the power for each fuel assembly segment. The mathematical basis for this step is not the focus of this paper. We have used a 3-D homogeneous reactor model to demonstrate the data analytic methodology; and have found that the data analytic methodology operates as intended.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

DOE OSTI2024

Adaptive laboratory evolution and genetic engineering improved terephthalate utilization in Pseudomonas putida KT2440

Poly(ethylene terephthalate) (PET) is one of the most ubiquitous plastics and can be depolymerized through biological and chemo-catalytic routes to its constituent monomers, terephthalic acid (TPA) and ethylene glycol (EG). TPA and EG can be re-synthesized into PET for closed-loop recycling or microbially converted into higher-value products for open-loop recycling. Here, in this study, we expand on our previous efforts engineering and applying Pseudomonas putida KT2440 for PET conversion by employing adaptive laboratory evolution (ALE) to improve TPA catabolism. Three P. putida strains with varying degrees of metabolic engineering for EG catabolism underwent an automation-enabled ALE campaign on TPA, a TPA and EG mixture, and glucose as a control. ALE increased the growth rate on TPA and TPA-EG mixtures by 4.1- and 3.5-fold, respectively, in approximately 350 generations. Evolved isolates were collected at the midpoints and endpoints of 39 independent ALE experiments, and growth rates were increased by 0.15 and 0.20 h -1 on TPA and a TPA-EG, respectively, in the best performing isolates. Whole-genome re-sequencing identified multiple converged mutations, including loss-of-function mutations to global regulators gacS, gacA, and turA along with large duplication and intergenic deletion events that impacted the heterologously-expressed tphAB II catabolic genes. Reverse engineering of these targets confirmed causality, and a strain with all three regulators deleted and second copies of tphAB II and tpaK displayed improved TPA utilization compared to the base strain. Taken together, an iterative strain engineering process involving heterologous pathway engineering, ALE, whole genome sequencing, and genome editing identified five genetic interventions that improve P. putida growth on TPA, aimed at developing enhanced whole-cell biocatalysts for PET upcycling.

36 MATERIALS SCIENCE↗

NASA NTRS1981

Software models: A Bayesian approach to parameter estimation in the Jelenski-Moranda software reliability model

Maximum likelihood estimation procedures for the Jelinski-Moranda software reliability model often give misleading answers. A reparameterization and a Bayesian analysis eliminate some of the problems incurred by MLE methods and often give better predictions on sets of real and simulated data. Practical difficulties in estimating the initial number of errors N and the failure rate of each error phi by the method of maximum likelihood are: N, the MLE of N, is occasionally infinite (i.e., the routines for calculating N and phi do not converge). It is shown that N is finite sub i only if the regression line of the interevent times t sub i vs. i has positive slope. A serious problem is that often N approximates n, the sample size, and sometimes N = n. Thus the MLE predicts that the program is perfect even when it is far from being so. Only when almost all failures have been removed can N and phi be trusted near the end of debugging.

Littlewood, B.↗

NASA NTRS1988

Investigation of supersonic chemically reacting and radiating channel flow

The 2-D time-dependent Navier-Stokes equations are used to investigate supersonic flows undergoing finite rate chemical reaction and radiation interaction for a hydrogen-air system. The explicit multistage finite volume technique of Jameson is used to advance the governing equations in time until convergence is achieved. The chemistry source term in the species equation is treated implicitly to alleviate the stiffness associated with fast reactions. The multidimensional radiative transfer equations for a nongray model are provided for a general configuration and then reduced for a planar geometry. Both pseudo-gray and nongray models are used to represent the absorption-emission characteristics of the participating species. The supersonic inviscid and viscous, nonreacting flows are solved by employing the finite volume technique of Jameson and the unsplit finite difference scheme of MacCormack. The specified problem considered is of the flow in a channel with a 10 deg compression-expansion ramp. The calculated results are compared with those of an upwind scheme. The problem of chemically reacting and radiating flows are solved for the flow of premixed hydrogen-air through a channel with parallel boundaries, and a channel with a compression corner. Results obtained for specific conditions indicate that the radiative interaction can have a significant influence on the entire flow field.

Mani, Mortaza↗

NASA NTRS1989

Extension during continental convergence, with application to the Tibetan Plateau

An explanation in terms of the thermal evolution of thickened continental lithosphere is offered to explain the transition, in the late Tertiary to Quaternary, from north-south compression to east-west extension in the strain rate field of the Tibetan Plateau. The lower part of the lithosphere consists of a thermal boundary layer which, when thickened by horizontal shortening, is colder and denser than its surroundings. Convective instability of the thickened thermal boundary layer and its replacement by hot asthenosphere would rapidly raise the surface elevation and gravitational potential energy of the overlying part of the lithosphere. The convective instability would happen in a time brief compared with the collision time scale but would only occur after there had already been substantial thickening of the lithosphere. The increase in surface height and of potential energy are sufficient for east-west extension to replace north-south compression as the dominant feature of the stress field.

England, Philip↗

DOE OSTI2023

Unraveling Anisotropic and Pulsating Etching of ZnO Nanorods in Hydrochloric Acid via Correlative Electron Microscopy

Despite much technical progress achieved so far, the exact surface and shape evolution during wet chemical etching is less unraveled, especially in ionically bonded ceramics. Herein, by using in situ liquid cell transmission electron microscopy, a repeated two-stage anisotropic and pulsating periodic etching dynamic is discovered during the pencil shape evolution of a single crystal ZnO nanorod in aqueous hydrochloric acid. Specifically, the nanopencil tip shrinks at a slower rate along [0001̅] than that along the $\langle$101̅0$\rangle$ directions, resulting in a sharper ZnO pencil tip. Afterward, rapid tip dissolution happens due to accelerated etching rates along various crystal directions. Concurrently, the vicinal base region of the original nanopencil tip emerges as a new tip followed by the repeated sequence of tip shrinking and removal. The high-index surfaces, such as {101̅m} (m = 0, 1, 2, or 3) and {21̅1̅n} (n = 0, 1, 2, or 3), are found to preferentially expose in different ratios. Our 3D electron tomography, convergent beam electron diffraction, middle-angle bright-field STEM, and XPS results indicate the dissociative Cl – species were bound to the Zn-terminated tip surfaces. Furthermore, DFT calculation suggests the preferential Cl – passivation over the {101̅1} and (0001) surfaces of lower energy than others, leading to preferential surface exposures and the oscillatory variation of different facet etching rates. The boosted reactivity due to high-index nanoscale surface exposures is confirmed by comparatively enhanced chemical sensing and CO 2 hydrogenation activity. In conclusion, these findings provide an in-depth understanding of anisotropic wet chemical etching of ionic nanocrystals and offer a design strategy for advanced functional materials.

36 MATERIALS SCIENCE↗