Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “multiple solutions”

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 73 records · Page 4

Quantum protocol for decision making and verifying truthfulness among N ‐quantum parties: Solution and extension of the quantum coin flipping game

Abstract The authors devised a protocol that allows two parties, who may malfunction or intentionally convey incorrect information in communication through a quantum channel, to verify each other's measurements and agree on each other's results. This has particular relevance in a modified version of the quantum coin flipping game. The key innovation of the authors’ work includes the new design of a quantum coin that excludes any advantage of cheating, by which the long‐standing problem of the fair design of the game is, affirmatively, solved. Furthermore, the analysis is extended to N ‐parties communicating with each other, where multiple solutions for the verification of each player's measurement is proposed. The results in the N ‐party scenario could have particular relevance for the implementation of future quantum networks, where verification of quantum information is a necessity.

Ikeda, Kazuki↗

Optimizing and Exploring Untapped Micro-Hydro Hybrid Systems: a Multi-Objective Approach for Crystal Lake as a Large-Scale Energy Storage Solution

Increasing electricity demand and concerns about climate change and fossil fuel consumption have highlighted the importance of renewable energy resources and storage systems. This paper proposes a method for exploring untapped pumped hydro storage potentials to accommodate intermittent renewable energy generation profiles. Hourly measured data from 2022 in Benzie County, Michigan, United States, were gathered for system sizing and a thorough, realistic analysis. By employing the multi-objective grey wolf optimization algorithm, we formulated optimal sizing and energy-management strategies for three different scenarios. Unlike similar studies, the 3rd with triple objective functions (OFs) scenario aims to maximize both reliability and ecological OFs while minimizing the cost OF. It has shown promising results with multiple solutions, considering economic, environmental, and reliability factors. A case study conducted in Crystal Lake, Michigan, revealed that although Crystal Lake would function only as a micro-hydro power facility, it is a promising and huge storage unit with a substantial storage capacity of around 14.9734GWh. The system investigated is significant in the USA due to its rapid deployment capabilities, minimal construction requirements, and ease of integration with the distribution grid. The fuzzy logic method was employed to identify the best non-dominant solution among the other solutions. Furthermore, these outcomes include a notably low levelized cost of energy at 0.046147$/kWh, a robust index of reliability of 99.705%, and a significant reduction in CO₂ emissions amounting to 7.9142×10 3 tons/year, when considering the triple OFs. The paper’s methodology provides valuable insights for regions aiming to utilize renewable energy from untapped storage sources.

13 HYDRO ENERGY↗

Check-Cases for Verification of 6-Degree-of-Freedom Flight Vehicle Simulations

The rise of innovative unmanned aeronautical systems and the emergence of commercial space activities have resulted in a number of relatively new aerospace organizations that are designing innovative systems and solutions. These organizations use a variety of commercial off-the-shelf and in-house-developed simulation and analysis tools including 6-degree-of-freedom (6-DOF) flight simulation tools. The increased affordability of computing capability has made highfidelity flight simulation practical for all participants. Verification of the tools' equations-of-motion and environment models (e.g., atmosphere, gravitation, and geodesy) is desirable to assure accuracy of results. However, aside from simple textbook examples, minimal verification data exists in open literature for 6-DOF flight simulation problems. This assessment compared multiple solution trajectories to a set of verification check-cases that covered atmospheric and exo-atmospheric (i.e., orbital) flight. Each scenario consisted of predefined flight vehicles, initial conditions, and maneuvers. These scenarios were implemented and executed in a variety of analytical and real-time simulation tools. This tool-set included simulation tools in a variety of programming languages based on modified flat-Earth, round- Earth, and rotating oblate spheroidal Earth geodesy and gravitation models, and independently derived equations-of-motion and propagation techniques. The resulting simulated parameter trajectories were compared by over-plotting and difference-plotting to yield a family of solutions. In total, seven simulation tools were exercised.

Murri, Daniel G.↗

Qualitative Headspace GCMS Analysis of Lunar Regolith and Volatile Simulant Mixtures

