Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “weak form”

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 OSTI2023

Electron scattering and neutrino physics

Abstract A thorough understanding of neutrino–nucleus scattering physics is crucial for the successful execution of the entire US neutrino physics program. Neutrino–nucleus interaction constitutes one of the biggest systematic uncertainties in neutrino experiments—both at intermediate energies affecting long-baseline deep underground neutrino experiment, as well as at low energies affecting coherent scattering neutrino program—and could well be the difference between achieving or missing discovery level precision. To this end, electron–nucleus scattering experiments provide vital information to test, assess and validate different nuclear models and event generators intended to test, assess and validate different nuclear models and event generators intended to be used in neutrino experiments. Similarly, for the low-energy neutrino program revolving around the coherent elastic neutrino–nucleus scattering (CEvNS) physics at stopped pion sources, such as at ORNL, the main source of uncertainty in the evaluation of the CEvNS cross section is driven by the underlying nuclear structure, embedded in the weak form factor, of the target nucleus. To this end, parity-violating electron scattering (PVES) experiments, utilizing polarized electron beams, provide vital model-independent information in determining weak form factors. This information is vital in achieving a percent level precision needed to disentangle new physics signals from the standard model expected CEvNS rate. In this white paper, we highlight connections between electron- and neutrino–nucleus scattering physics at energies ranging from 10 s of MeV to a few GeV, review the status of ongoing and planned electron scattering experiments, identify gaps, and lay out a path forward that benefits the neutrino community. We also highlight the systemic challenges with respect to the divide between the nuclear and high-energy physics communities and funding that presents additional hurdles in mobilizing these connections to the benefit of neutrino programs.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

A unified framework of stabilized finite element method for solving the Boltzmann transport equation

This paper presents a unified framework of stabilized finite element method for solving the Boltzmann transport equation. Unlike the traditional Petrov-Galerkin finite element method which modifies the test function to construct the stabilization term, we derive the stabilization methods from the standard Galerkin weak form with Sub-grid scale model. The basic idea of this method is to decompose the unknowns into its numerical solution and residual, with an approximation for the residual and embeds it in the Galerkin weak form to yield a stabilized variational formula. Different approximations of the residual lead to different stabilization methods, all the frequently used stabilized methods, including the Streamline Upwinding Petrov-Galerkin (SUPG) method, Galerkin/Least-Square (GLS) method, and Algebraic Sub-Grid Scale (ASGS) method can be obtained from this framework. The similarities and differences of the different approximations are compared in this paper. The numerical results show that the behaviors of the different methods area similar with the same stabilization parameter, and all these stabilized techniques can obtain a correct and stable solution. (authors)

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1997

Numerical Computation of Sensitivities and the Adjoint Approach

We discuss the numerical computation of sensitivities via the adjoint approach in optimization problems governed by differential equations. We focus on the adjoint problem in its weak form. We show how one can avoid some of the problems with the adjoint approach, such as deriving suitable boundary conditions for the adjoint equation. We discuss the convergence of numerical approximations of the costate computed via the weak form of the adjoint problem and show the significance for the discrete adjoint problem.

Lewis, Robert Michael↗
DOE OSTI2021

A finite element formulation for deformation twinning induced strain localization in polycrystal magnesium alloys

Deformation twinning induces shear strain localization in hexagonal close-packed crystals and is critical for the material’s ductility and failure. Cracks often occur at twin-twin or twin-grain boundary intersections and propagate along twin bands. However, most crystal plasticity models for deformation twinning are based on a “pseudo-slip” approach and do not capture the localized deformation associated with the formation of each discrete twin band. The few exceptions are discrete twin models that involve very complex numerical algorithms and are often compromised in accuracy due to the numerical convergence. These factors make the discrete twin models hard to adopt. This paper proposes a modification to the conventional finite element weak form, to fully incorporate a twin-induced heterogeneous deformation that does not depend on the “pseudo-slip” assumption. The model starts by splitting the deformation gradient into elastic-slip-twinning components. The twin-induced deformation gradient component is computed separately by solving a microstructural evolution problem and then implemented into finite element weak form by constructing a global “twin-force” vector. The constitutive update (e.g., in the user-defined material subroutine, or UMAT, for ABAQUS) therefore avoids dealing with the twinning and recovers to the form of a regular slip-based crystal plasticity model. The results presented here indicate that the twin-induced strain localization and the associated stress-reversal phenomena near the twin band were naturally captured in the model, which was validated against an in-situ synchrotron X-ray micro-diffraction experiment.

36 MATERIALS SCIENCE↗
NASA NTRS2002

Development of MLPG and LBIE Methods for Nonlinear Problems of Fracture

