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 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
NASA NTRS1989

A weak Hamiltonian finite element method for optimal control problems

A temporal finite element method based on a mixed form of the Hamiltonian weak principle is developed for dynamics and optimal control problems. The mixed form of Hamilton's weak principle contains both displacements and momenta as primary variables that are expanded in terms of nodal values and simple polynomial shape functions. Unlike other forms of Hamilton's principle, however, time derivatives of the momenta and displacements do not appear therein; instead, only the virtual momenta and virtual displacements are differentiated with respect to time. Based on the duality that is observed to exist between the mixed form of Hamilton's weak principle and variational principles governing classical optimal control problems, a temporal finite element formulation of the latter can be developed in a rather straightforward manner. Several well-known problems in dynamics and optimal control are illustrated. The example dynamics problem involves a time-marching problem. As optimal control examples, elementary trajectory optimization problems are treated.

Hodges, Dewey H.↗
NASA NTRS1991

Weak Hamiltonian finite element method for optimal control problems

A temporal finite element method based on a mixed form of the Hamiltonian weak principle is developed for dynamics and optimal control problems. The mixed form of Hamilton's weak principle contains both displacements and momenta as primary variables that are expanded in terms of nodal values and simple polynomial shape functions. Unlike other forms of Hamilton's principle, however, time derivatives of the momenta and displacements do not appear therein; instead, only the virtual momenta and virtual displacements are differentiated with respect to time. Based on the duality that is observed to exist between the mixed form of Hamilton's weak principle and variational principles governing classical optimal control problems, a temporal finite element formulation of the latter can be developed in a rather straightforward manner. Several well-known problems in dynamics and optimal control are illustrated. The example dynamics problem involves a time-marching problem. As optimal control examples, elementary trajectory optimization problems are treated.

Hodges, Dewey H.↗
NASA NTRS1990

A weak Hamiltonian finite element method for optimal control problems

A temporal finite element method based on a mixed form of the Hamiltonian weak principle is developed for dynamics and optimal control problems. The mixed form of Hamilton's weak principle contains both displacements and momenta as primary variables that are expanded in terms of nodal values and simple polynomial shape functions. Unlike other forms of Hamilton's principle, however, time derivatives of the momenta and displacements do not appear therein; instead, only the virtual momenta and virtual displacements are differentiated with respect to time. Based on the duality that is observed to exist between the mixed form of Hamilton's weak principle and variational principles governing classical optimal control problems, a temporal finite element formulation of the latter can be developed in a rather straightforward manner. Several well-known problems in dynamics and optimal control are illustrated. The example dynamics problem involves a time-marching problem. As optimal control examples, elementary trajectory optimization problems are treated.

Hodges, Dewey H.↗
DOE OSTI2021

Grid Forming Inverter: Laboratory-Scale Hardware Test Bed Setup and Weak Grid Operation

In recent times, the concept of grid forming inverters has gained popularity. Grid forming inverters have proven to be a promising alternative to the grid following inverters. In this paper, we present the hardware test bed implementation of grid forming inverter in islanding mode as well as in grid connected mode. The control structure includes inverter-level inner current control and outer voltage control, and plant-level P - f and Q - V droop control. The Q - V droop provides the voltage command, whereas P - f droop provides the frequency command and is also key for the synchronization process. While it is straightforward to set up the test bed for the islanded operation, it is a challenging process to set up the grid-forming inverter in the grid-connected mode. We present the key technologies and demonstrate the start-up process and weak grid operation using the hardware experiment results. Results from EMT model developed in MATLAB/SimPowerSystems are also presented for side-by-side comparison.

Mittal, Ratik↗
DOE OSTI2020

An inelastic homogenization framework for layered materials with planes of weakness

