Search NASA⌕ Search

SEARCH · Search NASA

Results for “coupling methodology”

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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Bacterial community dynamics as a result of growth-yield trade-off and multispecies metabolic interactions toward understanding the gut biofilm niche

Abstract Bacterial communities are ubiquitous, found in natural ecosystems, such as soil, and within living organisms, like the human microbiome. The dynamics of these communities in diverse environments depend on factors such as spatial features of the microbial niche, biochemical kinetics, and interactions among bacteria. Moreover, in many systems, bacterial communities are influenced by multiple physical mechanisms, such as mass transport and detachment forces. One example is gut mucosal communities, where dense, closely packed communities develop under the concurrent influence of nutrient transport from the lumen and fluid-mediated detachment of bacteria. In this study, we model a mucosal niche through a coupled agent-based and finite-volume modeling approach. This methodology enables us to model bacterial interactions affected by nutrient release from various sources while adjusting individual bacterial kinetics. We explored how the dispersion and abundance of bacteria are influenced by biochemical kinetics in different types of metabolic interactions, with a particular focus on the trade-off between growth rate and yield. Our findings demonstrate that in competitive scenarios, higher growth rates result in a larger share of the niche space. In contrast, growth yield plays a critical role in neutralism, commensalism, and mutualism interactions. When bacteria are introduced sequentially, they cause distinct spatiotemporal effects, such as deeper niche colonization in commensalism and mutualism scenarios driven by species intermixing effects, which are enhanced by high growth yields. Moreover, sub-ecosystem interactions dictate the dynamics of three-species communities, sometimes yielding unexpected outcomes. Competitive, fast-growing bacteria demonstrate robust colonization abilities, yet they face challenges in displacing established mutualistic systems. Bacteria that develop a cooperative relationship with existing species typically obtain niche residence, regardless of their growth rates, although higher growth yields significantly enhance their abundance. Our results underscore the importance of bacterial niche dynamics in shaping community properties and succession, highlighting a new approach to manipulating microbial systems.

Microbiology↗
DOE OSTI2020

A Mercury Model of the Molly-G Fast Burst Reactor (FBR)

In order to design neutron irradiation experiments at a given nuclear reactor, one requires a representative neutron spectrum for the configuration of the reactor that will exist during the experiment. This usually implies modeling the operation of the reactor using a neutron transport code. Since the modeling of pulsed neutron reactors is quite complex, requiring coupled modeling of neutronics and thermo-structural response, the typical methodology employs analysis of the neutron spectrum which is obtained from a static keff calculation. While this time-independent Eigenfunction is not truly representative of that during a reactor pulse, we postulate that it is a reasonable representation for our purposes. The Nuclear Survivability (NS) Program at LLNL is considering fielding such experiments at the Molly-G (‘Molybdenum Godiva-II’) Fast Burst Reactor (FBR) at the White Sands Missile Range (WSMR). The goal of these experiments is to (a) test / proof / calibrate diagnostics, and (b) test the efficacy of neutron shield configurations, before fielding at other nuclear reactors, such as the Annular Core Research Reactor (ACRR) at SNL/NM. After attempting to secure either drawings or a copy of an existing MCNP model of the Molly-G FBR from various institutions, it was decided to develop a new model for use with the Mercury Monte Carlo transport code. The model described herein will be used to design these experiments.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2022

Modeling Methods for Capturing System Interactions of Combined Technologies: A Study of PV + Battery

The costs of solar photovoltaics (PV) have been dropping in recent years, leading to increasing installations of solar PV systems and growing interest in how the technology will impact the electric grid at higher penetrations. Additionally, as battery costs decline it becomes important to understand the benefits and limitations of battery storage for the grid as well as potential benefits of co-locating battery storage and PV, particularly within the realm of future system planning. However, it is non-trivial to capture these potential benefits within a linearized capacity expansion model (CEM). This paper presents methodological developments to more fully represent the value and limitations of coupled PV and battery systems (PV + Battery) in CEMs using the Resource Planning Model (RPM), which co-optimizes capacity investments, transmission investments, and reduced-order dispatch in the Western Interconnection of North America through 2045. We use the model to simulate the evolution of the generation and transmission system under two core scenarios - a baseline scenario and a high renewable penetration scenario - coupled with sensitivities assuming low and midline PV and battery cost projections. When incorporating PV + battery, we find that it is important for CEMs to capture the ability of the coupled technology to provide firm capacity and reduce expected curtailment, compared to independent systems. These interactions can have even more dramatic impacts at higher solar penetrations.