The truly meshless local Petrov-Galerkin (MLPG) method holds a great promise in solving boundary value problems, using a local symmetric weak form as a natural approach. In the present paper, in the context of MLPG and the meshless interpolation of a moving least squares (MLS) type, a method which uses primary and secondary nodes in the domain and on the global boundary is introduced, in order to improve the accuracy of solution. The secondary nodes can be placed at any location where one needs to obtain a better resolution. The sub-domains for the shape functions in the MLS approximation are defined only from the primary nodes, and the secondary nodes use the same sub-domains. The shape functions based on the MLS approximation, in an integration domain, have a single type of a rational function, which reduces the difficulty of numerical integration to evaluate the weak form. The present method is very useful in an adaptive calculation, because the secondary nodes can be easily added and/or moved without an additional mesh. The essential boundary conditions can be imposed exactly, and non-convex boundaries can be treated without special techniques. Several numerical examples are presented to illustrate the performance of the present method.

Atluri, Satya N.↗
DOE OSTI2024

Coarse-graining Hamiltonian systems using WSINDy

Abstract Weak form equation learning and surrogate modeling has proven to be computationally efficient and robust to measurement noise in a wide range of applications including ODE, PDE, and SDE discovery, as well as in coarse-graining applications, such as homogenization and mean-field descriptions of interacting particle systems. In this work we extend this coarse-graining capability to the setting of Hamiltonian dynamics which possess approximate symmetries associated with timescale separation. A smooth $$\varepsilon$$ ε -dependent Hamiltonian vector field $$X_\varepsilon$$ X ε possesses an approximate symmetry if the limiting vector field $$X_0=\lim _{\varepsilon \rightarrow 0}X_\varepsilon$$ X 0 = lim ε → 0 X ε possesses an exact symmetry. Such approximate symmetries often lead to the existence of a Hamiltonian system of reduced dimension that may be used to efficiently capture the dynamics of the symmetry-invariant dependent variables. Deriving such reduced systems, or approximating them numerically, is an ongoing challenge. We demonstrate that WSINDy can successfully identify this reduced Hamiltonian system in the presence of large perturbations imparted in the $$\varepsilon >0$$ ε > 0 regime, while remaining robust to extrinsic noise. This is significant in part due to the nontrivial means by which such systems are derived analytically. WSINDy naturally preserves the Hamiltonian structure by restricting to a trial basis of Hamiltonian vector fields. The methodology is computationally efficient, often requiring only a single trajectory to learn the global reduced Hamiltonian, and avoiding forward solves in the learning process. In this way, we argue that weak-form equation learning is particularly well-suited for Hamiltonian coarse-graining. Using nearly-periodic Hamiltonian systems as a prototypical class of systems with approximate symmetries, we show that WSINDy robustly identifies the correct leading-order system, with dimension reduced by at least two, upon observation of the relevant degrees of freedom. While our main contribution is computational, we also provide a contribution to the literature on averaging theory by proving that first-order averaging at the level of vector fields preserves Hamiltonian structure in nearly-periodic Hamiltonian systems. This provides theoretical justification for our approach as WSINDy’s computations occur at the level of Hamiltonian vector fields. We illustrate the efficacy of our proposed method using physically relevant examples, including coupled oscillator dynamics, the Hénon–Heiles system for stellar motion within a galaxy, and the dynamics of charged particles.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2026

Weak-charge form-factor determination at the electron-ion collider

Determining the weak charge form factor, 𝐹 𝑊 ⁡(𝑄 2 ), of nuclei over a continuous range of momentum transfers, 0 ≲ 𝑄 2 ≲ 0.1 GeV 2 , is essential for mapping out the distribution of neutrons in nuclei. The neutron density distribution has significant implications for a broad range of areas, including studies of nuclear structure, neutron stars, and physics beyond the Standard Model. Currently, our knowledge of 𝐹 𝑊 ⁡(𝑄 2 ) comes primarily from fixed target experiments that measure the parity-violating asymmetry in coherent elastic electron-ion scattering. Fixed target experiments, such as CREX and PREX-1,2, have provided high-precision weak charge form factor extractions for the 48 Ca and 208 Pb nuclei, respectively. However, a major limitation of fixed target experiments is that they each provide data only at a single value of 𝑄 2 . With the proposed electron-ion collider (EIC) on the horizon, we explore its potential to impact the determination of the weak charge form factor. While it cannot compete with the precision of fixed target experiments, it can provide data over a wide and continuous range of 𝑄 2 values, and for a wide variety of nuclei. We show that with data corresponding to an integrated luminosity of ℒ ∼ 500/𝐴 fb −1 , where 𝐴 is the nucleus atomic weight, the EIC can significantly impact constraints by lifting degeneracies in theoretical models of the neutron density distribution. Ensuring EIC detector coverage at low 𝑄 2 and large negative pseudorapidities will be essential for such 𝐹 𝑊 ⁡(𝑄 2 ) measurements.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2019

