Search NASA⌕ Search

SEARCH · Search NASA

Results for “weak 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 199 records · Page 11
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

A fourth-order phase-field fracture model: Formulation and numerical solution using a continuous/discontinuous Galerkin method

Modeling crack initiation and propagation in brittle materials is of great importance to be able to predict sudden loss of load-carrying capacity and prevent catastrophic failure under severe dynamic loading conditions. Second-order phase-field fracture models have gained wide adoption given their ability to capture the formation of complex fracture patterns, e.g. via crack merging and branching, and their suitability for implementation within the context of the conventional finite element method. Higher-order phase-field models have also been proposed to increase the regularity of the exact solution and thus increase the spatial convergence rate of its numerical approximation. However, they require special numerical techniques to enforce the necessary continuity of the phase field solution. In this paper, we derive a fourth-order phase-field model of fracture in two independent ways; namely, from Hamilton’s principle and from a higher-order micromechanics-based approach. The latter approach is novel, and provides a physical interpretation of the higher-order terms in the model. In addition, we propose a continuous/discontinuous Galerkin (C/DG) method for use in computing the approximate phase-field solution. This method employs Lagrange polynomial shape functions to guarantee -continuity of the solution at inter-element boundaries, and enforces the required regularity with the aid of additional variational and interior penalty terms in the weak form. Finally, the phase-field equation is coupled with the momentum balance equation to model dynamic fracture problems in hyper-elastic materials. Two benchmark problems are presented to compare the numerical behavior of the C/DG method with mixed finite element methods.

42 ENGINEERING↗
DOE OSTI2021

Muon-neutrino disappearance with multiple liquid argon time projection chambers in the Fermilab Booster neutrino beam

The Short Baseline Neutrino (SBN) program consists of three liquid argon time projection chamber (LArTPC) experiments: SBND, MicroBooNE and ICARUS, with 110 m, 470 m and 600 m baselines respectively. The detectors are located in the Booster Neutrino Beam (BNB) at Fermilab which has a peak energy around 0.7 GeV and contains predominantly muon neutrinos. The baseline and energy range of the SBN program is conducive to measuring neutrino oscillation parameters under various sterile neutrino hypotheses. Sterile neutrinos have been proposed as a possible solution to the numerous short baseline anomalies. The proposed particles must be sterile in nature such that they do not interact via the weak force, however they may undergo oscillations with the active neutrino flavours. Their existence may consequently be confirmed through measurements of the appearance and disappearance of the active flavours. The analyses presented in this thesis aimed to calculate and understand the sensitivity of the SBN program to measuring the ?µ disappearance parameters under the (3+1) sterile neutrino oscillation hypothesis. The sensitivity of SBN to measuring the ?µ disappearance sterile oscillation parameters, sin2 2?µµ, ?m2 41, was calculated through semi-exclusive joint fits of the ?µ CC 0p and ?µ CC Other reconstructed neutrino energy spectra. The first iteration used truth-level Monte Carlo (MC) events, and determined that the 5s SBN sensitivity is comparable to the 90% MINOS/MINOS+ confidence level and supersedes the 90% MiniBooNE confidence level across entire phase space. Semi-exclusive joint fits of the aforementioned sample spectra were performed between the MC and multiple mock data sets in SBND. This analysis assessed the accuracy with which the near detector can disentangle systematic from physics effects in the oscillation analysis. The result was a 5.49% discrepancy between the ICARUS Monte Carlo and mock data event rates, when the systematic constraints from the near detector fit were extrapolated to the far detector. The second iteration of the SBN sensitivity analysis involved the application of an event selection procedure developed in SBND, following the full reconstruction chain. ?µ CC 0p events were selected from sample of neutrinos with 84.5% efficiency and 84.3% purity. The sterile neutrino sensitivity was determined once more at the near detector with these samples, and was shown to be consistent with the truth-level studies.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

Statistical Analysis of the Limiting Dynamics of Two dimensional Boussinesq Turbulent Systems

