SEARCH · Engineering Papers
Results for “Boundary condition”
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.
Experimental and numerical studies of seismic fluid-structure interaction in a base-supported cylindrical vessel
Seismic design and qualification of advanced reactors will rely heavily on the use of verified and validated numerical models capable of capturing the interaction of the vessel, its contained fluid, and the internal equipment: fluid-structure interaction (FSI) analysis. Analytical solutions can be used for preliminary sizing and design of such vessels but their application is limited to simple geometries and boundary conditions, and small amplitude, translational (and rotational) inputs. To validate numerical models for seismic FSI analysis in finite element codes, a comprehensive set of experiments was performed on a liquid-filled cylindrical vessel, using a 6 degree-of-freedom earthquake simulator. Results in terms of sloshing frequency, damping ratio in sloshing modes, and hydrodynamic responses (wave height, hydrodynamic pressure, base shear, and base moment) for multi-directional earthquake simulator inputs are reported and compared with analytical solutions for liquid-filled vessels. The impact of seismic (base) isolation on hydrodynamic responses was studied using earthquake simulator inputs generated using a virtual isolation system. Data from the experiments are used to validate a numerical model of the fluid-filled vessel using the Arbitrary Lagrangian Eulerian (ALE) solver in the commercial finite element program LS-DYNA. Validation studies are presented for multi-directional seismic inputs, including rocking motions. Lagrangian modeling approaches using an elastic material formulation for the fluid are also investigated and their limitations and possible applications are identified. Here, the results are broadly applicable to the seismic response of base supported, liquid-filled vessels.
Sensitivity analysis of a layered piezoelectric system using ZFEM
The complex variable finite element method (ZFEM) is a numerical technique which aims to find the partial derivatives of the independent variables with respect to variation in dependent parameters declared in the physics. This is done by combining the complex Taylor series expansion within the weak formulation of the governing equation in a coupled system of linear equations forming a complex valued block matrix given by the Cauchy–Riemann matrix representation. In this work, two-dimensional linear first-order elements have been implemented in ZFEM to predict the design derivatives of the mechanical displacement field and the voltage potential field for a layered piezoelectric system in a steady-state study with Dirichlet boundary condition applied at the top and bottom edges of the geometry. This approach allows the standard FEM solution to quantify the sensitivity of the mechanical displacement and voltage potential fields with respect to small variations in the material properties through the information obtained from the computation of the derivatives. The domain is formed by a layered body with PZT-4 and PZT-5 stacked together. For result verification, the numerical solution obtained with ZFEM was compared to results from a commercial FEM package and the solution from the imaginary part was compared to the exact solution of a well-known benchmark problem. In conclusion, comparison of the results showed good agreement for both the real and imaginary parts of the solution and the largest sensitivities were found in PZT-5 specifically in C 13 , C 33 , and ε 33 .
Strain localization criteria for viscoplastic geomaterials
This report presents a viscoplastic localization criterion to detect quasi-instantaneous (i.e., load-induced) and delayed (creep-induced) strain localization in rate-dependent solids. The study is based on the theory of controllability and a viscoplastic description of the mechanical response. Analytical precursors of unstable states are defined through systems of ordinary differential equations (OEDs). The use of the proposed criteria is illustrated at the material point level through a set of strain localization analyses simulating active strain localization of a porous rock. In addition, full-field finite element simulations of compression tests conducted under various pressures are reported to demonstrate the role of local unstable viscoplasticity in the spontaneous propagation of deformation bands under stationary boundary conditions. The study shows that the viscoplastic localization criterion maintains a negative sign as long as the behavior is unstable, that is, the rate of deformation is accelerating. The sign switch coincides with the transition to decelerating deformation. The analyses revealed that pulses of overstress always emerge in correspondence with the growth of unstable behavior, and the peak matches the transition to stable behavior. The local responses recovered from full-field analyses were consistent with those observed in analyses at material point level and the predictions of the presented theory.
A multiscale anisotropic polymer network model coupled with phase field fracture
The study of polymers has continued to gain substantial attention due to their expanding range of applications, spanning essential engineering fields to emerging domains like stretchable electronics, soft robotics, and implantable sensors. These materials exhibit remarkable properties, primarily stemming from their intricate polymer chain network, which, in turn, increases the complexity of precisely modeling their behavior. Especially for modeling elastomers and their fracture behavior, accurately accounting for the deformations of the polymer chains is vital for predicting the rupture in highly stretched chains. Despite the importance, many robust multiscale continuum frameworks for modeling elastomer fracture tend to simplify network deformations by assuming uniform behavior among chains in all directions. Recognizing this limitation, our study proposes a multiscale fracture model that accounts for the anisotropic nature of elastomer network responses. At the microscale, damage in the chains is assumed to be driven by both the chain's entropy and the internal energy due to molecular bond distortions. In order to bridge the stretching in the chains to the macroscale deformation, we employ the maximal advance path constraint network model, inherently accommodating anisotropic network responses. As a result, chains oriented differently can be predicted to exhibit varying stretch and, consequently, different damage levels. To drive macroscale fracture based on damages in these chains, we utilize the micromorphic regularization theory, which involves the introduction of dual local-global damage variables at the macroscale. The macroscale local damage variable is obtained through the homogenization of the chain damage values, resulting in the prediction of an isotropic material response. The macroscale global damage variable is subjected to nonlocal effects and boundary conditions in a thermodynamically consistent phase field continuum formulation. Moreover, the total dissipation in the system is considered to be mainly due to the breaking of the molecular bonds at the microscale. To validate our model, we employ the double-edge notched tensile test as a benchmark, comparing simulation predictions with existing experimental data. Additionally, to enhance our understanding of the fracturing process, we conduct uniaxial tensile experiments on a square film made up of polydimethylsiloxane (PDMS) rubber embedded with a hole and notches and then compare our simulation predictions with the experimental observations. Furthermore, we visualize the evolution of stretch and damage values in chains oriented along different directions to assess the predictive capacity of the model. In conclusion, the results are also compared with another existing model to evaluate the utility of our model in accurately simulating the fracture behavior of rubber-like materials.
Interface flux recovery framework for constructing partitioned heterogeneous time-integration methods
A common approach for the development of partitioned schemes employing different time integrators on different subdomains is to lag the coupling terms in time. This can lead to accuracy issues, especially in multistage methods. Here, in this article, we present a novel framework for partitioned heterogeneous time-integration methods, which allows the coupling of arbitrary multistage and multistep methods without reducing their order of accuracy. At the core of our approach are accurate estimates of the interface flux obtained from the Schur complement of an auxiliary monolithic system. We use these estimates to construct a polynomial-in-time approximation of the interface flux over the current time coupling window. This approximation provides the interface boundary conditions necessary to decouple the subdomain problems at any point within the coupling window. In so doing our framework enables a flexible choice of time-integrators for the individual subproblems without compromising the time-accuracy at the coupled problem level. This feature is the main distinction between our framework and other approaches. To demonstrate the framework, we construct a family of partitioned heterogeneous time-integration methods, combining multistage and multistep methods, for a simplified tracer transport component of the coupled air-sea system in Earth system models. We report numerical tests evaluating accuracy and flux conservation for different pairs of time-integrators from the explicit Runge-Kutta and Adams-Moulton families.
Imaging Bias‐Driven Domain Wall Motion With Scanning Oscillator Piezoresponse Force Microscopy
Understanding ferroelectric domain wall dynamics at the nanoscale across a broad range of timescales requires measuring domain wall position under different applied electric fields. The success of piezoresponse force microscopy (PFM) as a tool to apply local electric fields at different positions and imaging their changing position, together with the information obtained from associated switching spectroscopies has fueled numerous studies of the dynamics of ferroelectric domains to determine the impact of intrinsic parameters such as crystalline order, defects and pinning centers, as well as boundary conditions such as environment. However, the investigation of sub-coercive reversible domain wall vibrational modes requires the development of new tools that enable visualizing domain wall motion under varying applied fields with high temporal and spatial resolution while also accounting for spurious electrostatic effects. Here, scanning oscillator piezoresponse force microscopy extends the investigation of domain wall dynamics to new regimes, providing direct visualization of domain wall position as a function of an external electric field that varies in time and location. In conclusion, this enables studying the energetics of field-driven ferroelectric domain wall motion, which is shown to obey a thermally activated flow regime in the millisecond timescale.
Ensemble‐Based, Large‐Eddy Reconstruction of Wind Turbine Inflow in a Near‐Stationary Atmospheric Boundary Layer Through Generative Artificial Intelligence
ABSTRACT To validate the second‐by‐second dynamics of turbines in field experiments, it is necessary to accurately reconstruct the winds going into the turbine. Current time‐resolved inflow reconstruction techniques estimate wind behavior in unobserved regions using relatively simple spectral‐based models of the atmosphere. Here, we develop a technique for time‐resolved inflow reconstruction that is rooted in a large‐eddy simulation model of the atmosphere. Our “large‐eddy reconstruction” technique blends observations and atmospheric model information through a diffusion model machine learning algorithm, allowing us to generate probabilistic ensembles of reconstructions for a single 10‐min observational period. Our generated inflows can be used directly by aeroelastic codes or as inflow boundary conditions in a large‐eddy simulation. We verify the second‐by‐second reconstruction capability of our technique in three synthetic field campaigns, finding positive Pearson correlation coefficient values () between ground‐truth and reconstructed streamwise velocity, as well as smaller positive correlation coefficient values for unobserved fields (spanwise velocity, vertical velocity, and temperature). We validate our technique in three real‐world case studies by driving large‐eddy simulations with reconstructed inflows and comparing to independent inflow measurements. The reconstructions are visually similar to measurements, follow desired power spectra properties, and track second‐by‐second behavior ().
Entanglement island, miracle operators and the firewall
In this paper, we obtain some general results on information retrieval from the black hole interior, based on the recent progress on quantum extremal surface formula and entanglement island. We study an AdS black hole coupled to a bath with generic dynamics, and ask whether it is possible to retrieve information about a small perturbation in the interior from the bath system. We show that the one-norm distance between two reduced states in a bath region A is equal to the same quantity in the bulk quantum field theory for region AI where I is the entanglement island of A. This is a straightforward generalization of bulk-boundary correspondence in AdS/CFT. However, we show that a contradiction arises if we apply this result to a special situation when the bath dynamics includes a unitary operation that carries a particular measurement to a region A and send the result to another region W. Physically, the contradiction arises between transferability of classical information during the measurement, and non-transferability of quantum information which determines the entanglement island. We propose that the resolution of the contradiction is to realize that the state reconstruction formula does not apply to the special situation involving interior-information-retrieving measurements. This implies that the assumption of smooth replica AdS geometry with boundary condition set by the flat space bath has to break down when the particular measurement operator is applied to the bath. Using replica trick, we introduce an explicitly construction of such operator, which we name as “miracle operators”. From this construction we see that the smooth replica geometry assumption breaks down because we have to introduce extra replica wormholes connecting with the “simulated blackholes” introduced by the miracle operator. We study the implication of miracle operators in understanding the firewall paradox.
Heisenberg spin chain as a worldsheet coordinate for lightcone quantized string
Although the energy spectrum of the Heisenberg spin chain on a circle defined by H = 1 4 Σ k = 1 M σ k x σ k + 1 x + σ k y σ k + 1 y + Δ σ k z σ k + 1 z is well known for any fixed M, the boundary conditions vary according to whether Mϵ 4N + r, where r = -1,0,1,2, and also according to the parity of the number of overturned spins in the state, In string theory all these cases must be allowed because interactions involve a string with M spins breaking into strings with M 1 < M and M - M 1 spins (or vice versa). We organize the energy spectrum and degeneracies of H in the case Δ = 0 where the system is equivalent to a system of free fermions. In spite of the multiplicity of special cases, in the limit M → ∞ the spectrum is that of a free compactified worldsheet field. Such a field can be interpreted as a compact transverse string coordinate x(σ) ≡ x(σ) + R 0 . We construct the bosonization formulas explicitly in all separate cases, and for each sector give the Virasoro conformal generators in both fermionic and bosonic formulations. Furthermore from calculations in the literature for selected classes of excited states, there is strong evidence that the only change for Δ ≠ 0 is a change in the compactification radius R 0 → R Δ . As Δ→ -1 this radius goes to infinity, giving a concrete example of noncompact space emerging from a discrete dynamical system. Finally we apply our work to construct the three string vertex implied by a string whose bosonic coordinates emerge from this mechanism.
Inside the hologram: reconstructing the bulk observer’s experience
We develop a holographic framework for describing the experience of bulk observers in AdS/CFT, that allows us to compute the proper time and energy distribution measured along any bulk worldline. Our method is formulated directly in the CFT language and is universal: it does not require knowledge of the bulk geometry as an input. When used to propagate operators along the worldline of an observer falling into an eternal black hole, our proposal resolves a conceptual puzzle raised by Marolf and Wall. Notably, the prescription does not rely on an external dynamical Hamiltonian or the AdS boundary conditions and is, therefore, outlining a general framework for the emergence of time.
Trigonometric continuous-variable gates and hybrid quantum simulations of the sine-Gordon model
Hybrid qubit-qumode quantum computing platforms provide a natural setting for simulating interacting bosonic quantum field theories. However, existing continuous-variable gate constructions rely predominantly on polynomial functions of canonical quadratures. In this work, we introduce a complementary universality paradigm based on trigonometric continuous-variable gates, which enable a Fourier-like representation of bosonic operators and are particularly well suited for periodic and non-perturbative interactions. We present an ancilla-based framework for implementing trigonometric gates with arguments given by arbitrary Hermitian functions of qumode quadratures. The protocol yields unitary gates deterministically, and non-unitary gates through probabilistic post-selection. As a concrete application, we develop a hybrid qubit-qumode quantum simulation of the lattice sine-Gordon model. Using these gates, we prepare ground states via quantum imaginary-time evolution, simulate real-time dynamics, compute time-dependent vertex two-point correlation functions, and extract quantum kink profiles under topological boundary conditions. Our results demonstrate that trigonometric continuous-variable gates provide a physically natural framework for simulating interacting field theories on near-term hybrid quantum hardware, while establishing a parallel route to universality beyond polynomial gate constructions. We expect that the trigonometric gates introduced here to find broader applications, including quantum simulations of condensed matter systems, quantum chemistry, and biological models.
Gluing. Part I. Integrals and symmetries
We review some aspects of the cutting and gluing law in local quantum field theory (QFT) and study it from a new point of view. In particular, we emphasize the description of gluing by a path integral over a space of polarized boundary conditions, which are given by leaves of some Lagrangian foliation in the phase space. We think of this path integral as a non-local (d – 1)-dimensional gluing theory associated to the parent local d-dimensional QFT. This is a novel point of view paving the way for applications of the standard QFT techniques (that do not rely on locality) to the gluing theory. We describe various properties of this procedure and spell out conditions under which symmetries of the parent theory lead to symmetries of the gluing theory. The purpose of this paper is to set up a playground for the companion paper where these techniques are applied to obtain new results in supersymmetric theories.
Multi-instanton calculus in $c$ = 1 string theory
We formulate a strategy for computing the complete set of non-perturbative corrections to closed string scattering in c = 1 string theory from the worldsheet per spective. This requires taking into account the effect of multiple ZZ-instantons, including higher instantons constructed from ZZ boundary conditions of type (m, 1), with a careful treatment of the measure and contour in the integration over the instanton moduli space. The only a priori ambiguity in our prescription is a normalization constant $\mathcal{N}_m$ that ap pears in the integration measure for the (m, 1)-type ZZ instanton, at each positive integer m. We investigate leading corrections to the closed string reflection amplitude at the n instanton level, i.e. of order e -n/g s , and find striking agreement with our recent proposal on the non-perturbative completion of the dual matrix quantum mechanics, which in turn fixes $\mathcal{N}_m$ for all m.
Geometric entropies and their Hamiltonian flows
In holographic theories, the Hubeny-Rangamani-Takayanagi (HRT) area operator plays a key role in our understanding of the emergence of semiclassical Einstein-Hilbert gravity. When higher derivative corrections are included, the role of the area is instead played by a more general functional known as the geometric entropy. It is thus of interest to understand the flow generated by the geometric entropy on the classical phase space. In particular, the fact that the associated flow in Einstein-Hilbert or Jackiw-Teitelboim (JT) gravity induces a relative boost between the left and right entanglement wedges is deeply related to the fact that gravitational dressing promotes the von Neumann algebra of local fields in each wedge to type II. This relative boost is known as a boundary-condition-preserving (BCP) kink-transformation. In a general theory of gravity (with arbitrary higher-derivative terms), it is straightforward to show that the flow continues to take the above geometric form when acting on a spacetime where the HRT surface is the bifurcation surface of a Killing horizon. However, the form of the flow on other spacetimes is less clear. In this paper, we use the manifestly-covariant Peierls bracket to explore such flows in two-dimensional theories of JT gravity coupled to matter fields with higher derivative interactions. The results no longer take a purely geometric form and, instead, demonstrate new features that should be expected of such flows in general higher derivative theories. We also show how to obtain the above flows using Poisson brackets.
Fishnet four-point integrals: integrable representations and thermodynamic limits
In this work, we consider four-point integrals arising in the planar limit of the conformal “fishnet” theory in four dimensions. They define a two-parameter family of higher-loop Feynman integrals, which extend the series of ladder integrals and were argued, based on integrability and analyticity, to admit matrix-model-like integral and determinantal representations. In this paper, we prove the equivalence of all these representations using exact summation and integration techniques. We then analyze the large-order behaviour, corresponding to the thermodynamic limit of a large fishnet graph. The saddle-point equations are found to match known two-cut singular equations arising in matrix models, enabling us to obtain a concise parametric expression for the free-energy density in terms of complete elliptic integrals. Interestingly, the latter depends non-trivially on the fishnet aspect ratio and differs from a scaling formula due to Zamolodchikov for large periodic fishnets, suggesting a strong sensitivity to the boundary conditions. We also find an intriguing connection between the saddle-point equation and the equation describing the Frolov-Tseytlin spinning string in AdS 3 × S 1 , in a generalized scaling combining the thermodynamic and short-distance limits.
Investigating two-dimensional adjoint QCD on the lattice
We present our investigations of SU(N) adjoint QCD in two dimensions with one Majorana fermion on the lattice. We determine the relevant parameter range for the simulations with Wilson fermions and present results for Polyakov loop, chiral condensate, and string tension. In the theory with massive fermions, all observables we checked show qualitative agreement between numerical lattice data and theory, while the massless limit is more subtle since chiral and non-invertible symmetry of the continuum theory are explicitly broken by lattice regularization. In thermal compactification, we observe N perturbative vacua for the holonomy potential at high-T with instanton events connecting them, and a unique vacuum at low-T. At finite-N, this is a cross-over and it turns to a phase transition at large-N thermodynamic limit. In circle compactification with periodic boundary conditions, we observe a unique center-symmetric minimum at any radius. In continuum, the instantons in the thermal case carry zero modes (for even N) and indeed, in the lattice simulations, we observe that chiral condensate is dominated by instanton centers, where zero modes are localized. We present lattice results on the issue of confinement vs. screening in the theory and comment on the roles of chiral symmetry and non-invertible symmetry.
On the symmetry TFT of Yang-Mills-Chern-Simons theory
Three-dimensional Yang-Mills-Chern-Simons theory has the peculiar property that its one-form symmetry defects have nontrivial braiding, namely they are charged under the same symmetry they generate, which is then anomalous. This poses a few puzzles in describing the corresponding Symmetry TFT in a four-dimensional bulk. First, the braiding between lines at the boundary seems to be ill-defined when such lines are pulled into the bulk. Second, the Symmetry TFT appears to be too trivial to allow for topological boundary conditions encoding all the different global variants. We show that both of these puzzles can be solved by including endable (tubular) surfaces in the class of bulk topological operators one has to consider. In this way, we are able to reproduce all global variants of the theory, with their symmetries and their anomalies. We check the validity of our proposal also against a top-down holographic realization of the same class of theories.