Many geologic materials have a composite structure, in which macroscopic mechanical behavior is determined by the properties, shape, and heterogeneous distribution of individual constituents. In particular, sedimentary rocks commonly exhibit a layered microstructure, with distinct bedding planes that can also form planes of weakness. In this work, we present a homogenization framework for modeling inelastic layered media. The proposed constitutive model allows for distinct micro-constitutive laws for each layer, explicit representation of layer distributions, as well as incorporation of imperfect bonding at the interface between adjacent layers. No a priori assumptions are needed regarding the specific constitutive models used for the layers and interfaces, providing significant modeling flexibility. The overall framework provides a simple and physically-motivated way of defining anisotropic material behavior as an emergent property of the layered microstructure. The model is calibrated using triaxial and true-triaxial experimental data to demonstrate its ability to describe anisotropic deformation and multiple modes of failure.

58 GEOSCIENCES↗
NASA NTRS1987

A Taylor weak-statement algorithm for hyperbolic conservation laws

Finite element analysis, applied to computational fluid dynamics (CFD) problem classes, presents a formal procedure for establishing the ingredients of a discrete approximation numerical solution algorithm. A classical Galerkin weak-statement formulation, formed on a Taylor series extension of the conservation law system, is developed herein that embeds a set of parameters eligible for constraint according to specification of suitable norms. The derived family of Taylor weak statements is shown to contain, as special cases, over one dozen independently derived CFD algorithms published over the past several decades for the high speed flow problem class. A theoretical analysis is completed that facilitates direct qualitative comparisons. Numerical results for definitive linear and nonlinear test problems permit direct quantitative performance comparisons.

Baker, A. J.↗
DOE OSTI2022

Closures for multicomponent reacting flows based on dispersion analysis

This work presents algebraic closure models associated with advective transport and nonlinear reactions in a Reynolds-averaged Navier-Stokes context for a system of species subject to binary reactions and transport by advection and diffusion. Expanding upon analysis originally developed for non-reactive transport in the context of Taylor dispersion of scalars, this work extends the modified gradient diffusion model explicated by Peters [N. Peters, Turbulent Combustion, Cambridge Monographs on Mechanics (Cambridge University Press, Cambridge, 2000)] and based on work by Corrsin [S. Corrsin, The reactant concentration spectrum in turbulent mixing with a first-order reaction, J. Fluid Mech. 11, 407 (1961)] beyond single-component transport phenomena and involving nonlinear reactions. The presented model forms, from this weakly nonlinear extension of the original dispersion theory, lead to an analytic expression for the eddy diffusivity matrix that explicitly captures the influence of the reaction kinetics on the closure operators. Furthermore, we demonstrate that the derived model form directly translates between flow topologies through a priori and a posteriori testing of a binary species system subject to homogeneous isotropic turbulence. Using two- and three-dimensional direct numerical simulations involving laminar and turbulent flows, it is shown that this framework improves prediction of mean quantities compared to previous results. Lastly, the presented model form, collapses to the earlier gradient diffusion and its modified version derived by Corrsin in the limits of nonreactive species and linear reactions, respectively.

42 ENGINEERING↗
DOE OSTI2026

Beyond Mass and Multiscale Environments: What Shapes Low Surface Brightness Galaxies? Evidence from MaNGA

The origin of low surface brightness (LSB) galaxies remains a key open question in galaxy formation, reflecting the balance between internal mechanisms and environmental influence. Using MaNGA integral-field spectroscopy, we investigate whether LSB and high surface brightness (HSB) galaxies of comparable stellar mass (9 < logM * < 10) occupy distinct environments or differ primarily through internal evolution. Our late-type sample comprises 113 central and 29 satellite LSB galaxies, and 374 central and 142 satellite HSB galaxies. We characterize environments on scales from ∼100 kpc to 10 Mpc, analyzing radial profiles of stellar mass surface density (Σ * ), star formation activity, and gas-phase metallicity. Central LSB and HSB galaxies inhabit similarly low-density large-scale (>200 kpc) environments, but LSB galaxies are more isolated on small scales (∼100 kpc). Even after matching in stellar mass and environment, LSB galaxies show systematically lower Σ * , Σ SFR , and metallicities, often hosting diffuse, weakly star-forming bulges embedded in extended disks. These results indicate that LSB structure and star formation are not primarily governed by a large-scale environment or halo mass. While secondary halo properties such as spin, concentration, or gas accretion history are often invoked, their environmental dependence appears weak. Instead, LSB–HSB differences for centrals likely reflect divergent assembly or interaction histories and internal processes—such as angular momentum-driven disk evolution or inefficient gas conversion—largely decoupled from a large-scale environment. Nonetheless, the environment still influences the observed star formation and chemical differences between central and satellite LSB galaxies.

