Search NASA⌕ Search

SEARCH · Search NASA

Results for “unbounded domain”

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 19 records
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Peridynamic elastic waves in two-dimensional unbounded domains: Construction of nonlocal Dirichlet-type absorbing boundary conditions

The focus of this paper is on application of peridynamics (PD) to propagation of elastic waves in unbounded domains. We construct absorbing boundary conditions (ABCs) derived from a semi-analytical solution of the PD governing equation at the exterior region. This solution is made up of a finite series of plane waves, as fundamental solutions (modes), which satisfy the PD dispersion relations. The modes are adjusted to transmit the energy from the interior region (near field) to the exterior region (far field). The corresponding unknown coefficients of the series are found in terms of the displacement field at a layer of points adjacent to the absorbing boundary. This is accomplished through a collocation procedure at subregions (clouds) around each absorbing point. The proposed ABCs offer appealing advantages, which facilitate their application to PD. They are of Dirichlet-type, hence their implementation is relatively simple as no derivatives of the field variables are required. They are constructed in the time and space domains and thus application of Fourier and Laplace transforms, cumbersome for nonlocal models, is not required. At the discrete level, the modes satisfy the same numerical dispersion relations of the near field, which makes the far-field solution compatible with that of the near field. We scrutinize the performance of the proposed ABCs through several examples. So our investigation shows that the proposed ABCs perform stably in time with an appropriate level of accuracy even in problems characterized by highly-dispersive propagating waves, including crack propagation in semi-unbounded brittle solids.

42 ENGINEERING↗
DOE OSTI2020

Dirichlet absorbing boundary conditions for classical and peridynamic diffusion-type models

Diffusion-type problems in (nearly) unbounded domains play important roles in various fields of fluid dynamics, biology, and materials science. The aim of this paper is to construct accurate absorbing boundary conditions (ABCs) suitable for classical (local) as well as nonlocal peridynamic (PD) diffusion models. The main focus of the present study is on the PD diffusion formulation. The majority of the PD diffusion models proposed so far are applied to bounded domains only. In this study, we propose an effective way to handle unbounded domains both with PD and classical diffusion models. For the former, we employ a meshfree discretization, whereas for the latter the finite element method (FEM) is employed. The proposed ABCs are time-dependent and Dirichlet-type, making the approach easy to implement in the available models. The performance of the approach, in terms of accuracy and stability, is illustrated by numerical examples in 1D, 2D, and 3D.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Dirichlet-type absorbing boundary conditions for peridynamic scalar waves in two-dimensional viscous media

Construction of absorbing boundary conditions (ABCs) for nonlocal models is generally challenging, primarily due to the fact that nonlocal operators are commonly associated with volume constrained boundary conditions. Moreover, application of Fourier and Laplace transforms, which are essential for the majority of available methods for ABCs, to nonlocal models is complicated. In this paper, we propose a simple method to construct accurate ABCs for peridynamic scalar wave-type problems in viscous media. The proposed ABCs are constructed in the time and space domains and are of Dirichlet type. Consequently, their implementation is relatively simple, since no derivatives of the wave field are required. The proposed ABCs are derived at the continuum level, from a semi-analytical solution of the exterior domain using harmonic exponential basis functions in space and time (plane-wave modes). The numerical implementation is done using a meshfree collocation approach employed within a boundary layer adjacent to the interior domain boundary. The modes satisfy the peridynamic numerical dispersion relation, resulting in a compatible solution of the interior region (near-field) with that of the exterior region (far-field). The accuracy and stability of the proposed ABCs are demonstrated with several numerical examples in two-dimensional unbounded domains.

42 ENGINEERING↗
DOE OSTI2021

Conservative numerical schemes with optimal dispersive wave relations: Part I. Derivation and analysis

