Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Linear Parameter Varying”

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 469 records · Page 26

Command of active and responsive elastomers by topological defects and patterns

The project resulted in the development of stimuli-responsive elastomer coatings formed by photopolymerized liquid crystal molecules. The molecular orientation of the liquid crystal elastomers is coupled to rubber-like elasticity. The project developed an approach to produce complex patterns of molecular orientation in elastomers by employing photopatterning technique based on plasmonic metamasks that convert unpolarized incoming light into a transmitted light beam with spatially-varying linear polarization. The polarization-modulated light beam aligns the film of azodye molecules, which serves as the substrate for liquid crystal elastomer coating. The alignment of the azodye molecules follows the polarization pattern with a phase shift of the alignment direction by 90 degrees. The azodye substrate imposes orientation of molecules in the adjacent liquid crystal layer which is preserved after the monomers are polymerized into a liquid crystal elastomer. By using differently aligned substrates one could produce various patterns of molecular orientation in the elastomer, which dictate the mechanical properties of the elastomers and their response to external stimuli such as temperature, humidity, and ultraviolet irradiation. In particular, the predesigned molecular orientation produces deterministically defined topography response of the liquid crystal coatings, in which their free surface changes from flat to either locally elevated or locally compressed, depending on whether the pre-inscribed molecular orientation is predesigned with a splay or bend deformation. The research established that the mechanism of this deterministic relation between the dynamic topography and pre-inscribed molecular orientation is rooted in the change of the tensorial order parameter that characterizes the degree of orientational order, and in the contraction/expansion of polymer molecules in response to the order parameter changes. In particular, lowering the order causes contraction of polymer globules. The work demonstrated that the dynamic topography of elastomer coatings with three-dimensional pattern of orientational order includes up and down motion of the material (along the normal to the coating) and also shifts in lateral directions. We demonstrate that the deterministic relationship between the orientational patterns and variations of coatings’ profile is caused by forces that are mathematically equivalent to active forces in the system of “pullers” or “pushers” in active matter, demonstrating universality of the out-of-equilibrium behavior. The project resulted in theoretical models that predictively describe the dynamic response of coatings. The topography of liquid crystal elastomers is sensitive to stimuli such as temperature, ultraviolet irradiation, and humidity. The applicability of the patterned approach to produce dynamic coatings was expanded from nematic elastomers to their smectic counterparts. The project combined a battery of imaging techniques to unveil the dynamic properties of the coatings, ranging from the state-of-the-art dynamic holographic microscopy to fluorescent confocal polarizing microscopy and polscope microscopy. The research uncovered potential applications of dynamic elastomer coatings, such as control placement of colloidal particles, engineering of living tissues with alignment patterns of cells that follow the orientational order of elastomer coatings. The approach developed in the project can be used to design programmable dynamic coatings with functionalities that mimic biological tissues such as skin.

36 MATERIALS SCIENCE↗

Assimilation of MODIS Snow Cover Through the Data Assimilation Research Testbed and the Community Land Model Version 4

To improve snowpack estimates in Community Land Model version 4 (CLM4), the Moderate Resolution Imaging Spectroradiometer (MODIS) snow cover fraction (SCF) was assimilated into the Community Land Model version 4 (CLM4) via the Data Assimilation Research Testbed (DART). The interface between CLM4 and DART is a flexible, extensible approach to land surface data assimilation. This data assimilation system has a large ensemble (80-member) atmospheric forcing that facilitates ensemble-based land data assimilation. We use 40 randomly chosen forcing members to drive 40 CLM members as a compromise between computational cost and the data assimilation performance. The localization distance, a parameter in DART, was tuned to optimize the data assimilation performance at the global scale. Snow water equivalent (SWE) and snow depth are adjusted via the ensemble adjustment Kalman filter, particularly in regions with large SCF variability. The root-mean-square error of the forecast SCF against MODIS SCF is largely reduced. In DJF (December-January-February), the discrepancy between MODIS and CLM4 is broadly ameliorated in the lower-middle latitudes (2345N). Only minimal modifications are made in the higher-middle (4566N) and high latitudes, part of which is due to the agreement between model and observation when snow cover is nearly 100. In some regions it also reveals that CLM4-modeled snow cover lacks heterogeneous features compared to MODIS. In MAM (March-April-May), adjustments to snowmove poleward mainly due to the northward movement of the snowline (i.e., where largest SCF uncertainty is and SCF assimilation has the greatest impact). The effectiveness of data assimilation also varies with vegetation types, with mixed performance over forest regions and consistently good performance over grass, which can partly be explained by the linearity of the relationship between SCF and SWE in the model ensembles. The updated snow depth was compared to the Canadian Meteorological Center (CMC) data. Differences between CMC and CLM4 are generally reduced in densely monitored regions.