14 SOLAR ENERGY↗
DOE OSTI2020

Pin-pointing Node Failures in HPC Systems

Automated fault prediction and diagnosis in HPC systems needs to be efficient for better system resilience. With increasing scalability required for exascale, accurate fault prediction aiding in quick remedy is hard. With changing supercomputer architectures, distilling fault data from the noisy raw logs requires substantial efforts. Predicting node failures in such voluminous system logs is challenging. To this end, we investigate an interesting way to pin-point node failures in such supercomputing systems. Our study on Cray system data with automated machine learning tools suggests that specific patterns of event messages on node unavailability can be indicator to node failures. This data extraction coupled with system and job data correlation helps in devising a methodology to predict node failures and their location over a specific time frame. This work aims to enable broader applicability for a generic fault prediction framework.

Roman, E↗
DOE OSTI2025

Direct Study of Changes in Catalyst Structure-Kinetic Properties During Redox Transitions

The Temporal Analysis of Products (TAP) pulse response methodology is a transient technique that provides the time resolution needed to deconvolve reaction steps from the complex networks typical in industrial catalytic processes. Traditionally, TAP measurements observe gas phase dynamics at the reactor exit but lack direct measurements of changes in the catalyst itself. Recently, a new operando technique was developed that couples gas phase transients to dynamic changes in metal centers with the precise TAP methodology for nanomole titration. Using an industrial CrOx/Al2O3 catalyst used for propane dehydrogenation, we demonstrate the capabilities of this unique device to reveal key catalytic processes: 1) total propane oxidation not associated with chromia centers, 2) reduction of Cr6+ to Cr3+ correlated with selective product formation, and 3) subsequent carbon accumulation. By utilizing incremental pulsing in a diffusion-only transport regime, the spectrokinetic device allows us to resolve detailed changes in catalyst structure, composition, and kinetic function that are otherwise indistinguishable in conventional operando devices. The unification of the TAP methodology with time-resolved spectroscopic measurements offers new and unique insights into the complex kinetic phenomena regulated by solid catalyst surfaces.

03 - NATURAL GAS↗
DOE OSTI2021

A Simple Criterion for Feasibility of Heat Integration between Distillation Streams Based on Relative Volatilities

In a multicomponent distillation configuration, there are numerous sources and sinks of heat, and a potential way to reduce the heat duty requirement is to perform heat integration. Unfortunately, an algorithmic search of the optimal heat integration opportunities is intractable when the required temperatures of intermediate streams are computed via complex models. Instead, in this work, we introduce pressure-scaled pseudo relative volatility, a new metric to compare stream temperatures. Here, we justify the use of pseudo relative volatility by proving that this variable is a monotonically increasing function of the liquid fraction in a saturated mixture stream. Using this metric, we derive a shortcut criterion to check the feasibility of various heat integration opportunities, such as thermal coupling via heat transfer (TCH). The advantage of this approach is that it circumvents the need for explicit temperatures and instead relies on composition, component relative volatilities, and pressure—quantities that are readily available in shortcut models for optimization of distillation configurations. Leveraging this fact, we propose a new optimization framework that identifies feasible TCHs that we consider within the formulation while minimizing the total heat duty of a distillation configuration. We demonstrate, on a few examples, that our formulation can identify heat duty efficient configurations, some of which are multieffect configurations. Using this methodology, we discover configurations that are not only simpler than the fully thermally coupled (FTC) configuration but also have a much lower heat duty.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

A Transferable Force Field for Predicting Adsorption and Diffusion of Water in Cationic Zeolites with Coupled Cluster Accuracy