An energy-conserving and an energy-and-enstrophy conserving numerical schemes are derived by approximating the Hamiltonian formulation of the inviscid shallow water flows based on the vorticity-divergence variables. These schemes also conserve the first-order moments such as mass and vorticity, as usual. The conservative properties of the schemes stem from the skew-symmetry and singularities of the Poisson brackets, which are carefully retained in the discrete approximations. Here, the schemes operate on unstructured orthogonal dual meshes, over bounded or unbounded domains, and they are also shown to possess the same optimal dispersive wave relations as those of the Z-grid scheme, which is a consequence of the use of the vorticity and divergence variables.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2022

A Decision-Making Machine Learning Approach in Hermite Spectral Approximations of Partial Differential Equations

The accuracy and effectiveness of Hermite spectral methods for the numerical discretization of partial differential equations on unbounded domains are strongly affected by the amplitude of the Gaussian weight function employed to describe the approximation space. This is particularly true if the problem is under-resolved, i.e., there are no enough degrees of freedom. The issue becomes even more crucial when the equation under study is time-dependent, forcing in this way the choice of Hermite functions where the corresponding weight depends on time. In order to adapt dynamically the approximation space, it is here proposed an automatic decision-making process that relies on machine learning techniques, such as deep neural networks and support vector machines. The algorithm is numerically tested with success on a simple 1D problem, but the main goal is its exportability in the context of more serious applications. Here we also show at the end an application in the framework of plasma physics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

Lattice Green’s Functions for High-Order Finite Difference Stencils

Lattice Green's Functions (LGFs) are fundamental solutions to discretized linear operators, and as such they are a useful tool for solving discretized elliptic PDEs on domains that are unbounded in one or more directions. The majority of existing numerical solvers that make use of LGFs rely on a second-order discretization and operate on domains with free-space boundary conditions in all directions. Under these conditions, fast expansion methods are available that enable precomputation of 2D or 3D LGFs in linear time, avoiding the need for brute-force multi-dimensional quadrature of numerically unstable integrals. Here we focus on higher-order discretizations of the Laplace operator on domains with more general boundary conditions, by (1) providing an algorithm for fast and accurate evaluation of the LGFs associated with high-order dimension-split centered finite differences on unbounded domains, and (2) deriving closed-form expressions for the LGFs associated with both dimension-split and Mehrstellen discretizations on domains with one unbounded dimension. Through numerical experiments we demonstrate that these techniques provide LGF evaluations with near machine-precision accuracy, and that the resulting LGFs allow for numerically consistent solutions to high-order discretizations of the Poisson's equation on fully or partially unbounded 3D domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Structure-function studies of Vibrio cholerae quorum-sensing receptor CqsR signal recognition

Ethanolamine signaling through the transmembrane quorum-sensing receptor CqsR influencesVibrio choleraeniche recognition and host colonization. In this study, we present a comprehensive structure-function analysis of CqsR. Specifically, we have determined X-ray crystal structures of the CqsR periplasmic domain bound to the signaling agonist ethanolamine and its analogs, serinol and L-alaninol, as well as the ligand-free (apo) form of CqsR. The periplasmic ligand-binding domain of CqsR is a Cache domain, the most prevalent extracellular sensory module in prokaryotes. Our findings provide a rare structural comparison of ligand-bound and unbound states of a Cache domain receptor. Coupled with thermodynamic binding assays and genetic analyses, these structures elucidate the molecular basis of CqsR ligand specificity. This study not only advances the understanding of Cache domain function but also informs the identification of ligands for orphan Cache receptors and the rational design of signaling agonists and antagonists. Lastly, we discuss ligand-induced conformational changes in the CqsR Cache domains and explore the potential for the existence of additional regulatory ligands.

Microbiology↗
DOE OSTI2021

Theoretical analysis of key factors achieving reversed magnetic shear q-profiles sustained with lower hybrid waves on EAST

