Search NASASearch

SEARCH · Search NASA

Results for “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.

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

At least 55 records · Page 3

DUNE Phase II: scientific opportunities, detector concepts, technological solutions

The international collaboration designing and constructing the Deep Underground Neutrino Experiment (DUNE) at the Long-Baseline Neutrino Facility (LBNF) has developed a two-phase strategy toward the implementation of this leading-edge, large-scale science project. The 2023 report of the US Particle Physics Project Prioritization Panel (P5) reaffirmed this vision and strongly endorsed DUNE Phase I and Phase II, as did the European Strategy for Particle Physics. While the construction of the DUNE Phase I is well underway, this White Paper focuses on DUNE Phase II planning. DUNE Phase-II consists of a third and fourth far detector (FD) module, an upgraded near detector complex, and an enhanced 2.1 MW beam. The fourth FD module is conceived as a “Module of Opportunity”, aimed at expanding the physics opportunities, in addition to supporting the core DUNE science program, with more advanced technologies. This document highlights the increased science opportunities offered by the DUNE Phase II near and far detectors, including long-baseline neutrino oscillation physics, neutrino astrophysics, and physics beyond the standard model. It describes the DUNE Phase II near and far detector technologies and detector design concepts that are currently under consideration. A summary of key R&D goals and prototyping phases needed to realize the Phase II detector technical designs is also provided. DUNE's Phase II detectors, along with the increased beam power, will complete the full scope of DUNE, enabling a multi-decadal program of groundbreaking science with neutrinos.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS

Solute rejection by porous glass membranes. I - Hyperfiltration of sodium chloride and urea feed solutions.

Hyperfiltration of sodium chloride and urea was studied with porous glass membranes in closed-end capillary form, to determine the effect of pressure, temperature, and concentration variations, and lifetime rejection and flux characteristics. Rejection data for sodium chloride were consistent with the functioning of the porous glass as a low-capacity ion-exchange membrane.

Ballou, E. V.

Self adaptive solution strategies: Locally bound constrained Newton Raphson solution algorithms

A summary is given of strategies which enable the automatic adjustment of the constraint surfaces recently used to extend the range and numerical stability/efficiency of nonlinear finite element equation solvers. In addition to handling kinematic and material induced nonlinearity, both pre-and postbuckling behavior can be treated. The scheme employs localized bounds on various hierarchical partitions of the field variables. These are used to resize, shape, and orient the global constraint surface, thereby enabling essentially automatic load/deflection incrementation. Due to the generality of the approach taken, it can be implemented in conjunction with the constraints of an arbitrary functional type. To benchmark the method, several numerical experiments are presented. These include problems involving kinematic and material nonlinearity, as well as pre- and postbuckling characteristics. Also included is a list of papers published in the course of the work.

Padovan, Joe

The Pursuit of K: Reflections on the Current State-of-the-Art in Stress Intensity Factor Solutions for Practical Aerospace Applications

