Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “local dynamic operator”

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 91 records · Page 5

Dynamics of Current-Sharing Within a REBCO Tape-Stack Cable

Current redistribution in a non-insulating REBCO tape-stack cable may enable the use of such a cable in a high field dipole without provision for transposition of the constituent tapes. As current is increased in the REBCO cable a dynamic resistance arises in the superconducting layer near operation at critical current. The critical current of each tape in the stack is determined by the local magnetic field variation within the cable. The superconducting-normal transition in REBCO occurs over a working range of current and, within that range, the longitudinal Ohmic electric field produces a transverse electric field between adjacent tapes with different critical currents. Circuit models indicate that, as a given tape approaches critical operation, current will naturally redistribute to neighboring tapes with higher current carrying capacity due to the dynamic rise in resistivity thus preventing premature quench. Here, a multi-scale model is being developed to study the dynamics of current-sharing, the limits of stability, and the impact of fluctuations in critical-current density along each tape.

43 PARTICLE ACCELERATORS↗

Interfaces in Heterogeneous Catalysts: Advancing Mechanistic Understanding through Atomic-Scale Measurements

Developing novel catalysts with high efficiency and selectivity is critical for enabling future clean energy conversion technologies. Interfaces in catalyst systems have long been considered the most critical factor in controlling catalytic reaction mechanisms. Interfaces include not only the catalyst surface but also interfaces within catalyst particles and those formed by constructing heterogeneous catalysts. The atomic and electronic structures of catalytic surfaces govern the kinetics of binding and release of reactant molecules from surface atoms. Interfaces within catalysts are introduced to enhance the intrinsic activity and stability of the catalyst by tuning the surface atomic and chemical structures. Examples include interfaces between the core and shell, twin or domain boundaries, or phase boundaries within single catalyst particles. In supported catalyst nanoparticles (NPs), the interface between the metallic NP and support serves as a critical tuning factor for enhancing catalytic activity. Surface electronic structure can be indirectly tuned and catalytically active sites can be increased through the use of supporting oxides. Additionally, tuning interfaces in catalyst systems has been identified as an important strategy in the design of novel catalysts. However, the governing principle of how interfaces contribute to catalyst behavior, especially in terms of interactions with intermediates and their stability during electrochemical operation, are largely unknown. This is mainly due to the evolving nature of such interfaces. Small changes in the structural and chemical configuration of these interfaces may result in altering the catalytic performance. These interfacial arrangements evolve continuously during synthesis, processing, use, and even static operation. A technique that can probe the local atomic and electronic interfacial structures with high precision while monitoring the dynamic interfacial behavior in situ is essential for elucidating the role of interfaces and providing deeper insight for fine-tuning and optimizing catalyst properties.Scanning transmission electron microscopy (STEM) has long been a primary characterization technique used for studying nanomaterials because of its exceptional imaging resolution and simultaneous chemical analysis. Over the past decade, advances in STEM, that is, the commercialization of both aberration correctors and monochromators, have significantly improved the spatial and energy resolution. Imaging atomic structures with subangstrom resolution and identifying chemical species with single-atom sensitivity are now routine for STEM. These advancements have greatly benefitted catalytic research. For example, the roles of lattice strain and surface elemental distribution and their effect on catalytic stability and reactivity have been well documented in bimetallic catalysts. In addition, three-dimensional atomic structures revealed by STEM tomography have been integrated in theoretical modeling for predictive catalyst NP design. Recent developments in stable electronic and mechanical devices have opened opportunities to monitor the evolution of catalysts in operando under synthesis and reaction conditions; high-speed direct electron detectors have achieved sub-millisecond time resolutions and allow for rapid structural and chemical changes to be captured. Investigations of catalysts using these latest microscopy techniques have provided new insights into atomic-level catalytic mechanisms. Further integration of new microscopy methods is expected to provide multidimensional descriptions of interfaces under relevant synthesis and reaction conditions.In this Account, we discuss recent insights on understanding catalyst activity, selectivity, and stability using advanced STEM techniques, with an emphasis on how critical interfaces dictate the performance of precious metal-based heterogeneous catalysts. The role of extended interfacial structures, including those between core and shell, between separate phases and twinned grains, between the catalyst surface and gas, and between metal and support are discussed. We also provide an outlook on how emerging electron microscopy techniques, such as vibrational spectroscopy and electron ptychography, will impact future catalysis research.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Investigation of the Impact of Flow on MHD Perturbations with the NIMROD Code (Final Report)