In the current study, we investigate rotational and stratified turbulent systems within the framework of the Boussinesq equations. The effects of rotation and stratification give rise to a variety of dynamical regimes in atmospheric and oceanic turbulence. To develop a better understanding of the fundamental flow properties through a more straightforward mathematical framework, we commence our research with a simplified two-dimensional Boussinesq model that incorporates rotation and stratification effects. We explore two distinct limiting dynamics, considering pure rotation and pure stratification separately. We examine these intriguing dynamics through numerical investigations, applying the exact solution theory of Boussinesq equations and examining small-scale perturbations to the exact solution under rotation and stratification. We conduct an analysis of statistical quantities associated with the turbulent state variables. The results indicate that, under weak rotation at the equilibrium statistical state, the flow field exhibits vortex flow. However, strong rotation transforms the flow into vertical shear flow. We observe both downscale and upscale energy transfers with a decay rate of k −3 in the absence of external forcing and dissipation. Furthermore, we delve into the study of internal gravity wave mode interactions in the simplified two-dimensional Boussinesq system, which is a crucial aspect of geophysical turbulence due to gravity. Our findings reveal that the time series of the mode coefficients exhibit wave-like interactions between wave modes.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

A provably stable numerical method for the anisotropic diffusion equation in confined magnetic fields

We present a novel numerical method for solving the anisotropic diffusion equation in magnetic fields confined to a periodic box which is accurate and provably stable. We derive energy estimates of the solution of the continuous initial boundary value problem. A discrete formulation is presented using operator splitting in time with the summation by parts finite difference approximation of spatial derivatives for the perpendicular diffusion operator. Weak penalty procedures are derived for implementing both boundary conditions and parallel diffusion operator obtained by field line tracing. We prove that the fully-discrete approximation is unconditionally stable. Discrete energy estimates are shown to match the continuous energy estimate given the correct choice of penalty parameters. A nonlinear penalty parameter is shown to provide an effective method for tuning the parallel diffusion penalty and significantly minimises rounding errors. Several numerical experiments, using manufactured solutions, the “NIMROD benchmark” problem and a single island problem, are presented to verify numerical accuracy, convergence, and asymptotic preserving properties of the method. Finally, we present a magnetic field with chaotic regions and islands and show the contours of the anisotropic diffusion equation reproduce key features in the field.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Foldable semi-ladder polymers: novel aggregation behavior and high-performance solution-processed organic light-emitting transistors

A critical issue in developing high-performance organic light-emitting transistors (OLETs) is to balance the trade-off between charge transport and light emission in a semiconducting material. Although traditional materials for organic light-emitting diodes (OLEDs) or organic field-effect transistors (OFETs) have shown modest performance in OLET devices, design strategies towards high-performance OLET materials and the crucial structure–performance relationship remain unclear. Our research effort in developing cross-conjugated weak acceptor-weak donor copolymers for luminescent properties lead us to an unintentional discovery that these copolymers form coiled foldamers with intramolecular H-aggregation, leading to their exceptional OLET properties. An impressive external quantum efficiency (EQE) of 6.9% in solution-processed multi-layer OLET devices was achieved.

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

Invisible axion search methods

In the late 1970s, the axion was proposed as a solution to the strong CP problem, i.e., the puzzle why the strong interactions conserve parity P and the product CP of charge conjugation and parity in spite of the fact that the standard model of elementary particles as a whole violates those symmetries. The original axion was soon ruled out by laboratory experiments and astrophysical considerations, but a new version was invented that is much more weakly coupled and that evades the laboratory and astrophysical constraints. It was named the “invisible” axion. However, the axion cannot be arbitrarily weakly coupled because it is overproduced in the early Universe by vacuum realignment in the limit of vanishing coupling. The axions produced by vacuum realignment are a form of cold dark matter today. The axion provides a solution then not only to the strong CP problem but also to the dark matter problem. Various methods have been proposed to search for dark matter axions and for axions emitted by the Sun. Their implementation and improvement has led to significant constraints on the notion of an invisible axion. Even purely laboratory methods may place significant constraints on invisible axions or axionlike particles. This review discusses the various methods that have been proposed and provides theoretical derivations of their signals.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Convergence of sum-up rounding schemes for cloaking problems governed by the Helmholtz equation