The elastic stress intensity factor (SIF, commonly denoted as K) is the foundation of practical fracture mechanics (FM) analysis for aircraft structures. This single parameter describes the first-order effects of stress magnitude and distribution as well as the geometry of both structure/component and crack. Hence, the calculation of K is often the most significant step in fatigue analysis based on FM. This presentation will provide several reflections on the current state-of-the-art in SIF solution methods used for practical aerospace applications, including a brief historical perspective, descriptions of some recent and ongoing advances, and comments on some remaining challenges. Newman and Raju made significant early contributions to practical structural analysis by developing closed-form SIF equations for surface and corner cracks in simplified geometries, often based on empirical fits of finite element (FE) solutions. Those solutions (and others like them) were sometimes revised as new analyses were conducted or limitations discovered. The foundational solutions have exhibited striking longevity, despite the relatively "coarse" FE models employed many decades ago. However, in recent years, the accumulation of different generations of solutions for the same nominal geometry has led to some confusion (which solution is correct?), and steady increases in computational capabilities have facilitated the discovery of inaccuracies in some (not all!) of the legacy solutions. Some examples of problems and solutions are presented and discussed, including the challenge of maintaining consistency with legacy design applications. As computational power has increased, the prospect of calculating large numbers of SIF solutions for specific complex geometries with advanced numerical methods has grown more attractive. Fawaz and Andersson, for example, have been generating literally millions of new SIF solutions for different combinations of multiple cracks under simplified loading schemes using p-version FE methods. These data are invaluable, but questions remain about their practical use, because the tabular databases of key results needed to support practical life analysis can occupy gigabytes of storage for only a few classes of geometries. The prospect of using such advanced numerical methods to calculate in real time only those K solutions actually needed to support a specific crack growth analysis is also tempting, but the stark reality is that the computational cost is still so high that the approach is not practical except for specific, critical application problems. Some thoughts are offered about alternative paradigms. Compounding approaches are some of the earliest building blocks of SIF development for more complex geometries. These approaches are especially attractive because of their very low computational cost and their conceptual robustness; they are, in some ways, an intriguing contrast and complement to the brute-force numerical methods. In recent years, researchers at NRC-Canada have published remarkable results showing how compounding approaches can be used to generate accurate solutions for very difficult problems. Examples are provided of some successes--and some limitations--using this approach. These closed-form, tabulated numerical, and compounding approaches have typically been used for simple remote loading with simple load paths to the crack. However, many significant cracks occur in complex stress gradient fields. This is a job for weight function (WF) methods, where the arbitrary stress distribution on the crack plane in the corresponding uncracked body (typically determined using FE methods) is used to determine K. Several significant recent advances in WF methods and solutions are highlighted here. Fueled by advanced 3D numerical methods, many new solutions have been generated for classic geometries such as surface and corner cracks with wide ranges of geometrical validity. A new WF formulation has also be developed for part-through cracks considering the arbitrary stress gradients in all directions in the crack plane (so-called bivariant solutions). Basic WF methods have recently been combined with analytical expressions for crack plane stresses to develop a large family of accurate SIF solutions for corner, surface, and through cracks at internal or external notches with very wide ranges of shapes, sizes, acuities, and offsets. Finally, WF solutions are much faster than FE or boundary element solutions, but can still be much slower than simple closed-form solutions, especially for bivariant solutions that can require 2D numerical integration. Novel pre-integration and dynamic tabular methods have been developed that substantially increase the speed of these advanced WF solutions. The practical utility of advanced SIF methods, including both WF and direct numerical methods, is greatly enhanced if the FM life analysis can be directly and efficiently linked with digital models of the actual structure or component (e.g., FE models for stress analysis). Two recent advances of this type will be described. One approach directly interfaces the FM life analysis with the FE model of the uncracked component (including stress results). Through a powerful graphical user interface, simplified FM life models can be constructed (and visualized) directly on the component model, with the computer collecting the geometry and stress gradient information needed for the life calculation. An even more powerful paradigm uses expert logic to automatically build an optimum simple fracture model at any and every desired location in the component model, perform the life calculation, and even generate fatigue crack growth life contour maps, all with minimal user intervention. This paradigm has also been extended to the automatic calculation of fracture risk, considering uncertainty or variability in key input parameters such as initial crack size or location. Another new integrated approach links the engineering life analysis, the component model, and a 3D numerical fracture analysis built with the same component model to generate a table of SIF values at a specific location that can then be employed efficiently to perform the life calculation. Some attention must be given to verification and validation (V&V) issues and challenges: how good are these SIF solutions, how good is good enough, and does anyone believe the life answer? It is important to think critically about the different sources of error or uncertainty and to perform V&V in a hierarchal, building-block manner. Some accuracy issues for SIF solutions, for example, may actually involve independent material behavior issues, such as constraint loss effects for crack fronts near component surfaces, and can be a source of confusion. Recommendations are proposed for improved V&V approaches. This presentation will briefly but critically survey the range of issues and advances mentioned above, with a particular view towards assembling an integrated approach that combines different methods to create practical tools for real-world design and analysis problems. Examples will be selectively drawn from the recent literature, from recent enhancements in the NASGRO and DARWIN computer codes, and from previously unpublished research

CraigMcClung, R.

Gravity Induced Formation of Concentration Gradients in Supersaturated Binary Solutions