This work will enable improved understanding of tokamak edge flows and impurity-species dynamics and their impact on MHD stability as it relates to ELMs, RMPs and QH-mode by leveraging computational developments of the NIMROD code. Experimentally, it is established that the tokamak flows and impurities species can have a substantial effect on MHD stability. This is particularly true for the edge plasma where, for example, large flow shear is correlated with the occurrence of the Quiescent H-mode (QH- mode) state as opposed to operation with edge localized modes (ELMs). Understanding plasma flows is challenging because of the many physical effects that come into play: particle orbits loss, neutrals dynamics, and interaction between multiple ion species with decoupled motion. We will develop a model that incorporates multiple ion species into the MHD framework that includes self-consistent magnetic-field evolution. With this model we will understand how the presence of impurities impacts the tokamak edge-pedestal flows. Finally, we will quantify the impact of the new multiple species modeling on 3D QH-mode simulations. This work will have broader impacts on all MHD studies when multiple species are present (as is the case in modern tokamaks), as well as contributing to the edge and transport communities who also are interested in the studies of multiple-ion-species flows in the edge pedestal region.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The effect of local chemical ordering on Frank-Read source activation in a refractory multi-principal element alloy

In this work, we investigate the operation of Frank-Read (FR) sources in a refractory multi-principal element alloy (MPEA). Simulations of discrete dislocation motion in MPEAs is enabled by the development of a phase field dislocation dynamics model that treats the atomic-scale fluctuations in lattice energies across the glide plane, present due to local ordering in the chemical composition within the nominally random MPEA atomic structure. We consider, through simulation, a range of length scales over which ordering occurs, varying from short-range lengths, a few times dislocation core width, to long-range lengths, an order of magnitude longer than the core. Characteristic of this body-centered cubic MPEA, the simulations also include screw/edge character dependence in glide resistance, as informed by atomic scale simulation. The critical stresses to activate the source for the same source size are found to be statistically distributed, as a direct consequence of the underlying variation in lattice energy. Furthermore, analysis of critical state for activating edge and screw FR sources in the MPEA reveals that FR source operation occurs via a two-step mechanism, involving athermal kink-pair formation, unlike the conventional FR source operation in a material with no composition fluctuations. This mechanism lowers the average critical stress required to activate the FR source and causes the statistical dispersion in critical stress to depend on the range of composition ordering. More importantly, it leads to a more severe dependence of source strength on FR source length than predicted by line tension alone.

36 MATERIALS SCIENCE↗

Spatial dynamics of lithiation and lithium plating during high-rate operation of graphite electrodes

The principal inhibitor of fast charging lithium ion cells is the graphite negative electrode, where favorable conditions for lithium plating occur at high charge rates, causing accelerated degradation and safety concerns. The local response of graphite, both at the electrode and particle level, when exposed to fast charging conditions of around 6C is not well understood. Consequently, the conditions that lead to the onset of lithium plating, as well as the local dynamics of lithium plating and stripping, have also remained elusive. Here, we use high-speed (100 Hz) pencil-beam X-ray diffraction to repeatedly raster along the depth of a 101 μm thick graphite electrode in 3 μm steps during fast (up to 6C) charge and discharge conditions. Consecutive depth profiles from separator to current collector were each captured in 0.5 seconds, giving an unprecedented spatial and temporal description of the state of the electrode and graphite's staging dynamics during high rate conditions. The electrode is preferentially activated near the separator, and the non-uniformity increases with rate and is influenced by free-energy barriers between graphite's lithiation stages. The onset of lithium plating and stripping was quantified, occurring only within the first 15 μm from the separator. The presence of lithium plating changed the behavior of the underlying graphite, such as causing co-existence of LiC 6 and graphite in the fully discharged state. Finally, the staging behavior of graphite at different rates was quantified, revealing a high dependency on rate and drastic hysteresis between lithiation and delithiation.

25 ENERGY STORAGE↗

On the validity and limitations of 1D model for heat and mass transfer performance evaluation in a multilayer binder-free desiccant dehumidifier: isothermal dehumidification with internal cooling

