Search NASA⌕ Search

SEARCH · Search NASA

Results for “continuum”

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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Voltage-Dependent First-Principles Simulation of Insertion of Chloride Ions into Al/Al 2 O 3 Interfaces Using the Quantum Continuum Approximation

Experiments have shown that pitting corrosion can develop in aluminum surfaces at potentials > − 0.5 V relative to the standard hydrogen electrode (SHE). Until recently, the onset of pitting corrosion in aluminum has not been rigorously explored at an atomistic scale because of the difficulty of incorporating a voltage into density functional theory (DFT) calculations. We introduce the Quantum Continuum Approximation (QCA) which self-consistently couples explicit DFT calculations of the metal-insulator and insulator-solution interfaces to continuum Poisson-Boltzmann electrostatic distributions describing the bulk of the insulating region. By decreasing the number of atoms necessary to explicitly simulate with DFT by an order of magnitude, QCA makes the first-principles prediction of the voltage of realistic electrochemical interfaces feasible. After developing this technique, we apply QCA to predict the formation energy of chloride atoms inserting into oxygen vacancies in Al(111)/α-Al 2 O 3 (0001) interfaces as a function of applied voltage. We predict that chloride insertion is only favorable in systems with a grain boundary in the Al 2 O 3 for voltages > − 0.2 V (SHE). Here our results roughly agree with the experimentally demonstrated onset of corrosion, demonstrating QCA's utility in modeling realistic electrochemical systems at reasonable computational cost.

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

Achieving extreme light confinement in low-index dielectric resonators through quasi-bound states in the continuum

Obtaining large field enhancement in low-refractive-index dielectric materials is highly relevant to many photonic and quantum optics applications. However, confining light in these materials is challenging, owing to light leakage through coupling to continuum modes in the surrounding environment. Here, we investigate the possibility of achieving high quality factors in low-index dielectric resonators through the bound states in the continuum (BIC). Our simulations demonstrate that destructive interference between leaky modes can be achieved by tuning the geometrical parameters of the resonator arrays, leading to the emergence of quasi-BIC in resonators that have a small index contrast to the underlying substrates. The resultant large field enhancement gives rise to giant quality factors and Purcell effects. By introducing vertical mirror symmetry, the quasi-BIC can be tuned into an ideal BIC. In addition, the quasi-BIC can modify the emission patterns of the coupled emitters, rendering highly directional and focused far-field emission. These findings may provide a path for the practical implementation of photonic and quantum devices based on low-index dielectric materials.

42 ENGINEERING↗
DOE OSTI2022

The Low-redshift Lyman Continuum Survey. II. New Insights into LyC Diagnostics

The Lyman continuum (LyC) cannot be observed at the epoch of reionization (z ≳ 6) owing to intergalactic H I absorption. To identify LyC emitters (LCEs) and infer the fraction of escaping LyC, astronomers have developed various indirect diagnostics of LyC escape. Using measurements of the LyC from the Low-redshift Lyman Continuum Survey (LzLCS), we present the first statistical test of these diagnostics. While optical depth indicators based on Lyα, such as peak velocity separation and equivalent width, perform well, we also find that other diagnostics, such as the [O III ]/[O II ] flux ratio and star formation rate surface density, predict whether a galaxy is an LCE. The relationship between these galaxy properties and the fraction of escaping LyC flux suggests that LyC escape depends strongly on H I column density, ionization parameter, and stellar feedback. We find that LCEs occupy a range of stellar masses, metallicities, star formation histories, and ionization parameters, which may indicate episodic and/or different physical causes of LyC escape.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Application of the Continuum Damage Mechanics Wilshire-Cano-Stewart (WCS) Model