Introduction: Future Artemis missions aim to return the volatile-bearing samples collected near lunar polar craters. We, as advanced curation scientists, are responsible for developing techniques and methodologies for preserving the integrity of returned samples and the science value those samples contain. The extent to which that preservation is possible, and the trade-offs preservation requires (e.g. monetary costs, sample volume limitations) all must be considered. Even less-than pristine volatile-bearing samples will be of tremendous value to the scientific community seeking to unravel the history of lunar surface volatiles and, more broadly, volatiles in the solar system. A sample collected on the lunar surface will experience at least five distinctive periods during which any changes here referred to as “alteration” will certainly occur at some scale: collection on the lunar surface; transportation back to Earth; long-term storage; curatorial processing; and allocation/distribution. The Planetary Exploration and Astromaterials Research Lab (PEARL) seeks to understand temperature and pressure effects on returned volatile samples by working with high-fidelity volatile-containing regolith simulants, setting the foundation for the future of cold curation. This abstract is focused on gas-surface interactions between LCROSS volatiles and readily available lunar regolith simulants. Experiments involved analyzing differences in headspace gas composition for various combinations of volatile and regolith simulants using gas chromatography/mass spectrometry (GC/MS). Experimental Procedure: The volatile simulants were chosen based on the molecules detected during the LCROSS mission.1 Stock solutions of condensed lunar volatile analytes were: methanol, 7 N ammonia in methanol, and 0.4% hydrogen sulfide in water. The regolith simulants used were primarily JSC-1A and NU-LHT-4M. Additional regolith simulant control studies were conducted with <150 μm sieved sand and KBr. A nested vial sample preparation approach separated the liquid stock solutions from the regolith, eliminating potential matrix effects between liquid and solid phases. Fifteen microliter aliquots of liquid volatile simulants were added to a 2 mL liquid GC vial and capped in atmosphere. An 18G needle punctured the 2 mL GC vial immediately before being transferred and sealed in a 20 mL GC vial containing 0-0.3 g of regolith simulant, see Figure 1. Separating the analytes ensures any changes observed in the total headspace gases is a result of gas-surface and/or gas-gas interactions. Equipment and Method: Initial GC/MS method development for the separation and identification of relevant headspace gases can be found in Amick, et al. 2023.2 The only hardware change is a different column: a TG-1701MS 30 m × 0.25 mm × 1.00 μm column. The vials were sampled at 10°C, room temperature (~25°C) and 50°C. Low temperature samples were kept in a chilled autosampler stage for at least 1 hour prior to sampling. High temperature samples were agitated at 50°C for 5 minutes immediately prior to injection onto the column. Headspace chromatograms were collected for each combination of temperature, regolith, and volatile simulant, including controls without one or both types of analytes, in triplicate. The chromatogram elution window for each analyte or significant atmospheric gas was identified using the peak mass spectrum cross-referenced with Figure 1.Picture of the nested vial set-up. a NIST MS library search. Each analyte peak was integrated after filtering the mass spectrum trace for the parent or most unique mass-to-charge ratio. For example, the mass-to-charge ratio used to identify, filter for, and integrate the carbon dioxide peak was centered around 44 m/z. Results: Figure 2 shows the integrated peak area for hydrogen sulfide in all combinations of regolith simulants, temperatures, and the addition of ammonia in methanol solution. Multiple repeat experiments with H2S and regolith simulants have confirmed hydrogen sulfide is removed from the headspace within 1 hour when exposed to JSC-1A, NU-LHT-4M and sand (not pictured). The consumption of H2S by lunar regolith simulants at different temperatures indicates surface chemistry will be an integral component in sample integrity and preservation. It is important to note that while the effect of surface chemistry on gaseous hydrogen sulfide is intriguing, further investigation into more chemically accurate regolith simulants is necessary and ongoing. The sulfur in hydrogen sulfide is in its most reduced state while off-the-shelf, terrestrially sourced JSC-1A and NU-LHT-4M are more oxidized than most lunar materials,3which likely leads to different oxidation-reduction reactions than would be expected in lunar regolith. Figure 3 shows the integrated peak area of carbon dioxide for each sample combination that contained ammonia in methanol solution. The addition of ammonia to the GC vials results in a consistent and reproducible decrease in carbon dioxide gas, even at 50°C. The effect became more pronounced when JSC-1A or NU-LHT-4M were present. This set of experiments demonstrated that the sample composition will affect the chemical and physical state of each component present. Future spectroscopic and microscopy experiments will be geared towards identifying the cause for the CO2(g) concentration decrease and the consumption of hydrogen sulfide. As cold and volatile curation scientists, this information provides necessary insight on how to appropriately handle and analyze volatile bearing returned samples, as well as predict the effect chemical composition has on the various sample phases we will analyze upon return to Earth. Unlike traditional curation of geologic materials, the molecules in a sample cannot be identified or processed using the naked eye or even an optical microscope. Volatile curation will require a combination of analytical techniques, including but not limited to highly sensitive gas and solid/condensed phase spectroscopy. This set of experiments has demonstrated the need for more detailed studies of volatile mixtures with mineralogically and geochemically analogous lunar regolith simulants to prepare for the curation of volatile-rich lunar samples from the south polar region of the Moon. References: [1] Colaprete, A., et al. (2010) Science, 330, (463-468). [2] Amick, C. L., et al.(2023) Houston, Texas, [3] Heiken, G. H., et al.(1991) (778-778)