We present a transferable force field for water in proton-exchanged, alkali (Li, Na, K, Rb, and Cs) metal-exchanged, and alkaline-earth (Mg, Ca, Sr, and Ba) metal-exchanged zeolites. The fitting methodology is based on adsorbate–adsorbent interaction energies obtained from periodic density functional theory calculations and corrected using the coupled-cluster method applied to small model clusters. To ensure an accurate prediction of both adsorption and diffusion properties of water, sets of configurations that sample both adsorption sites and intracrystalline hopping transition states were used in the fitting. The quality of the force field is assessed for a wide range of zeolites with different topologies and chemical compositions, demonstrating good agreement between theoretical predictions and experimental measurements of water adsorption and diffusion.

Adsorption↗
DOE OSTI2024

Forecasting fluid-injection induced seismicity to choose the best injection strategy for safety and efficiency

Induced seismicity poses a challenge to the development of Enhanced Geothermal Systems (EGS). Improving monitoring and forecasting techniques is essential to mitigate induced seismicity and thereby fostering a positive perception of EGS projects among local authorities and population. Induced seismicity is the result of complex and coupled thermo-hydro-mechanical-chemical mechanisms. Injection flux and pressure are crucial controlling parameters for both hydraulic stimulation and circulation protocols. We develop a methodology combining a hydro-mechanical model with a seismicity rate model to estimate the magnitude and frequency of mainshocks and aftershocks induced by fluid injection. We apply the methodology to the case of the Basel EGS (2006, Switzerland) to compare the effects of progressive, cyclic and constant injections on the mechanical response of discrete faults. Results from the coupled hydro-mechanical models show that the pore pressure diffusion and consequent enhancement of fault permeability are limited to the vicinity of the injection well during cyclic injection. Additionally, constant injection induces seismicity from the start of the injection but enhances the permeability of most of the faults within a shorter duration, inducing less post-injection seismicity. The methodology can be adapted to any numerical model and allows new projects to be developed by anticipating the safest injection protocol. This article is part of the theme issue ‘Induced seismicity in coupled subsurface systems’.

58 GEOSCIENCES↗
DOE OSTI2022

Modeling the viscoplastic behavior of a semicrystalline polymer

In this study, a complex constitutive relation is identified using inverse modeling with the nominal mechanical response as sole experimental input. The methodology is illustrated for a semicrystalline thermoplastic in the presence of strain localization at finite deformations. The experimental database includes cylindrical tensile bars, compression pins and round notched bars loaded at strain rates spanning up to five decades and temperatures below and above T g . The data is organized into a calibration set and a validation set. The response of tensile specimens is determined using finite element analyses and a two-phase constitutive relation for semicrystalline polymers that accounts for temperature- and rate-sensitive plastic flow, pressure-sensitivity, small-strain softening and large-strain orientational hardening of the amorphous phase, along with the evolution of crystallinity. The large number of constitutive parameters is identified using an optimization tool coupled with the finite element solver and the calibration set from experiments. The methodology is shown to be successful in predicting the response of round notched bars and replicating the effects of temperature and strain rate on the severity of necking in tensile bars. The proposed model identification strategy is both simple and effective in comparison with other elaborate methods that attempt to access intrinsic behavior directly from high-fidelity experimental measurements.

36 MATERIALS SCIENCE↗
DOE OSTI2022

The Schwarz Alternating Method for the Seamless Coupling of Nonlinear Reduced Order Models and Full Order Models

Projection-based model order reduction allows for the parsimonious representation of full order models (FOMs), typically obtained through the discretization of a set of partial differential equations (PDEs) using conventional techniques (e.g., finite element, finite volume, finite difference methods) where the discretization may contain a very large number of degrees of freedom. As a result of this more compact representation, the resulting projection-based reduced order models (ROMs) can achieve considerable computational speedups, which are especially useful in real-time or multi-query analyses. One known deficiency of projection-based ROMs is that they can suffer from a lack of robustness, stability and accuracy, especially in the predictive regime, which ultimately limits their useful application. Another research gap that has prevented the widespread adoption of ROMs within the modeling and simulation community is the lack of theoretical and algorithmic foundations necessary for the “plug-and-play” integration of these models into existing multi-scale and multi-physics frameworks. This paper describes a new methodology that has the potential to address both of the aforementioned deficiencies by coupling projection-based ROMs with each other as well as with conventional FOMs by means of the Schwarz alternating method [41]. Leveraging recent work that adapted the Schwarz alternating method to enable consistent and concurrent multiscale coupling of finite element FOMs in solid mechanics [35, 36], we present a new extension of the Schwarz framework that enables FOM-ROM and ROM-ROM coupling, following a domain decomposition of the physical geometry on which a PDE is posed. In order to maintain efficiency and achieve computation speed-ups, we employ hyper-reduction via the Energy-Conserving Sampling and Weighting (ECSW) approach [13]. We evaluate the proposed coupling approach in the reproductive as well as in the predictive regime on a canonical test case that involves the dynamic propagation of a traveling wave in a nonlinear hyper-elastic material.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Simulation-Based Recovery Action Analysis Using the EMRALD Dynamic Risk Assessment Tool