In this study, the applications of the continuum damage mechanics-based Wilshire-Cano-Stewart (WCS) model are explored to predict rupture time, minimum-creep-strain-rate (MCSR), damage, damage evolution, and creep deformation. Increase knowledge in manufacturing methods has pushed the limit of material science and the development of new materials. Conventional testing is required to qualify materials against creep which according to the ASME B&PV III code, 10,000+ hours of experiments are necessary for each heat before materials are put into service. This process is costly and not feasible for new materials. As an alternative, models have been employed to predict creep behaviors and reduce the amount of time necessary for material qualification. Many models have been developed to predict distinct creep behaviors and the question of which model is best remains. Amongst current models, the WCS model has emerge with the ability to predict multiple behaviors using an explicit analytical approach with the ability to predict long-term creep. In this study the novel continuum damage mechanics WCS is employed in multiple applications. The goals of the study are (a) to discuss and determine the framework of the WCS model and validated it mathematically and using parametric simulations, (b) applied the model to accelerated creep data to show the capabilities of the model with non-conventional data, and (c) applied a novel numerical method, the datum temperature method (DTM) to show the model extrapolations and interpolations capabilities with limited and reduce data sets. To accomplish these goals, data is gathered for alloy P91 and Inconel 718 to develop and post-audit validate the model. The benefits of using the WCS model is that it provides an explicit stress and temperature dependency ideal for extrapolations, the coupled equations are suitable for finite element analysis (FEA) implementation, and it follows an explicit calibration approach. The model also proves that it can be applied to accelerated testing data and using the DTM.

Cano, Jaime A↗
DOE OSTI2026

Continuum contribution to charged-current absorption of low-energy $ν_e$ on $^{40}$Ar

Accurate modeling of the absorption of tens-of-MeV $ν_e$ on $^{40}$Ar is needed to enable measurements of astrophysical neutrinos using large liquid argon time projection chamber (LArTPC) detectors, such as those planned for the Deep Underground Neutrino Experiment (DUNE). We revisit the MARLEY neutrino interaction model used in present estimates of DUNE sensitivity to supernova and solar neutrino signals. Multiple theoretical refinements are pursued, especially in the unbound continuum region of nuclear excitation energy. Inclusive charged-current neutrino-argon cross sections are calculated using a hybrid strategy. Nuclear transitions to unbound states are treated using a Hartree-Fock Continuum Random Phase Approximation (HF-CRPA) model, including forbidden contributions. Allowed transitions to low-lying discrete levels are also included using indirect measurements and approximate corrections for the momentum transfer dependence. Exclusive predictions are obtained by coupling these calculations with a statistical nuclear de-excitation model. The impact on observables of interest for DUNE and similar experiments is examined in terms of both total and differential cross sections. Our refined calculations predict a lower allowed portion of the cross section relative to the prior MARLEY model. At neutrino energies appreciably below 100 MeV, the inclusion of forbidden transitions does not fully compensate for the loss of allowed strength. For a representative neutrino burst from a galactic core-collapse supernova, our results suggest that MARLEY 1.2.0 overestimates the event yield in a DUNE-like detector by approximately 20%. However, because this overestimation is more severe at backwards angles, use of the charged-current $ν_e$-$^{40}$Ar reaction for supernova pointing may be more feasible than previously expected.

Gardiner, Steven [Fermilab]↗
DOE OSTI2023

Half‐Integer Topological Charge Polarization of Quasi‐Dirac Bound States in the Continuum

Abstract The non‐trivial polarization topology of bound states in the continuum (BICs) provides new strategies in nanophotonics. The polarization topology depends on the geometric parameters and energy‐momentum dispersion of the system and can be engineered to add specific functionalities for light molding. Herein, such a possibility is investigated by studying the topology of the polarization states associated with the optical field radiated by BICs when Dirac‐cone‐degeneracy is lifted. The opening of a pseudogap in the Dirac cone dispersion of square‐lattice dielectric photonic crystal slabs is achieved by tuning the slab thickness. First, the emergence of half‐integer topological charges without the requirement of BIC annihilation is theoretically shown, which instead occurs when in‐plane inversion symmetry is broken. Then, using spin‐to‐orbital angular momentum conversion, the theory of half‐integer topological charges mediated by BICs is demonstrated and experimentally proved. The same device is able to give rise to vortices with different orbital angular momentum depending on the way it is illuminated, thus improving the potential of optical multiplexing. In addition, the additive character of the topology‐induced phase‐vortex generation is finally demonstrated for both integer and half‐integer charges using also vortex states as input beams, which is of relevance for information delivery.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Molecularly Imprinted Polymer Sensor Empowered by Bound States in the Continuum for Selective Trace-Detection of TGF-beta