Efficient humidity control is essential for maintaining indoor thermal comfort, yet conventional vapor-compression-based dehumidifiers are energy-intensive. Employing separate sensible and latent cooling through desiccant-coated heat exchangers (DCHEs) combined with evaporative coolers offers energy savings of up to 80 % compared to conventional systems. However, the dehumidification performance of DCHEs remains limited due to the use of polymer binders for coating desiccant materials onto heat exchange surfaces. In our previous study, we developed a multilayer fixed-bed binder-free desiccant dehumidifier (MFBDD) that demonstrated high dehumidification capacity and low pressure drop compared to rotary desiccant wheels. Nevertheless, its potential for further enhancement through internal cooling and the use of step-shaped adsorption isotherms has not been explored. In this study, a physics-based one-dimensional (1D) transient model is developed and validated to capture the coupled heat and mass transfer processes in the MFBDD and extended to simulate internal cooling using a high-capacity composite metal–organic framework, MIL-101/GO-6 (water uptake ≈1.6 g/g within 35–47 % RH). The model enables detailed analysis of local air and bed temperature dynamics and quantifies how internal cooling affects the dehumidification performance under a wide range of operating conditions. Results show that integrating internal cooling and using MIL-101/GO-6 enhance mass adsorbed, moisture removal capacity, and dehumidification effectiveness by 50 %–99 % compared with the M.S. Gel baseline. The study further reveals that achieving near-isothermal operation requires simultaneous enhancement of the convective heat transfer coefficient and heat exchange surface area. In conclusion, this work provides the first detailed physical insight into the interplay between internal cooling and step-shaped isotherms in a binder-free desiccant device and establishes a validated modeling framework for scaling up and system-level performance evaluation of next-generation energy-efficient dehumidification systems.

Heat and mass transfer↗

Global-, local-, and intermediate-scale structures in prototype spiral galaxies

The relationship between galactic spiral structure and the matter in the underlying disk constitutes one of the central problems in galactic dynamics. In Bertin et al. (1989), disk matter characterized by a low-dispersive speed is shown to be capable of playing a key role in the generation of large-scale spiral structure. In Roberts et al. (1992), this self-gravitating, low-dispersion disk matter is shown to be capable of playing an essential role in the formation of structure on local and intermediate scales. Both in computed cases where large-scale spiral structure is present and in those where it is not, the same dominant physical processes and fundamental dynamical mechanisms are active on local scales. The new perception, in which large-scale and small-scale phenomena operate somewhat independently as evidenced in the computational studies, permits a range of flocculent, multiarmed, and grand design spiral types to be simulated. In particular, grand design galaxies with ragged appearances exhibiting spurs, arm branchings, and interarm bridges in addition to the major spiral arms, similar to those often observed, can be generated.

Roberts, William W., Jr.↗

An Overview of the Waste-to-Energy System Simulation (WESyS) Model

The leveraging of waste streams for energy and chemical production could add revenue to waste disposal operations, and it presents opportunities for addressing a variety of economic and environmental objectives at the local, state, and national levels. The Waste-to-Energy System Simulation (WESyS) model is a system dynamics model that was created to simulate the development of the U.S. waste-to-energy industry over time. For each of the three primary waste resources modeled (landfills, concentrated animal feeding operations, and publicly owned treatment works), WESyS simulates technically feasible scenarios for use of the waste, including direct conversion to fuels, and anaerobic digestion followed by flaring, electricity generation, combined heat and power, cleanup and compression to compressed natural gas, and cleanup and injection into an existing pipeline. The model allows users to explore numerous plausible future scenarios for the development of the U.S. waste-to-energy industry. This report provides an overview of the WESyS model and documents the key assumptions, equations, and data sources used to create the model.

09 BIOMASS FUELS↗

XPRESS: eXascale PRogramming Environment and System Software

