Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “numerical optimal control”

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

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

At least 1,063 records · Page 59
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Incentivizing Cooperative Merging Control: Insights from Multi-Agent Deep Reinforcement Learning

Cooperative driving automation enables connected and automated vehicles (CAVs) to devise cooperative merging control, introducing great potentials to alleviate traffic congestion, reduce energy consumption, and enhance safety for highway on-ramp operations. Although numerous CAV cooperative merging algorithms have been developed to improve energy and traffic performance, the agreement-seeking among CAV users and their local benefits have been understudied. This can lead to rejections of cooperative merging plans and jeopardizing CAV performance, as a cooperation may entail certain CAVs to sacrifice their local benefits to achieve a system optimum. To address this issue, the study first leverages multi-agent deep reinforcement learning (MADRL) factoring both local reward and regional reward to demonstrate the discrepancies between CAV users’ local benefits and system optimum. Next, the existence of a correlated equilibrium is proved to characterize the convergence of MADRL training. This further facilitates the incorporation of incentives (computed based on reward discrepancies) to compensate for CAV users’ local benefits and facilitate system-optimal agreements in cooperative merging operations.

Zhou, Anye [ORNL] (ORCID:0000000301455579)↗
DOE OSTI2023

Driving Next-Generation Workflows from the Data Plane

We observe the emergence of a new generation of scientific workflows that process data produced at a sustained rate by scientific instruments and large scale numerical simulations. This data is consumed by multiple analysis, visualization, or Machine Learning components not only to enable inference and justify the scientific program, but also to monitor and steer the evolution of these experiments. In such workflows, moving intermediate data efficiently is key to performance, more than efficiently scheduling computational tasks. However, most traditional workflow management systems focus on optimizing task scheduling and then deal with data management, assuming a “move little, compute for long” model, which makes them unfit to the efficient management of this new generation of workflows. Therefore, we advocate for a new way to manage scientific workflows. We propose to consider an efficiently and independently managed data plane that can store and stream data. Workflows compute components, in the application plane can then interact with the data plane, abstracted from complexities of data management. Then, the role of a workflow management system would become that of a control plane that allows users to connect services together to execute the workflow and manages connections between the application and data planes. In this position paper, we characterize several next-generation workflow motifs and describe how their interaction with the data plane is a challenge to traditional workflow management systems. Then, we express a set of requirements that a workflow management system should meet to efficiently manage next-generation workflows at different scales. Based on these requirements, we expose our vision of driving next-generation workflows from the data plane and list remaining open challenges.

Suter, Fred↗
DOE OSTI2021

The NREL Sensor Laboratory: Status and Future Directions for Hydrogen Detection

The NREL Hydrogen Sensor Laboratory was commissioned in 2010 as a resource for the national and international hydrogen community to ensure the availability and proper use of hydrogen sensors. Since then, the Sensor Laboratory has provided unbiased verification of hydrogen sensor performance for sensor developers, end-users, and regulatory agencies and has also provided active support for numerous code and standards development organizations. Although sensor performance assessment remains a core capability, the mission of the NREL Sensor Laboratory has expanded toward a more holistic approach regarding the role of hydrogen detection and its implementation strategy for both assurance of facility safety as well as for process control applications. Active monitoring for detection of unintended releases has been identified by the NFPA 2 Storage Task Group as a viable approach for improving facility safety and lowering setbacks. The current research program for the Sensor Laboratory addresses both conventional and advanced developing detection strategies in response to the emerging large-scale markets, such as those envisioned by H2@Scale. These emerging hydrogen applications may require alternative detection strategies that supplement and may ultimately supplant the use of traditional sensors for monitoring hydrogen releases. Research focus areas for the NREL Sensor Laboratory now encompasses the characterization of released hydrogen behavior to optimize detection strategies for both indoor and outdoor applications, advanced methods of hydrogen leak detection such as hydrogen wide area monitoring approaches for large scale and distributed applications, the incorporation of active monitoring as a risk reduction strategy to improve safety at hydrogen facilities, and continuing support of hydrogen safety codes and standards. In addition to assurance of safety, detection will be critical for process control applications, such as hydrogen fuel quality verification for fuel cell vehicle applications and for monitoring and controlling of hydrogen-natural gas blend composition.

