SEARCH · Search NASA

Results for “monotone solution”

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.

SourceResearch connectionFrom yearThrough year

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

At least 91 records · Page 5

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Lyapunov exponent as a signature of dissipative many-body quantum chaos

A distinct feature of Hermitian quantum chaotic dynamics is the exponential increase of certain out-of-time-order correlation (OTOC) functions around the Ehrenfest time with a rate given by a Lyapunov exponent. Physically, the OTOCs describe the growth of quantum uncertainty that crucially depends on the nature of the quantum motion. Here, we employ the OTOC in order to provide a precise definition of dissipative quantum chaos. For this purpose, we compute analytically the Lyapunov exponent for the vectorized formulation of the large- q limit of a q -body Sachdev-Ye-Kitaev model coupled to a Markovian bath. These analytic results are confirmed by an explicit numerical calculation of the Lyapunov exponent for several values of q ≥ 4 based on the solutions of the Schwinger-Dyson and Bethe-Salpeter equations. We show that the Lyapunov exponent decreases monotonically as the coupling to the bath increases and eventually becomes negative at a critical value of the coupling signaling a transition to a dynamics which is no longer quantum chaotic. Therefore, a positive Lyapunov exponent is a defining feature of dissipative many-body quantum chaos. The observation of the breaking of the exponential growth for sufficiently strong coupling suggests that dissipative quantum chaos may require in certain cases a sufficiently weak coupling to the environment. Published by the American Physical Society 2024

García-García, Antonio M.↗

NASA NTRS1988

Galerkin approximation for inverse problems for nonautonomous nonlinear distributed systems

An abstract framework and convergence theory is developed for Galerkin approximation for inverse problems involving the identification of nonautonomous nonlinear distributed parameter systems. A set of relatively easily verified conditions is provided which are sufficient to guarantee the existence of optimal solutions and their approximation by a sequence of solutions to a sequence of approximating finite dimensional identification problems. The approach is based on the theory of monotone operators in Banach spaces and is applicable to a reasonably broad class of nonlinear distributed systems. Operator theoretic and variational techniques are used to establish a fundamental convergence result. An example involving evolution systems with dynamics described by nonstationary quasilinear elliptic operators along with some applications are presented and discussed.

Banks, H. T.↗

NASA NTRS2017

Two-Pendulum Model of Propellant Slosh in Europa Clipper PMD Tank

The objective of this fluids analysis is to model propellant slosh for the Europa Clipper mission using a two-pendulum model, such that controls engineers can predict slosh behavior during the mission. Propellant slosh causes shifts in center of mass and exerts forces and torques on the spacecraft which, if not adequately controlled, can lead to mission failure. The two-pendulum model provides a computationally simple model that can be used to predict slosh for the Europa Clipper tank geometry. The Europa Clipper tank is cylindrical with a domed top and bottom and includes a propellant management device (PMD). Due to the lack of experimental data in low gravity environments, computational fluid dynamics (CFD) simulation results were used as 'real' slosh behavior for two propellants at three fill fractions. Key pendulum parameters were derived that allow the pendulum model's center of mass, forces, and moments to closely match the CFD data. The parameter trends were examined as a function of tank fill fraction and compared with solutions to analytic equations that describe the frequency of slosh in tanks with simple geometries. The trends were monotonic as expected, and parameters resembled analytical predictions; any differences could be explained by the specific differences in the geometry of the tank. This paper summarizes the new method developed at Goddard Space Flight Center (GSFC) for deriving pendulum parameters for two-pendulum equivalent sloshing models. It presents the results of this method and discusses the validity of the results. This analysis is at a completed stage and will be applied in the immediate future to the evolving tank geometry as Europa Clipper moves past its preliminary design review (PDR) phase.

Propellant Slosh↗

NASA NTRS2018

Parallel Monotonic Basin Hopping for Low Thrust Trajectory Optimization