Cecilia L Amick↗

Analysis of nonlinear features associated with radio-frequency sheaths

An analytic study is conducted of nonlinear features that arise when radio-frequency (RF) sheath boundary conditions are solved self-consistently with the electromagnetic fields in a magnetized plasma. An instability associated with RF sheath–plasma interactions is identified and shown to be closely linked with sheath–plasma resonance and the occurrence of multiple solutions. It is proposed that this instability serves as a physical mechanism that enables jumping from one solution branch to another, providing a physical basis for the prediction of an RF sheath hysteresis in the time domain. It is further shown that this instability can lead to non-intuitive features such as spontaneous symmetry breaking in an otherwise completely symmetric problem. In plane-stratified geometry, two analytically solvable examples are considered for the cases of an unbounded half-space and a bounded domain

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Robust wind farm layout optimization

Wake interactions in wind farms cause losses in annual energy production (AEP) on the order of 10%. Wind farm designers optimize the layout of the farm to mitigate wake losses, especially in the dominant site-specific wind directions. As wind turbines and wind farms grow in scale, optimization becomes more complex. Offshore wind farms regularly comprise more than 100 wind turbines and are characterized by complex boundaries due to shipping lanes, neighboring wind farms, and other constraints. Layout optimization methods are broadly split between gradient-based and gradient-free approaches. Gradient-based approaches can converge quickly and perform well for smaller, academic problems but are often sensitive to initial conditions and tuning parameters and require expert knowledge to use. On the other hand, gradient-free approaches can be more robust to problem complexities. We present a robust layout optimization approach based on a random search algorithm. The algorithm is intended for those who are not optimization experts and has few tuning parameters that need specification to achieve satisfactory results. Unlike off-the-shelf methods, which use generally available, non-domain-specific optimization routines that accept as inputs an optimization function and constraint definitions, this approach takes advantage of the relative computational costs of the different evaluations by evaluating cheaper computations first (boundary and minimum distance constraints) and running expensive AEP evaluations only if all other checks pass. Moreover, an outer genetic algorithm allows multiple solutions to evolve in parallel, enabling rapid solution development on high-performance computers. We discuss the relative ease of selecting necessary tuning parameters and demonstrate the efficacy of the genetic random search on a complex layout problem consisting of placing 70 turbines in a nonconvex and unconnected boundary region.

17 WIND ENERGY↗

Shedding light on shadow generalized parton distributions

The feasibility of extracting generalized parton distributions (GPDs) from deeply virtual Compton scattering (DVCS) data has recently been questioned because of the existence of an infinite set of so-called “shadow GPDs” (SGPDs). These SGPDs depend on the process and manifest as multiple solutions (at a fixed scale 𝑄 2 ) to the inverse problem that needs to be solved to infer GPDs from DVCS data. SGPDs therefore pose a significant challenge for extracting GPDs from DVCS data. With this motivation we study the extent to which QCD evolution can provide constraints on SGPDs. This is possible because the known classes of SGPDs begin to contribute to observables after evolution, and can then be constrained (at the input scale 𝑄$^{2}_{0}$) by data that has a finite 𝑄 2 range. The impact that SGPDs could have on determining the total angular momentum, pressure and sheer force distributions, and tomography is also discussed. We find that 𝑄 2 evolution, given empirical data over a wide range of skewness (𝜉 ) and 𝑄 2 scales, constrains the SGPDs studied at low 𝑥 and low 𝜉. If this finding generalizes to all SGPDs—and we have not been able to find any SGPDs with qualitatively different behavior—then this potentially makes extraction of GPDs from DVCS possible over a limited range of 𝑥 and 𝜉.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

A Cross-Domain Optimization Framework of PMU and Communication Placement for Multidomain Resiliency and Cost Reduction