codes and standards↗
NASA NTRS2026

The Future of in-Situ Sequencing-Based Microbial Monitoring: Development of a Shelf-Stable Method for Artemis and Beyond

Microbial monitoring onboard the International Space Station (ISS) is essential for assessing the efficiency of the Environmental Control and Life Support Systems (ECLSS) and providing insight into potential risk to both crew and spacecraft. Historically, this monitoring required the need to culture organisms onboard, return these cultures to Earth, and then complete the identifications, a process that would take months. Over the past decade, and through numerous payloads, advances in molecular biology have enabled in-flight microbial identifications using nanopore sequencing. The swab-to-sequencer method resulting from these efforts was transitioned from research to operations for microbial monitoring under the Crew Health Care Systems (CHeCS) BioMole. Collectively, these accomplishments have propelled the swab-to-sequencer method to be selected as the Microbial Surface Monitor (MSM) for Gateway, as well as a payload on Artemis IV. However, the lack of cold stowage availability for Artemis requires modifications to the entire method due to the thermal instability of the reagents required for sample preparation. To achieve this, new development, optimization, and validations were undertaken. Key considerations included enzyme concentration, buffer compatibility, and equal or enhanced sensitivity and specificity. At each step, thorough side-by-side comparisons with the current ISS method were performed. The development of a robust shelf-stable method will ensure continued sequencing-based microbial monitoring for Artemis and beyond, providing data in near real-time, enhancing risk response time, and yielding clear insight into the microbiome of spacecraft.

Christian G Mena↗
DOE OSTI2021

Atomic Layer Deposition for Materials-Based H2 Storage: Opportunities and Limitations

The transportation demands in our growing hydrogen (H2) economy requires robust storage systems. The current commercially implemented technology in fuel cell cars relies on the well-established technology of hydrogen gas compressed to 350-700 bar, depending on the application. The compressed gas tanks in use today are bulky and cost intensive. To address this challenge, material-based storage is one of the long-term alternatives considered and constitutes the focus of this talk. Material-based storage is broadly defined as hydrogen bound to solid materials, with its binding strength varying from physisorption to porous materials, such as zeolites and metal organic frameworks, to chemisorption in (complex) metal hydrides. The ultimate targets set by the U.S. Department of Energy for this technology include a system gravimetric capacity of 6.5 wt% and volumetric capacity of 40 g/L at 100 bar, operating temperatures ranging between -40 C and +40 C and adsorption/desorption timescales of < 5 min. Storage in the form of physisorbed or chemisorbed hydrogen has guided the materials research, which metal- organic framework and (complex) metal hydrides being the most promising materials classes. A variety of these materials have met one or more of the targets, but it has remained elusive for a single material system to meet all these stringent requirements. Nano-encapsulation and low-concentration chemical additives have previously been employed separately to overcome such challenges. Functionalization via atomic layer deposition (ALD), however, offers unique characteristics that make it suitable for both, nano-encapsulation and "doping" with low-concentration additives. This deposition technique has sub-monolayer thickness control, is highly conformal in high-surface area materials and is self-limiting, i.e., once the gas-phase precursor reacts with the available surface sites, the surface reactions stop. In this presentation, we will showcase examples where ALD, more generally vapor-phase functionalization, on (complex) metal hydrides and organic frameworks has improved the material properties for H2 storage. In our first study, Al2O3 was deposited on magnesium borohydride, Mg(BH4)2, at room temperature using trimethylaluminum (TMA) and water. From our findings, encouraging initial results were obtained: the H2 desorption temperature was lowered by 60-120 degrees C, the desorbed gravimetric H2 capacity at temperatures < 250 degrees C was doubled, and the desorption kinetics increased by a factor of ~6 compared to uncoated Mg(BH4)2. However, hydrolysis reactions caused by residual surface -OH groups from the ALD water-pulse degraded the sample substantially. Through this study, the use of TMA was observed to be highly reactive with the Mg(BH4)2 surface species, leading to the strategy of exposing the sample only to TMA. The relative mole fraction of the vapor-phase additive can be precisely controlled with the number of TMA pulses and pulse duration. By tuning these parameters, we show in our second study that 10 pulses of TMA at ambient conditions were able to decrease the H2 desorption temperature by ~100 degrees C while retaining <95 % of its H2 capacity. We applied this approach with other additives such as BBr3, TiCl4 and tetrahydrofuran, and demonstrate that this unique approach opens the door to a new class of molecular additives and catalysts which cannot easily be introduced with conventional mechano-chemical or solvent-based techniques for (complex) metal-hydrides. In the case of metal organic frameworks, ALD is a promising technique to functionalize the pores with metal atoms or functional groups able to tune the gas selectivity and binding energy, for which ~15 kJ/mol has been established as the optimal value for H2 storage in sorbent materials. Vapor-phase techniques such as ALD have numerous benefits over other functionalization tools for H2 storage materials opening innumerable opportunities to the field. To exploit these opportunities, the current limitations on room temperature and water-less ALD processes needs to be overcome which will greatly expand the possibilities of encapsulation and incorporation of additives for organic frameworks and (complex) metal hydrides.

