Search NASA⌕ Search

SEARCH · Search NASA

Results for “Bounded domains”

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 37 records · Page 2
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Boundary Corrections for Kernel Approximation to Differential Operators

The kernel-based approach to operator approximation for partial differential equations has been shown to be unconditionally stable for linear PDEs and numerically exhibit unconditional stability for non-linear PDEs. These methods have the same computational cost as an explicit finite difference scheme but can exhibit order reduction at boundaries. In previous work on periodic domains, order reduction was addressed, yielding high-order accuracy. The issue addressed in this work is the elimination of order reduction of the kernel-based approach for a more general set of boundary conditions. Further, we consider the case of both first and second order operators. To demonstrate the theory, we provide not only the mathematical proofs but also experimental results by applying various boundary conditions to different types of equations. The results agree with the theory, demonstrating a systematic path to high order for kernel-based methods on bounded domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Conservative Numerical Schemes with Optimal Dispersive Wave Relations: Part II. Numerical Evaluations

A new energy and enstrophy conserving scheme (EEC) for the shallow water equations is proposed and evaluated using a suite of test cases over the global spherical or bounded domain. The evaluation is organized around a set of pre-defined properties: accuracy of individual operators, accuracy of the whole scheme, conservation of key quantities, control of the divergence variable, representation of the energy and enstrophy spectra, and simulation of nonlinear dynamics. The results confirm that the scheme is between the first and second order accurate, and conserves the total energy and potential enstrophy up to the time truncation errors. Here, the scheme is capable of producing more physically realistic energy and enstrophy spectra, indicating that it can help prevent the unphysical energy cascade towards the finest resolvable scales. With an optimal representation of the dispersive wave relations, the scheme is able to keep the flow close to being non-divergent, and maintain the geostrophically balanced structures with large-scale geophysical flows over long-term simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Boundary-consistent B-spline filtering schemes and application to high-fidelity simulations of turbulence

A filtering operation, based on B-spline discretizations, is introduced to target weakly growing mesh-scale oscillations that can arise in high-fidelity turbulence simulations. This is a spectral regularization that can be described using the singular values of a banded matrix operator, with the filtering strength set by a scalar- or vector-valued penalty parameter. The penalty parameter can be specified though it can also be advantageously selected to minimize the generalized cross validation (GCV) measure of distance between the pre- and post-filtered solutions. Efficient algorithms are developed to compute both the scalar and vector penalty parameters. The B-spline filter has a sharper localization to high-wavenumber than compact or explicit filters of the same stencil width and is demonstrated for solutions of the Burgers' equation, decaying Burgers' turbulence, and compressible Navier–Stokes turbulent channel flow. Furthermore, these simulations confirm the scheme's numerical stability and ability to narrowly target the high wavenumber components of numerical solutions. An advantage over finite-difference filters is that these B-spline filters are stable on bounded domains and even preserve formal order of accuracy.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Onset of cavitation and vapor bubble development over hydrophilic and hydrophobic surfaces

Cavitation, the formation of vapor bubbles as the liquid pressure is reduced below the saturated vapor pressure, often requires a substantial negative relative pressure in a pure liquid. Classical nucleation theory (CNT) provides an estimate for the rate of cavitation but there is often a disconnect between the predictions at the molecular scale compared to observations at the macroscale. We report on mesoscale simulations of cavitation based on many-body dissipative particle dynamics (mDPD), a coarse-grained molecular dynamics (MD), which bridges the two scales. A liquid layer is confined between smooth planar walls at a constant temperature, while the pressure is reduced slowly by expanding the wall-bounded domain. The wetting properties of the liquid are determined by the parameters of the interaction potentials. With hydrophilic walls, homogeneous nucleation is observed in the liquid bulk. As a bubble forms and grows, it creates a strong pressure pulse and oscillations that cause other bubbles that may have formed slightly later to collapse. For a nearly neutral wall with a contact angle close to 90 ° , heterogeneous nucleation occurs at the walls at a smaller negative pressure and generates weaker pressure oscillations. With hydrophobic walls or seed particles, heterogeneous nucleation readily occurs, where fluctuations and the merger of transient surface bubbles are significant.