Phasor measurement units (PMUs) play a crucial role in real-time monitoring and control of power grids. They rely on a communication network to transfer measurement data to the phasor data concentrator (PDC) for further processing and analysis. In this paper, a resilient cross-domain PMU and communication link placement method for minimizing the overall installation cost of the wide-area measurement system (WAMS) is proposed. Here, the main idea is to break down the barrier between the power grid domain and the communication domain, and consider the impact of one when design the other. The PMU placement in the power grid domain takes into account the cost of communication links by generating multiple solutions with equally minimum PMU costs for communication link placement evaluation. On the other hand, the communication link placement problem reduces the cost by customizing the routing policies based on the different roles of PMUs in grid observability. The proposed WAMS design is capable of withstanding any single component failure in the power domain (PMU failure or power branch failure) or in the communication domain (communication link failure or PDC failure). Numerical study on the IEEE 57-bus system reveals that the developed cross-domain optimization framework can significantly reduce the overall installation cost of WAMS while attaining multi-domain resiliency.

24 POWER TRANSMISSION AND DISTRIBUTION↗

On the Viability of Quantum Annealers to Solve Fluid Flows

This paper explores the suitability of upcoming novel computing technologies, particularly adiabatic annealing based quantum computers, to solve fluid dynamics problems that form a critical component of several science and engineering applications. For our experiments, we start with a well-studied one-dimensional simple flow problem, and provide a framework to convert such problems in continuum to a form amenable for deployment on such quantum annealers. Since the DWave annealer returns multiple states sampling the energy landscape of the problem, we explore multiple solution selection strategies to approximate the solution of the problem. We analyze the continuum solutions obtained both qualitatively and quantitatively as well as their sensitivities to the particular solution selection scheme.

42 ENGINEERING↗

Stability of absorption phenomena in laser-thermal propulsion

The mean flow and stability characteristics of laser absorption phenomena in a choked converging-diverging nozzle are considered. Calculations are presented for a given nozzle geometry and a series of laser intensities. Gas absorptivities corresponding to a hydrogen-cesium mixture are used with different initial temperatures being selected to investigate the effects of changes in the shape of the k-T curve. Both stability and mean flow calculations are limited to the one-dimensional case. The mean flow results show a decrease in mass flow as laser power is increased, along with increasingly steep temperature profiles. Calculations span regions of partial and complete absorption. One region is found where multiple solutions exist. Local stability results indicate the u-c characteristic is the only unstable mode in the unheated case. Laser heat addition makes this mode more unstable and also destabilizes the u-characteristic. Numerical calculations of disturbance propagation show that the instability of the u-c disturbances is counteracted by their reflection to u + c disturbances at the upstream end. The growth of the u-disturbances is localized in regions where the temperature profile is steep and they are damped in other regions. The increasing destabilization that is observed with increased laser power is probably the reason for difficulty in obtaining converged mean flow solutions at high laser intensities.

Merkle, C. L.↗

Squeeze-film dampers for turbomachinery stabilization

A technique for investigating the stability and damping present in centrally preloaded radially symmetric multi-mass flexible rotor bearing systems is presented. In general, one needs to find the eigenvalues of the linearized perturbation equations, though zero frequency stability maps may be found by solving as many simultaneous non-linear equations as there are dampers; and in the case of a single damper, such maps may be found directly, regardless of the number of degrees of freedom. The technique is illustrated for a simple symmetric four degree of freedom flexible rotor with an unpressurized damper. This example shows that whereas zero frequency stability maps are likely to prove to be a simple way to delineate multiple solution possibilities, they do not provide full stability information. Further, particularly for low bearing parameters, the introduction of an unpressurized squeeze film damper may promote instability in an otherwise stable system.

Mclean, L. J.↗

Breakdown of the conservative potential equation

The conservative full-potential equation is used to study transonic flow over five airfoil sections. The results of the study indicate that once shock waves are present in the flow, the qualitative approximation is different from that observed with the Euler equations. The difference in behavior of the potential eventually leads to multiple solutions.

Salas, M. D.↗

Breakdown of the conservative potential equation

The conservative full-potential equation is used to study transonic flow over five airfoil sections. The results of the study indicate that once shock are present in the flow, the qualitative approximation is different from that observed with the Euler equations. The difference in behavior of the potential eventually leads to multiple solutions.

Salas, M. D.↗

Calculation of symmetric and asymmetric vortex seperation on cones and tangent ogives based on discrete vortex models