atomic layer deposition↗
NASA NTRS1993

Spray combustion experiments and numerical predictions

The next generation of commercial aircraft will include turbofan engines with performance significantly better than those in the current fleet. Control of particulate and gaseous emissions will also be an integral part of the engine design criteria. These performance and emission requirements present a technical challenge for the combustor: control of the fuel and air mixing and control of the local stoichiometry will have to be maintained much more rigorously than with combustors in current production. A better understanding of the flow physics of liquid fuel spray combustion is necessary. This paper describes recent experiments on spray combustion where detailed measurements of the spray characteristics were made, including local drop-size distributions and velocities. Also, an advanced combustor CFD code has been under development and predictions from this code are compared with experimental results. Studies such as these will provide information to the advanced combustor designer on fuel spray quality and mixing effectiveness. Validation of new fast, robust, and efficient CFD codes will also enable the combustor designer to use them as additional design tools for optimization of combustor concepts for the next generation of aircraft engines.

Mularz, Edward J.↗
DOE OSTI2024

Batch Extraction Studies to Evaluate Trace Element Behavior in PUREX Conditions

The multilab Intentional Forensics Venture is working to identify which stable elements (i.e., taggants) at trace concentrations relative to U would persist throughout the nuclear fuel cycle in a voluntary fuel tagging scheme. A taggant would provide the nuclear forensics community with a “barcode” to help identify nuclear materials found outside of regulatory control. A portion of this project was focused on reprocessing effects and determining which, if any, elements would coextract with U(VI) in standard Pu–U reduction extraction (PUREX) conditions. Elements with a propensity to coextract could, in theory, be used as taggants from a PUREX perspective. Although retention is not a performance requirement, the taggant signature would need to partition predictably from the U stream after the PUREX process to maintain forensic utility. This report documents results from several batch extraction studies with numerous trace elements from HNO 3 (1.5–5 M), with and without U(VI), into 30% tri-n-butyl phosphate (TBP) in kerosene. Extraction and back-extraction tests were used to evaluate nearly 60 elements in surrogate conditions for PUREX, and distribution coefficients (i.e., D-values) for most species were <0.1, indicating few species are likely to co-extract with U through PUREX. Additional studies are needed to optimize sample volumes and dilutions to dial in these low D-values. The D-values (D) were determined for several of the more promising elements, including Re and Se. Ultimately, we conclude that only a limited number of the ~ 60 elements investigated are extractable in the U stream of PUREX, based on measured D values, meaning most candidate elemental taggants would likely be lost at this stage of the nuclear fuel cycle, even when considering a range of acid concentrations.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗
NASA NTRS1993