Experimental and theoretical studies of the formation of solute concentration gradient in supersaturated binary solutions in a gravitational field were carried out. The formation of solute concentration gradient was associated with the gravity induced redistribution of subcritical solute clusters. The birth-death process of the new solute-rich phase domains (subcritical solute clusters) was described in terms of the time-dependent Ginzburg Landau model developed for metastable state relaxation in binary (solute + solvent) non-critical solutions in the presence of a gravitational field. A new mathematical Ansatz was developed for solution of the model equations. This Ansatz has allowed to approach for the first time the following important problems: (1) Microstructure of solute distribution inside of the subcritical solute clusters. The analytical results obtained demonstrate that solute inside of the subcritical solute clusters is heterogeneously distributed with a spatially periodic structure. (2) Macrostructure of the solute subcritical clusters distribution in a gravitational field. The subcritical solute clusters are found to be distributed heterogeneously in a gravitational field. This heterogeneity, which is due to the heterogeneous birth-death process of the subcritical solute clusters in a gravitational field, initiates a noticeable solute concentration gradient in vertical columns of supersaturated binary solutions. An analysis and comparison of theoretical results and experimental data related to the solute concentration gradient formation in a gravitational field are presented. It is also demonstrated that the critical radius of solute clusters (radius of nucleation) and induction time are gravity-dependent.

Izmailov, Alexander F.

Signatures of Thermoreversible Associations in X-ray and Neutron Scattering from Dilute Polyzwitterion Solutions

In aqueous solutions of polyzwitterions (PZs), an interplay between dipole–dipole interactions and hydration of zwitterionic groups can lead to thermoreversible associations, which have been difficult to detect in experiments. Here, in this study, we investigated dilute aqueous solutions of poly(1-(3-sulphonatopropyl)-2-vinylpyridinium) (P2VPPS) using small-angle X-ray and neutron scattering (SAXS and SANS) to probe the structure and neutron spin-echo (NSE) spectroscopy to probe dynamics. The SAXS and SANS data show that the correlation length increases with an increase in the concentration of P2VPPS for three different molecular weights. The addition of 0.1 M NaCl to one of the solutions led to almost no dependence of the correlation length on the concentration. Such a concentration dependence of the correlation length suggests the formation of clusters driven by thermoreversible associations in the solutions. The NSE measurements show that the solutions with a larger correlation length display slower relaxation, reflecting the reduced mobility of larger clusters. These results should be considered as signatures of thermoreversible associations in the solutions of P2VPPS. To establish a quantitative link between local structure (clusters) and dynamics in dilute solutions of P2VPPS, we combined a thermoreversible gelation theory for the structure of PZ solutions (Li, S.-F.; Muthukumar, M. Theory of Thermoreversible Gelation and Anomalous Concentration Fluctuations in Polyzwitterion Solutions. J. Chem. Phys. 2024, 161, 024903) with a model for the dynamics of the clusters (generalized Zimm model), developed in this work. Using such a theoretical framework, we have predicted the distribution of clusters in the solutions probed with SANS and SAXS. With the distributions, the generalized Zimm model has been used to extract the diffusion constant of the clusters and their characteristic size from the NSE data, where the latter agrees with the values estimated from the SAXS/SANS data. These findings confirm the presence of thermoreversible associations and establish a quantitative link between local structure (clusters) and dynamics in dilute solutions of P2VPPS. Furthermore, with growing interest in technological applications, this work can provide useful insights into the structural and dynamical properties of other PZ solutions.

field theory

Improved solution accuracy for TDRSS-based TOPEX/Poseidon orbit determination