Science & Technology - Other Topics↗
DOE OSTI2023

Enabling attractive-repulsive potentials in binary-collision-approximation monte-carlo codes for ion-surface interactions

Abstract Binary Collision Approximation (BCA) codes for ion-material interactions, such as SRIM, Tridyn, F-TRIDYN, and SDtrimSP, have historically been limited to screened Coulomb potentials even at low energies due to the difficulty in numerically solving the Distance of Closest Approach (DOCA) problem for attractive-repulsive potentials. Techniques such as direct n-body simulation or modifications to Newton’s method are either prohibitively costly or not guaranteed to work for all potentials. Advanced rootfinding techniques, such as companion matrix solvers, offer a solution. For many attractive-repulsive potentials, however, a companion matrix cannot be used directly, because there is no way to put the associated functions into a monomial basis form. A complementary technique is proxy rootfinding—by finding the best-fit polynomial approximant of a function, the zeros of the approximant can be guaranteed to be close to the zeros of the function. Using the Chebyshev basis and grid offers additional guarantees with regards to the quality of the approximation, the speed of convergence, and the avoidance of Runge’s phenomenon. By finding Chebyshev interpolants and using the Chebyshev-Frobenius companion matrix, the zeros of any real function on a bounded domain can be found. Here we show that using an Adaptive Chebyshev Proxy Rootfinder with Automatic Subdivision (ACPRAS) with appropriate scaling functions, numerical issues presented by attractive-repulsive potentials, including those of scale, can be handled. Using these techniques, we show that it is possible to include any physically reasonable interatomic potential in a BCA code, and to guarantee correctness of the resulting scattering angle calculations.

Materials Science↗
DOE OSTI2021

A Parallel Cut-Cell Algorithm for the Free-Boundary Grad--Shafranov Problem

A parallel cut-cell algorithm is described to solve the free-boundary problem of the Grad--Shafranov equation. The algorithm reformulates the free-boundary problem in an irregular bounded domain and its important aspects include a searching algorithm for the magnetic axis and separatrix, a surface integral along the irregular boundary to determine the boundary values, an approach to optimize the coil current based on a targeting plasma shape, Picard iterations with Aitken's acceleration for the resulting nonlinear problem, and a Cartesian grid embedded boundary method to handle the complex geometry. Here the algorithm is implemented in parallel using a standard domain-decomposition approach and a good parallel scaling is observed. Numerical results verify the accuracy and efficiency of the free-boundary Grad--Shafranov solver.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

Optimal Experimental Design for Inverse Problems in the Presence of Observation Correlations

Optimal experimental design (OED) is the general formalism of sensor placement and decisions on the data collection strategy for engineered or natural experiments. This problem is prevalent in many critical fields such as battery design, numerical weather prediction, geosciences, environmental and urban studies. State-of-the-art computational methods for experimental design do not accommodate correlation structure in observational errors produced by many expensive-to-operate devices such as X-ray machines, radars, and satellites. Discarding evident data correlations leads to biased results, higher expenses, and waste of valuable resources. We present a general formulation of the OED formalism for model-constrained large-scale Bayesian linear inverse problems, where measurement errors are generally correlated. The proposed approach utilizes the Hadamard product of matrices to formulate the weighted-likelihood, and is valid for both finite as well as infinite-dimensional Bayesian inverse problems. Extensive numerical experiments are carried out for empirical verification of the proposed approach using an advection-diffusion model, where the objective is to optimally place a small set of sensors, under a limited budget, to predict the concentration of a contaminant in a closed and bounded domain.

Correlated observations↗
DOE OSTI2025

Efficient Training of Deep Neural Operator Networks via Randomized Sampling