Spray Combustion Experiments and Numerical Predictions

The next generation of commercial aircraft will include turbofan engines with performances significantly better than those in the current fleet. Control of particulate and gaseous emissions will also be an integral part of the engine design criteria. These performance and emission requirements present a technical challenge for the combustor: control of the fuel and air mixing and control of the local stoichiometry will have to be maintained much more rigorously than combustors in current production. A better understanding of the flow physics of liquid fuel spray combustion is necessary. This presentation describes recent experiments on spray combustion where detailed measurements of the spray characteristics were made, including local drop-size distributions and velocities. In addition, an advanced combustor CFD code has been under development and predictions from this code are presented and compared with measurements. Studies such as these will provide information to the advanced combustor designer on fuel spray quality and mixing effectiveness. Validation of new fast, robust, and efficient CFD codes will also enable the combustor designer to use them as valuable design tools for optimization of combustor concepts for the next generation of aircraft engines.

Mularz, Edward J.↗
DOE OSTI2025

A dynamic 2D Borehole Thermal Energy Storage (BTES) model for enhanced computational efficiency

Progressing toward a future increasingly reliant on renewable energy sources, the development of effective, durable energy storage solutions becomes essential to balance supply and demand fluctuations. Borehole Thermal Energy Storage (BTES) is a long-duration thermal energy storage technology that captures excess heat generated from renewable energy sources and stores it underground for later use, enabling the efficient utilization of sustainable energy. This approach is particularly valuable in district energy networks when integrated with Ground Source Heat Pumps (GSHP) to provide stable heating and cooling. However, traditional three-dimensional (3D) numerical models of BTES systems demand extensive computational resources, limiting their practicality for real-time and large-scale applications. This study introduces a novel two-dimensional (2D) modeling approach that reduces computational costs while maintaining high accuracy. By employing a radial ring-based discretization method, the model simulates heat injection, retention, and retrieval dynamics over seasonal cycles. A new thermal-mass weighted-average temperature parameter is introduced to evaluate the performance of BTES systems. Model validation against FEFLOW simulations demonstrates a 17-fold improvement in computational speed compared to traditional Computational Fluid Dynamics (CFD) models while achieving a mean absolute percentage error (MAPE) of 2 % during charging and 4 % during discharging. Additionally, a trade-off analysis between computational efficiency and accuracy is conducted, ensuring the model's applicability for real-world scenarios. The findings of this research contribute to the development of computationally efficient BTES models, facilitating better optimization, control, and integration into renewable energy systems. This work provides a foundation for further studies in techno-economic analysis, multi-year performance evaluation, and real-time operational strategies for BTES applications, supporting a more sustainable energy future.

2D modeling↗
DOE OSTI2025

Tree tensor network hierarchical equations of motion based on time-dependent variational principle for efficient open quantum dynamics in structured thermal environments

