Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Extended Lagrangian dynamics”

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 55 records · Page 3

Physics-informed machine learning of the Lagrangian dynamics of velocity gradient tensor

Reduced models describing the Lagrangian dynamics of the velocity gradient tensor (VGT) in homogeneous isotropic turbulence (HIT) are developed under the physics-informed machine learning (PIML) framework. We consider the VGT at both Kolmogorov scale and coarse-grained scale within the inertial range of HIT. Building reduced models requires resolving the pressure Hessian and subfilter contributions, which is accomplished by constructing them using the integrity bases and invariants of the VGT. The developed models can be expressed using the extended tensor basis neural network (TBNN) introduced by Ling et al. [J. Fluid Mech. 807, 155 (2016)]. Physical constraints, such as Galilean invariance, rotational invariance, and incompressibility condition, are thus embedded in the models explicitly. Our PIML models are trained on the Lagrangian data from a high-Reynolds number direct numerical simulation (DNS). To validate the results, we perform a comprehensive out-of-sample test. We observe that the PIML model provides an improved representation for the magnitude and orientation of the small-scale pressure Hessian contributions. Statistics of the flow, as indicated by the joint PDF of second and third invariants of the VGT, show good agreement with the “ground-truth” DNS data. A number of other important features describing the structure of HIT are reproduced by the model successfully. We have also identified challenges in modeling inertial range dynamics, which indicates that a richer modeling strategy is required. This helps us identify important directions for future research, in particular towards including inertial range geometry into the TBNN.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Dynamic Analysis and Control of Lightweight Manipulators with Flexible Parallel Link Mechanisms

The objective is the theoretical analysis and the experimental verification of dynamics and control of a two link flexible manipulator with a flexible parallel link mechanism. Nonlinear equations of motion of the lightweight manipulator are derived by the Lagrangian method in symbolic form to better understand the structure of the dynamic model. The resulting equation of motion have a structure which is useful to reduce the number of terms calculated, to check correctness, or to extend the model to higher order. A manipulator with a flexible parallel link mechanism is a constrained dynamic system whose equations are sensitive to numerical integration error. This constrained system is solved using singular value decomposition of the constraint Jacobian matrix. Elastic motion is expressed by the assumed mode method. Mode shape functions of each link are chosen using the load interfaced component mode synthesis. The discrepancies between the analytical model and the experiment are explained using a simplified and a detailed finite element model.

Lee, Jeh Won↗

Application of a novel finite difference method to dynamic crack problems

A versatile finite difference method (HEMP and HEMP 3D computer programs) was developed originally for solving dynamic problems in continuum mechanics. It was extended to analyze the stress field around cracks in a solid with finite geometry subjected to dynamic loads and to simulate numerically the dynamic fracture phenomena with success. This method is an explicit finite difference method applied to the Lagrangian formulation of the equations of continuum mechanics in two and three space dimensions and time. The calculational grid moves with the material and in this way it gives a more detailed description of the physics of the problem than the Eulerian formulation.

Chen, Y. M.↗

Reduced scale stripline platform to extend accessible pressures on the Z machine: Shockless compression of platinum to 650 GPa

Reaching astrophysically relevant high energy density (HED) material states in the laboratory is an ongoing effort at multiple experimental facilities. We have developed a new dynamic compression platform for the Z Pulsed Power Facility that allows for sample sizes 100s of [Formula: see text]m in thickness that accommodate multiple grains in order to fully capture bulk properties, such as material strength. A pair of experiments compressed platinum (Pt) to HED conditions and conventional inverse Lagrangian analysis as well as a recent Bayesian calibration technique were used to determine the principal isentrope to 650 GPa with density uncertainties of <2%. These low uncertainties are calculated for single sample experiments, presenting the possibility of even smaller experimental uncertainties with multiple samples the platform allows. Our new platform extends the accessible Pt ramp pressures achievable on the Z machine to over 80% of the pressure recently achieved using the National Ignition Facility planar Hohlraum platform. This new capability, the next generation evolution of the stripline platform, was made possible by advancements in both our understanding of the Z pulsed power driver and our overall magnetohydrodynamic modeling capabilities.

Porwitzky, Andrew↗

Population balance modeling of polyurethane foam formation with pressure‐dependent growth kernel