data assimilation↗

Driven resonant current amplification in self-organized plasma configurations with uniform λ and plasma pressure confinement

Driven resonant current amplification in uniform λ=$^{μoJ}_{B}$, self-organized plasma configurations with plasma pressure confinement is studied. This work extends previous analyses of force-free states (J ∥ B) to equilibria with plasma pressure gradients (∇P = J × B ≠ 0). An analytical treatment shows that in equilibria with uniform λ and a plasma pressure gradient that varies linearly with poloidal magnetic flux, driven resonances are either downshifted or upshifted from their force-free values. Additionally, a modification of off resonance current amplification is observed in pressure confining equilibria. Computational results add support to analytical predictions, and also include driven resonance behavior with equilibrium profiles not considered analytically. Furthermore, these results suggest not only persistence, but also modification of driven resonant current amplification in sustained spheromak configurations with plasma pressure confinement when compared to force-free states. By providing a basis for the existence of driven, high current amplification, pressure confining equilibria, this work provides insight into the accessibility of sustained spheromak states of most interest for potential fusion energy applications.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Reciprocal swimming at intermediate Reynolds number

In Stokes flow, Purcell's scallop theorem forbids objects with time-reversible (reciprocal) swimming strokes from moving. In the presence of inertia, this restriction is eased and reciprocally deforming bodies can swim. A number of recent works have investigated dimer models that swim reciprocally at intermediate Reynolds numbers Re≈1–1000. These show interesting results (e.g. switches of the swim direction as a function of inertia) but the results vary and seem to be case specific. Here, in this work, we introduce a general model and investigate the behaviour of an asymmetric spherical dimer of oscillating length for small-amplitude motion at intermediate Re. In our analysis we make the important distinction between particle and fluid inertia, both of which need to be considered separately. We asymptotically expand the Navier–Stokes equations in the small-amplitude limit to obtain a system of linear partial differential equations. Using a combination of numerical (finite element) and analytical (reciprocal theorem, method of reflections) methods we solve the system to obtain the dimer's swim speed and show that there are two mechanisms that give rise to motion: boundary conditions (an effective slip velocity) and Reynolds stresses. Each mechanism is driven by two classes of sphere–sphere interactions, between one sphere's motion and (1) the oscillating background flow induced by the other's motion, and (2) a geometric asymmetry induced by the other's presence. We can thus unify and explain behaviours observed in other works. Our results show how sensitive, counterintuitive and rich motility is in the parameter space of finite inertia of particles and fluid.

42 ENGINEERING↗

Theoretical investigation of the triggering of neoclassical tearing modes by transient resonant magnetic perturbations in NSTX

The extended perturbed equilibrium code asymptotic matching code is used to simulate the triggering of n = 1 neoclassical tearing modes (NTMs) by a pulsed, rotating, n = 1, resonant magnetic perturbation (RMP) in two example NSTX discharges. Although the two discharges are significantly different from one another, the results of the two sets of simulations are quite similar. The critical n = 1 RMP pulse amplitude required to trigger an n = 1 NTM is minimized when the RMP pulse rotation frequency matches the linear natural frequency of an n = 1 tearing mode, resonant within the plasma, that is metastable to an NTM. However, if there is a frequency mismatch, then the seed magnetic island chain driven at the relevant resonant surface is forced to rotate with respect to the RMP, because the RMP pulse amplitude is nowhere near sufficient to lock the island chain to the RMP. This rotation causes the critical RMP pulse amplitude required to trigger an NTM to oscillate as the RMP pulse duration is varied. Furthermore, the critical amplitude is minimized when the RMP pulse duration is such that seed island chain executes a half-integer number of rotations with respect to the pulse. All of the minima have the same value.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ignition of metals and alloys in gaseous oxygen by frictional heating