In this work, we introduce an efficient method, TTN-HEOM, for exactly calculating the open quantum dynamics for driven quantum systems interacting with highly structured bosonic baths by combining the tree tensor network (TTN) decomposition scheme with the bexcitonic generalization of the numerically exact hierarchical equations of motion (HEOM). The method yields a series of quantum master equations for all core tensors in the TTN that efficiently and accurately capture the open quantum dynamics for non-Markovian environments to all orders in the system–bath interaction. These master equations are constructed based on the time-dependent Dirac–Frenkel variational principle, which isolates the optimal dynamics for the core tensors given the TTN ansatz. The dynamics converges to the HEOM when increasing the rank of the core tensors, a limit in which the TTN ansatz becomes exact. We introduce TENSO, tensor equations for non-Markovian structured open systems, as a general-purpose Python code to propagate the TTN-HEOM dynamics. We implement three general propagators for the coupled master equations: two fixed-rank methods that require a constant memory footprint during the dynamics and one adaptive-rank method with a variable memory footprint controlled by the target level of computational error. We exemplify the utility of these methods by simulating a two-level system coupled to a structured bath containing one Drude–Lorentz component and eight Brownian oscillators, which is beyond what can presently be computed using the standard HEOM. Our results show that the TTN-HEOM is capable of simulating both dephasing and relaxation dynamics of driven quantum systems interacting with structured baths, even those of chemical complexity, with an affordable computational cost.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2023

Operator Relaxation and the Optimal Depth of Classical Shadows

Classical shadows are a powerful method for learning many properties of quantum states in a sample-efficient manner, by making use of randomized measurements. Here we study the sample complexity of learning the expectation value of Pauli operators via “shallow shadows,” a recently proposed version of classical shadows in which the randomization step is effected by a local unitary circuit of variable depth t. Here we show that the shadow norm (the quantity controlling the sample complexity) is expressed in terms of properties of the Heisenberg time evolution of operators under the randomizing (“twirling”) circuit—namely the evolution of the weight distribution characterizing the number of sites on which an operator acts nontrivially. For spatially contiguous Pauli operators of weight k, this entails a competition between two processes: operator spreading (whereby the support of an operator grows over time, increasing its weight) and operator relaxation (whereby the bulk of the operator develops an equilibrium density of identity operators, decreasing its weight). From this simple picture we derive (i) an upper bound on the shadow norm which, for depth t~log⁡(k), guarantees an exponential gain in sample complexity over the t=0 protocol in any spatial dimension, and (ii) quantitative results in one dimension within a mean-field approximation, including a universal subleading correction to the optimal depth, found to be in excellent agreement with infinite matrix product state numerical simulations. Our Letter connects fundamental ideas in quantum many-body dynamics to applications in quantum information science, and paves the way to highly optimized protocols for learning different properties of quantum states.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2020

Underlying Fundamentals of Kalman Filtering for River Network Modeling

The grand challenge of producing hydrometeorological estimates every time and everywhere has motivated the fusion of sparse observations with dense numerical models, with a particular interest on discharge in river modeling. Ensemble methods are largely preferred as they enable the estimation of error properties, but at the expense of computational load and generally with underestimations. These imperfect stochastic estimates motivate the use of correction methods, that is, error localization and inflation, although the physical justifications for their optimality are limited. The purpose of this study is to use one of the simplest forms of data assimilation when applied to river modeling and reveal the underlying mechanisms impacting its performance. Our framework based on assimilating daily averaged in situ discharge measurements to correct daily averaged runoff was tested over a 4-yr case study of two rivers in Texas. Results show that under optimal conditions of inflation and localization, discharge simulations are consistently improved such that the mean values of Nash–Sutcliffe efficiency are enhanced from211.32 to 0.55 at observed gauges and from212.24 to21.10 at validation gauges. Yet, parameters controlling the inflation and the localization have a large impact on the performance. Further investigations of these sensitivities showed that optimal inflation occurs when compensating exactly for discrepancies in the magnitude of errors while optimal localization matches the distance traveled during one assimilation window. These results may be applicable to more advanced data assimilation methods as well as for larger applications motivated by upcoming river-observing satellite missions, such as NASA’s Surface Water and Ocean Topography mission.

Streamflow↗
NASA NTRS1993

Spray combustion experiments and numerical predictions