Neural operators (NOs) employ deep neural networks to learn the mappings between infinitedimensional function spaces. Deep operator network (DeepONet), a popular NO architecture, has demonstrated success in the real-time prediction of complex dynamics across various scientific and engineering applications. In this work, we introduce a random sampling technique to be adopted during the training of DeepONet, aimed at improving the generalization ability of the model, while significantly reducing the computational time. The proposed approach targets the trunk network of the DeepONet model that outputs the basis functions corresponding to the spatiotemporal locations of the bounded domain on which the physical system is defined. While constructing the loss function, DeepONet training traditionally considers a uniform grid of spatiotemporal points at which all the output functions are evaluated for each iteration. This approach leads to a larger batch size, resulting in poor generalization and increased memory demands, due to the limitations of the stochastic gradient descent (SGD) optimizer. The proposed random sampling over the inputs of the trunk net mitigates these challenges, improving generalization and reducing the memory requirements during training, resulting in significant computational gains. We validate our hypothesis through three benchmark examples, demonstrating substantial reductions in training time while achieving comparable or lower overall test errors relative to the traditional training approach. Our results indicate that incorporating randomization in the trunk network inputs during training enhances the efficiency and robustness of DeepONet, offering a promising avenue for improving the framework’s performance in modeling complex physical systems.