The integration of advanced materials and photonic nanostructures can lead to enhanced biodetection capabilities, crucial in clinical scenarios and point-of-care diagnostics, where simplified strategies are essential. Herein, a molecularly imprinted polymer (MIP) photonic nanostructure is demonstrated, which selectively binding to transforming growth factor-beta (TGF-β), in which the sensing transduction is enhanced by bound states in the continuum (BICs). The MIP operating as a synthetic antibody matrix and coupled with BIC resonance, enhances the optical response to TGF-β at imprinted sites, leading to an augmented detection capability, thoroughly evaluated through spectral shift and optical lever analogue readout. The validation underscores the MIP-BIC sensor capability to detect TGF-β in spiked saliva, achieving a limit of detection of 10 fM and a resolution of 0.5 pM at physiological concentrations, with a precision of two orders of magnitude above discrimination threshold in patients. The MIP tailored selectivity is highlighted by an imprinting factor of 52, showcasing the sensor resistance to interference from other analytes. The MIP-BIC sensor architecture streamlines the detection process eliminating the need for complex sandwich immunoassays and demonstrates the potential for high-precision quantification. This positions the system as a robust tool for biomarker detection, especially in real-world diagnostic scenarios.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2025

Lines of Bound States in the Continuum in Monodisperse Microsphere Structures

Incorporating epitaxially grown active semiconductor media into self-assembled microsphere crystals has remained a major challenge as the curved geometry of the spheres is incompatible with conventional planar growth and patterning techniques. As a result, active photonic devices based on such self-assembled structures, which have other desirable properties such as enabling the rapid assembly of large area systems, have been elusive. Here, we numerically demonstrate an alternative route to creating high-Q resonances in microsphere-based nanophotonic assemblies compatible with standard optical gain materials. We show that close-packed layers of monodisperse microspheres can exhibit symmetry bandgaps, regardless of their layer stacking order. Thus, this microsphere environment can be used to create lines of bound states in the continuum (BICs) of an embedded slab made of active material. Using a realistic ZnS@SiO 2 (core-shell) architecture and sphere diameters compatible with scalable colloidal assembly, we predict sharp spectral features in angle-resolved transmission that trace these lines of BICs and persist in the presence of some fabrication disorder. Altogether, in enabling the integration of unpatterned epitaxial gain slabs with large-area colloidal microsphere assemblies, this approach opens a practical path toward active, dynamically reconfigurable photonic devices that leverage BIC physics for narrow-linewidth lasing and enhanced light–matter interactions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Resolving SPARC–HSA binding kinetics with an ultrasensitive photonic sensor based on bound states in the continuum

Secreted protein acidic and rich in cysteine (SPARC) is critical in cell-matrix interactions and tissue remodeling. It influences tumor progression through its affinity for human serum albumin (HSA) - the most abundant plasma protein, which also plays a crucial role in drug delivery. Strong molecular binding leads to a dissociation constant KD in the nanomolar range. Thus, determining KD requires detecting sub-nanomolar concentrations with ultrasensitive methods. This may be crucial for elucidating the nature of SPARC-HSA binding, as their interaction remains a subject of debate. Capturing these interactions accurately requires a platform capable of resolving rapid binding kinetics at extremely low analyte concentrations. In this work, we report on a microfluidics-integrated photonic nanostructure that supports bound states in the continuum (BICs) and is optimized for studying the fast kinetics of high-affinity protein-protein interactions. The unprecedented capability of detecting sub-nanomolar concentrations allows quantifying KD between SPARC and HSA beyond the state of the art. We leverage an all-dielectric photonic crystal slab (PhCS) sustaining two BIC branches arising from gapped Dirac cone dispersion. HSA is covalently immobilized on the PhCS bonded to a PDMS microfluidic chamber. SPARC dissociation is carried out using PBS buffer (pH 7.4), ensuring complete protein release through precise control of the flow rate and continuous spectral monitoring of the BICs. The measured KD=8.2±0.8 nM confirms the strong affinity of SPARC for HSA. This study highlights the potential of BIC-based sensing as a versatile tool for investigating protein interactions. These results also have implications for the optimization of drug delivery systems and cancer treatment strategies.