Monotonic Basin Hopping has been shown to be an effective method of solving low thrust trajectory optimization problems. This paper outlines an extension to the common serial implementation by parallelizing it over any number of available compute cores. The Parallel Monotonic Basin Hopping algorithm described herein is shown to be an effective way to more quickly locate feasible solutions, and improve locally optimal solutions in an automated way without requiring a feasible initial guess. The increased speed achieved through parallelization enables the algorithm to be applied to more complex problems that would otherwise be impractical for a serial implementation. Low thrust cislunar transfers and a hybrid Mars example case demonstrate the effectiveness of the algorithm. Finally, a preliminary scaling study quantifies the expected decrease in solve time compared to a serial implementation.,

McCarty, Steven L.↗

NASA NTRS2018

Parallel Monotonic Basin Hopping for Low Thrust Trajectory Optimization

McCarty, Steven L.↗

NASA NTRS2014

Tuning Monotonic Basin Hopping: Improving the Efficiency of Stochastic Search as Applied to Low-Thrust Trajectory Optimization

Trajectory optimization methods using monotonic basin hopping (MBH) have become well developed during the past decade [1, 2, 3, 4, 5, 6]. An essential component of MBH is a controlled random search through the multi-dimensional space of possible solutions. Historically, the randomness has been generated by drawing random variable (RV)s from a uniform probability distribution. Here, we investigate the generating the randomness by drawing the RVs from Cauchy and Pareto distributions, chosen because of their characteristic long tails. We demonstrate that using Cauchy distributions (as first suggested by J. Englander [3, 6]) significantly improves monotonic basin hopping (MBH) performance, and that Pareto distributions provide even greater improvements. Improved performance is defined in terms of efficiency and robustness. Efficiency is finding better solutions in less time. Robustness is efficiency that is undiminished by (a) the boundary conditions and internal constraints of the optimization problem being solved, and (b) by variations in the parameters of the probability distribution. Robustness is important for achieving performance improvements that are not problem specific. In this work we show that the performance improvements are the result of how these long-tailed distributions enable MBH to search the solution space faster and more thoroughly. In developing this explanation, we use the concepts of sub-diffusive, normally-diffusive, and super-diffusive random walks (RWs) originally developed in the field of statistical physics.

autonomous↗

DOE OSTI2025

Concentration-Driven Slowdown of Chain-Exchange Kinetics in Block Copolymer Micelle Solutions

Here, we investigate the dependence of chain-exchange kinetics on polymer concentration for spherical micelles prepared from a polystyrene-b-poly(ethylene-alt-propylene) (SEP) diblock copolymer ( M n = 75 kDa, Ð = 1.05, and f PS = 0.18) in the PEP-selective solvent squalane. Solutions with concentrations of 15 vol% or more adopt a body-centered-cubic (BCC) packing of micelles, or a dense liquid-like packing (LLP), as confirmed by small-angle X-ray scattering. Concurrently, the mean aggregation number of the micelles grows steadily from about 60 in dilute solution to nearly 200 by 50 vol%. Time-resolved small-angle neutron scattering (TR-SANS) analysis reveals a monotonic and substantial (i.e., more than four orders of magnitude) decrease in the rate of chain exchange as the polymer concentration increases from 1 to 50 vol%. The TR-SANS relaxation functions are analyzed, taking into account both changes in contrast (due to chain exchange) and evolution in the structure factor (due to refinement of micelle packing). The activation energy for chain pullout is extracted using an established model and can be expressed as a linear function of polymer concentration, indicating that, in addition to the enthalpic barrier for extracting the core block, the evolution of micelle characteristics with increasing concentration imposes further constraints on chain exchange. Consistent with a theory of Halperin, these results suggest that the dominant factor reflects steric penalties during core block pullout arising from a reduced interfacial area per chain and associated increased corona crowding, due to the increase in aggregation number. The rate of exchange is independent of whether the micelles adopt a BCC or LLP arrangement, confirming that the rate-limiting step is extraction of the core block.

chemical structure↗

DOE OSTI2025

A simple centrifuge cell method for ex situ quantification of electrical conductivity of slurry electrode materials