The phase space analysis technique is applied to provide new insights into the fully non-inductive lower hybrid current driven (LHCD) discharges on EAST (Garofalo et al 2017 Nucl. Fusion 57 076037). The analysis shows that there are bounded and unbounded topologies of lower hybrid (LH) waves in phase space. For typical parameters on EAST, the propagation domain for 4.6 GHz LH waves is bounded, while unbounded for 2.45 GHz LH waves and one of the conditions is recognized to be good for achieving an off-axis current profile driven by 4.6 GHz LH waves on EAST. The parametric analysis on the potential power deposition (PPD) region for those experiments demonstrates that the reversed magnetic shear dominates the wave behavior, and confines the LH power absorption to the far off-axis region as long as it occurs in the 4.6 GHz dominated LHCD discharge. GENRAY/CQL3D simulations also confirm this effect of reversed shear on the power absorption of 4.6 GHz waves. A bounded propagation domain and a positive feedback loop between magnetic shear reversal and off-axis LHCD profile could explain the LH-only-sustained strong reversed magnetic shear observed in experiments on multiple machines. In contrast, there is no such effect for 2.45 GHz waves on EAST since the PPD region is less sensitive to the change of reversed shear due to the unbounded propagation domain.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Analysis of nonlinear features associated with radio-frequency sheaths

An analytic study is conducted of nonlinear features that arise when radio-frequency (RF) sheath boundary conditions are solved self-consistently with the electromagnetic fields in a magnetized plasma. An instability associated with RF sheath–plasma interactions is identified and shown to be closely linked with sheath–plasma resonance and the occurrence of multiple solutions. It is proposed that this instability serves as a physical mechanism that enables jumping from one solution branch to another, providing a physical basis for the prediction of an RF sheath hysteresis in the time domain. It is further shown that this instability can lead to non-intuitive features such as spontaneous symmetry breaking in an otherwise completely symmetric problem. In plane-stratified geometry, two analytically solvable examples are considered for the cases of an unbounded half-space and a bounded domain

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

Inverse point source location with the Helmholtz equation on a bounded domain

The problem of recovering acoustic sources, more specifically monopoles, from point-wise measurements of the corresponding acoustic pressure at a limited number of frequencies is addressed. To this purpose, a family of sparse optimization problems in measure space in combination with the Helmholtz equation on a bounded domain is considered in this study. A weighted norm with unbounded weight near the observation points is incorporated into the formulation. Optimality conditions and conditions for recovery in the small noise case are discussed, which motivates concrete choices of the weight. The numerical realization is based on an accelerated conditional gradient method in measure space and a finite element discretization.

97 MATHEMATICS AND COMPUTING↗
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 OSTI2026

On the Trotter Error in Many-body Quantum Dynamics with Coulomb Potentials

Efficient simulation of many-body quantum systems is central to advances in physics, chemistry, and quantum computing, with a key question being whether the simulation cost scales polynomially with the system size. Here, in this work, we analyze many-body quantum systems with Coulomb interactions, which are fundamental to electronic and molecular systems. We prove that Trotterization for such unbounded Hamiltonians achieves a 1/4-order convergence rate, with explicit polynomial dependence on the number of particles. The result holds for all initial wavefunctions in the domain of the Hamiltonian, and the 1/4-order convergence rate is optimal, as previous work has numerically demonstrated that it can be saturated by a specific initial ground state. The main challenges arise from the many-body structure and the singular nature of the Coulomb potential. Our proof strategy differs from prior state-of-the-art Trotter analyses, addressing both difficulties in a unified framework. Our analysis treats the Coulomb potential as an unbounded operator without modification or regularization, and does not rely on spatial discretization, making it compatible with both first- and second-quantized circuit constructions.

Fang, Di [Duke Univ., Durham, NC (United States)]↗
DOE OSTI2024

Probing Supercritical Phase Transition using Ultrafast X-ray Diagnostics (Final Technical Report)