Coherent elastic neutrino-nucleus scattering on Ar 40 from first principles

Coherent elastic neutrino scattering on the 40 Ar nucleus is computed with coupled-cluster theory based on nuclear Hamiltonians inspired by effective field theories of quantum chromodynamics. Our approach is validated by calculating the charge form factor and comparing it to data from electron scattering. We make predictions for the weak form factor, the neutron radius, and the neutron skin and estimate systematic uncertainties. The neutron-skin thickness of 40 Ar is consistent with results from density-functional theory. Precision measurements from coherent elastic neutrino-nucleus scattering could potentially be used to extract these observables and help to constrain nuclear models.

39 ≤ A ≤ 58↗
DOE OSTI2021

Generalized Boozer coordinates: A natural coordinate system for quasisymmetry

We prove the existence of a straight-field-line coordinate system we call generalized Boozer coordinates. This coordinate system exists for magnetic fields with nested toroidal flux surfaces provided ∮ d l / B   ( j · ∇ ψ ) = 0 , where symbols have their usual meaning, and the integral is taken along closed magnetic field lines. All quasisymmetric fields, regardless of their associated form of equilibria, must satisfy this condition. This coordinate system presents itself as a convenient form to describe general quasisymmetric configurations and their properties. Insight can be gained analytically into the difference between strong and weak forms of quasisymmetry, as well as axisymmetry, and the interaction of quasisymmetry with different forms of equilibria.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Formation of Water Networks on Anionic Perylene

We present infrared photodissociation spectra of hydrated perylene anion clusters with up to four water molecules, as well as electronic structure calculations based on density functional theory. Water molecules form weak hydrogen bonds to the π system of the perylene anion. For clusters with more than one water molecule, water−water hydrogen bonds are formed, which generally appear to be stronger than water−π hydrogen bonds, especially for the trihydrate and tetrahydrate. The resulting water networks exist as water subclusters on the surface of the carbon frame of perylene. Here, we observe temperature-dependent dynamic effects, which highlight large amplitude motions of the water network and the shallowness of the potential energy surfaces governing the structures of these clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS1993

New results on the realizability of Reynolds stress turbulence closures

The realizability of Reynolds stress models in homogeneous turbulence is critically assessed from a theoretical standpoint. It is proven that a well known second-order closure formulated by Shih and Lumley using the strong realizability constraints of Schumann is, in fact, not a realizable model. The problem arises from the failure to properly satisfy the necessary positive second time derivative constraint when a principal Reynolds stress vanishes - a fatal flaw that becomes apparent when the non-analytic terms in their model are made single-valued as required on physical grounds. It is furthermore shown that the centrifugal acceleration generated by rotations of the principal axes of the Reynolds stress tensor can make the second derivative singular at the most extreme limits of realizable turbulence. This previously overlooked effect appears to make it impossible to identically satisfy the strong form of realizability in any version of the present generation of second-order closures. On the other hand, models properly formulated to satisfy the weak form of realizability - wherein states of one or two component turbulence are not accessible in finite time are found to be realizable. However, unlike the simpler and more commonly used second order closures, these models can be ill-behaved near the extreme limits of realizable turbulence due to the way that higher-degree nonlinearities are often unnecessarily introduced to satisfy realizability. Illustrative computations of homogeneous shear flows are presented to demonstrate these points which can have important implications for turbulence modeling.

Speziale, Charles G.↗
NASA NTRS1989

Tau approximation techniques for identification of coefficients in parabolic PDE

A variant of the Tau method, called the weak Tau method, is developed on the basis of the weak form of the PDE for use in least-squares parameter estimation; also presented is a suitable abstract convergence framework. The emphasis is on the theoretical framework that allows treatment of the weak Tau method when it is applied to a wide class of inverse problems, including those for diffusion-advection equations, the Fokker-Planck model for population dynamics, and damped beam equations. Extensive numerical testing of the weak Tau method has demonstrated that it compares quite favorably with existing methods.

Banks, H. T.↗
DOE OSTI2021

Accurate Determination of the Neutron Skin Thickness of Pb 208 through Parity-Violation in Electron Scattering