The XPRESS project is the research and development of a complete software architecture for Exascale computing. The four major R&D thrusts include: 1) an Exascale lightweight kernel operating system (LXK), based on the Kitten operating system (OS), to manage billion way hardware parallelism, management of faults and power, management of global virtual name space, and other features of future system architectures; 2) a runtime system (HPX-5), co-designed with LXK, based on the ParalleX execution model that supports dynamic resource management and task scheduling; 3) system interfaces for interoperability between the runtime system and both the OS and the runtime system (RIOS); 4) compilation strategies and systems to translate MPI and OpenMP legacy codes to a form that can be run with performance at least as good as a native code implementation. The software architecture development has been driven by performance evaluation of mini-applications developed at Department of Energy Co-Design Centers as well as ancillary in-house developed applications and benchmarks. Research on the ParalleX execution model has been driven by the development of formal semantics and model checkers for the execution model. The major components of the software architecture are: ParalleX execution model – guiding principles for co-design of components of the system stack with advances in locality management and task prioritization through introspection; HPX runtime system – support of application dynamic adaptive resource management, task scheduling, and introspective control policies; LXK operating system – lightweight kernel operating system for order-constant scalability and low/no noise to manage resources; RIOS – a realm of Exascale system operation unique in the XPRESS project that includes the relationship between the new generation of lightweight kernel operating systems and runtime system software; XPI advanced programming model – intermediate form and low-level (readable) programming interface reflecting the ParalleX model, providing a target for source-to- source high level parallel language translation, and supporting early direct programming experimentation and measurement; Performance models & metrics – providing parameters and their mutual sensitivities to guide co-design and quantify operational behavior; Legacy migration – ensuring seamless transition of legacy codes and programming methods to the future generation of ParalleX-based Exascale systems; . Experiments and evaluation – critical to determining degree of effectiveness and likelihood of ultimate success as well as guiding corrective design changes to achieve DOE objectives; and Applications – collaborations with Co-Design Centers and other mission critical codes. The exemplar implementation of the ParalleX execution model, HPX-5, has been released along with five fully ported applications (LULESH, HPCG, SSSP, WAMR, and DASHMM) and multiple microbenchmarks. Using these applications, the system architecture demonstrated dynamic load balancing through active migration, legacy code migration, and data-driven approaches to for graph-like applications. Crucially for Exascale and billion way parallelism, the needs and patterns of irregular applications in HPX-5 were addressed using local flow control rather than relying on globally synchronized data structures.

97 MATHEMATICS AND COMPUTING↗

The symbolic computation and automatic analysis of trajectories

Research was generally done on computation of trajectories of dynamical systems, especially control systems. Algorithms were further developed for rewriting expressions involving differential operators. The differential operators involved arise in the local analysis of nonlinear control systems. An initial design was completed of the system architecture for software to analyze nonlinear control systems using data base computing.

Grossman, Robert↗

Fault Detection and Diagnosis In Hall–Héroult Cells Based on Individual Anode Current Measurements Using Dynamic Kernel PCA

Individual anode current signals in aluminum reduction cells provide localized cell conditions in the vicinity of each anode, which contain more information than the conventionally measured cell voltage and line current. One common use of this measurement is to identify process faults that can cause significant changes in the anode current signals. While this method is simple and direct, it ignores the interactions between anode currents and other important process variables. This paper presents an approach that applies multivariate statistical analysis techniques to individual anode currents and other process operating data, for the detection and diagnosis of local process abnormalities in aluminum reduction cells. Specifically, since the Hall–Héroult process is time-varying with its process variables dynamically and nonlinearly correlated, dynamic kernel principal component analysis with moving windows is used. The cell is discretized into a number of subsystems, with each subsystem representing one anode and cell conditions in its vicinity. The fault associated with each subsystem is identified based on multivariate statistical control charts. The results show that the proposed approach is able to not only effectively pinpoint the problematic areas in the cell, but also assess the effect of the fault on different parts of the cell.

36 MATERIALS SCIENCE↗

In situ catalyst activation and regeneration enable energy-efficient high-current CO 2 reduction to ethanol-rich C 2+ mixtures