The next generation of commercial aircraft will include turbofan engines with performance levels significantly better than those in the current fleet. Control of particulate and gaseous emissions will also be an integral part of the engine design criteria. These performance and emission requirements present a technical challenge for the combustor: control of the fuel and air mixing and control of the local stoichiometry will have to be maintained much more rigorously than with combustors in current production. A better understanding of the flow physics of liquid fuel spray combustion is necessary. This paper describes recent experiments on spray combustion where detailed measurements of the spray characteristics were made, including local drop-size distributions and velocities. Also, an advanced combustor CFD code has been under development and predictions from this code are compared with experimental results. Studies such as these will provide information to the advanced combustor designer on fuel spray quality and mixing effectiveness. Validation of new fast, robust, and efficient CFD codes will also enable the combustor designer to use than as valuable additional design tools for optimization of combustor concepts for the next generation of aircraft engines.

Mularz, Edward J.↗
NASA NTRS1993

Next generation initiation techniques

Four-dimensional data assimilation strategies can generally be classified as either current or next generation, depending upon whether they are used operationally or not. Current-generation data-assimilation techniques are those that are presently used routinely in operational-forecasting or research applications. They can be classified into the following categories: intermittent assimilation, Newtonian relaxation, and physical initialization. It should be noted that these techniques are the subject of continued research, and their improvement will parallel the development of next generation techniques described by the other speakers. Next generation assimilation techniques are those that are under development but are not yet used operationally. Most of these procedures are derived from control theory or variational methods and primarily represent continuous assimilation approaches, in which the data and model dynamics are 'fitted' to each other in an optimal way. Another 'next generation' category is the initialization of convective-scale models. Intermittent assimilation systems use an objective analysis to combine all observations within a time window that is centered on the analysis time. Continuous first-generation assimilation systems are usually based on the Newtonian-relaxation or 'nudging' techniques. Physical initialization procedures generally involve the use of standard or nonstandard data to force some physical process in the model during an assimilation period. Under the topic of next-generation assimilation techniques, variational approaches are currently being actively developed. Variational approaches seek to minimize a cost or penalty function which measures a model's fit to observations, background fields and other imposed constraints. Alternatively, the Kalman filter technique, which is also under investigation as a data assimilation procedure for numerical weather prediction, can yield acceptable initial conditions for mesoscale models. The third kind of next-generation technique involves strategies to initialize convective scale (non-hydrostatic) models.

Warner, Tom↗
NASA NTRS1996

Application of pneumatic lift and control surface technology to advanced transport aircraft

The application of pneumatic (blown) aerodynamic technology to both the lifting and the control surfaces of advanced transport aircraft can provide revolutionary changes in the performance and operation of these vehicles, ranging in speed regime from Advanced Subsonic Transports to the High Speed Civil Transport, and beyond. This technology, much of it based on the Circulation Control Wing blown concepts, can provide aerodynamic force augmentations of 80 to 100 (i.e., return of 80-100 pounds of force per pound of input momentum from the blowing jet). This can be achieved without use of external mechanical surfaces. Clever application of this technology can provide no-moving-part lifting surfaces (wings/tails) integrated into the control system to greatly simplify aircraft designs while improving their aerodynamic performance. Lift/drag ratio may be pneumatically tailored to fit the current phase of the flight, and takeoff/landing performance can be greatly improved by reducing ground roll distances and liftoff/touchdown speeds. Alternatively, great increases in liftoff weights and payloads are possible, as are great reductions in wing and tail planform size, resulting in optimized cruise wing designs. Furthermore, lift generation independent of angle of attack provides much promise for increased safety of flight in the severe updrafts/downdrafts of microbursts and windshears, which is further augmented by the ability to sustain flight at greatly reduced airspeeds. Load-tailored blown wings can also reduce tip vorticity during highlift operations and the resulting vortex wake hazards near terminal areas. Reduced noise may also be possible as these jets can be made to operate at low pressures. The planned presentation will support the above statements through discussions of recent experimental and numerical (CFD) research and development of these advanced blown aerodynamic surfaces, portions of which have been conducted for NASA. Also to be presented will be predicted performance of advanced transports resulting from these devices. Suggestions will be presented for additional innovative high-payoff research leading to further confirmation of these concepts and their application to advanced efficient commercial transport aircraft.