Abstract Polyurethane foams are widely used materials often chosen for their useful characteristics such as low thermal conductivity, ease of application, and high strength‐to‐weight ratios. Computational models are needed to predict the dynamics of the flow and expansion, and the resulting material properties, to improve manufacturing processes. In this paper, a model for PMDI, a water‐blown polyurethane foam, is presented. By extending a kinetics‐based approach by adding bubble‐scale information via a population balance equation (PBE) using the quadrature method of moments, we can track bubble size distributions during foaming. We present results from a three‐dimensional computational fluid dynamics model using arbitrary Lagrangian–Eulerian interface tracking implemented in finite element software. The model compares favorably with experimental data, including dynamics, bubble distributions measured by both camera and diffusion wave spectroscopy, and post‐test bubble size from scanning electron microscopy and density measurements from x‐ray computed tomography.

Ortiz, Weston↗

Bulk viscosity of two-flavor color superconducting quark matter in neutron star mergers

This work investigates the bulk viscosity of warm, dense, neutrino-transparent, color-superconducting quark matter, where damping of density oscillations in the kHz frequency range arises from weak-interaction-driven direct Urca processes involving quarks. We study the two-flavor red-green paired color-superconducting (2SC) phase, while allowing for the presence of unpaired strange quarks and blue color light quarks of all flavors. Our calculations are based on the SU(3) Nambu–Jona–Lasinio model, extended to include both vector interactions and the ‘t Hooft determinant term. The primary focus is on how variations in the NJL Lagrangian parameters—specifically, the diquark and vector coupling strengths—affect both the static properties of quark matter, such as its equation of state and composition, and its dynamical behavior, including bulk viscosity and associated damping timescales. We find that the bulk viscosity and corresponding damping timescale can change by more than an order of magnitude upon varying the vector coupling by a factor of two at high densities and by a lesser degree at lower densities. This sensitivity primarily arises from the susceptibility of 2SC matter, with a smaller contribution from modifications to the weak interaction rates. In comparison, changes in the diquark coupling have a more limited impact. The damping of density oscillations in 2SC matter is similar quantitatively to nucleonic matter and can be a leading mechanism of dissipation in merging hybrid stars containing color superconducting cores.

Gravitational Waves↗

Intermittent Frontogenesis in the Alboran Sea

Abstract We present a phenomenological description and dynamical analysis of the Alboran fronts using a realistic simulation at submesoscale resolution. The study is focused on east Alboran fronts emerging within relatively strong flows that separate from the Spanish coast into the basin interior. Despite modest lateral shifting associated with the position of the Alboran anticyclonic gyres and variations in intensity, these fronts present a similar structure and dynamical configuration as the climatological Almeria–Oran front. The statistical analysis of our solution shows that strained-induced frontogenesis is a recurrent submesoscale mechanism associated with these fronts, and the process is assessed in terms of the advective Lagrangian frontogenetic tendencies associated with buoyancy and velocity horizontal gradients. Intermittency in their strength and patterns is indicative of high variability in the occurrence of active frontogenesis in association with the secondary (overturning) circulation across the frontal gradient. As a result, we find many episodes with strong surface fronts that do not have much associated downwelling. Frontogenesis and the associated secondary circulation are further explored during two particular frontal events, both showing strong downwelling of (1) cm s −1 extending down into the pycnocline. A frontogenetic contribution of turbulent vertical momentum mixing to the secondary circulation is identified in the easternmost region during the cold season, when the dynamics are strongly influenced by the intrusion of the salty Northern Current. The background vertical velocity fields observed during the analyzed events indicate other currents in the submesoscale range, including tidal and topographic internal waves.

Capó, Esther↗

Stability of trapped solutions of a nonlinear Schrödinger equation with a nonlocal nonlinear self-interaction potential

This work focuses on the study of the stability of trapped soliton-like solutions of a (1 + 1)-dimensional nonlinear Schrödinger equation (NLSE) in a nonlocal, nonlinear, self-interaction potential of the form [|Ψ(x,t)| 2 +|Ψ(-x,t)| 2 ] κ where κ is an arbitrary nonlinearity parameter. Although the system with κ = 1 (i.e. fully integrable case) was first reported by Yang (2018 Phys. Rev. E 98 042202), here in the present work, we extend this model to the one in which κ is arbitrary. This allows us to compare the stability properties of the now trapped solutions to previously found solutions of the more usual NLSE with κ ≠ 1 which are moving soliton solutions. We show that there is a simple, one-component, nonlocal Lagrangian and corresponding action governing the dynamics of the system. Using a collective coordinate method derived from the action as well as assuming the validity of Derrick's theorem, we find that these trapped solutions are stable for 0 < κ < 2 and unstable when κ > 2. At the critical value of κ, i.e. κ = 2, the solution can either collapse or blowup linearly in time when q 0 = 0, where q 0 is the center of the initial density ρ(x, t = 0) = ψ*ψ of the solution. For q 0 ≠ 0 the displaced solution collapses. When κ > 2 initial small displacements from the origin also lead to collapse of the wave function. This phenomenon is not seen in the usual NLSE.