We present the design, experimental procedure, and experimental evaluation of a system for fast, simple, and ex situ characterization of electrical conductivity of slurry electrode materials. The system uses a custom-designed electrochemical cell compatible with centrifugation in a swing-bucket centrifuge. The cell features cylindrical graphite electrodes that are partially sheathed so as to expose only 2 mm of the electrode surface to the bottom region of the packed particulate pellet. Also presented is a conduction model that provides a shape factor for estimating effective conductivity. We tested aqueous solutions of carbon black (CB), activated carbon (AC), and mixtures thereof. These particles were dispersed in 0.0 and 0.5 M NaCl solutions. Measurements show that the effective conductivity initially increases linearly with pellet mass and then saturates at higher masses. Notably, CB exhibited a fivefold increase in conductivity than AC at equal pellet masses. CB/AC mixtures at a fixed pellet mass were tested with CB mass fractions of 0 to 100%. Interestingly, the mixture conductivity was found to be a non-monotonic function of CB mass fraction, with a maximum conductivity at about 60 % CB mass fraction. At this maximum, the mixture conductivity is approximately 30 % higher than pure CB. NaCl concentration in the slurry solution had no effects on conductivity. These results highlight the interactions between slurry electrode composition and compaction, offering insights for optimizing slurry electrodes. Furthermore, the system and results may also be applicable to evaluation of particulate materials (including slurries) used for Li-ion batteries, capacitive deionization, fuel cells, and flow electrodes.

Capacitive deionization↗

DOE OSTI2023

Dissolution and reprecipitation of amorphous silica in silica Rich shales induces Non-Monotonic evolution of porosity in acidic reactive environments

Advances in sustainable subsurface energy technologies are crucial for meeting our energy and resource needs for a climate-resilient future. Novel strategies to harness subsurface shale reservoirs for recovering valuable metals and for enabling CO 2 storage are influenced by the morphological and mineralogical heterogeneities of these materials. In this context, delineating the interactions of highly acidic solutions such as wet supercritical CO 2 on shales with varying mineralogy is crucial to inform the stability of caprock seal for CO 2 storage and enhancements in permeability for fluid transport, reactivity, and storage. The feedback chemical effects associated with the interactions of acidic solutions on the morphologies and mineralogies of shales have not been extensively investigated. These insights are crucial for assessing temporal changes in the reactivity and the fate of the fluids in subsurface environments. In this study, we investigate the effect of 1M HCl solution on the chemistry and morphology of three different shale samples with varying carbonate, clay and silica contents. An increase in the amorphous content, from 37% to 41.3%, of silica-rich and carbonate/clay lean shale is noted due to reactions with an acidic solution which is attributed to the dissolution of Si-bearing phases such as clays, accompanied by SiO 2 precipitation. In shales bearing high content of clays and carbonates, significant increase in the pore volumes and surface areas are noted. Non-monotonic changes in the micron-scale porosity of silica rich–carbonate/clay lean (e.g., Mowry shale) are noted using in-situ X-ray microtomography experiments. Due to the initial mobilization of silica and dissolution of carbonate/clay phases, the total porosity slightly increases from 6.7% to 10.7% followed by a decrease to similar to 4% caused by SiO 2 reprecipitation. In conclusion, these findings suggest that even though silica is less reactive in acidic environments, the changes in the amorphous and crystalline content due to dissolution and reprecipitation alter the porosity and fluid flow paths.

04 OIL SHALES AND TAR SANDS↗

DOE OSTI2025

Improved Guarantees for Optimal Nash Equilibrium Seeking and Bilevel Variational Inequalities