Albumin↗
DOE OSTI2020

A physics-informed operator regression framework for extracting data-driven continuum models

The application of deep learning toward discovery of data-driven models requires careful application of inductive biases to obtain a description of physics which is both accurate and robust. We present here a framework for discovering continuum models from high fidelity molecular simulation data. Our approach applies a neural network parameterization of governing physics in modal space, allowing a characterization of differential operators while providing structure which may be used to impose biases related to symmetry, isotropy, and conservation form. Here, we demonstrate the effectiveness of our framework for a variety of physics, including local and nonlocal diffusion processes and single and multiphase flows. For the flow physics we demonstrate this approach leads to a learned operator that generalizes to system characteristics not included in the training sets, such as variable particle sizes, densities, and concentration.

42 ENGINEERING↗
DOE OSTI2021

Multiscale Concurrent Atomistic-Continuum (CAC) modeling of multicomponent alloys

We report strengthening in complex multicomponent systems such as solid solution alloys is controlled primarily by the dynamic interactions between dislocation lines and heterogeneously distributed solute species. Modeling of extended defect length scales in such multicomponent systems becomes prohibitively expensive, motivating the development of reduced order approaches. This work explores the application of the Concurrent Atomistic-Continuum (CAC) method to model dislocation mobility in random alloys at extended length scales. By employing recently developed average-atom interatomic potentials, the average “bulk” material response in coarse-grained regions interacts with true random solute species in the atomistic-scale domain. We demonstrate that spurious stresses in domain resolution transition regions are eliminated entirely due to the CAC formulation. Simultaneously, the key details of local stress fluctuation due to randomness in the dislocation core region are captured, and fluctuating stress smoothly decays to the long-range dislocation stress field response. Dislocation mobility calculations, for line lengths over 400 nm, are computed as a function of alloy composition in the model FeNiCr system and compared to full molecular dynamics (MD). The results capture the composition-dependent trends, while reducing degrees of freedom by nearly 40%. This approach can be readily extended to any system described by an EAM potential and facilitates the study of large-scale defect dynamics in complex solute environments to support computational alloy design.

36 MATERIALS SCIENCE↗
DOE OSTI2021

On the computational solution of vector-density based continuum dislocation dynamics models: A comparison of two plastic distortion and stress update algorithms

Continuum dislocation dynamics models of mesoscale plasticity consist of dislocation transport-reaction equations coupled with crystal mechanics equations. The coupling between these two sets of equations is such that dislocation transport gives rise to the evolution of plastic distortion (strain), while the evolution of the latter fixes the stress from which the dislocation velocity field is found via a mobility law. Earlier solutions of these equations employed a staggered solution scheme for the two sets of equations in which the plastic distortion was updated via time integration of its rate, as found from Orowan's law. In this work, we show that such a direct time integration scheme can suffer from accumulation of numerical errors. We introduce an alternative scheme based on field dislocation mechanics that ensures consistency between the plastic distortion and the dislocation content in the crystal. The new scheme is based on calculating the compatible and incompatible parts of the plastic distortion separately, and the incompatible part is calculated from the current dislocation density field. Stress field and dislocation transport calculations were implemented within a finite element based discretization of the governing equations, with the crystal mechanics part solved by a conventional Galerkin method and the dislocation transport equations by the least squares method. A simple test was first performed to show the accuracy of the two schemes for updating the plastic distortion, which shows that the solution method based on field dislocation mechanics is more accurate. This method then was used to simulate an austenitic steel crystal under uniaxial loading and multiple slip conditions. By considering dislocation interactions caused by junctions, a hardening rate similar to discrete dislocation dynamics simulation results was obtained. Finally, the simulations show that dislocations exhibit some self-organized structures as the strain is increased.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Accelerating high-order continuum kinetic plasma simulations using multiple GPUs