The ignition of metals and alloys have been investigated by rotating the end of a hollow cylinder against an identical stationary cylinder on a common axis (frictional heating) in gaseous oxygen. A ranking criterion that measures the resistance of metals and alloys to ignition is discussed. It consists of the power per unit area required for ignition conveniently expressed as the product of contact pressure P and average linear velocity ALV. Data are presented that demonstrate that materials that are high in nickel and copper require greater P-ALV products for ignition (more resistant to ignition) than materials that are high in iron. Aluminum and titanium alloys are shown to require the lowest P-ALV products for ignition. The effects of varying surface velocity, contact pressure, coefficient of friction, and oxygen pressure on P-ALV products required for ignition are discussed. The results indicate that P-ALV product required for ignition increases as surface velocity increases, as contact pressure decreases, and as the coefficient of friction decreases. Increasing oxygen pressure will cause the P-ALV product required for ignition to decrease at low pressures and increase at high pressures. In some cases, the relative ranking of materials based on P-ALV products will also change as conditions are varied. The data presented in this paper are from a development program to determine the merits of using such a dynamic test method for ranking materials for oxygen service. The results to date have indicated that more testing is required before this method can be made into a standard test. The data presented do show some of the major parameters affecting the ignition of materials when exposed to frictional heating. The repeatability of the data at any set of conditions have not been thoroughly investigated.

Benz, F. J.↗

Nonlinear kinetic simulation study of the ion–ion streaming instability in single- and multi-ion species plasmas

The nonlinear evolution of the ion–ion streaming instability (IISI) is studied using numerical techniques novel to this problem that afford direct insight into the evolution of the particle distributions of each species. Here, during the linear phase of the instability, we demonstrate quantitative agreement with linear kinetic theory. Subsequently, the electrostatic field generated by the IISI causes ring-like velocity distributions of ions to form that are both heated and slowed to varying degrees relative to their initial flows. Due to variation in the trapping conditions for ion species of differing charge-to-mass ratio, when flows of multiple species interact, the nonlinear evolution of each species can be starkly different: we show a case where a lighter ion species is completely stopped by a heavier ion species via the IISI alone (i.e., without collisions) and, for the first time, demonstrate how the IISI can introduce a relative flow between ion species that initially have the same flow velocities, thereby separating them.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Plasma gradient effets on doubel-probe measurements in the magnetosphere

The effects on double-probe electric field measurements induced by electron density and temperature gradients are investigated. We show that on some occasions such gradients may lead to marked spurious electric fields if the probes are assumed to lie at the same probe potential with repect to the plasma. The use of a proper bias current will decrease the magnitude of such an error. When the probes are near the plasma potential, the magnitude of these error signals, delta Epsilon, can vary as delta Epsilon approx. T(sub e)(Delta n(sub e)/n(sub e)) + 0.5 Delta T(sub e), where T(sub e) is the electron temperature, Delta n(sub e)/ n(sub e) the relative electron density variation between the two sensors, and Delta T(sub e) the electron temperature difference between the two sensors. This not only implies that the error signals will increase linearly with the density variations but also that such signatures grow with Delta T(sub e) i.e., such effects are 10 times larger in a 10-eV plasma than in a 1-eV plasma. This type of error is independent of the probe separation distance provided the gradient scale length is much larger than the distance. The largest errors occur when the probes are near to the plasma potential. During the crossing of a small structure (e.g, a double layer) the error signal appears as a bipolar signature. Our analysis shows that errors in double-probe measurements caused by plasma gradients are not significant at large scale (much greater than 1 km) plasma boundaries, and may only be important in cases where small-scale (less than 1 km), internal gradient structures exist. Bias currents tailored for each plasma parameter regime (i.e., variable bias current) would improve the double-probe response to gradient effects considerably.

Laakso, Harri↗

Description, Usage, and Validation of the MVL-15 Modified Vortex Lattice Analysis Capability