We consider a class of hierarchical variational inequality (VI) problems that subsumes VI-constrained optimization and several other problem classes, including the optimal solution selection problem and the optimal Nash equilibrium (NE) seeking problem. Our main contribution is threefold. (i) We consider bilevel VIs with monotone and Lipschitz continuous mappings and devise a single-timescale iteratively regularized extragradient method, named IR-EG 𝚖,𝚖 . We improve the existing iteration complexity results for addressing both bilevel VI and VI-constrained convex optimization problems. (ii) Under the strong monotonicity of the outer-level mapping, we develop a method named IR-EG 𝚜,𝚖 and derive faster guarantees than those in (i). We also study the iteration complexity of this method under a constant regularization parameter. These results appear to be new for both bilevel VIs and VI-constrained optimization. (iii) To our knowledge, complexity guarantees for computing the optimal NE in nonconvex settings do not exist. Motivated by this lacuna, we consider VI-constrained nonconvex optimization problems and devise an inexactly projected gradient method, named IPR-EG, where the projection onto the unknown set of equilibria is performed using IR-EG 𝚜,𝚖 with a prescribed termination criterion and an adaptive regularization parameter. We obtain new complexity guarantees in terms of a residual map and an infeasibility metric for computing a stationary point. Here, we validate the theoretical findings using preliminary numerical experiments for computing the best and the worst NEs.

bilevel optimization↗

NASA NTRS1993

Interactive displays for trajectory planning and proximity operations

Rendezvous, docking, and other Space Station proximity operations (PROX OPS) will be conducted routinely in space. Real-time interactive visual aids and planning tools will be helpful, if not necessary, for future missions both in preflight training and on orbit. Two such displays, eivaN and Navie, are currently available for examination and human factors testing. A study was conducted in which data were collected from eight test subjects. Solution times for both devices decreased rapidly with experience. Neither fuel usage nor the number of waypoints (burns) decreased with experience. With Navie, medians of solution time and fuel consumption totaled over all subjects peaked at one of two starting points above the V-bar with monotonically decreasing values in both directions. This pattern did not appear with eivaN values. Since the docking tasks were fundamentally different with each device, and because Navie imposed more constraints on the users than eivaN did, the orbital mechanics effects had a more pronounced effect on the Navie results than on the eivaN data.

Brody, Adam R.↗

NASA NTRS2013

Silicon Heat Pipe Array

Improved methods of heat dissipation are required for modern, high-power density electronic systems. As increased functionality is progressively compacted into decreasing volumes, this need will be exacerbated. High-performance chip power is predicted to increase monotonically and rapidly with time. Systems utilizing these chips are currently reliant upon decades of old cooling technology. Heat pipes offer a solution to this problem. Heat pipes are passive, self-contained, two-phase heat dissipation devices. Heat conducted into the device through a wick structure converts the working fluid into a vapor, which then releases the heat via condensation after being transported away from the heat source. Heat pipes have high thermal conductivities, are inexpensive, and have been utilized in previous space missions. However, the cylindrical geometry of commercial heat pipes is a poor fit to the planar geometries of microelectronic assemblies, the copper that commercial heat pipes are typically constructed of is a poor CTE (coefficient of thermal expansion) match to the semiconductor die utilized in these assemblies, and the functionality and reliability of heat pipes in general is strongly dependent on the orientation of the assembly with respect to the gravity vector. What is needed is a planar, semiconductor-based heat pipe array that can be used for cooling of generic MCM (multichip module) assemblies that can also function in all orientations. Such a structure would not only have applications in the cooling of space electronics, but would have commercial applications as well (e.g. cooling of microprocessors and high-power laser diodes). This technology is an improvement over existing heat pipe designs due to the finer porosity of the wick, which enhances capillary pumping pressure, resulting in greater effective thermal conductivity and performance in any orientation with respect to the gravity vector. In addition, it is constructed of silicon, and thus is better suited for the cooling of semiconductor devices.

Yee, Karl Y.↗

NASA NTRS2014

Influence of Mixed Mode I-Mode II Loading on Fatigue Delamination Growth Characteristics of a Graphite Epoxy Tape Laminate