A recovery action is defined as the action that prevents deviant conditions from producing unwanted effects. It generally indicates a kind of countermeasure performed in response to a failure of human action. The recovery actions especially play an important role in complex systems like nuclear power plants (NPPs), which consist of highly sophisticated controllers to ensure that desired performance and safety must be achieved and maintained. This is because a combination of human error and its recovery failure may be able to cause a catastrophic effect on a system. Analyzing recovery actions has been a critical part of HRA, which is a technique to evaluate human errors and provide human error probabilities (HEPs) for application in probabilistic safety assessment (PSA). If recovery actions are not adequately analyzed and applied to PSA models, the PSA results may be under-estimated or be not able to reasonably account for the failure of human actions in the context of PSA. For this reason, some regulatory documents such as ASME/ANS RA-Sb-2013 by the American Society for Mechanical Engineers and the American Nuclear Society and NUREG-1792 by U.S. Nuclear Regulatory Commission have emphasized the importance of recovery analysis within the HRA. A couple of existing HRA methods, such as the Technique for Human Error-Rate Prediction (THERP), the Cause-Based Decision Tree (CBDT), and the Korean Standard HRA (K-HRA), have respectively suggested their own approaches to the HRA recovery analysis. However, there are a couple of limitations to treating recovery actions using only the current HRA methods available. The biggest limitation is that the existing recovery analysis does not explicitly consider a variety of recovery action types and recovery sequences as they occur in actual NPPs. To handle the limitations of existing recovery analysis, this study proposes a simulation-based recovery analysis method using the Event Modeling Risk Assessment Using Linked Diagram (EMRALD) software. The EMRALD software is a dynamic simulation tool for PSA. It supports realistic and dynamic modeling of human actions as they would be performed at NPPs. It is also favorable to simultaneously model the specific moment at which an action is performed, the time it takes to perform the action, and the failure probability of that action. In this paper, a detailed methodology for modeling recovery actions in the simulation platform is proposed with a couple of examples. Then, outputs from the simulation are discussed as reviewing if this novel approach can complement the challenges of existing recovery analyses.

99 GENERAL AND MISCELLANEOUS↗
DOE OSTI2024

Energetic and Spectroscopic Properties of Astrophysically Relevant MgC 4 H Radicals Using High-Level Ab Initio Calculations

Considering the importance of magnesium-bearing hydrocarbon molecules (MgC n H; n = 2, 4, and 6) in the carbon-rich circumstellar envelopes (e.g., IRC+10216), a total of 28 constitutional isomers of MgC 4 H have been theoretically investigated using density functional theory (DFT) and coupled-cluster methods. The zero-point vibrational energy corrected relative energies at the ROCCSD(T)/cc-pCVTZ level of theory reveal that the linear isomer, 1-magnesapent-2,4-diyn-1-yl (1, 2 Σ + ), is the global minimum geometry on the MgC 4 H potential energy surface. The latter has been detected both in the laboratory and in the evolved carbon star, IRC+10216. The calculated spectroscopic data for 1 match well with the experimental observations (error ~ 0.78%) which validates our theoretical methodology. Plausible isomerization processes happening among different isomers are examined using DFT and coupled-cluster methods. CASPT2 calculations have been performed for a few isomers exhibiting multireference characteristics. The second most stable isomer, 1-ethynyl-1λ 3 -magnesacycloprop-2-ene-2,3-diyl (2, 2 A 1 , μ = 2.54 D), is 146 kJ mol –1 higher in energy than 1 and possibly the next promising candidate to be detected in the laboratory or in the interstellar medium in future.

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