MVL-15 is the most recent version of the Modified Vortex-Lattice (MVL) code developed within the Aerodynamics Systems Analysis Branch (ASAB) at NASA LaRC. The term "modified" refers to the primary modification of the core vortex-lattice methodology: inclusion of viscous aerodynamics tables that are linked to the linear solution via iterative processes. The inclusion of the viscous aerodynamics inherently converts the MVL-15 from a purely analytic linearized method to a semi-empirical blend which retains the rapid execution speed of the linearized method while empirically characterizing the section aerodynamics at all spanwise lattice points. The modification provides a means to assess non-linear effects on lift that occur at angles of attack near stall, and provides a means to determine the drag associated with the application of design strategies for lift augmentation such as the use of flaps or blowing. The MVL-15 code is applicable to the analyses of aircraft aerodynamics during cruise, but it is most advantageously applied to the analysis of aircraft operating in various high-lift configurations. The MVL methodology has been previously conceived and implemented; the initial concept version was delivered to the ASAB in 2001 (van Dam, C.), subsequently revised (Gelhausen, P. and Ozoroski, T. 2002 / AVID Inc., Gelhausen, P., and Roberts, M. 2004), and then overhauled (Ozoroski, T., Hahn, A. 2008). The latest version, MVL-15 has been refined to provide analysis transparency and enhanced to meet the analysis requirements of the Environmentally Responsible Aviation (ERA) Project. Each revision has been implemented with reasonable success. Separate applications of the methodology are in use, including a similar in-house capability, developed by Olson, E. that is tailored for structural and acoustics analyses. A central premise of the methodology is that viscous aerodynamic data can be associated with analytic inviscid aerodynamic results at each spanwise wing section, thereby providing a pathway to map viscous data to the inviscid results. However, a number of factors can sidetrack the analysis consistency during various stages of this process. For example, it should be expected that the final airplane lift curve and drag polar results depend strongly on the geometry and aerodynamics of the airfoil section; however, flap deflections and flap chord extensions change the local reference geometry of the input airfoil, the airplane wing, the tabulated non-dimensional viscous aerodynamics, and the spanwise links between the linear and the viscous aerodynamics. These changes also affect the bound circulation and therefore, calculation and integration of the induced angle of attack and induced drag. MVL-15 is configured to ensure these types of challenges are properly addressed. This report is a comprehensive manual describing the theory, use, and validation of the MVL-15 analysis tool. Section 3 summarizes theoretical, procedural, and characteristic features of MVL-15, and includes a list of the files required to setup, execute, and summarize an analysis. Section 4, Section 5, Section 6, and Section 7 combine to comprise the User's Guide portions of this report. The MVL-15 input and output files are described in Section 4 and Section 5, respectively; the descriptions are supplemented with example files and information about the file formats, parameter definitions, and typical parameter values. Section 6 describes the Wing Geometry Setup Utility and the 2d-Variants Utility files that simplify and assist setting up a consistent set of MVL-15 geometry and aerodynamics input parameters and input files. Section 7 describes the use of the 3d-Results Presentation Utility file that can be used to automatically create summary tables and charts from the MVL-15 output files. Section 8 documents the Validation Results of an extensive and varied validation test matrix, including results of an airplane analysis representative of the ERA Program. A start-to-finish example of the airplane analysis procedure is described in Section 7.

Ozoroski, Thomas A.↗

Design and Evaluation of a Dynamic Programming Flight Routing Algorithm Using the Convective Weather Avoidance Model

The optimization of traffic flows in congested airspace with varying convective weather is a challenging problem. One approach is to generate shortest routes between origins and destinations while meeting airspace capacity constraint in the presence of uncertainties, such as weather and airspace demand. This study focuses on development of an optimal flight path search algorithm that optimizes national airspace system throughput and efficiency in the presence of uncertainties. The algorithm is based on dynamic programming and utilizes the predicted probability that an aircraft will deviate around convective weather. It is shown that the running time of the algorithm increases linearly with the total number of links between all stages. The optimal routes minimize a combination of fuel cost and expected cost of route deviation due to convective weather. They are considered as alternatives to the set of coded departure routes which are predefined by FAA to reroute pre-departure flights around weather or air traffic constraints. A formula, which calculates predicted probability of deviation from a given flight path, is also derived. The predicted probability of deviation is calculated for all path candidates. Routes with the best probability are selected as optimal. The predicted probability of deviation serves as a computable measure of reliability in pre-departure rerouting. The algorithm can also be extended to automatically adjust its design parameters to satisfy the desired level of reliability.