Over recent years, several hypotheses have been put forward to elucidate supercritical transition states. Despite significant progress on the fundamental understanding of fluids at these supercritical conditions, important questions concerning the microstructure and dynamical processes remain. A particular research issue is hereby the fundamental understanding of the morphology of the molecular microstructure and its effect on the macroscopic behaviour and thermodynamic response functions. It is widely believed that the supercritical state is homogeneous without structural and dynamic observables to distinguish between a liquid and a vapor. However, recent investigations have identified regions of distinct liquid-like or vapor-like properties even under supercritical conditions. Specifically, it was shown that the transition between liquid-like and vapor-like states occurs across an extension to the coexistence line, marked by almost discontinuous changes in fluid properties. This transition was first identified experimentally by Nishikawa and Tanaka and stands in contrast to the classical presentation of the supercritical state space as a featureless, homogeneous domain. Structurally, the most important properties of supercritical fluids are the dynamic heterogeneities and local density fluctuations that are present within the fluid state at the microscopic level. These inhomogeneities are associated with the formation of molecular clusters of various sizes with liquid-like properties separated by voids of unbound gas-phase molecules which continuously restructure itself at picosecond timescales. By utilizing ultrafast X-ray Photon Correlation Spectroscopy (XPCS) at the Linac Coherent Light Source (LCLS) and X-ray synchrotron measurements as an experimental probe of these structural changes at the molecular level we seek to elucidate the higher-order phase transition from a liquid-like to a vapor-like state in the region extending the critical point. Open questions we seek to address particularly with XPCS are: (i) the underlying mechanisms responsible for the density fluctuations at supercritical conditions and how these dynamical processes at the molecular level affect the thermodynamic response functions and (ii) on what timescales do dynamical processes evolve that are associated with thermalization, molecular diffusion, and intermolecular cluster transfer.

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

Structures of the mycobacterial membrane protein MmpL3 reveal its mechanism of lipid transport

The mycobacterial membrane protein large 3 (MmpL3) transporter is essential and required for shuttling the lipid trehalose monomycolate (TMM), a precursor of mycolic acid (MA)-containing trehalose dimycolate (TDM) and mycolyl arabinogalactan peptidoglycan (mAGP), in Mycobacterium species, including Mycobacterium tuberculosis and Mycobacterium smegmatis. However, the mechanism that MmpL3 uses to facilitate the transport of fatty acids and lipidic elements to the mycobacterial cell wall remains elusive. Here, we report 7 structures of the M. smegmatis MmpL3 transporter in its unbound state and in complex with trehalose 6-decanoate (T6D) or TMM using single-particle cryo-electron microscopy (cryo-EM) and X-ray crystallography. Combined with calculated results from molecular dynamics (MD) and target MD simulations, we reveal a lipid transport mechanism that involves a coupled movement of the periplasmic domain and transmembrane helices of the MmpL3 transporter that facilitates the shuttling of lipids to the mycobacterial cell wall.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2026

Structural and mechanistic insights into protective non-neutralizing antibodies targeting Crimean-Congo hemorrhagic fever virus nucleocapsid protein

Abstract Crimean-Congo Hemorrhagic Fever Virus (CCHFV) is a tick-borne virus endemic to Africa, Asia, and expanding regions within Europe. With mortality rates approaching 40%, rising incidence, and no currently approved countermeasures, CCHFV is recognized as a priority public health threat. CCHFV nucleocapsid protein (NP) has long been a key target for diagnostics. Recently, NP-specific humoral responses have also been correlated with protection conferred by protective vaccines candidates. Additionally, the first non-neutralizing monoclonal antibody (mAb) 9D5 demonstrated protective efficacy against CCHFV challenge, underscoring NP as a viable antiviral target. Here, nine anti-NP mAb were utilized to identify four antigenic sites on NP as well as localize these sites to the head or stalk domains. These mAb also revealed variable levels of in vivo protection, independent from whether the epitope site is located in the head or stalk regions. Additionally, three X-ray crystallography structures were obtained that included CCHFV NP from strain Afg09-2990 in complex with the most potent mAb (9D5). This, along with additional structures of two unbound NPs, revealed structural elements critical for mAb-9D5 broad-spectrum protective characteristics. These findings provide a path towards the rapid identification of broadly protective anti-NP mAb countermeasures.

Moresco, Vanessa↗
DOE OSTI2023

Less Is More: Oligomer Extraction and Hydrothermal Annealing Increase PDMS Adhesion Forces for Materials Studies and for Biology-Focused Microfluidic Applications