Shen, Mengting [Xiamen University (China); SDSS Co↗
NASA NTRS1989

A finite element based method for solution of optimal control problems

A temporal finite element based on a mixed form of the Hamiltonian weak principle is presented for optimal control problems. The mixed form of this principle contains both states and costates as primary variables that are expanded in terms of elemental values and simple shape functions. Unlike other variational approaches to optimal control problems, however, time derivatives of the states and costates do not appear in the governing variational equation. Instead, the only quantities whose time derivatives appear therein are virtual states and virtual costates. Also noteworthy among characteristics of the finite element formulation is the fact that in the algebraic equations which contain costates, they appear linearly. Thus, the remaining equations can be solved iteratively without initial guesses for the costates; this reduces the size of the problem by about a factor of two. Numerical results are presented herein for an elementary trajectory optimization problem which show very good agreement with the exact solution along with excellent computational efficiency and self-starting capability. The goal is to evaluate the feasibility of this approach for real-time guidance applications. To this end, a simplified two-stage, four-state model for an advanced launch vehicle application is presented which is suitable for finite element solution.

Bless, Robert R.↗
DOE OSTI2024

Superconductivity in a topological lattice model with strong repulsion

The highly tunable nature of synthetic quantum materials-both in the solid-state and cold atom contexts-invites examining which microscopic ingredients aid in the realization of correlated phases of matter such as superconductors. Recent experimental advances in moiré materials suggest that unifying the features of the Fermi-Hubbard model and quantum Hall systems creates a fertile ground for the emergence of such phases. Here, we introduce the "double Hofstadter"model, a minimal 2D lattice model that incorporates exactly these features: Time-reversal symmetry, band topology, and strong repulsive interactions. By using infinite cylinder density matrix renormalization group methods (cylinder iDMRG), we investigate the ground state phase diagram of this model. We find that it hosts an interaction-induced quantum spin Hall insulator and demonstrate that weakly hole-doping this state gives rise to a superconductor at a finite circumference, with indications that this behavior persists on larger cylinders. Further, at the aforementioned circumference, the superconducting phase is surprisingly robust to perturbations including additional repulsive interactions in the pairing channel. By developing a technique to probe the superconducting gap function in iDMRG, we phenomenologically characterize the superconductor. Namely, we demonstrate that it is formed from the weak pairing of holes atop the quantum spin Hall insulator. Furthermore, we determine the pairing symmetry of the superconductor, finding it to be p-wave-reminiscent of the unconventional superconductivity reported in experiments on twisted bilayer graphene (TBG). Motivated by this, we elucidate structural similarities and differences between our model and those of TBG in its chiral limit. Finally, to provide a more direct experimental realization, we detail an implementation of our Hamiltonian in a system of cold fermionic alkaline-earth atoms in an optical lattice.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
NASA NTRS2002

Fundamental Mixing and Combustion Experiments for Propelled Hypersonic Flight

The first experiment is a study of a coaxial jet discharging into stagnant laboratory air, with center jet of a mixture of 5% oxygen and 95% helium by volume and coflow jet of air. The exit flow pressure of both center-jet and coflow nozzles is 1 atmosphere. The presence of oxygen in the center jet is to allow the use of an oxygen flow-tagging technique (RELIEF4) to obtain non-intrusive velocity measurements. Both jets are nominally Mach 1.8, but, because of the greater speed of sound, the center jet velocity is more than twice that of the coflow. The mixing layer which forms between the center jet and the coflow near the nozzle exit is compressible, with a calculated convective Mach number of approximately 0.7. This geometry has several advantages: The streamwise development of the flow is generally dominated by turbulent stresses (rather than pressure forces), and thus calculations are sensitive to turbulence modeling. It includes features present in supersonic combustors, including a compressible mixing layer near the nozzle exit and a light-gas/air plume downstream. Since it is a free jet, it provides easy access for both optical instrumentation and probes. Since it is axisymmetric, it requires fewer experimental measurements to fully characterize, and calculations can be performed with more modest computer resources. However, weak shock waves formed at the nozzle exit strengthen and turn normal as they approach the axis, complicating the flow. Care is thus taken in the design of the facility to provide as near as possible to 1-D flow at the exit of both center and coflow nozzles, and to minimize the strength of waves generated at the nozzle exit. Results from this experiment are compared to CFD solutions obtained by VULCAN, a previously developed code used in engine analysis. The second experiment is a study of a supersonic combustor consisting of a diverging duct with single downstream-angled wall injector. Thus, the geometry is relatively simple and large regions of subsonic recirculating flow are avoided. The nominal entrance Mach number is 2 and the enthalpy of the test gas (hot air "simulant") is nominally that of Mach 7 flight. It was believed, on the basis of calculations performed that this would produce mixing-limited flow, that is to say, one for which chemical reaction to equilibrium proceeds at a much greater rate than mixing. It later proved that this was not the case. The primary experimental technique employed is coherent anti-Stokes Raman spectroscopy, known by its acronym CARS. The species probed is molecular nitrogen and the quantity measured is temperature. Intrusive probes, such as Pitot, total temperature, hot-wire, etc., are not used due to access difficulty and high heat flux in the combustor, and because they may alter the flow. CARS has several advantages over other optical methods. It is a relatively mature and well-understood technique. Signal levels are relatively high and the signal is in the form of a coherent (laser) beam that can be collected through small windows. Incoherent (non-CARS) interferences are rejected by spatial filtering.

Diskin, G. S.↗
NASA NTRS1989

A weak Hamiltonian finite element method for optimal guidance of an advanced launch vehicle

A temporal finite-element method based on a mixed form of the Hamiltonian weak principle is presented for optimal control problems. The mixed form of this principle contains both states and costates as primary variables, which are expanded in terms of nodal values and simple shape functions. Time derivatives of the states and costates do not appear in the governing variational equation; the only quantities whose time derivatives appear therein are virtual states and virtual costates. Numerical results are presented for an elementary trajectory optimization problem; they show very good agreement with the exact solution along with excellent computational efficiency and self-starting capability. The feasibility of this approach for real-time guidance applications is evaluated. A simplified model for an advanced launch vehicle application that is suitable for finite-element solution is presented.

Hodges, Dewey H.↗
DOE OSTI2022

From rotating to charged black holes and back again

The mild form of the Weak Gravity Conjecture (WGC) requires higher derivative corrections to extremal charged black holes to increase their charge-to-mass ratio. This allows decay via emission of a smaller extremal black hole. In this paper, we investigate if similar constraints hold for extremal rotating black holes. We do so by considering the leading higher derivative corrections to the four-dimensional Kerr black hole and five-dimensional Myers-Perry black hole. We use a known mapping of these rotating solutions to a four-dimensional non-rotating dyonic Kaluza-Klein black hole and impose the WGC on this charged solution. Going back again to the rotating solutions, this fixes the sign of the corrections to the rotating extremality bounds. The sign of the corrections is non-universal, depending on the black hole under consideration. We argue that this is not at odds with black hole decay, because of the presence of a superradiant instability that persists in the extremal limit. When this instability is present, the WGC is implied for the four-dimensional charged black hole.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Constant-potential molecular dynamics simulations of molten salt double layers for FLiBe and FLiNaK

We report the results of constant-potential molecular dynamics simulations of the double layer interface between molten 2LiF–BeF 2 (FLiBe) and 23LiF–6NaF–21KF (FLiNaK) fluoride mixtures and idealized solid electrodes. Employing methods similar to those used in studies of chloride double layers, we compute the structure and differential capacitance of molten fluoride electric double layers as a function of applied voltage. The role of molten salt structure is probed through comparisons between FLiBe and FLiNaK, which serve as models for strong and weak associate-forming salts, respectively. In FLiBe, screening involves changes in Be–F–Be angles and alignment of the oligomers parallel to the electrode, while in FLiNaK, the electric field is screened mainly by rearrangement of individual ions, predominantly the polarizable potassium cation.

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

Comparison of stellar populations in simulated and real post-starburst galaxies in MaNGA

ABSTRACT Recent integral field spectroscopic (IFS) surveys have revealed radial gradients in the optical spectral indices of post-starburst (PSB) galaxies, which can be used to constrain their formation histories. We study the spectral indices of post-processed mock IFS datacubes of binary merger simulations, carefully matched to the properties of the MaNGA IFS survey, with a variety of black hole (BH) feedback models, progenitor galaxies, orbits, and mass ratios. Based on our simulation sample, we find that only major mergers on prograde–prograde or retrograde–prograde orbits in combination with a mechanical BH feedback model can form galaxies with weak enough ongoing star formation, and therefore absent H α emission, to be selected by traditional PSB selection methods. We find strong fluctuations in nebular emission line strengths, even within the PSB phase, suggesting that H α selected PSBs are only a subsample of the underlying population. The global PSB population can be more robustly identified using stellar continuum-based approaches. The difficulty in reproducing the very young PSBs in simulations potentially indicates that new sub-resolution star formation recipes are required to properly model the process of star formation quenching. In our simulations, we find that the starburst peaks at the same time at all radii, but is stronger and more prolonged in the inner regions. This results in a strong time evolution in the radial gradients of the spectral indices that can be used to estimate the age of the starburst without reliance on detailed star formation histories from spectral synthesis models.

Zheng, Yirui↗
DOE OSTI2025

Robust preparation of ground state phases under noisy imaginary time evolution

Nonunitary state preparation protocols such as imaginary time evolution (ITE) offer substantial advantages relative to unitary ones, including the ability to prepare certain long-range correlated states more efficiently. Here, we ask whether such protocols are also robust to noise arising due to coupling to the environment. We consider a nonunitary ITE “circuit” subjected to a variety of noise models and investigate whether the resulting steady state remains in the same phase as the target state of the ITE at finite noise strength. Taking the one-dimensional quantum Ising model as a concrete example, we find that the ground-state order and associated phase transition persist in the presence of noise, provided the noise does not explicitly break the symmetry that protects the phase transition. That is, the noise must possess the protecting symmetry in a weak (or average) form. Our analysis is facilitated by a mapping to an effective Hamiltonian picture in a doubled Hilbert space. We discuss possible implications of these findings for quantum simulation on noisy quantum hardware. Published by the American Physical Society 2025

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2023

PW2 Pumping: March 18 thru March 22, 2020

Processed Coherence-length-gated Microwave Photonics Interferometry (CMPI) distributed vertical strain on March 18 through March 22, 2020 at the hydromechanical test facility located at the Clemson University Simpson Station research facility in Central, South Carolina, USA. This includes ambient data in addition to a pumping test performed between 14:31-20:31 EST on 3/22/2020 from Pumping Well 2 (PW2). These data are from a strain ribbon that was deployed in the vadose zone at approximately 34.670130°N, -82.729524°E from the ground surface to a depth of 8m in saprolite containing weak reflectors pairs were placed every 0.85 m along the fiber. Each weak reflector pair forms a 15 cm cavity strain meter. Strain unit: meters/meter (strain). Sampling rate 0.1 Hz. Each column corresponds to the strain data taken at a depth d from the ground surface. c1: d = 1.00m c2: d = 1.85m c3: d = 2.70m c4: d = 3.55m (large demodulation error) c5: d = 4.40m c6: d = 5.25m c7: d = 6.10m (large demodulation error) c8: d = 6.95m

Distributed Strain Sensing↗
NASA NTRS1972

Solar large-scale positive polarity magnetic fields and geomagnetic disturbances

Unlike the negative polarity solar magnetic field large-scale regular features that correlate with enhanced solar activity regions, the positive polarity regular formations formed in the weak and old background magnetic fields seem to correlate well with geomagnetically enhanced periods of time (shifted for 4 days), which means that they seem to be the source of the quiet solar wind. This behavior of the large intervals of heliographic longitude with prevailing positive polarity fields may be followed to the end of the 18th cycle, during the declining part of the 19th cycle, and during the first half of the present 20th cycle of solar activity.

Bumba, V.↗