Electrochemical conversion of dissolved CO 2 in bicarbonate electrolytes, i.e., bicarbonate electrolysis, offers distinct advantages over gas diffusion electrode systems by enabling direct utilization of the CO 2 capture electrolyte while bypassing the energy-intensive CO 2 release step. However, bicarbonate electrolysis faces challenges such as CO 2 mass-transfer limitation, local pH-driven CO 2 depletion, and high cathodic potentials. The higher potential often causes catalyst surface reorganization, leading to a gradual loss of active sites and variations in selectivity during CO 2 reduction. Here, we report a directed, in situ activation and regeneration method that allows precatalysts to equilibrate under dynamic (pulsed) electrolysis conditions. We demonstrate in situ activation of a scalable Cu 2 O/Cu mesh that, under short-width (t = 4 s) pulsed electrolysis, provides stable mixed oxidation states of Cu, favoring the formation of an ethanol-rich crude mixture. The pulsed electrolysis waveform, consisting of six distinct segments, is tuned to form Cu + oxides, which are then reduced to generate local alkaline conditions favoring C–C coupling. This synergistic effect results in FEs of 73% for C2+ products and 39% for ethanol at an applied current density of −150 mA cm −2 and a cathodic potential of −1.45 V (vs. RHE). The overall half-cell energy efficiency is ∼30% for C 2+ products. The in situ Raman experiments confirm the role of pCO 2 R in dynamically regenerating Cu+-containing surface species during pulsed operation, thereby steering selectivity towards C 2+ products. A comprehensive multiscale, multiphysics model is developed to investigate the dynamic behavior of copper surface species (Cu, Cu + , and Cu 2+ ) and local microenvironmental conditions during the pCO 2 R. The results reveal that the coexistence of different copper oxidation states, especially the Cu+ intermediate, is critical in steering selectivity towards multicarbon (C 2+ ) products. The dynamic modulation of surface redox states via tailored pulsing strategies favors C–C coupling pathways by inducing localized alkaline conditions and stabilizing reactive intermediates. This work establishes a predictive modeling platform that links pulse waveform design with mechanistic insights into catalyst state evolution and product selectivity. Overall, this study provides valuable insights into the synergistic effect of in situ activation of pre-catalysts and pulsed electrolysis for higher selectivity towards C 2+ products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

End-to-end microgrid protection using distributed data-driven methods

This paper introduces an end-to-end microgrid protection framework that offers real-time system monitoring, fault-related decision making, and circuit breaker control. This is achieved through the design of distributed data-driven techniques based on the support vector machine method, where each relay is responsible for distributed data collection, fault detection, fault localization, and fault isolation. Local communication is established among neighboring relays, fostering cooperative fault localization and isolation. This decentralized design not only reduces the computational and communication requirements but also enables the adaptability of each relay under varying operational dynamics. The proposed end-to-end protection framework was validated using MATLAB/Simulink simulations on a 100% renewable microgrid, achieving an accuracy of 93.1% with response time of 0.0523 s, in protecting against a range of fault scenarios that are characterized by various types, locations, impedances, load conditions, photovoltaic power levels, and microgrid operating modes.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Plasma Imaging, LOcal Measurement, and Tomographic Experiment (PILOT): A Mission Concept for Transformational Multi-Scale Observations of Mass and Energy Flow Dynamics in Earth’s Magnetosphere

We currently do not understand the fundamental physical processes that govern mass and energy flow through the Earth’s magnetosphere. Knowledge of these processes is critical to understanding the mass loss rate of Earth’s atmosphere, as well as for determining the role that a planetary magnetic field plays in atmospheric retention, and therefore habitability, for Earth-like planets beyond the solar system. Mass and energy flow processes are challenging to determine at Earth in part because Earth’s planetary magnetic field creates a complex “system of systems” composed of interdependent plasma populations and overlapping spatial regions that perpetually exchange mass and energy across a broad range of temporal and spatial scales. Further, the primary mass carrier in the magnetosphere is cold plasma (as cold as ~0.1 eV), which is invisible to many space-borne instruments that operate in the inner magnetosphere. The Plasma Imaging LOcal and Tomographic experiment (PILOT) mission concept, described here, provides the transformational multi-scale observations required to answer fundamental open questions about mass and energy flow dynamics in the Earth’s magnetosphere. PILOT uses a constellation of spacecraft to make radio tomographic, remote sensing, and in-situ measurements simultaneously, fully capturing cold plasma mass dynamics and its impact on magnetospheric systems over an unprecedented range of spatial and temporal scales. This article details the scientific motivation for the PILOT mission concept as well as a potential mission implementation.

79 ASTRONOMY AND ASTROPHYSICS↗

Gain dynamics of quantum dot devices for dual-state operation