Orbit determination results are obtained by the Goddard Space Flight Center (GSFC) Flight Dynamics Division (FDD) using a batch-least-squares estimator available in the Goddard Trajectory Determination System (GTDS) and an extended Kalman filter estimation system to process Tracking and Data Relay Satellite (TDRS) System (TDRSS) measurements. GTDS is the operational orbit determination system used by the FDD in support of the Ocean Topography Experiment (TOPEX)/Poseidon spacecraft navigation and health and safety operations. The extended Kalman filter was implemented in an orbit determination analysis prototype system, closely related to the Real-Time Orbit Determination System/Enhanced (RTOD/E) system. In addition, the Precision Orbit Determination (POD) team within the GSFC Space Geodesy Branch generated an independent set of high-accuracy trajectories to support the TOPEX/Poseidon scientific data. These latter solutions use the geodynamics (GEODYN) orbit determination system with laser ranging and Doppler Orbitography and Radiopositioning integrated by satellite (DORIS) tracking measurements. The TOPEX/Poseidon trajectories were estimated for November 7 through November 11, 1992, the timeframe under study. Independent assessments were made of the consistencies of solutions produced by the batch and sequential methods. The batch-least-squares solutions were assessed based on the solution residuals, while the sequential solutions were assessed based on primarily the estimated covariances. The batch-least-squares and sequential orbit solutions were compared with the definitive POD orbit solutions. The solution differences were generally less than 2 meters for the batch-least-squares and less than 13 meters for the sequential estimation solutions. After the sequential estimation solutions were processed with a smoother algorithm, position differences with POD orbit solutions of less than 7 meters were obtained. The differences among the POD, GTDS, and filter/smoother solutions can be traced to differences in modeling and tracking data types, which are being analyzed in detail.

Doll, C. E.

Comparisons of the Solutions for the Category 3--Problem 2: Cascade-Gust Interaction

Six different solutions were submitted for this benchmark problem. These were obtained using a variety of methods that can be conveniently categorized in two main groups, a nonlinear time-domain group and a linearized frequency-domain group. The first includes solutions submitted by (1) Hixon, (2) Nallasamy et. al, (3) Shieh et. al, and (4) Wang et. al, and the second includes solutions submitted by (5) Coupland and (6) Serrano et. al. Methods (1) and (2) use sixth order compact differencing schemes and the rest are essentially second order in space. With the exception of the solution submitted by Shieh et. al, all are individually discussed in great detail in the workshop proceedings. Comparisons of the submitted solutions with the benchmark solution are presented below. Due to differences in the level of solution detail provided to the author by the participants, the comparisons do not always include results from all submissions. It should be noted at the outset that, since the benchmark solution itself was numerically computed, the comparisons are somewhat subjective. In order to provide maximum latitude for the participants of the workshop, no restrictions were placed on the type of method that could be used to solve the problem. Neither was there were any stipulations to use a particular grid topology or grid density. Therefore, without a detailed study of the critical features of the computed solutions, it is not possible to make concrete statements about the relative merits of one method over another. Such a study is beyond the scope of the current exercise, especially since complete flowfield details were not provided to the author by all participants. Instead a package, containing the information about the benchmark solution (both the steady and unsteady parts of it), is included on the proceedings CD should the authors who submitted solutions for this problem wish to examine in detail the benchmark solution and compare their results to it.

Envia, Edmane

Learning and discovering multiple solutions using physics-informed neural networks with random initialization and deep ensemble

In this work we explore the capability of physics-informed neural networks (PINNs) to discover multiple solutions. Many real-world phenomena governed by nonlinear differential equations (DEs), such as fluid flow, exhibit multiple solutions under the same conditions, yet capturing this solution multiplicity remains a significant challenge. A key difficulty lies in providing appropriate initial conditions or guesses, as widely used time-marching schemes and Newton’s method are highly sensitive to these choices when solving complex computational problems. While machine learning models, particularly PINNs, have shown promise in solving DEs, their ability to capture multiple solutions remains underexplored. In this work, we propose a simple and practical approach using PINNs to learn and discover multiple solutions. We first demonstrate that PINNs, when combined with random initialization and deep ensemble method—originally developed for uncertainty quantification—can effectively uncover multiple solutions to nonlinear ordinary and partial DEs. Although training large ensembles of PINNs may appear computationally demanding, this can be done efficiently using vectorization techniques supported by modern deep learning frameworks, allowing many networks to be trained simultaneously. Our approach highlights the critical role of initialization in shaping solution diversity, addressing an often-overlooked aspect of machine learning for scientific computing. Furthermore, we propose utilizing PINN-generated solutions as initial conditions or initial guesses for conventional numerical solvers to enhance accuracy and efficiency in capturing multiple solutions. Extensive numerical experiments, including the Allen–Cahn equation and cavity flow, where our approach successfully identifies both stable and unstable solutions, validate the effectiveness of our method. These findings establish a general and efficient framework for addressing solution multiplicity in nonlinear DEs.