Cues in the micro-environment are key determinants in the emergence of complex cellular morphologies and functions. Primary among these is the presence of neighboring cells that form networks. For high-resolution analysis, it is crucial to develop micro-environments that permit exquisite control of network formation. This is especially true in cell science, tissue engineering, and clinical biology. We introduce a new approach for assembling polydimethylsiloxane (PDMS)-based microfluidic environments that enhances cell network formation and analyses. We report that the combined processes of PDMS solvent-extraction and hydrothermal annealing create unique conditions that produce high-strength bonds between solvent-extracted PDMS (E-PDMS) and glass—properties not associated with conventional PDMS. Extraction followed by hydrothermal annealing removes unbound oligomers, promotes polymer cross-linking, facilitates covalent bond formation with glass, and retains the highest biocompatibility. Herein, our extraction protocol accelerates oligomer removal from 5 to 2 days. Resulting microfluidic platforms are uniquely suited for cell-network studies owing to high adhesion forces, effectively corralling cellular extensions and eliminating harmful oligomers. We demonstrate the simple, simultaneous actuation of multiple microfluidic domains for invoking ATP- and glutamate-induced Ca 2+ signaling in glial-cell networks. These E-PDMS modifications and flow manipulations further enable microfluidic technologies for cell-signaling and network studies as well as novel applications.

36 MATERIALS SCIENCE↗
DOE OSTI2022

All-Atom Simulations Uncover Structural and Dynamical Properties of STING Proteins in the Membrane System

Recent studies have shown that the stimulator of interferon gene (STING) protein plays a central role in the immune system by facilitating the production of type I interferons in cells. The STING signaling pathway is also a prominent activator of cancer-killing T cells that initiate a powerful adaptive immune response. Since biomolecular signaling pathways are complicated and not easily identified through traditional experiments, molecular dynamics (MD) has often been used to study structural and dynamical responses of biological pathways. In this work, we carried out MD simulations for full-length chicken and human STING (chSTING and hSTING) proteins. Specifically, we investigated ligand-bound closed (holo) and ligand-unbound open (apo) forms of STING in the membrane system by comparing their conformational and dynamical differences. Our research provides clues for understanding the mechanism of the STING signaling pathway by uncovering detailed insights for the examined systems: the residues from each chain in the binding pocket are strongly correlated to one another in the open STING structure compared with those in the closed STING structure. Ligand-bound closed STING displays ~174° rotation of the ligand-binding domain (LBD) relative to the open STING structure. The dynamical analysis of residue Cys148 located in the linker region of hSTING does not support the earlier hypothesis that Cys148 can form disulfide bonds between adjacent STING dimers. We also demonstrate that using the full-length proteins is critical, since the MD simulations of the LBD portion alone cannot properly describe the global conformational properties of STING.

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

A new method for obtaining model-free viscoelastic material properties from atomic force microscopy experiments using discrete integral transform techniques

Viscoelastic characterization of materials at the micro- and the nanoscale is commonly performed with the aid of force–distance relationships acquired using atomic force microscopy (AFM). The general strategy for existing methods is to fit the observed material behavior to specific viscoelastic models, such as generalized viscoelastic models or power-law rheology models, among others. Here we propose a new method to invert and obtain the viscoelastic properties of a material through the use of the Z-transform, without using a model. We present the rheological viscoelastic relations in their classical derivation and their z-domain correspondence. We illustrate the proposed technique on a model experiment involving a traditional ramp-shaped force–distance AFM curve, demonstrating good agreement between the viscoelastic characteristics extracted from the simulated experiment and the theoretical expectations. We also provide a path for calculating standard viscoelastic responses from the extracted material characteristics. The new technique based on the Z-transform is complementary to previous model-based viscoelastic analyses and can be advantageous with respect to Fourier techniques due to its generality. Additionally, it can handle the unbounded inputs traditionally used to acquire force–distance relationships in AFM, such as ramp functions, in which the cantilever position is displaced linearly with time for a finite period of time.

36 MATERIALS SCIENCE↗