Ground state gain dynamics of In(Ga)As-quantum dot excited state lasers are investigated via single-color ultrafast pump-probe spectroscopy below and above lasing threshold. Two-color pump-probe experiments are used to localize lasing and non-lasing quantum dots within the inhomogeneously broadened ground state. Single-color results yield similar gain recovery rates of the ground state for lasing and non-lasing quantum dots decreasing from 6 ps to 2 ps with increasing injection current. We find that ground state gain dynamics are influenced solely by the injection current and unaffected by laser operation of the excited state. This independence is promising for dual-state operation schemes in quantum dot based optoelectronic devices.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

MPC-Based Local Voltage Control Strategy of DGs in Active Distribution Networks

The increasing penetration of distributed generators (DGs) exacerbates the voltage violations in active distribution networks (ADNs). The commonly used centralized voltage control may aggravate the burden of computation and communication. Local voltage control can realize fast response to frequent voltage fluctuations. However, effective local control for DGs still needs to be further investigated. In this paper, a model predictive control (MPC)-based local voltage control strategy of DGs is proposed to cope with voltage violation problems and enhance the adaptation to the DG volatility. DGs are regulated with local Q-V control curves to respond to voltage fluctuations rapidly. Based on the dynamic network partition, a unified control curve for DGs is set in the same area to realize efficient operation. To further improve the performance of voltage control, the parameter tuning of control curves is regularly conducted based on decentralized inter-area coordination with the rolling of the control horizon. Finally, case study is verified on the modified IEEE 33-node distribution system. The proposed local control strategy of DGs can effectively improve the voltage profile and reduce power losses of the whole network without high computation and communication burden.

active distribution network (ADN)↗

Fluid System Stability Analysis Techniques

Fluid systems, or networks, consist of multiple components that work together to achieve some desired thermofluid state. For any generic application, this fluid state can be a combination of the fluid pressure, flow rate, enthalpy, or species concentration. Fluid system components, such as pumps and valves, are often governed by nonlinear differential equations, resulting in complex component-to-component interactions. System-level fluid network stability occurs when the flow through the system can maintain a steady-state solution in the presence of small perturbations, which depends on these component interactions. System instability, however, can go undetected until issues arise during integrated system testing. This presentation explores a method for system designers to think of the fluid network as an assembly of components, each with their own thermofluid surfaces of partial stability, called nullclines. The intersections of all nullclines yields system-level solutions, called equilibrium points. When designers define operating points, they are tuning system parameters so that these equilibrium points move to the desired location in the thermofluid state plane. However, linearization theory shows us that the dynamic behavior around these points can be unstable. The local stability of these equilibrium points can be assessed analytically with eigen-analysis, or numerically by propagating state-plane samples to construct a phase portrait. Investigating a phase portrait can help designers gain a qualitative understanding of a system’s dynamic performance. This understanding can then help inform requirement definitions, component selection, and operational procedures. This presentation includes an example of the phase portrait technique on a system featuring a centrifugal pump and a back-pressure regulator (BPR). Numerical modeling of this system suggests that equilibrium points on the left-hand side of the pump curve are dynamically unstable.

Fluid System↗

Fluid System Stability Analysis Techniques

Fluid systems, or networks, consist of multiple components that work together to achieve some desired thermofluid state. For any generic application, this fluid state can be a combination of the fluid pressure, flow rate, enthalpy, or species concentration. Fluid system components, such as pumps and valves, are often governed by nonlinear differential equations, resulting in complex component-to-component interactions. System-level fluid network stability occurs when the flow through the system can maintain a steady-state solution in the presence of small perturbations, which depends on these component interactions. System instability, however, can go undetected until issues arise during integrated system testing. This presentation explores a method for system designers to think of the fluid network as an assembly of components, each with their own thermofluid surfaces of partial stability, called nullclines. The intersections of all nullclines yields system-level solutions, called equilibrium points. When designers define operating points, they are tuning system parameters so that these equilibrium points move to the desired location in the thermofluid state plane. However, linearization theory shows us that the dynamic behavior around these points can be unstable. The local stability of these equilibrium points can be assessed analytically with eigen-analysis, or numerically by propagating state-plane samples to construct a phase portrait. Investigating a phase portrait can help designers gain a qualitative understanding of a system’s dynamic performance. This understanding can then help inform requirement definitions, component selection, and operational procedures. This presentation includes an example of the phase portrait technique on a system featuring a centrifugal pump and a back-pressure regulator (BPR). Numerical modeling of this system suggests that equilibrium points on the left-hand side of the pump curve are dynamically unstable.

Fluid System↗