97 MATHEMATICS AND COMPUTING

Evaluation of the Potential for Precipitation of Solids during Storage of Non-Aluminum SNF Solutions

Non-aluminum clad spent nuclear fuels (NASNF) stored in the L-Area basin will be dissolved in H-Canyon using the 6.3D electrolytic dissolver. The solutions will be stored in either the hot or warm canyon until the preparation of a sludge batch for the Defense Waste Processing Facility. Spent nuclear fuel solutions could be stored for 1-2 years before transfer to the H-Area Tank Farm depending on the interval between sludge batches. The solution level in the storage tanks will be maintained; therefore, precipitation of solids due to evaporation is not an issue. However, the precipitation of solids from completely dissolved SNF due to solution instabilities has been observed during intermediate storage of solutions generating hydrated oxides.The presence of fissile material in these solids is generally associated with zirconium molybdate, which is known to act as a host lattice for Pu and can carry the actinides upon precipitation. The formation of zirconium molybdate solids which carry fissile material is a potential concern for the storage of NASNF solutions. To address this concern, the Savannah River National Laboratory performed a literature review to identify knowledge gaps which may require experimental work to determine if the formation of solids is a concern during storage of these solutions. Based on the literature review, the precipitation of zirconium molybdate solids from the Campaign 1 NASNF solutions during intermediatestorage is expected. This conclusion is supported by the identification of zirconium molybdate solids found on the H-Canyon 6.1D Dissolver MK-12 insert spacer. The formation of the zirconium molybdate solids is attributed to hydrolysis and radiolytic processes in the nitric acid solution. As the molybdate solids form, U and Pu can substitute for Zr in the crystal lattice resulting in co-precipitation. Generally, the Pu substitutes directly into the crystal lattice during precipitation while the U associated with the molybdate solids more likely absorbs from the solution. The U in the NASNF solutions is present as uranyl nitrate, a 2+ cation which will not substitute as easily into the molybdate crystal lattice for the Zr 4+ ion.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W

Uncertainty estimation of bifurcated solutions in the Rayleigh–Bénard problem for advanced nuclear reactors applications

Multiphysics models of nuclear reactors frequently comprise nonlinear systems of equations. The nonlinear nature of these models could lead to solution bifurcations, where a small change in a certain parameter, e.g., the thermophysical properties of the coolant, can lead to a sudden change in the system’s behavior. At the point in parameter space where this happens, called a critical point, the Jacobian matrix of the model’s nonlinear operator becomes singular potentially permitting multiple solutions to coexist. In this paper, we perform uncertainty estimation (UE) in a parameter range that includes bifurcated solutions within the context of Rayleigh–Bénard problem. We perform this analysis assuming uncertain temperature difference, and tilt angle for the iterative solution algorithm with a unit Prandtl number (Pr = 1). Also, we perform this analysis under uncertain thermophysical properties for both FLiBe molten salt and liquid sodium as working fluid. We deploy two approaches to compute statistical moments for the resulting distributions of selected flow-field variables. The first approach is the blind computation of the mean and the standard deviation without any consideration of solution bifurcation, while the second approach utilizes k-means clustering to cluster each branch’s solutions together and compute separate statistical moments for each branch. The statistical distributions are obtained by perturbing the selected parameters about nominal values that correspond to a solution on one of the valid branches, and that solution is used as initial guess for the iterative solution algorithm. We found that perturbation of any parameter when its nominal value is close to its critical point always leads to branch jumping, i.e., the iterations converge to a solution on a branch different from the branch of the initial guess. This produces a statistical ensemble comprised of fundamentally different solutions leading to wrong mean values and uncertainty estimates, whereas clustering provides an efficient way to deal with this type of computation. This work is important for developing Gen IV nuclear systems because many of these systems rely on natural convection for cooling especially in accident conditions.

97 - MATHEMATICS AND COMPUTING

Supersaturated Electrolyte Solutions: Theory and Experiment