Mixed mode I-mode II interlaminar tests were conducted on IM7/8552 tape laminates using the mixed-mode bending test. Three mixed mode ratios, G(sub II)/G(sub T) = 0.2, 0.5, and 0.8, were considered. Tests were performed at all three mixed-mode ratios under quasi-static and cyclic loading conditions, where the former static tests were used to determine initial loading levels for the latter fatigue tests. Fatigue tests at each mixed-mode ratio were performed at four loading levels, Gmax, equal to 0.5G(sub c), 0.4G(sub c), 0.3G(sub c), and 0.2G(sub c), where G(sub c) is the interlaminar fracture toughness of the corresponding mixed-mode ratio at which a test was performed. All fatigue tests were performed using constant-amplitude load control and delamination growth was automatically documented using compliance solutions obtained from the corresponding quasi-static tests. Static fracture toughness data yielded a mixed-mode delamination criterion that exhibited monotonic increase in Gc with mixed-mode ratio, G(sub II)/G(sub T). Fatigue delamination onset parameters varied monotonically with G(sub II)/G(sub T), which was expected based on the fracture toughness data. Analysis of non-normalized data yielded a monotonic change in Paris law exponent with mode ratio. This was not the case when normalized data were analyzed. Fatigue data normalized by the static R-curve were most affected in specimens tested at G(sub II)/G(sub T)=0.2 (this process has little influence on the other data). In this case, the normalized data yielded a higher delamination growth rate compared to the raw data for a given loading level. Overall, fiber bridging appeared to be the dominant mechanism, affecting delamination growth rates in specimens tested at different load levels and differing mixed-mode ratios.

Ratcliffe, James G.↗

NASA NTRS1977

The artificial compression method for computation of shocks and contact discontinuities. I - Single conservation laws

The paper discusses the use of the artificial compression method for the computation of discontinuous solutions of a single conservation law by finite difference methods. The single conservation law has either a shock or a contact discontinuity. Any monotone finite difference scheme applied to the original equation smears the discontinuity, while the same scheme applied to the equation modified by an artificial compression flux produces steady progressing profiles. If L is any finite difference scheme in conservation form and C is an artificial compressor, the split flux artificial compression method CL is a corrective scheme: L smears the discontinuity while propagating it; C compresses the smeared transition toward a sharp discontinuity. Numerical implementation of artificial compression is described.

Harten, A.↗

DOE OSTI2022

Understanding effects of chemical complexity on helium bubble formation in Ni-based concentrated solid solution alloys based on elemental segregation measurements

Here, helium bubble formation and swelling were systematically studied in Ni-based concentrated solid solution alloys containing different numbers and types of elements. Our microscopy analysis showed that although increasing the alloy chemical complexity helps suppress bubble formation in general, there is no monotonic relationship between the bubble growth rate and the number of alloying elements. Certain elements (e.g., Fe and Pd) are more effective in suppressing bubble growth than others (e.g., Cr and Mn). Atom probe tomography was applied to accurately measure elemental segregation around bubbles, revealing unique effects of certain alloying elements on vacancy migration towards bubbles. More specifically, the high vacancy mobility via Cr sites leads to a large vacancy flux and an increased bubble size, while the high degree of atomic size mismatch introduced by Pd helps deflect vacancy flow away from bubbles and decrease the amount of swelling. The effects identified in this study provide new strategies to design concentrated solid solutions with superior resistance to swelling.

36 MATERIALS SCIENCE↗

NASA NTRS1989

Estimation on nonlinear damping in second order distributed parameter systems

An approximation and convergence theory for the identification of nonlinear damping in abstract wave equations is developed. It is assumed that the unknown dissipation mechanism to be identified can be described by a maximal monotone operator acting on the generalized velocity. The stiffness is assumed to be linear and symmetric. Functional analytic techniques are used to establish that solutions to a sequence of finite dimensional (Galerkin) approximating identification problems in some sense approximate a solution to the original infinite dimensional inverse problem.

Banks, H. T.↗

NASA NTRS1995

Nonlinear evolution of magnetic flux ropes. 2: Finite beta plasma