In this work, we consider the problem of designing a cloak for waves described by the Helmholtz equation from an integer programming point of view. The problem can be modeled as a PDE-constrained optimization problem with integer-valued control inputs that are distributed in the computational domain. A first-discretize-then-optimize approach results in a large-scale mixed-integer nonlinear program that is in general intractable because of the large number of integer variables that arise from the discretization of the domain. Instead, we propose an efficient algorithm that is able to approximate the local infima of the underlying nonconvex infinite-dimensional problem arbitrarily close without the need to solve the discretized finite-dimensional integer programs to optimality. We optimize only the continuous relaxations of the approximations for local minima and then apply the sum-up rounding methodology to obtain integer-valued controls. If the solutions of the discretized continuous relaxations converge to a local minimizer of the continuous relaxation, then the resulting discrete-valued control sequence converges weakly \(^*\) in \(L^\infty\) to the same local minimizer. These approximation properties follow under suitable refinements of the involved discretization grids. Our results use familiar concepts arising from the analytical properties of the underlying PDE and complement previous results, derived from a topology optimization point of view.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A mixed formulation of the plane-stress problem to facilitate reuse of constitutive models in finite-element programs

Here, the plane-stress assumption can be challenging to support in a finite element program because it traditionally requires separate implementations of constitutive models than those intended for three-dimensional or two-dimensional plane-strain simulations. As a solution to this issue, this paper presents a method to solve the plane-stress problem using a mixed formulation. In this formulation, the out-of-plane strain is treated as a field variable that is solved for in addition to the standard in-plane displacement variables, in a manner that weakly enforces the condition that the out-of-plane stress is zero. The proposed formulation is non-intrusive, requiring no modifications to the constitutive models in contrast to the conventional plane-stress formulation. The proposed mixed formulation has been benchmarked against analytical solutions and numerical solutions, with good performance and accuracy.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Selective Sampling for Sensor Type Classification in Buildings

A key barrier to applying any smart technology to a building is the requirement of locating and connecting to the necessary resources among the thousands of sensing and control points, i.e., the metadata mapping problem. Existing solutions depend on exhaustive manual annotation of sensor metadata --- a laborious, costly, and hardly scalable process. To reduce the amount of manual effort required, this paper presents a multi-oracle selective sampling framework to leverage noisy labels from information sources with unknown reliability such as existing buildings, which we refer to as weak oracles, for metadata mapping. This framework involves an interactive process, where a small set of sensor instances are progressively selected and labeled for it to learn how to aggregate the noisy labels as well as to predict sensor types. Two key challenges arise in designing the framework, namely, weak oracle reliability estimation and instance selection for querying. To address the first challenge, we develop a clustering-based approach for weak oracle reliability estimation to capitalize on the observation that weak oracles perform differently in different groups of instances. For the second challenge, we propose a disagreement-based query selection strategy to combine the potential effect of a labeled instance on both reducing classifier uncertainty and improving the quality of label aggregation. We evaluate our solution on a large collection of real-world building sensor data from 5 buildings with more than 11,000 sensors of 18 different types. The experiment results validate the effectiveness of our solution, which outperforms a set of state-of-the-art baselines.

96 KNOWLEDGE MANAGEMENT AND PRESERVATION↗
DOE OSTI2024

Accelerating Multivariate Functional Approximation Computation with Domain Decomposition Techniques⋆

Modeling large datasets through Multivariate Functional Approximations (MFA) provide an elegant way to handle many visualization and scientific analysis workflows. The process necessitates scalable data partitioning methods to compute MFA representations efficiently without compromising the accuracy or continuity of the reconstructed solution. We propose a domain -decomposed method for computing the MFA with B -spline bases, which reduces the total work per task and uses a restricted Additive Schwarz (RAS) method to converge the control point data degrees -of -freedom along subdomain boundaries. We provide an in-depth analysis of the parallel approach with domain decomposition solvers, aiming to minimize local subdomain error residuals and recover high -order continuity at subdomain interfaces with appropriate choices of knot overlaps. The communication cost, determined by the overlap regions in the RAS implementation, is optimized to recover the numerical error profile of the single subdomain case. Our proposed method stands in contrast to previous methods, which typically only recover either C 0 or at best C 1 continuity for arbitrary B -spline degree expansions, or those that require post -processing to blend discontinuities in the reconstructed data. We demonstrate the effectiveness of our approach using analytical and real -world datasets in 1D, 2D, and 3D through both strong and weak scaling studies. The performance results indicate that the overall cost of computing the approximation is directly proportional to the underlying nearest -neighbor communication implementation, and is only weakly dependent on the overlap region size that determines the size of the messages. This finding underscores the efficiency and scalability of our proposed method, making it a promising solution for handling large datasets in scientific workflows.