Kinetic plasma simulations solve the Vlasov-Poisson or Vlasov-Maxwell equations to evolve scalar-variable distribution functions in position-velocity phase space and vector-variable electromagnetic fields in configuration space. The immense computational cost of evolving high-dimensional variables, and their large number of degrees of freedom, often limits the utility of continuum kinetic simulations and presents a challenge when it comes to accurately simulating real-world physical phenomena. To address this challenge, we present techniques that accelerate and minimize the computational work required for a scalable Vlasov-Poisson solver. We show theoretical hardware compute and communication bounds for solving a fourth-order finite-volume Vlasov-Poisson system. These bounds are then used to inform and evaluate the design of performance portable algorithms for a multiple graphics processing unit (GPU) accelerated version of the Vlasov-Poisson solver VCK-CPU [1]. We demonstrate that the multi-GPU Vlasov solver implementation, VCK-GPU, simultaneously minimizes required inter-process data transfer while also being bounded by the machine network performance limits. This results in an overall strong scaling speedup per timestep of up to 40x in three-dimensional phase space (one position, two velocity coordinates) and 54x in four dimensional phase space (two position, two velocity coordinates) and a 341x increase in simulation throughput of the GPU accelerated code over the existing CPU code. The GPU code is also able to weak scale up to 256 compute nodes and 1024 GPUs. In conclusion, we demonstrate that the improved compute performance enables exploring configurations which were previously computationally infeasible, including resolving fine-scale distribution function filamentation and multi-species dynamics with realistic electron-proton mass ratios.

Continuum kinetics↗
DOE OSTI2021

Incorporating point defect generation due to jog formation into the vector density-based continuum dislocation dynamics approach

During plastic deformation of crystalline materials, point defects such as vacancies and interstitials are generated by jogs on moving dislocations. A detailed model for jog formation and transport during plastic deformation was developed within the vector density-based continuum dislocation dynamics framework. Here as a part of this model, point defect generation associated with jog transport was formulated in terms of the volume change due to the non-conservative motion of jogs. Balance equations for the vacancies and interstitials including their rate of generation due to jog transport were also formulated. A two-way coupling between point defects and dislocation dynamics was then completed by including the stress contributed by the eigen-strain of point defects. A jog drag stress was further introduced into the mobility law of dislocations to account for the energy dissipation during point defects generation. A number of test problems and a fully coupled simulation of dislocation dynamics and point defects generation and diffusion were performed. The results show that there is an asymmetry of vacancy and interstitial generation due to the different formation energies of the two types of defects. The results also show that a higher hardening rate and a higher dislocation density are obtained when the point defect generation mechanism is coupled to dislocation dynamics.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Quantifying the Influence of Charge Rate and Cathode-Particle Architectures on Degradation of Li-Ion Cells Through 3D Continuum-Level Damage Models

In this article, we develop a 3D, continuum-level damage model implemented on statistically generated LiNi 0.5 Mn 0.3 Co 0.2 O (NMC 532) secondary cathode particles. The primary motivation of the particle-level model is to inform cathode-particle design through detailed exploration of the influence of secondary and primary particle sizes on the damage predicted during operation, and determine charging profiles that reduce cathode fracture. The model considers NMC 532 secondary particles containing an agglomeration of anisotropic, randomly oriented grains. These brittle, Ni-based cathodes are prone to mechanical degradation, which reduces overall battery cycle life. The model predicts that secondary-particle fracture is primarily due to non-ideal grain interactions and high-rate charge demands. The model predicts that small secondary-particles with large grains develop significantly less damage than larger secondary particles with small grains. The model predicts most of the chemo-mechanical damage accumulates in the first few cycles. The chemo-mechanical model predicts monotonically increasing capacity fade with cycling and rate. Comparing to experimental results, the model is well suited for capturing initial capacity fade mechanisms, but additional physics is required to capture long-term capacity fade effects.