collective coordinates↗

Multivariable control of the Space Shuttle Remote Manipulator System using linearization by state feedback

This paper demonstrates an approach to nonlinear control system design that uses linearization by state feedback to allow faster maneuvering of payloads by the Shuttle Remote Manipulator System (SRMS). A nonlinear feedback law is defined to cancel the nonlinear plant dynamics so that a linear controller can be designed for the SRMS. First a nonlinear design model was generated via SIMULINK. This design model included nonlinear arm dynamics derived from the Lagrangian approach, linearized servo model, and linearized gearbox model. The current SRMS position hold controller was implemented on this system. Next, a trajectory was defined using a rigid body kinematics SRMS tool, KRMS. The maneuver was simulated. Finally, higher bandwidth controllers were developed. Results of the new controllers were compared with the existing SRMS automatic control modes for the Space Station Freedom Mission Build 4 Payload extended on the SRMS.

Gettman, Chang-Ching L.↗

Coupling of Sph and Finite Element Codes for Multi-Layer Orbital Debris Shield Design

Particle-based hydrodynamics models offer distinct advantages over Eulerian and Lagrangian hydrocodes in particular shock physics applications. Particle models are designed to avoid the mesh distortion and state variable diffusion problems which can hinder the effective use of Lagrangian and Eulerian codes respectively. However conventional particle-in-cell and smooth particle hydrodynamics methods employ particles which are actually moving interpolation points. A new particle-based modeling methodology, termed Hamiltonian particle hydrodynamics, was developed by Fahrenthold and Koo (1997) to provide an alternative, fully Lagrangian, energy-based approach to shock physics simulations. This alternative formulation avoids the tensile and boundary instabilities associated with standard smooth particle hydrodynamics formulations and the diffusive grid- to-particle mapping schemes characteristic of particle-in-cell methods. In the work described herein, the method of Fahrenthold and Koo has been extended, by coupling the aforementioned hydrodynamic particle model to a hexahedral finite element based description of the continuum dynamics. The resulting continuum model retains all of the features (including general contact-impact effects) of Hamiltonian particle hydrodynamics, while in addition accounting for tensile strength, plasticity, and damage effects important in the simulation of hypervelocity impact on orbital debris shielding. A three dimensional, vectorized, and autotasked implementation of the extended particle method described here has been coded for application to orbital debris shielding design. Source code for the pre-processor (PREP), analysis code (EXOS), post-processor (POST), and rezoner (ZONE), have been delivered separately, along with a User's Guide describing installation and application of the software.

Fahrenthold, Eric P.↗

Distributed Constrained Optimization with Semicoordinate Transformations

Recent work has shown how information theory extends conventional full-rationality game theory to allow bounded rational agents. The associated mathematical framework can be used to solve constrained optimization problems. This is done by translating the problem into an iterated game, where each agent controls a different variable of the problem, so that the joint probability distribution across the agents moves gives an expected value of the objective function. The dynamics of the agents is designed to minimize a Lagrangian function of that joint distribution. Here we illustrate how the updating of the Lagrange parameters in the Lagrangian is a form of automated annealing, which focuses the joint distribution more and more tightly about the joint moves that optimize the objective function. We then investigate the use of "semicoordinate" variable transformations. These separate the joint state of the agents from the variables of the optimization problem, with the two connected by an onto mapping. We present experiments illustrating the ability of such transformations to facilitate optimization. We focus on the special kind of transformation in which the statistically independent states of the agents induces a mixture distribution over the optimization variables. Computer experiment illustrate this for &sat constraint satisfaction problems and for unconstrained minimization of NK functions.

Macready, William↗

A Lagrangian Study of Interfaces at the Edges of Cumulus Clouds