additive Schwarz solvers↗
DOE OSTI2020

Leaching of Rare Earth Elements from Central Appalachian Coal Seam Underclays

Rare earth elements (REE) are necessary for advanced technological and energy applications. To support the emerging need, it is necessary to identify new domestic sources of REE and technologies to separate and recover saleable REE product in a safe and economical manner. Underclay rock associated with Central Appalachian coal seams and prevalent in coal utilization waste products is an alternative source of REE to hard rock ores that are mainly composed of highly refractory REE-bearing minerals. This study utilizes a suite of analytical techniques and benchtop leaching tests to characterize the properties and leachability of the coal seam underclays sampled. Laboratory bench-top and flow-through reactor leaching experiments were conducted on underclay rock powders to produce a pregnant leach solution (PLS) that has relatively low concentrations of gangue elements Al, Si, Fe, and Th and is amenable to further processing steps to recover and produce purified REE product. The leaching method described here uses a chelating agent, the citrate anion, to solubilize elements that are adsorbed, or weakly bonded to the surface of clay minerals or other mineral solid phases in the rock. The citrate PLS produced from leaching specific underclay powders contains relatively higher concentrations of REE and lower concentrations of gangue elements compared to PLS produced from sequential digestion using ammonium sulfate and mineral acids. Citrate solution leaching of underclay produces a PLS with lower concentrations of gangue elements and higher concentrations of REE than achieved with hydrochloric acid or sulfuric acid. The results provide a preliminary assessment of the types of REE-bearing minerals and potential leachability of coal seam underclays from the Central Appalachian basin.

58 GEOSCIENCES↗
DOE OSTI2025

Overcoming efficiency and cost barriers for large-area quantum dot photovoltaics through stable ink engineering

The bottom-up construction of electronics from colloidal quantum dots (CQDs) could innovate nanotechnology manufacturing through printing. However, the unstable and expensive semiconductive CQD inks make the scaling up of CQD electronics challenging. Here, in this study, we develop a strategy for engineering the solution chemistry of lead sulfide (PbS) CQD inks prepared from a low-cost direct synthesis method. By creating an iodine-rich environment in weakly coordinating solvents, we convert the iodoplumbates into functional anions, which condense into a robust surface shell. The fully charged electrostatic surface layer prevents aggregation and epitaxial fusion of CQDs, yielding stable inks. By eliminating the fusion-induced inter-band states, we print a compact CQD film with uniformity in three dimensions, flattened energy landscape and improved carrier transport. We achieved a certified efficiency of 13.40% on 0.04 cm 2 cells, with a 300-fold increase in active area, scaling up to a 12.60 cm 2 module with a certified efficiency of 10%.

14 SOLAR ENERGY↗
DOE OSTI2024

Stress-controlled medium-amplitude oscillatory shear (MAOStress) of PVA–Borax

We report the first-ever complete measurement of MAOStress material functions, which reveal that stress can be more fundamental than strain or strain rate for understanding linearity limits as a function of Deborah number. The material used is a canonical viscoelastic liquid with a single dominant relaxation time: polyvinyl alcohol (PVA) polymer solution cross-linked with tetrahydroborate (Borax) solution. We outline experimental limit lines and their dependence on geometry and test conditions. These MAOStress measurements enable us to observe the frequency dependence of the weakly nonlinear deviation as a function of stress amplitude. The observed features of MAOStress material functions are distinctly simpler than MAOStrain, where the frequency dependence is much more dramatic. The strain-stiffening transient network model was used to derive a model-informed normalization of the nonlinear material functions that accounts for their scaling with linear material properties. Moreover, we compare the frequency dependence of the critical stress, strain, and strain-rate for the linearity limit, which are rigorously computed from the MAOStress and MAOStrain material functions. While critical strain and strain-rate change by orders of magnitude throughout the Deborah number range, critical stress changes by a factor of about 2, showing that stress is a more fundamental measure of nonlinearity strength. This work extends the experimental accessibility of the weakly nonlinear regime to stress-controlled instruments and deformations, which reveal material physics beyond linear viscoelasticity but at conditions that are accessible to theory and detailed simulation.