25 ENERGY STORAGE↗
DOE OSTI2021

A continuum meshfree method for sandbox-style numerical modeling of accretionary and doubly vergent wedges

Here, in this study, we introduce and utilize GEOSPH, a Lagrangian particle-based continuum Smoothed Particle Hydrodynamics (SPH) code to simulate the deformation of geologic materials in sandbox-style experiments using an elastoplastic constitutive model. Due to the meshfree nature of SPH, our proposed framework is capable of handling the large deformations characteristic of sandbox-style experiments. SPH is also able to capture strain localization and discontinuities, allowing us to resolve faulting patterns emerging in the sandbox experiments. We successfully validate our SPH implementation against the benchmark experiments of Buiter and co-workers, and then apply SPH to understand the deformational processes occurring in doubly vergent orogens resultant from the S-point experimental setup and from finitely high wedge-like indenter backstops of different inclinations. We find that the topographic slope of the prowedge results from the competing processes of vertical uplift due to the indenter and forward thrusting, while the topographic slope of the retrowedge is a product of repeated shallow slope failure. Our findings show that SPH is a promising numerical method which should be added to the current palette of computational tools used by structural geologists for solving problems relating to tectonics and crustal deformation.

58 GEOSCIENCES↗
DOE OSTI2021

Unraveling the Ion Adsorption Kinetics in Microporous Carbon Electrodes: A Multiscale Quantum-Continuum Simulation and Experimental Approach

With the understanding that sorption in porous carbon electrodes is crucial to many environmental and energy technologies, such as capacitive deionization (CDI), supercapacitor energy storage, and activated carbon filters. In each of these examples, a practical model that can describe ion electrosorption kinetics is highly desirable for accelerating material design. Here, we proposed a multiscale model to study the ion electrosorption kinetics in porous carbon electrodes by combining quantum mechanical simulations with continuum approaches. Our model integrates the Butler-Volmer (BV) equation for sorption kinetics and a continuously stirred tank reactor (CSTR) formulation with atomistic calculations of ion hydration and ion-pore interactions based on density functional theory (DFT). We validated our model experimentally by using ion mixtures in a flow-through electrode CDI device and developed an in-line UV absorption system to provide unprecedented resolution of individual ions in the separation process. We showed that the multiscale model captures unexpected experimental phenomena that cannot be explained by the traditional ion electrosorption theory. The proposed multiscale framework provides a viable approach for modeling separation processes in systems where pore sizes and ion hydration effects strongly influence the sorption kinetics, which can be leveraged to explore possible strategies for improving carbon-based and, more broadly, pore-based technologies

36 MATERIALS SCIENCE↗
DOE OSTI2026

Enhanced Second and Third Harmonic Generations in Metasurfaces Enabled through Point Group Design and Accidental Bound States in the Continuum

Here, we investigate enhanced second‑ and third‑harmonic generation (SHG and THG) in GaAs metasurfaces by combining electromagnetic point‑group designs with accidental bound states in the continuum (BICs). Using group theory, we predict the existence, suppression, and polarization selection rules of SHG and THG for metasurfaces belonging to the D 2h and C 2v point groups. In D 2h metasurfaces, inversion symmetry suppresses SHG while enabling polarization‑selective THG, which is strongly enhanced near accidental BICs with ultrahigh Q‑factors approaching 107. By breaking inversion symmetry to form C 2v metasurfaces through the introduction of nanogaps, SHG is now allowed and is dramatically enhanced near accidental BICs. Importantly, accidental BICs preserve polarization selection rules while providing a robust parameter range for high‑Q resonances, in contrast to quasi‑BICs induced by symmetry breaking. Our results demonstrate that the combination of point‑group theory and accidental‑BIC engineering offers a predictive and robust route toward highly efficient nonlinear metasurfaces.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