We report a precision measurement of the parity-violating asymmetry A PV in the elastic scattering of longitudinally polarized electrons from 208 Pb. We measure A PV = 550 ± 16 (stat) ±8 (syst) parts per billion, leading to an extraction of the neutral weak form factor F W (Q 2 = 0.00616 GeV 2 ) = 0.368 ± 0.013. Combined with our previous measurement, the extracted neutron skin thickness is R n -R p = 0.283 ± 0.071 fm. The result also yields the first significant direct measurement of the interior weak density of 208 Pb: ρ$^0_W$ = -0.0796 ± 0.0036(exp) ± 0.0013(theo) fm -3 leading to the interior baryon density ρ$^0_b$ = 0.1480 ± 0.0036(exp) ± 0.0013(theo) fm -3 . Finally, the measurement accurately constrains the density dependence of the symmetry energy of nuclear matter near saturation density, with implications for the size and composition of neutron stars.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2023

Precise Measurement of the Neutron Skin Thicknesses of 208Pb and 48Ca

Despite great leaps of development in nuclear physics over the past century, we still do not have a comprehensive understanding of nuclear structure. This is mainly due to the lack of precise knowledge of the neutron distribution inside nucleus. Electromagnetic probes are ineffective in probing neutrons, unlike their charged partners, the protons. Therefore, it is of great importance to constrain the neutron distributions experimentally. Heavy nuclei have more neutrons than protons in order to balance the repulsion between protons. In such neutron-rich nuclei, the extra neutrons are pushed out to the surface by the nuclear symmetry energy, forming the so-called ?neutron skin?. The neutron skin can be probed with a well-established experimental technique ? parity-violating electron scattering (PVES). Using the scattering of longitudinally polarized electrons by an unpolarized target, PREX-II and CREX measure the small parity-violating asymmetry in cross sections. Employing electrons with opposite helicities the weak form factor, the neutron distribution and the neutron skin thickness of the target nucleus are extracted. With excellent beam qualities and dedicated instrumentation at Jefferson Lab, the asymmetry measurements are statistics-limited. We report the results of these two high-precision measurements and their implications on broad topics, from the nuclear structures to the neutron stars.

Zhang, Weibin↗
NASA NTRS2020

Model-Reference Adaptive Control of Distributed Lagrangian Infinite-Dimensional Systems Using Hamilton’s Principle

This paper presents a Hamilton's principle for distributed control of infinite-dimensional systems modeled by a distributed form of the Euler-Lagrange method. The distributed systems are governed by a system of linear partial differential equations in space and time. A generalized potential energy expression is developed that can capture most physical systems including those systems that have no spatial distribution. The Hamilton's principle is applied to derive distributed feedback control methods without resorting to the standard weak-form discretization approach to convert an infinite-dimensional systems to a finite-dimensional systems. It can be shown by the principle of least action that the distributed control synthesized by the Hamilton's principle is a minimum-norm control. A model-reference adaptive control framework is developed for distributed Lagrangian systems in the presence of uncertainty. The theory is demonstrated by an application of adaptive flutter suppression control of a flexible aircraft wing.

Nguyen, Nhan T.↗
DOE OSTI2025

A higher-order finite-element implementation of the nonlinear Fokker–Planck collision operator for charged particle collisions in a low density plasma

Collisions between particles in a low density plasma are described by the Fokker–Planck collision operator. In applications, this nonlinear integro-differential operator is often approximated by linearised or ad-hoc model operators due to computational cost and complexity. In this work, we present an implementation of the nonlinear Fokker–Planck collision operator written in terms of Rosenbluth potentials in the Rosenbluth–MacDonald–Judd (RMJ) form. The Rosenbluth potentials may be obtained either by direct integration or by solving partial differential equations (PDEs) similar to Poisson's equation: we optimise for performance and scalability by using sparse matrices to solve the relevant PDEs. We represent the distribution function using a tensor-product continuous-Galerkin finite-element representation and we derive and describe the implementation of the weak form of the collision operator. We present tests demonstrating a successful implementation using an explicit time integrator and we comment on the speed and accuracy of the operator. Finally, we speculate on the potential for applications in the current and next generation of kinetic plasma models.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Numerical integration in the virtual element method with the scaled boundary cubature scheme

Abstract The virtual element method (VEM) is a stabilized Galerkin method on meshes that consist of arbitrary (convex and nonconvex) polygonal and polyhedral elements. A crucial ingredient in the implementation of low‐ and high‐order VEM is the numerical integration of monomials and nonpolynomial functions over such elements. In this article, we apply the recently proposed scaled boundary cubature (SBC) scheme to compute the weak form integrals in various virtual element formulations over polygonal and polyhedral meshes. In doing so, we demonstrate the flexibility of the approach and the accuracy that it delivers on a broad suite of boundary‐value problems in 2D and 3D over polytopes with affine faces as well as on elements with curved boundaries. In addition, the use of the SBC scheme is exemplified in an enriched Poisson formulation of the VEM in which weakly singular functions are required to be integrated. This study establishes the SBC method as a simple, accurate and efficient integration scheme for use in the VEM.

Chin, Eric B.↗