Englar, Robert J.↗
NASA NTRS1992

A generalized reusable guidance algorithm for optimal aerobraking

A practical real-time guidance algorithm was developed for guiding aerobraking vehicles in such a way that the maximum heating rate, the maximum structural loads, and the post-aeropass delta-V requirements (for post-aeropass orbit insertion) are all minimized. The algorithm is general and reusable in the sense that a minimum of assumptions are made, thus minimizing the number of gains and mission-dependent parameters that must be laboriously determined prior to a particular mission. A particularly interesting feature is that inplane guidance performance is tuned by simply adjusting one mission-dependent parameter, the bank margin; similarly, the out-of-plane guidance performance is turned by simply adjusting a plane controller time constant. Other objectives in the algorithm development are simplicity, efficiency, and ease of use. The algorithm is developed for, but not necessarily restricted to, a single pass mission and a trimmed vehicle with a bank angle modulation as the method of trajectory control. Guidance performance is demonstrated via results obtained using this algorithm integrated into an aerobraking test-bed program. Comparisons are made with numerical results from a version of the aerobraking guidance algorithm that was to be flown onboard NASA's aeroassist flight experiment (AFE) vehicle. Promising results are obtained with a minimum of development effort.

Dukeman, G. A.↗
DOE OSTI2025

Degrees of rate control and AutoDiff-driven direct sensitivity analysis in heterogeneous catalysis

Despite the wide application and benefits of the degree of rate control (DRC) analysis, several details remain argued, particularly about the conservation of DRCs at transient (TR) and steady-state (SS) conditions, especially for complex reaction networks. This work argues that previous proofs about the conservation properties of DRCs have been incomplete, and we provide new mathematical proofs at TR and SS conditions. In addition, we use both analytical (automatic differentiation) and numerical (finite difference) approaches to compute DRCs for the case study of ethane hydrogenolysis (EH) over Pt(111). This work confirms that at both TR and SS conditions, the sum of all DRCs, i.e., sum of the degrees of kinetic (DKRC) and thermodynamic rate control (DTRC), is conserved at zero. At SS conditions, the sum of DKRC is conserved at 1 while the sum of DTRC is conserved at −1. In corroboration of previous works, we show that the DTRC for any adsorbate at SS is equal to the product of the species coverage and a constant. In contrast, at TR conditions, the individual sums of both DTRC and DKRC are not conserved and can be any real number, with potential implications for the novel field of dynamic catalysis. Finally, we show that the conventional finite difference (FD) approach, only useful at SS, is prone to inaccuracy and very sensitive to the value of the differential change applied. The optimal differential value also varies significantly with system and rate definition. Consequently, we describe and illustrate in this work the application of the automatic differentiation (AD) approach for the more accurate determination of DRCs at both TR and SS conditions.

Automatic differentiation↗
NASA NTRS1991

Comparison of ice-phase microphysical parameterization schemes using numerical simulations of tropical convection

The performance of several ice parameterizations has been evaluated through a numerical cloud model. Ice effects using different schemes are contrasted with each other and with an ice-free control by incorporating them into the cloud model and by applying them to simulations of tropical squall systems. The latter are simulated in 2D so that a large domain can be used to incorporate a complete anvil. Nonsquall-type convective lines are simulated in 3D owing to their smaller horizontal scale. It is concluded that inclusion of ice microphysics in the cloud model enhanced the agreement of the simulated convection with some features of observed convection, including the proportion of surface rainfall in the anvil region and the intensity and structure of the radar brightband near the melting level in the anvil. In the experiments with bulk microphysics, three ice categories produced much better results than two ice categories, which in turn was better than no ice. For the tropical squall-type and nonsquall-type systems the optimal mix was ice, snow, and graupel.

Mccumber, Michael↗