Highly supersaturated electrolyte solutions can be prepared and studied employing an electrodynamic levitator trap (ELT) technique. The ELT technique involves containerless suspension of a microdroplet thus eliminating dust, dirt, and container walls which normally cause heterogeneous nucleation. This allows very high supersaturations to be achieved. A theoretical study of the experimental results obtained for the water activity in microdroplets of various electrolyte solutions is based on the development of the Cahn-Hilliard formalism for electrolyte solutions. In the approach suggested the metastable state for electrolyte solutions is described in terms of the conserved order parameter omega(r,t) associated with fluctuations of the mean solute concentration n(sub 0). Parameters of the corresponding Ginzburg-Landau free energy functional which defines the dynamics of metastable state relaxation are determined and expressed through the experimentally measured quantities. A correspondence of 96-99 % between theory and experiment for all solutions studied was achieved and allowed the determination of an analytical expression for the spinodal concentration n(sub spin), and its calculation for various electrolyte solutions at 298 K. The assumption that subcritical solute clusters consist of the electrically neutral Bjerrum pairs has allowed both analytical and numerical investigation of the number-size N(sub c) of nucleation monomers (aggregates of the Bjerrum pairs) which are elementary units of the solute critical clusters. This has also allowed estimations for the surface tension Alpha, and equilibrium bulk energy Beta per solute molecule in the nucleation monomers. The dependence of these properties on the temperature T and on the solute concentration n(sub 0) through the entire metastable zone (from saturation concentration n(sub sat) to spinodal n(sub spin) is examined. It has been demonstrated that there are the following asymptotics: N(sub c), = I at spinodal concentration and N(sub c) = infinity at saturation.

Izmailov, Alexander F.

Gravitational Effects on the Morphology and Kinetics of Photodeposition of Polydiacetylene Thin Films From Monomer Solutions

The goal of this proposed work is to study gravitational effects on the photodeposition of polydiacetylene thin films from monomer solutions onto transparent substrates. Polydiacetylenes have been an extensively studied class of organic polymers because they exhibit many unusual and interesting properties, including electrical conductivity and optical nonlinearity. Their long polymeric chains render polydiacetylenes readily conducive to thin film formation, which is necessary for many applications. These applications require thin polydiacetylene films possessing uniform thicknesses, high purity, minimal inhomogeneities and defects (such as scattering centers), etc. Also, understanding and controlling the microstructure and morphology of the films is important for optimizing their electronic and optical properties. The lack of techniques for processing polydiacetylenes into such films has been the primary limitation to their commercial use. We have recently discovered a novel method for the formation of polydiacetylene thin films using photo-deposition from monomer solutions onto transparent substrates with UV light. This technique is very simple to carry out, and can yield films with superior quality to those produced by conventional methods. Furthermore, these films exhibit good third-order properties and are capable of waveguiding. We have been actively studying the chemistry of diacetylene polymerization in solution and the photo-deposition of polydiacetylene thin films from solution. It is well-known that gravitational factors such as buoyancy-driven convection and sedimentation can affect chemical and mass transport processes in solution. One important aspect of polydiacetylene thin film photodeposition in solution, relevant to microgravity science, is that heat generated by absorption of UV radiation induces thermal density gradients that under the influence of gravity, can cause fluid flows (buoyancy-driven convection). Additionally, changes in the chemical composition of the solution during polymerization may cause solutal convection. These fluid flows affect transport of material to and from the film surface and thereby affect the kinetics of the growth process. This manifests itself in the morphology of the resulting films; films grown under the influence of convection tend to have less uniform thicknesses, and can possess greater inhomogeneities and defects. Specifically, polydiacetylene films photodeposited from solution, when viewed under a microscope, exhibit very small particles of solid polymer which get transported by convection from the bulk solution to the surface of the growing film and become embedded. Even when carried out under conditions designed to minimize unstable density gradients (i.e., irradiating the solution from the top), some fluid flow still takes place (particles remain present in the films). It is also possible that defect nucleation may be occurring within the films or on the surface of the substrate; this, too, can be affected by convection (as is the case with crystal growth). Hence films grown in 1-g will, at best, still possess some defects. The objective of this proposal is to investigate, both in 1-g and in low-g, the effects of gravitational factors (primarily convection) on the dynamics of these processes, and on the quality, morphology, and properties of the films obtained.

Paley, Mark S.