Abstract Interfaces at the edge of an idealized, nonprecipitating, warm cloud are studied using direct numerical simulation (DNS) complemented with a Lagrangian particle tracking routine. Once a shell has formed, four zones can be distinguished: the cloud core, visible shell, invisible shell, and the environment. The union of the visible and invisible regions is the shell commonly referred to in literature. The boundary between the invisible shell and the environment is the turbulent–nonturbulent interface (TNTI), which is typically not considered in cloud studies. Three million particles were seeded homogeneously across the domain and properties were recorded along individual trajectories. The results demonstrate that the traditional cloud boundary (separating cloudy and noncloudy regions using thresholds applied on liquid condensate or updraft velocity) are some distance away from the TNTI. Furthermore, there is no dynamic difference between the traditional liquid-condensate boundary and the region extending to the TNTI. However, particles crossing the TNTI exhibit a sharp jump in enstrophy and a smooth increase in buoyancy. The traditional cloud boundary coincides with the location of minimum buoyancy in the shell. The shell premixes the entraining and detraining air and analysis reveals a highly skewed picture of entrainment and detrainment at the traditional cloud boundary. A preferential entrainment of particles with velocity and specific humidity higher than the mean values in the shell is observed. Large-eddy simulation of a more realistic setup detects an interface with similar properties using the same thresholds as in the DNS, indicating that the DNS results extrapolate beyond their idealized conditions.

54 ENVIRONMENTAL SCIENCES↗

Background-field Method and QCD Factorization

One method for deriving a factorization for QCD processes is to use successive integration over fields in the functional integral. In this approach, we separate the fields into two categories: dynamical fields with momenta above a relevant cutoff, and background fields with momenta below the cutoff. The dynamical fields are then integrated out in the background of the low-momentum background fields. This strategy works well at tree level, allowing us to quickly derive QCD factorization formulas at leading order. However, to extend the approach to higher loops, it is necessary to rigorously define the functional integral over dynamical fields in an arbitrary background field. This framework was carefully developed for the calculation of the effective action in a background field at the two-loop level in the classic paper by Abbott «The Background Field Method Beyond One Loop», Nucl. Phys. B 185 , 189 (1981). Building on this work, I specify the renormalized background-field Lagrangian and define the notion of the quantum average of an operator in a background field, consistent with the “separation of scales” scheme mentioned earlier. As examples, I discuss the evolution of the twist-2 gluon light-ray operator and the one-loop gluon propagator in a background field near the light cone.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

The completed SDSS-IV extended Baryon Oscillation Spectroscopic Survey: N-body mock challenge for the eBOSS emission line galaxy sample

ABSTRACT Cosmological growth can be measured in the redshift space clustering of galaxies targeted by spectroscopic surveys. Accurate prediction of clustering of galaxies will require understanding galaxy physics, which is a very hard and highly non-linear problem. Approximate models of redshift space distortion (RSD) take a perturbative approach to solve the evolution of dark matter and galaxies in the universe. In this paper, we focus on extended Baryon Oscillation Spectroscopic (eBOSS) emission line galaxies (ELGs) that live in intermediate mass haloes. We create a series of mock catalogues using haloes from the Multidark and outer rim dark matter only N-body simulations. Our mock catalogues include various effects inspired by baryonic physics such as assembly bias and the characteristics of satellite galaxies kinematics, dynamics, and statistics deviating from dark matter particles. We analyse these mocks using the TNS RSD model in Fourier space and the convolution Lagrangian perturbation theory (CLPT) in configuration space. We conclude that these two RSD models provide an unbiased measurement of RSD within the statistical error of our mocks. We obtain the conservative theoretical systematic uncertainty of $3.3{{\ \rm per\ cent}}$, $1.8{{\ \rm per\ cent}}$, and $1.5{{\ \rm per\ cent}}$ in fσ8, α∥, and α⊥, respectively, for the TNS and CLPT models. We note that the estimated theoretical systematic error is an order of magnitude smaller than the statistical error of the eBOSS ELG sample and hence are negligible for the purpose of the current eBOSS ELG analysis.

Alam, Shadab↗

Froggatt-Nielsen meets the SMEFT

We study the matching of Froggatt-Nielsen theories of flavour onto the Standard Model Effective Field Theory (SMEFT), upon integrating out a heavy Beyond-the-Standard-Model (BSM) scalar ‘flavon’ whose vacuum expectation value breaks an Abelian flavour symmetry at energies Λ FN well above the electroweak scale, Λ FN > Λ SM . We include matching contributions to the infrared d SM = 6 (Warsaw basis) SMEFT sourced from ultraviolet contact terms suppressed up to order 1/${\Lambda}_{\textrm{UV}}^2$ in the Froggatt-Nielsen Lagrangian, where Λ UV > Λ FN is an arbitrary ultraviolet scale where further unspecified BSM particles are dynamical. This includes tree-level (one-loop) ultraviolet diagrams with d FN = 6 (5) effective vertices. We first do so with a toy model, but then generalize our findings to arbitrary Frogatt-Nielsen charges. Our results indicate a rich and non-trivial signature of Froggatt-Nielsen theories on the (otherwise) model-independent operators of the SMEFT, and we briefly speculate on extending our analysis to broader classes of BSM flavour models, e.g. non-Abelian and/or gauged theories. We thus take an important step towards determining how to use rapidly developing theoretical and experimental SMEFT technologies to gain unambiguous insight into the SM’s longstanding fermion flavour puzzle.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Contributions of biofilm-induced flow heterogeneities to solute retention and anomalous transport features in porous media