Ng, Hok K.↗

Performance Analysis and Limiting Parameters of Cross-flow Membrane-based Liquid-desiccant Air Dehumidifiers

We report that to dehumidify a humid air stream, existing air conditioning (AC) systems substantially overcool the outdoor humid air below its dew point, thereby significantly reducing energy efficiency. Directly capturing humidity, membrane-based liquid-desiccant dehumidification systems separate sensible and latent cooling (SSLC) loads and thus offer a promising pathway for a high-performance AC solution. Design of an energy-efficient SSLC-AC system, however, rests largely on detailed understating of the dehumidification process. While some studies have identified the dehumidification process mainly depends on membrane characteristics, other studies have argued that the process is limited by desiccant liquid or alternatively air thermo-hydraulic physics for typical humid climate conditions. The present study examines performance and physics of the membrane-based liquid-desiccant dehumidification process over a wide range of climate conditions through a novel 3D, two-phase, multi-species CFD model. Decoupling the thermodynamic and hydraulic effects, the study reveals that the dehumidification rate is a linear function of the water vapor pressure potential ($J=α ΔP$) summarizing the system's thermodynamic state. The slope of the curve (i.e., α) depends on hydraulic transport characteristics of the membrane pores, air stream, and desiccant solution. More importantly, it was found that the air dehumidification process is mainly limited by the air-side transport physics for thin liquid-desiccant films and commonly used porous superhydrophobic membranes. Additionally, results show that, depending on ambient/desiccant conditions and physical dehumidifier characteristics, energy effectiveness and dehumidification rate vary from 13 to 34% and from 0.13 to 1.4 g m -2 s -1 , respectively. Therefore, the present study allows to efficiently design future SSLC-based AC systems exhibiting high performance energy metrics.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Photo-z outlier self-calibration in weak lensing surveys