In this second paper on the evolution of magnetic flux ropes we study the effects of gas pressure. We assume that the energy transport is described by a polytropic relationship and reduce the set of ideal MHD equations to a single, second-order, nonlinear, ordinary differential equation for the evolution function. For this conservative system we obtain a first integral of motion. To analyze the possible motions, we use a mechanical analogue -- a one-dimensional, nonlinear oscillator. We find that the effective potential for such an oscillator depends on two parameters: the polytropic index gamma and a dimensionless quantity kappa the latter being a function of the plasma beta, the strength of the azimuthal magnetic field relative to the axial field of the flux rope, and gamma. Through a study of this effective potential we classify all possible modes of evolution of the system. In the main body of the paper, we focus on magnetic flux ropes whose field and gas pressure increase steadily towards the symmetry axis. In this case, for gamma greater than 1 and all values of kappa, only oscillations are possible. For gamma less than 1, however, both oscillations and expansion are allowed. For gamma less than 1 and kappa below a critical value, the energy of the nonlinear oscillator determines whether the flux rope will oscillate or expand to infinity. For gamma less than 1 and kappa above critical, however, only expansion occurs. Thus by increasing kappa while keeping gamma fixed (less than 1), a phase transition occurs at kappa = kappa(sub critical) and the oscillatory mode disappears. We illustrate the above theoretical considerations by the example of a flux rope of constant field line twist evolving self-similarly. For this example, we present the full numerical MHD solution. In an appendix to the paper we catalogue all possible evolutions when (1) either the magnetic field or (2) the gas pressure decreases monotonically toward the axis. We find that in these cases critical conditions can occur for gamma greater than 1. While in most cases the flux rope collapses, there are notable exceptions when, for certain ranges of kappa and gamma, collapse may be averted.

Osherovich, V. A.↗

DOE OSTI2022

Superconductivity and phase diagrams of CaK(Fe 1- x Mn x ) 4 As 4 single crystals

In order to study the effects of Mn substitution on the superconducting and magnetic ground state of CaKFe 4 As 4 (T c = 35 K), members of the CaK (Fe 1–x Mn x ) 4 As 4 series have been synthesized by high-temperature solution growth in single-crystalline form and characterized by elemental analysis, thermodynamic, and transport measurements. These measurements show that the superconducting transition temperature decreases monotonically and is finally suppressed below 1.8 K as x is increased from 0 to 0.036. For x values greater than 0.016, signatures of a magnetic transition can be detected in both thermodynamic and transport measurements in which kinklike features allow for the determination of the transition temperature T* that increases as Mn substitution increases. Here, a temperature-composition (T–x) phase diagram is constructed, revealing a half-dome of superconductivity with the magnetic transition temperature T* appearing near 26 K for x ~ 0.017 and rising slowly up to 33 K for x ~ 0.036. In addition to the creation of the T–x phase diagram for CaK(Fe 1–x Mn x ) 4 As 4 , specific-heat data are used to track the jump in specific heat at T c ; the CaK(Fe 1–x Mn x ) 4 As 4 data do not follow the scaling of ΔC p with $T$$^{3}_{c}$ as many of the other Fe-based superconducting systems do. These data suggest that, as magnetic pair breaking is present, the jump in C p for a given T c is reduced. Elastoresistivity coefficients 2m 66 and m 11 – m 12 as a function of temperature are also measured. 2m 66 and m 11 – m 12 are qualitatively similar to CaK(Fe 1–x Ni x ) 4 As 4 . This may indicate that the magnetic order in Mn-substituted system may be still the same as CaK(Fe 1–x Mn x ) 4 As 4 . Superconductivity of CaK(Fe 1–x Mn x ) 4 As 4 is also studied as a function of magnetic field. A clear change in H' c2 (T)/T c , where H' c2 (T) is dH' c2 (T)/dT, at x ~ 0.015 is observed and probably is related to change of the Fermi surface due to magnetic order. Coherence lengths and the London penetration depths are also calculated based on H c1 and H c2 data. Coherence lengths as the function of x also show the changes near x = 0.015, again consistent with Fermi-surface changes associated with the magnetic ordering seen for higher-x values.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