Karumuri, Sharmila [Department of Civil & Systems ↗
DOE OSTI2024

Learning Only on Boundaries: A Physics-Informed Neural Operator for Solving Parametric Partial Differential Equations in Complex Geometries

Recently, deep learning surrogates and neural operators have shown promise in solving partial differential equations (PDEs). However, they often require a large amount of training data and are limited to bounded domains. In this work, we present a novel physics-informed neural operator method to solve parameterized boundary value problems without labeled data. By reformulating the PDEs into boundary integral equations (BIEs), we can train the operator network solely on the boundary of the domain. This approach reduces the number of required sample points from $O(N^d)$ to $O(N^{d-1}$), where $d$ is the domain’s dimension, leading to a significant acceleration of the training process. Additionally, our method can handle unbounded problems, which are unattainable for existing physics-informed neural networks (PINNs) and neural operators. Finally, our numerical experiments show the effectiveness of parameterized complex geometries and unbounded problems.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Canonical Transformation of Potential Model Hamiltonian Mechanics to Geometrical Form I

Using the methods of symplectic geometry, we establish the existence of a canonical transformation from potential model Hamiltonians of standard form in a Euclidean space to an equivalent geometrical form on a manifold, where the corresponding motions are along geodesic curves. The advantage of this representation is that it admits the computation of geodesic deviation as a test for local stability, shown in recent previous studies to be a very effective criterion for the stability of the orbits generated by the potential model Hamiltonian. We describe here an algorithm for finding the generating function for the canonical transformation and describe some of the properties of this mapping under local diffeomorphisms. We give a convergence proof for this algorithm for the one-dimensional case, and provide a precise geometric formulation of geodesic deviation which relates the stability of the motion in the geometric form to that of the Hamiltonian standard form. We apply our methods to a simple one-dimensional harmonic oscillator and conclude with a discussion of the relation of bounded domains in the two representations for which Morse theory would be applicable.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Peptidoglycan Contribution to the B Cell Superantigen Activity of Staphylococcal Protein A

Staphylococcus aureus causes reiterative and chronic persistent infections. This can be explained by the formidable ability of this pathogen to escape immune surveillance mechanisms. Cells of S. aureus display the abundant staphylococcal protein A (SpA). SpA binds to immunoglobulin (Ig) molecules and coats the bacterial surface to prevent phagocytic uptake. SpA also binds and cross-links variable heavy 3 (V H 3) idiotype (IgM) B cell receptors, promoting B cell expansion and the secretion of nonspecific V H 3-IgM via a mechanism requiring CD4 + T cell help. SpA binding to antibodies is mediated by the N-terminal Ig-binding domains (IgBDs). The so-called region X, uncharacterized LysM domain, and C-terminal LPXTG sorting signal for peptidoglycan attachment complete the linear structure of the protein. Here, we report that both the LysM domain and the LPXTG motif sorting signal are required for the B cell superantigen activity of SpA in a mouse model of infection. SpA molecules purified from staphylococcal cultures are sufficient to exert B cell superantigen activity and promote immunoglobulin secretion as long as they carry intact LysM and LPXTG motif domains with bound peptidoglycan fragments. The LysM domain binds the glycan chains of peptidoglycan fragments, whereas the LPXTG motif is covalently linked to wall peptides lacking glycan. These findings emphasize the complexity of SpA interactions with B cell receptors.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Cryo-EM structure of AAV2 Rep68 bound to integration site AAVS1: insights into the mechanism of DNA melting

Abstract The Rep68 protein from Adeno-Associated Virus (AAV) is a multifunctional SF3 helicase that performs most of the DNA transactions necessary for the viral life cycle. During AAV DNA replication, Rep68 assembles at the origin of replication, catalyzing the DNA melting and nicking reactions during the hairpin rolling replication process to complete the second-strand synthesis of the AAV genome. We report the cryo-electron microscopy structures of Rep68 bound to the adeno-associated virus integration site 1 in different nucleotide-bound states. In the nucleotide-free state, Rep68 forms a heptameric complex around DNA, with three origin-binding domains (OBDs) bound to the Rep-binding element sequence, while three remaining OBDs form transient dimers with them. The AAA+ domains form an open ring without interactions between subunits and DNA. We hypothesize that the heptameric structure is crucial for loading Rep68 onto double-stranded DNA. The ATPγS complex shows that only three subunits associate with the nucleotide, leading to a conformational change that promotes the formation of both intersubunit and DNA interactions. Moreover, three phenylalanine residues in the AAA+ domain induce a steric distortion in the DNA. Our study provides insights into how an SF3 helicase assembles on DNA and provides insights into the DNA melting process.

Jaiswal, Rahul↗
DOE OSTI2025

Structures of a synthetic antibody selected against and bound to the C‐terminal domain of Clostridium perfringens enterotoxin

Abstract Clostridium perfringensenterotoxin (CpE) causes cytotoxic gastrointestinal disease in mammalian epithelium by binding membrane protein receptors called claudins. Claudins direct the formation of cell/cell tight junctions through oligomerization and govern the transport of molecules between individual cells. CpE binds claudins through its C‐terminal domain (cCpE) and induces cytotoxicity through its N‐terminal domain. The non‐toxic cCpE is a useful tool to study claudins, tight junctions, and for translational applications, such as increasing the permeability of restrictive tissues like the blood–brain barrier or selective targeting of claudin overexpressing cancers. Conversely, there are no specialized molecular tools to study CpE or cCpE, or to modulate or inhibit their functions. We previously reported the development of synthetic antigen‐binding fragments (sFabs) that bind cCpE, and low‐resolution structures of them bound to claudin/cCpE complexes. Here, we determine high‐resolution structures of sFab COP‐2 bound to cCpE using X‐ray crystallography and cryogenic electron microscopy. The structures and biophysical findings provide the mechanism of COP‐2 binding to cCpE and the molecular determinants driving their interactions. These insights can advance the design of new antibody‐based tools from our COP‐2 scaffold to study or alter cCpE function and give rise to a “Trojan horse” strategy that exploits cCpE's tight junction barrier disrupting function to selectively deliver conjugated therapeutics through normally impermeable tissues.

Biochemistry & Molecular Biology↗
DOE OSTI2020

The effects of stress, temperature and facet structure on growth of { 10 1 ¯ 2 } twins in Mg: a molecular dynamics and phase field study

The present study focuses on the roles of structure, stress and temperature on the mobility of facets bounding twin domains and on the overall impact on the growth kinetics of { 10 1 ¯ 2 } twins in Mg. We first use molecular dynamics simulations (MD) to study the relationship between interface structures, stress, temperature and mobility. As facets can exhibit distinct atomic scale structures, we also quantify the effect of stress relaxation mediated by misfit dislocations on the mobility of the basal-prismatic interface. The information gathered for ten different facets is used to calibrate an anisotropic phase field model predicting the effects of temperature and stress on the kinetics and shape of a growing { 10 1 ¯ 2 } 3D twin, by comparison with MD data. The comparative analysis suggests that, as compared to simulations of individual facets, the effective mobility of facets can be significantly different during growth of a 3D twin. This can be attributed to the complex three-dimensional internal stress state of the twin and to coupling in facet motion triggered by facet junctions. As such, by directly fitting both the twin growth kinetics and morphology as predicted by MD for different stresses and temperatures, we can identify a limited set of rate-limiting facets that control twin growth and morphology, and measure their effective mobility.

36 MATERIALS SCIENCE↗
DOE OSTI2022

A mass conservative, well balanced, tangency preserving and energy decaying method for the shallow water equations on a sphere

Here a fully discrete surface finite element method is proposed for solving the viscous shallow water equations in a bounded Lipschitz domain on the sphere based on a general triangular mesh. The method consists of a modified Crank–Nicolson method in time and a Galerkin surface finite element method in space for the fluid thickness H and the fluid velocity u. A finite element space tangential to the sphere at all finite element nodes is proposed to approximate the fluid velocity u. The proposed method has second-order accuracy in time and first-order accuracy in space, and preserves mass conservation, well balancedness, tangency of velocity to the sphere, and energy decay. Numerical experiments are presented to illustrate the accuracy of the proposed method and the preservation of the physical properties, including mass conservation, well balancedness, and energy decay. A numerical simulation of ocean mesoscale activity on a circular basin with a continental shelf is provided.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Neutron skins: A perspective from dispersive optical models

An overview of neutron skin predictions obtained using an empirical nonlocal dispersive optical model (DOM) is presented. The DOM links both scattering and bound-state experimental data through a subtracted dispersion relation which allows for fully consistent, data-informed predictions for nuclei where such data exist. Large skins were predicted for both 48 Ca ( R$^{48}_{skin}$ = 0.25 ± 0.023 fm in 2017) and 208 Pb (R$^{208}_{skin}$) = 0.25 ± 0.05 fm in 2020). Whereas the DOM prediction in 208 Pb is within 1σ of the subsequent PREX-2 measurement, the DOM prediction in 48 Ca is over 2σ larger than the thin neutron skin resulting from CREX. From the moment it was revealed, the thin skin in 48 Ca has puzzled the nuclear-physics community as no adequate theories simultaneously predict both a large skin in 208 Pb and a small skin in 48 Ca. The DOM is unique in its ability to treat both structure and reaction data on the same footing, providing a unique perspective on this R skin puzzle. It appears vital that more neutron data be measured in both the scattering and bound-state domain for 48 Ca to clarify the situation.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