Calibrating photometric redshift errors in weak lensing surveys with external data is extremely challenging. In this work, we show that both Gaussian and outlier photo-z parameters can be self-calibrated from the data alone. This comes at no cost for the neutrino masses, curvature and dark energy equation of state w0, but with a 65% degradation when both w0 and wa are varied. We perform a realistic forecast for the Vera Rubin Observatory (VRO) Legacy Survey of Space and Time (LSST) 3× 2 analysis, combining cosmic shear, projected galaxy clustering and galaxy - galaxy lensing. We confirm the importance of marginalizing over photo-z outliers. We examine a subset of internal cross-correlations, dubbed "null correlations", which are usually ignored in 3× 2 analyses. Despite contributing only ~ 10% of the total signal-to-noise, these null correlations improve the constraints on photo-z parameters by up to an order of magnitude. Using the same galaxy sample as sources and lenses dramatically improves the photo-z uncertainties too. Together, these methods add robustness to any claim of detected new Physics, and reduce the statistical errors on cosmology by 15% and 10% respectively. Finally, including CMB lensing from an experiment like Simons Observatory or CMB-S4 improves the cosmological and photo-z posterior constraints by about 10%, and further improves the robustness to systematics. To give intuition on the Fisher forecasts, we examine in detail several toy models that explain the origin of the photo-z self-calibration. Our Fisher code LaSSI (Large-Scale Structure Information), which includes the effect of Gaussian and outlier photo-z, shear multiplicative bias, linear galaxy bias, and extensions to LCDM, is publicly available at \href{https://github.com/EmmanuelSchaan/LaSSI}{https://github.com/EmmanuelSchaan/LaSSI}.

79 ASTRONOMY AND ASTROPHYSICS↗

Collisionless adiabatic afterglow

In this work, we study, by numerical and analytical means, the evolution of a collisionless plasma initiated between absorbing walls. The ensuing flow is described by rarefaction waves that travel inward from the boundaries, interact, and eventually vanish after crossing through, leading up to the asymptotic stage of the decay. Particle simulations indicate that the kinetic evolution strongly resembles one found in isentropic gas dynamics. Namely, a very gradual density profile forms in the expanding central region where the rarefaction waves interact, with an accompanying linear velocity profile. Asymptotically, the density falls off as 1/t. The density and the flux at the boundary show little variation over the period when rarefaction waves still exist. Plasma potential, on the other hand, drops quite rapidly (on the underlying ion-acoustic timescale) to less than initial electron temperature Te when over 70% of the particles still remain in the system. This is due to electron kinetics being governed by conservation of adiabatic invariant in a slowly varying potential well. Analytical model of the velocity distribution is presented to explain the simulations. The results have implications for afterglow plasmas used in material processing and also for ion-extraction devices. One property of potential interest is good uniformity of the decaying plasma that occurs after approximately one ion-acoustic time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulation and Analysis of Technology and Operational Procedures to Reduce the Combined Effects of Emissions and Contrails

The development and evaluation of concepts and technology to support future air traffic management systems require a hierarchy of models ranging from real-time simulations to extensive field evaluations. Air traffic simulation models such as Airspace Concept Evaluation System, Center Tracon Automation System, Future Air traffic management Concept Evaluation Tool and others are used to design air traffic systems balancing the conflicting objectives of maximizing safety, meeting future demands for airports and airspace and increase efficiency of traffic flows in the presence of uncertain weather. The impact of aviation emissions and contrails on climate imposes another constraint on the design of aircraft and aviation operations. The understanding of the complex interaction between physical climate system, the carbon and other greenhouse gas emissions and aviation activity can be improved by the development of integrated assessment models that include emission and climate models together with air traffic simulations. The modeling of aircraft emissions and their interaction with each other to change the concentration levels of different gasses in the atmosphere and the resulting impact of the radiative forcing on the equilibrium of the Earth's atmosphere is complex and requires the use of coupled atmosphere-ocean general circulation models together with three-dimensional models of carbon cycle and chemistry of other non-CO2 greenhouse gases. These models are computationally intensive and unsuitable for studies involving the generation of multiple scenarios. Simple emission and climate models, based on the input-output relations of linear systems, capture the fundamental emission to climate impact behavior by careful selection of key variables and their dynamics. The impact of various greenhouse gases depends on the total concentration, effect per unit change in atmospheric concentration and the spatial distribution of the gas. All these quantities are influenced by the lifetime of the gas. The impact of a greenhouse gas depends on the interval of assessment, which may vary from a few decades to a few centuries. Climate metrics are aimed at providing a common scale to compare different greenhouse gases. If the metrics are to be used as a tool in developing and evaluating aviation operations, they should be transparent and easy to apply. Global Warming Potential and Aggregate Global Temperature Potential are some of the commonly used metrics. This paper integrates a national-level air traffic simulation and optimization capability with simple climate models and carbon cycle models, and climate metrics to assess the impact of aviation on climate. The capability brings together metrics, which are useful in aviation operations together with metrics used in climate studies. The capability can be used to make trade-offs between extra fuel cost and reduction in climate impact. There is considerable uncertainty in our understanding of the radiative forcing associated with emissions and contrails. The parameters in the simulation can be used to evaluate the effect of various uncertainties in emission models and contrails. It can also be used to evaluate the impact of different decision horizons. Alternatively, the optimization results from the simulation can be used as inputs to other tools that monetize global climate impacts like the FAA's Aviation Environmental Portfolio Management Tool for Impacts.

impact of aviation on the environment↗

Data-Informed Evaluation Framework for Integrated Energy Systems: Insights from Power, Process Heat, and Hydrogen Production Applications

The multi-criteria decision analysis (MCDA) framework provides a systematic evaluation of the diverse preferences and performance metrics associated with alternative solutions. This approach is advantageous over a single-criterion methodology, which are only valid under conditions that assume ceteris paribus or an "apples-to-apples" comparison. However, selecting suitable technologies for integrated energy systems (IES) can be likened to an "apples-to-oranges" comparison, given the heterogeneous factors at stake. These factors include economics and performance parameters, geological compatibility, and environmental impacts. Consequently, past research has often employed a mixture of qualitative and quantitative criteria tailored to the specific interests of each study. While the method proves effective in handling the intricate interplay of criteria, the resulting rankings and scores can vary from study to study. This inconsistency is introduced from the use of subjectively defined thresholds and weights. As a result, decision-makers frequently find it challenging to establish clear connections between specific criteria and the resulting scores, as the transformation of criteria into ordinal scores results in a substantial loss of information. To address this challenge, we introduce a data-informed IES evaluation framework that offers comprehensive, interpretable, and traceable evaluations backed by quantifiable rationale. First, we identified key IES evaluation criteria from a decade of literature, focusing on relevant IES applications in power, process heat, and hydrogen production. We leveraged state-of-the-art cost estimates from the Idaho National Laboratory (INL) and technical data from 78 reactor designs from the International Atomic Energy Agency (IAEA) and the Organization for Economic Co-operation and Development - Nuclear Energy Agency (OECD-NEA). Lastly, we established thresholds by analyzing the mean, variance, root mean square, and slope of values across alternatives, categorizing the preferences of decision-makers into distinct utility functions, such as linear, saturating, exponential, and stepwise. Our approach yielded two main outcomes: (1) it provided consistent assessments across different stakeholder groups and (2) it visualized uncertainties in the decision-making context via comprehensive sensitivity analysis. To demonstrate the impact of our framework, we conducted case studies on 6 reactor designs (AP1000, NuScale, BWRX-300, Xe-100, eVinci, iMSR) for the three applications. Our data-driven framework proved to be highly effective in addressing heterogenous uncertainties faced by varied decision-makers? preferences and IES applications, as well as cost and technical estimates of advanced reactors.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Data-Informed Evaluation Framework for Integrated Energy Systems: Insights from Power, Process Heat, and Hydrogen Production Applications

The multi-criteria decision analysis (MCDA) framework provides a systematic evaluation of the diverse preferences and performance metrics associated with alternative solutions. This approach is advantageous over a single-criterion methodology, which are only valid under conditions that assume ceteris paribus or an "apples-to-apples" comparison. However, selecting suitable technologies for integrated energy systems (IES) can be likened to an "apples-to-oranges" comparison, given the heterogeneous factors at stake. These factors include economics and performance parameters, geological compatibility, and environmental impacts. Consequently, past research has often employed a mixture of qualitative and quantitative criteria tailored to the specific interests of each study. While the method proves effective in handling the intricate interplay of criteria, the resulting rankings and scores can vary from study to study. This inconsistency is introduced from the use of subjectively defined thresholds and weights. As a result, decision-makers frequently find it challenging to establish clear connections between specific criteria and the resulting scores, as the transformation of criteria into ordinal scores results in a substantial loss of information. To address this challenge, we introduce a data-informed IES evaluation framework that offers comprehensive, interpretable, and traceable evaluations backed by quantifiable rationale. First, we identified key IES evaluation criteria from a decade of literature, focusing on relevant IES applications in power, process heat, and hydrogen production. We leveraged state-of-the-art cost estimates from the Idaho National Laboratory (INL) and technical data from 78 reactor designs from the International Atomic Energy Agency (IAEA) and the Organization for Economic Co-operation and Development - Nuclear Energy Agency (OECD-NEA). Lastly, we established thresholds by analyzing the mean, variance, root mean square, and slope of values across alternatives, categorizing the preferences of decision-makers into distinct utility functions, such as linear, saturating, exponential, and stepwise. Our approach yielded two main outcomes: (1) it provided consistent assessments across different stakeholder groups and (2) it visualized uncertainties in the decision-making context via comprehensive sensitivity analysis. To demonstrate the impact of our framework, we conducted case studies on 6 reactor designs (AP1000, NuScale, BWRX-300, Xe-100, eVinci, iMSR) for the three applications. Our data-driven framework proved to be highly effective in addressing heterogenous uncertainties faced by varied decision-makers? preferences and IES applications, as well as cost and technical estimates of advanced reactors.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Techniques to improve maneuver stability characteristics of a nonlinear wide-body transport airplane in cruise flight

The maneuver control stability characteristics of an aircraft are a flying qualities parameter of critical importance, to ensure structural protection as well as adequate predictability to the pilot. Currently, however, maneuver stability characteristics are not uniquely addressed in the Federal Aviation Regulations (FAR) Part 25, for transport aircraft. In past transport category certification programs, the Federal Aviation Administration (FAA) has used a combination of requirements (longitudinal control, vibration and buffeting, high-speed characteristics, and out-of-trim characteristics) to ensure safe and controllable maneuver stability characteristics over a range of flight conditions and airplane configurations. Controversies exist regarding each of these regulations, however, and considerable expenditures in terms of design studies and testing time have resulted from the requirements. It is also recognized that additional engineering guidance is needed for identifying acceptable nonlinear maneuver stability characteristics, particularly as they relate to relaxed stability, highly augmented transport configurations. The current trend in large aircraft design is toward relaxed, or even negative, static margins for improved fuel efficiency. The advanced flight control systems developed for these aircraft, in many instances, have rendered current aforementioned maneuver stability criteria either too stringent or of little practical use. Current design requirements do not account for these advanced designs. The objective was to evaluate a broad spectrum of linear and nonlinear longitudinal stability characteristics to generate data for defining satisfactory and unacceptable maneuver characteristics, as defined by pilot opinion. Primary emphasis was placed on two techniques of varying column force per normal acceleration. This study was a joint venture with four pilots participating; one from NASA, one from the FAA, and two from industry.

Grantham, William D.↗

Silent and Efficient Supersonic Bi-Directional Flying Wing

The supersonic bi-directional (SBiDir) flying wing (FW) concept has a great potential to achieve low sonic boom with high supersonic aerodynamic performance due to removal of performance conflict between high speed and low speed by rotating goo in flight. This NIAC Phase 1 research has achieved three objectives: 2) prove the concept based on simulation that it can achieve very low boom with smooth Sine wave ground over-pressure signature and excellent aerodynamic efficiency; 3) conduct trade study to correlate the geometric parameters with sonic boom and aerodynamic performance for further automated design optimization in Phase II. The design methodology developed in Phase I includes three parts: 1) an advanced geometry model, which can vary airfoil meanline angle distribution to control the expansion and shock waves on the airplane surface to mitigate sonic boom and improve aerodynamic efficiency. 2) a validated CFD procedure to resolve near field flow with accurate shock strength. The sonic boom propagation from near field to far field ground is simulated by NASA NF Boom code. The surface friction drag prediction is based on fiat plate correlation adopted by Seebass and supported by the experimental study of Winter and Smith, which is on the conservative side and is more reliable than CFD RANS simulation. 3) a mission analysis tool based on Corke's model that provides design requirements and constraints of supersonic airplanes for range, payload, volume, size, weight, etc. The design mission target is a supersonic transport with cruise Mach number 1.6, 100 passengers, and 4000nm range. The trade study has several very important findings: 1) The far field ground sonic boom signature is directly related to the smoothness of the flow on the airplane surface. The meanline angle distribution is a very effective control methodology to mitigate surface shock and expansion wave strength, and mitigating compression wave coalescing by achieving smooth loading distribution chord-wise. Compared with a linear meanline angle distribution, a design using nonlinear and non-monotonic meanline angle distribution is able to reduce the sonic boom ground loudness by over 20dBP1. The design achieves sonic boom ground loudness less than 70dBP1 and aerodynamic dynamic efficiency 1/D of 8.4. 2) Decreasing sweep angle within the Mach cone will increase 1/D as well as sonic boom. A design with variable sweep from 84° at the very leading edge to 68° at the tip achieves an extraordinarily high 1/D of 10.4 at Mach number 1.6 due to the low wave drag. If no sonic boom constraint is attached, SBiDir-FW concept still has a lot of room to increase the 1/D. 3) The round leading edge and trailing edge under high sweep angle are beneficial to improve aerodynamic performance, sonic boom, and to increase volume of the airplane. 4) Subsonic performance is benefited greatly from the high slenderness of supersonic configuration after rotating goo. A design with excellent supersonic aspect ratio of 0.44, 1/D of 8.g, gives an extraordinary subsonic aspect ration of 10 and 1/D of 1g.7. Two configurations are designed in details to install internal seats, landing gears, and engine installation to demonstrate the feasibility of SBiDir-FW configuration to accommodate all the required volume for realistic airplane. Here we emphasize that the qualitative findings in Phase I are very encouraging, more important than the quantitative results. Qualitative findings give the understanding of physics and provide the path to achieve the ultimate high performance design. The promising quantitative results achieved in Phase I need to be confirmed by wind tunnel testing in Phase II and ultimately proved by flight test. The other important step forward will be made to study the rotation transition from both CFD unsteady simulation and wind tunnel testing.

Zha, Gecheng↗