Direct isotopic analysis of solid uranium particulates on cotton swipes by microextraction-ICP-MS

Direct isotope ratio analysis of solid uranium particulates on cotton swipes was achieved using a solution-based microextraction technique, coupled to a quadrupole inductively coupled plasma - mass spectrometer (ICP-MS). This microextraction-ICP-MS methodology provides rapid isotopic analysis which could be applicable to nuclear safeguards measurements. Particulates of uranyl nitrate hexahydrate (UO 2 (NO 3 ) 2 ·6H 2 O) and uranyl fluoride (UO 2 F 2 ) ranging from 6 μm to 40 μm in length were transferred to cotton swipes with a particle manipulator. The microextraction probe then delivers a 5% nitric acid (HNO3) solvent onto the swipe surface to extract the uranium species. The extracted sample is then delivered to the ICP-MS for isotopic determination. The majority of uranium signal (~99% and ~94% for UO2(NO3)2·6H2O and UO 2 F 2 , respectively) was detected in the first 15 s extraction, while subsequent extractions on the same location had low or no U signal, suggesting near complete removal of the solid uranium compounds from the swipe surface. Ten samples (for each of the uranium compounds), were analyzed for their isotopic composition. For UO 2 (NO 3 ) 2 ·6H 2 O, the determined isotope ratios resulted in a % relative difference (% RD) from the referenced isotope ratios of 0.97, 1.0, and 7.3% for 234 U/ 238 U, 235 U/ 238 U, and 236 U/ 238 U, respectively. The % RD of the UO 2 F 2 isotope ratios were 1.9 and 0.60% for 234 U/ 238 U and 235 U/ 238 U, respectively. The preliminary limits of detection were determined to be 0.002, 0.4, and 60 pg for 234 U, 235 U and 238 U, respectively This work demonstrates that microextraction ICP-MS is a rapid and sensitive method that could directly determine uranium isotope ratios of UO 2 (NO 3 ) 2 ·6H 2 O and UO 2 F 2 particulates on cotton swipes.

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

Neutron transport methods for multiphysics heterogeneous reactor core simulation in Griffin

Griffin is a reactor physics application based on the Multiphysics Object-Oriented Simulation Environment (MOOSE). This work discloses the methods, algorithms, and implementation for simulating heterogeneous reactor dynamics models. Griffin utilizes a discontinuous finite-element method with discrete ordinates (DFEM-S ) to discretize the field variable of the multigroup neutron transport equation. Multiphysics feedback is handled using two-step tabulated cross-section methodology. Feedback quantities are evaluated using the MOOSE-MultiApp system to couple various engineering phenomena, such as heat conduction and thermal fluids. The multiphysics DFEM-S system is solved using fixed-point iteration with a fully asynchronous parallel sweeper, unstructured coarse-mesh finite difference acceleration, and a multi-timescale improved quasi-static method scheme. The implementation is applied to a multiphysics microreactor model, with two transients: one initiated by a single heat-pipe failure and another by control drum rotation. Importantly, these examples demonstrate the ability of Griffin to tractably solve the neutron transport equation considering seven independent variables and feedback.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Goal-oriented real-time Bayesian inference for linear autonomous dynamical systems with application to digital twins for tsunami early warning

We present a goal-oriented framework for constructing digital twins with the following properties: (1) they employ discretizations of high-fidelity partial differential equation (PDE) models governed by autonomous dynamical systems, leading to large-scale forward problems; (2) they solve a linear inverse problem to assimilate observational data to infer uncertain model components followed by a forward prediction of the evolving dynamics; and (3) the entire end-to-end, data-to-inference-to-prediction computation is carried out without approximation and in real time through a Bayesian framework that rigorously accounts for uncertainties. Several challenges must be overcome to realize this framework, including the large scale of the forward problem, the high dimensionality of the parameter space, and for a class of problems including those we target, the slow decay of the singular values of the parameter-to-observable map. Here we introduce a methodology to overcome these challenges by exploiting the autonomous structure of the forward model to decompose the solution of the inverse problem into a one-time-only offline phase in which the PDE model is solved a limited number of times (equal to the number of sensors), and an online phase that maps well onto GPUs and computes the parameter inference and prediction of quantities of interest in real time, given observational data. Our ultimate goal is to apply this framework to construct digital twins for subduction zones, including Cascadia, to provide early warning for tsunamis generated by megathrust earthquakes. To this end, we demonstrate how our methodology can be used to employ seafloor pressure observations, along with the coupled acoustic–gravity wave equations, to infer the earthquake-induced spatiotemporal seafloor motion (discretized with $\mathscr{O}$ (10 9 ) parameters) and forward predict the tsunami propagation. We present results of an end-to-end inference, prediction, and uncertainty quantification for a representative test problem with $\mathscr{O}$ (10 8 ) inversion parameters for which goal-oriented Bayesian inference is accomplished exactly and in real time, that is, in a matter of seconds.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Inverse-designed metasurfaces with facile fabrication parameters