Mechanics↗
DOE OSTI2024

Scalable Multiphysics Block Preconditioning for Low Mach Number Compressible Resistive MHD with Application to Magnetic Confinement Fusion

This study investigates multiphysics block preconditioners that are critical in devising scalable Newton–Krylov iterative solvers for longer time-scale fully implicit fluid plasma models. The specific model of interest is the visco-resistive, low Mach number, compressible magnetohydrodynamics (MHD) model. This model describes the dynamics of conducting fluids in the presence of electromagnetic fields and can be used to study aspects of astrophysical phenomena, important science and technology applications, and basic plasma physics. The specific application of interest that motivates this study is the macroscopic simulation of longer time-scale stability and disruptions of magnetic confinement fusion devices, specifically the ITER Tokamak. The computational solution of the governing balance equations for mass, momentum, heat transfer, and magnetic induction for resistive MHD systems can be extremely challenging. These difficulties arise from both the strong nonlinear, nonsymmetric coupling of fluid and electromagnetic phenomena as well as the significant range of time and length scales that the interactions of these physical mechanisms produce. To handle the range of time and spatial scales of interest, a fully implicit unstructured variational multiscale finite element formulation is employed. For the scalable solution of the Newton linearized systems, fully coupled block preconditioners are designed to leverage algebraic multigrid subsolves. In conclusion, results are presented for the strong and weak scaling of the method as well as the robustness of these techniques for a large range of Lundquist numbers.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Ultra-high SERS detection of consumable coloring agents using plasmonic gold nanostars with high aspect-ratio spikes

Solution-based SERS detection by using a portable Raman instrument has emerged as an important tool due to its simplicity, and flexibility for rapid and on-site screening of analyte molecules. However, this method has several shortcomings, including poor sensitivity especially for weak-affinity analyte molecules, where there is no close contact between the plasmonic metal surface and analyte molecule. Examples of weak-affinity molecules include pigment molecules that are commonly used as a consumable coloring agent, such as allura red (AR), and sunset yellow (SY). As high consumption of colorant agents has been shown to cause adverse effects on human health, there is a strong need to develop a simple and practical sensing system with high sensitivity for these agents. Here we present a novel, highly sensitive solution-based SERS detection method for AR, and SY by using CTAC capped gold nanostars (GNS) having different aspect ratios (GNS-2, GNS-4, and GNS-5) without utilizing any aggregating agents which can enhance SERS signal however it reduces batch to batch reproducibility. The influence of the aspect ratio of GNS on SERS properties was investigated. We have achieved a limit of detection (LOD) of AR and SY as low as 0.5 and 1 ppb, respectively by using GNS-5 with the advantages of minimal sample preparation by just mixing the analyte solution into a well plate containing GNS solution. In addition, excellent colloidal stability and reproducibility have further enhanced the applicability in real-world samples. Altogether, our results evidence that the solution-based SERS detection platform using high aspect-ratio GNS can be applied for practical application to detect pigment molecules in real samples with satisfactory results.

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

Error Estimates of Residual Minimization Using Neural Networks for Linear PDES

We propose an abstract framework for analyzing the convergence of least-squares methods based on residual minimization when feasible solutions are neural networks. With the norm relations and compactness arguments, we derive error estimates for both continuous and discrete formulations of residual minimization in strong and weak forms. The formulations cover recently developed physicsinformed neural networks based on strong and variational formulations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Predicting phase behavior in high entropy and chemically complex alloys

The interest in high entropy alloys and other metallic compounds with four or more elements at near-equiatomic ratios has drawn attention to the ability to rapidly predict phase behavior of these complex materials, particularly where existing thermodynamic data are lacking. In this paper, we discuss aspects of this from the point of view of predicting without utilizing (or fitting) experimental data. Of particular interest are heuristic approaches that provide prediction of single-phase compositions, more rigorous approaches that tackle the thermodynamics from a more fundamental point of view, and simulation approaches that provide further insight into the behaviors. This paper covers cases of all three of these, in order to examine the strengths and weaknesses of each approach, and to indicate directions where these may be utilized and improved upon. Of particular interest is moving beyond “which composition may form a solid solution,” to recognizing the importance of underlying thermodynamic realities that affect the temperature- and composition-dependent transformations of these materials.

36 MATERIALS SCIENCE↗