Functional characterization of a cellulose synthase, CtCESA1, from the marine red alga Calliarthron tuberculosum (Corallinales)

In land plants and algae, cellulose is important for strengthening cell walls and preventing breakage due to physical forces. Though our understanding of cellulose production by cellulose synthases (CESAs) has seen significant advances for several land plant and bacterial species, functional characterization of this fundamental protein is absent in red algae. Here we identify CESA gene candidates in the calcifying red alga Calliarthron tuberculosum using sequence similarity-based approaches, and elucidate their phylogenetic relationship with other CESAs from diverse taxa. One gene candidate, CtCESA1, was closely related to other putative red algal CESA genes. To test if CtCESA1 encoded a true cellulose synthase, CtCESA1 protein was expressed and purified from insect and yeast expression systems. CtCESA1 showed glucan synthase activity in glucose tracer assays. CtCESA1 activity was relatively low when compared with plant and bacterial CESA activity. In an in vitro assay, a predicted N-terminal starch-binding domain from CtCESA1 bound red algal floridean starch extracts, representing a unique domain in red algal CESAs not present in CESAs from other lineages. When the CtCESA1 gene was introduced into Arabidopsis thaliana cesa mutants, the red algal CtCESA1 partially rescued the growth defects of the primary cell wall cesa6 mutant, but not cesa3 or secondary cell wall cesa7 mutants. A fluorescently tagged CtCESA1 localized to the plasma membrane in the Arabidopsis cesa6 mutant background. This study presents functional evidence validating the sequence annotation of red algal CESAs. Furthermore, the relatively low activity of CtCESA1, partial complementation in Arabidopsis, and presence of unique protein domains suggest that there are probably functional differences between the algal and land plant CESAs.

59 BASIC BIOLOGICAL SCIENCES↗