Abstract Optical metasurfaces are planar nanostructured devices that are industrially attractive in part because they utilize high-throughput microelectronic fabrication techniques for implementation. It is therefore critical to develop design paradigms that can balance the realization of highly efficient wavefront responses together with device manufacturability. We introduce a gradient-based design framework for freeform metasurfaces in which nanoscale elements are explicitly constrained to feature basic shapes, nearly uniform feature sizes, and exceptionally low aspect ratios. In spite of the apparent uniformity of the metasurface geometric features, the devices are able to utilize nonlocal near-field optical coupling to achieve highly efficient and extreme wavefront scattering beyond conventional design methodologies. Utilizing this approach, we design facile high-numerical-aperture devices such as beam deflectors and large-area metalenses capable of diffraction-limited focusing. We anticipate that these concepts can facilitate the design and integration of metasurfaces into monolithic optical systems.

Optics↗
DOE OSTI2022

Federated Architecture for Secure and Transactive Distributed Energy Resource Management Solutions (FAST-DERMS)

This document provides system-level specifications for a federated architecture for secure and transactive distributed energy resource management solutions (FAST-DERMS), presents a solution, and describes operational concepts for the proposed solution. FAST-DERMS enables the provision of reliable, resilient, and secure transmission and distribution (T&D) grid services through the scalable aggregation and near-real-time management of utility-scale and small-scale distributed energy resources (DERs). We first present the principles and objectives of FAST-DERMS. Then, after discussing important system concepts, we present the specifications for FAST-DERMS and a solution that employs a distributed and federated control methodology in which the DERs connected to a single point of common coupling with the rest of the system, such as individual substations, are optimized coordinately to provide system-level grid services. FAST-DERMS aims to aggregate and coordinate the operations of DERs to support T&D grid operations. The key optimization and control component of this FAST-DERMS reference implementation is a flexible resource scheduler (FRS) that aggregates the DERs within a substation service area. These FRSs operate at the substation level and perform constrained economic dispatch of DERs, either directly or through a transactive market or aggregator, as shown in Figure ES-1. An FRS Coordinator at the distribution system operator (DSO) level aggregates distribution substations operated by FRSs and interfaces with the transmission system operator (TSO) to provide transmission services. FAST-DERMS also allows for the integration of the FRS Coordinator with an existing distribution utility management system that could be employed by the DSO to enhance distribution grid operations.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2026

Control Co-Design Framework for Joint Optimization of PID Governor and Hydraulic Dynamics for a Hydropower Plant

This paper presents a control co-design (CCD) framework tailored to the hydraulic and mechanical subsystems of hydropower systems, aimed at enhancing their dynamic performance and responsiveness to grid demands. Conventional sequential design approaches often fall short in capturing the coupled interactions between plant dynamics and control objectives. The proposed CCD methodology enables simultaneous optimization of key hydro plant parameters - such as governor settings, penstock characteristics, and valve control - with the control system architecture. Results indicate that both CCD and sequential design approaches achieve comparable transient and steady-state performance, with minor discrepancies arising from differences in their optimized parameters. However, CCD demonstrates a distinct advantage by tripling the water time constant (Tw) while maintaining similar dynamic performance to the sequential approach. This increase in Tw enables optimization of penstock dimensions, as it is directly influenced by penstock design, thereby contributing to overall cost optimization.

24 POWER TRANSMISSION AND DISTRIBUTION↗