An inviscid discrete vortex model, with newly derived expressions for the tangential velocity imposed at the separation points, is used to investigate the symmetric and asymmetric vortex separation on cones and tangent ogives. The circumferential locations of separation are taken from experimental data. Based on a slender body theory, the resulting simultaneous nonlinear algebraic equations in a cross-flow plane are solved with Broyden's modified Newton-Raphson method. Total force coefficients are obtained through momentum principle with new expressions for nonconical flow. It is shown through the method of function deflation that multiple solutions exist at large enough angles of attack, even with symmetric separation points. These additional solutions are asymmetric in vortex separation and produce side force coefficients which agree well with data for cones and tangent ogives.

Chin, S.↗

Effect of joint damping and joint nonlinearity on the dynamics of space structures

Analyses of the effect of linear joint characteristics on the vibrations of a free-free, three-joint beam model show that increasing joint damping increases resonant frequencies and increases modal damping but only to the point where the joint gets 'locked up' by damping. This behavior is different from that predicted by modeling joint damping as proportional damping. Nonlinear analyses of the three-joint model with cubic springs at the joints show all the classical single DOF nonlinear response behavior at each resonance of the multiple DOF system: nondoubling of response for a doubling of forcing amplitude, multiple solutions, jump behavior, and resonant frequency shifts. These properties can be concisely quantified by characteristic backbone curves, which show the locus of resonant peaks for increasing forcing amplitude.

Bowden, Mary↗

Warm O(+) polar wind and the DE-1 polar cap electron density profile

Theoretical steady state semikinetic polar wind density profiles, based on DE1/RIMS polar wind data (up to 3700 km), were obtained which agree very well with the power law electron density profile measured by the DE1/PWI for high altitudes. The polar wind is found to be O(+) dominated for the full altitude range considered (up to 8 R(E)). Multiple solutions are obtained for various combinations of base altitude ion temperatures and electron temperatures, such that the densities fit the Persoon et al. (1983) profile. For example, good fits to measured density profile are found for low base ion temperatures (5000 K) and high electron temperatures (9000 K), and also for unheated H(+) and O(+)(3000 K) with electron temperatures of 11,000 K. Below 2.8 R(E) the theoretical polar wind density deviates somewhat from the r exp -3.85 power law. It is concluded that this theoretical polar wind density profile, with a sum of base electron and ion temperatures of 14,000 K, yields a close match with the measured DE-1 electron density profile.

Ho, C. W.↗

Measuring attitude with a gradiometer

This paper explores using a gravity gradiometer to measure the attitude of a satellite, given that the gravity field is accurately known. Since gradiometers actually measure a combination of the gradient and attitude rate and acceleration terms, the answer is far from obvious. The paper demonstrates that it can be done and at microradian accuracy. The technique employed is dynamic estimation, based on the momentum biased Euler equations. The satellite is assumed nominally planet pointed, and subject to control, gravity gradient, and partly radom drag torques. The attitude estimator is unusual. While the standard method of feeding back measurement residuals is used, the feedback gain matrix isn't derived from Kalman theory. instead, it's chosen to minimize a measure of the terminal covariance of the error in the estimate. This depends on the gain matrix and the power spectra of all the process and measurement noises. An integration is required over multiple solutions of Lyapunov equations.

Sonnabend, David↗

Optimal Estimation of Rain-Rate Profiles from Single-Frequency Radar Echoes

The significant ambiguities inherent in the determination of a particular vertical rain intensity profile from a given time profile of radar echo powers measured by a downward-looking (spaceborne or airborne) radar at a single attenuating frequency are well documented. Indeed, one already knows that by appropriately varying the parameters of the reflectivity-rain rate (Z-R) and/or attenuation-rain rate (k- R) relationships one can produce several substantially different rain-rate profiles that would produce the same radar power profile. Imposing the additional constraint that the path-averaged rain rate be a given fixed number does reduce the ambiguities but falls far short of eliminating them. While formulas to generate all mutually ambiguous rain-rate profiles from a given profile of received radar reflectivities have already been derived, there remains to be produced a quantitative measure to assess how likely each of these profiles is, what the appropriate "average" profile should be, and what the "variance" of these multiple solutions is. To do this, one needs to spell out the stochastic constraints that can allow us to make sense of the words "average" and "variance" in a mathematically rigorous way. Such a quantitative approach would be particularly well suited for such systems as the planned precipitation radar of the Tropical Rainfall Measuring Mission (TRMM). Indeed, one would then be able to use the radar reflectivities measured by the TRMM radar to estimate the rain-rate profile that would most likely have produced the measurements, as well as the uncertainty in the estimated rain rates as a function of range. Such an optimal approach is described in this paper.

Haddad, Ziad S.↗