We report that microbial biofilms are ubiquitous within porous media and the dynamics of their growth influence surface and subsurface flow patterns which impacts the physical properties of porous media and large-scale transport of solutes. Here, a two-dimensional pore-scale numerical model was used to evaluate the impact of biofilm-induced flow heterogeneities on conservative transport. Our study integrates experimental biofilm images of Paenibacillus 300A strain in a microfluidic device packed with cylindrical grains in a hexagonal distribution, with mathematical modeling. Biofilm is represented as a synthetic porous structure with locally varying physical properties that honors the impact of biofilm on the porous medium. We find that biofilm plays a major role in shaping the observed conservative transport dynamics by enhancing anomalous characteristics. More specifically, when biofilm is present, the pore structure in our geometry becomes more spatially correlated. We observe intermittent behavior in the Lagrangian velocities that switches between fast transport periods and long trapping events. Our results suggest that intermittency enhances solute spreading in breakthrough curves which exhibit extreme anomalous slope at intermediate times and very marked late solute arrival due to solute retention. The efficiency of solute retention by the biofilm is controlled by a transport regime which can extend the tailing in the concentration breakthrough curves. These results indicate that solute retention by the biofilm exerts a strong control on conservative solute transport at pore-scale, a role that to date has not received enough attention.

54 ENVIRONMENTAL SCIENCES↗

Mature thunderstorm cloud top structure - Three-dimensional numerical simulation versus satellite observations

Preliminary results of four runs with a three-dimensional model of the effects of vertical wind shear on cloud top height/temperature structure and the internal properties of isolate midlatitude thunderstorms are reported. The model is being developed as an aid to analyses of GEO remote sensing satellite data. The grid is a 27 x 27 x 20 mesh with 2 km horizontal resolution and 0.9 vertical resolution. The total grid is 54 km on a side and 18 km deep. A second-order Crowley scheme for advecting momentum is extended with a third-order correction for spatial truncation error, and the earth-relative horizontal surface wind components are decreased to 50 percent of their values at 0.45 km. A temperature increase with height is included, together with an initial impulse consisting of a nonrotating cylindrical weak buoyant updraft 10 km in radius. The results of the runs are discussed in terms of the time variation of the vertical velocity extrema, the effects of strong and weak shear on a storm, the cloud top height, the Lagrangian dynamics of a thermal couplet, and data from a real storm.

Schlesinger, R. E.↗

Thermodynamically consistent semi-compressible fluids: a variational perspective

This paper presents (Lagrangian) variational formulations for single and multicomponent semi-compressible fluids with both reversible (entropy-conserving) and irreversible (entropy-generating) processes. Semi-compressible fluids are useful in describing low-Mach dynamics, since they are soundproof. These models find wide use in many areas of fluid dynamics, including both geophysical and astrophysical fluid dynamics. Specifically, the Boussinesq, anelastic and pseudoincompressible equations are developed through a unified treatment valid for arbitrary Riemannian manifolds, thermodynamic potentials and geopotentials. By design, these formulations obey the 1st and 2nd laws of thermodynamics, ensuring their thermodynamic consistency. This general approach extends and unifies existing work, and helps clarify the thermodynamics of semi-compressible fluids. To further this goal, evolution equations are presented for a wide range of thermodynamic variables: entropy density s, specific entropy η, buoyancy b, temperature T, potential temperature θ and a generic entropic variable χ; along with a general definition of buoyancy valid for all three semicompressible models and arbitrary geopotentials. Finally, the elliptic equation for the pressure perturbation (the Lagrange multiplier that enforces semi-compressibility) is developed for all three equation sets in the case of reversible dynamics, and for the Boussinesq/anelastic equations in the case of irreversible dynamics; and some discussion is given of the difficulty in formulating it for the pseudoincompressible equations with irreversible dynamics.

97 MATHEMATICS AND COMPUTING↗