Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Nonlinear models”

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 343 records · Page 19
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Estimating Leaf Area Index in Row Crops Using Wheel-Based and Airborne Discrete Return Light Detection and Ranging Data

Leaf area index (LAI) is an important variable for characterizing plant canopy in crop models. It is traditionally defined as the total one-sided leaf area per unit ground area and is estimated by both direct and indirect methods. This paper explores the effectiveness of using light detection and ranging (LiDAR) data to estimate LAI for sorghum and maize with different treatments at multiple times during the growing season from both a wheeled vehicle and Unmanned Aerial Vehicles. Linear and nonlinear regression models are investigated for prediction utilizing statistical and plant structure-based features extracted from the LiDAR point cloud data with ground reference obtained from an in-field plant canopy analyzer (indirect method). Results based on the value of the coefficient of determination ( R 2 ) and root mean squared error for predictive models ranged from ∼0.4 in the early season to ∼0.6 for sorghum and ∼0.5 to 0.80 for maize from 40 Days after Sowing to harvest.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2024

Analysis of DESI × DES using the Lagrangian effective theory of LSS

In this work we use Lagrangian perturbation theory to analyze the harmonic space galaxy clustering signal of the Bright Galaxy Survey (BGS) and luminous red galaxies (LRGs) targeted by the dark energy spectroscopic instrument (DESI), combined with the galaxy-galaxy lensing signal measured around these galaxies using Dark Energy Survey Year 3 source galaxies. The BGS and LRG galaxies are extremely well characterized by DESI spectroscopy and, as a result, lens galaxy redshift uncertainty and photometric systematics contribute negligibly to the error budget of our “2 × 2-point” analysis. On the modeling side, this work represents the first application of the SPINOSAURUS code, implementing an effective field theory model for galaxy intrinsic alignments, and we additionally introduce a new scheme ( MAIAR ) for marginalizing over the large uncertainties in the redshift evolution of the intrinsic alignment signal. Furthermore, this is the first application of a hybrid effective field theory model for galaxy bias based on the A emulus 𝜈 simulations. Our main result is a measurement of the amplitude of the lensing signal, 𝑆 8 = 𝜎 8 ⁢(Ω 𝑚 /0.3) 0.5 = 0.85⁢0$^{+0.042}_{−0.050}$, consistent with values of this parameter derived from the primary cosmic microwave background. This constraint is artificially improved by a factor of 51% if we assume a more standard, but restrictive parametrization for the redshift evolution and sample dependence of the intrinsic alignment signal, and 63% if we additionally assume the nonlinear alignment model. Furthermore, we show that when fixing the cosmological model to the best-fit values from Planck PR4 there is > 5⁢𝜎 evidence for a deviation of the evolution of the intrinsic alignment signal from the functional form that is usually assumed in cosmic shear and galaxy-galaxy lensing studies.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

Multilevel well modeling in aggregation-based nonlinear multigrid for multiphase flow in porous media

A full approximation scheme (FAS) nonlinear multigrid solver for two-phase flow and transport problems driven by wells with multiple perforations is developed here. It is an extension to our previous work on FAS solvers for diffusion and transport problems. The solver is applicable to discrete problems defined on unstructured grids as the coarsening algorithm is aggregation-based and algebraic. To construct coarse basis that can better capture the radial flow near wells, coarse grids in which perforated well cells are not near the coarse-element interface are desired. This is achieved by an aggregation algorithm proposed in this paper that makes use of the location of well cells in the cell-connectivity graph. Numerical examples in which the FAS solver is compared against Newton's method on benchmark problems are given. In particular, for a refined version of the SAIGUP model, the FAS solver is at least 35% faster than Newton's method for time steps with a CFL number greater than 10.

58 GEOSCIENCES↗
DOE OSTI2021

Convex Optimization of Integrated Power-Gas Energy Flow Model With Applications to Probabilistic Energy Flow

Energy flow calculation is a fundamental problem of the integrated power and gas system (IPGS) operation and planning. However, the nonlinear gas flow model introduces major challenges to the energy flow calculation. In this paper, we propose a tractably convex optimization model to solve the energy flow problem in IPGSs. It is demonstrated that the proposed optimization model has the same optimal solution as the original nonlinear steady energy flow model. Also, piecewise linearization is adopted to tightly linearize the nonlinear objective function of the model, which transforms the formulated convex optimization into a linear program one. Thus, the computation complexity of the proposed energy flow model is significantly reduced as compared with the existing methods. In addition, the proposed model can be extended to probabilistic energy flow estimation. Extensive case studies are conducted to validate the effectiveness of the proposed energy flow model using three IPGSs.

42 ENGINEERING↗
DOE OSTI2021

Full Vorticity Budget of the Arabian Sea from a 0.1° Ocean Model: Sverdrup Dynamics, Rossby Waves, and Nonlinear Eddy Effects

Abstract The Arabian Sea, influenced by the Indian monsoon, has many unique features, including its basin-scale seasonally reversing surface circulation and the Great Whirl, a seasonal anticyclonic system appearing during the southwest monsoon close to the western boundary. To establish a comprehensive dynamical picture of the Arabian Sea, we utilize numerical model output and design a full vorticity budget that includes a fully decomposed nonlinear term. The ocean general circulation model has 0.1° resolution and is mesoscale eddy-resolving in the region. In the western boundary current system, we highlight the role of nonlinear eddies in the life cycle of the Great Whirl. The nonlinear eddy term is of leading-order importance in this feature’s vorticity balance. Specifically, it contributes to the Great Whirl’s persistence in boreal fall after the weakening of the southwesterly winds. In the open ocean, Sverdrup dynamics and annual Rossby waves are found to dominate the vorticity balance; the latter is considered as a key factor in the formation of the Great Whirl and the seasonal reversal of the western boundary current. In addition, we discuss different forms of vertically integrated vorticity equations in the model and argue that the bottom pressure torque term can be interpreted analogously as friction in the western boundary and vortex stretching in the open ocean.

58 GEOSCIENCES↗
DOE OSTI2022

Extending the Heat Index

The Heat Index is a widely used measure of apparent temperature that accounts for the effects of humidity using Steadman’s model of human thermoregulation. Steadman’s model, however, gives unphysical results when the air is too hot and humid or too cold and dry, leading to an undefined Heat Index. For example, at a relative humidity of 80%, the Heat Index is only defined for temperatures in the range of 288-304 K (59-88 degrees Fahrenheit). Here, Steadman’s thermoregulation model is extended to define the Heat Index for all combinations of temperature and humidity, allowing for an assessment of Earth’s future habitability. The extended Heat Index can be mapped onto physiological responses of an idealized human, such as heat exhaustion, heat stroke, and even heat death, providing an indication of regional health outcomes for different degrees of global warming.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A sharp interface Lagrangian-Eulerian method for flexible-body fluid-structure interaction

This paper introduces a sharp-interface approach to simulating fluid-structure interaction (FSI) involving flexible bodies described by general nonlinear material models and across a broad range of mass density ratios. This new flexible-body immersed Lagrangian-Eulerian (ILE) scheme extends our prior work on integrating partitioned and immersed approaches to rigid-body FSI. Our numerical approach incorporates the geometrical and domain solution flexibility of the immersed boundary (IB) method with an accuracy comparable to body-fitted approaches that sharply resolve flows and stresses up to the fluid-structure interface. Unlike many IB methods, our ILE formulation uses distinct momentum equations for the fluid and solid subregions with a Dirichlet-Neumann coupling strategy that connects fluid and solid subproblems through simple interface conditions. As in earlier work, we use approximate Lagrange multiplier forces to treat the kinematic interface conditions along the fluid-structure interface. This penalty approach simplifies the linear solvers needed by our formulation by introducing two representations of the fluid-structure interface, one that moves with the fluid and another that moves with the structure, that are connected by stiff springs. This approach also enables the use of multi-rate time stepping, which allows us to use different time step sizes for the fluid and structure subproblems. Our fluid solver relies on an immersed interface method (IIM) for discrete surfaces to impose stress jump conditions along complex interfaces while enabling the use of fast structured-grid solvers for the incompressible Navier-Stokes equations. The dynamics of the volumetric structural mesh are determined using a standard finite element approach to large-deformation nonlinear elasticity via a nearly incompressible solid mechanics formulation. This formulation also readily accommodates compressible structures with a constant total volume, and it can handle fully compressible solid structures for cases in which at least part of the solid boundary does not contact the incompressible fluid. Selected grid convergence studies demonstrate second-order convergence in volume conservation and in the pointwise discrepancies between corresponding positions of the two interface representations as well as between first and second-order convergence in the structural displacements. The time stepping scheme is also demonstrated to yield second-order convergence. To assess and validate the robustness and accuracy of the new algorithm, comparisons are made with computational and experimental FSI benchmarks. Test cases include both smooth and sharp geometries in various flow conditions. Furthermore, we also demonstrate the capabilities of this methodology by applying it to model the transport and capture of a geometrically realistic, deformable blood clot in an inferior vena cava filter.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2021

Flux-driven algebraic damping of m = 2 diocotron mode

Experiments with pure electron plasmas in a Malmberg–Penning trap have observed linear in time, algebraic damping of m = 2 diocotron modes. Transport due to small field asymmetries produces a low-density halo of electrons moving radially outward from the plasma core, and the mode damping begins when the halo reaches the resonant radius of the mode. The damping rate is proportional to the flux of halo particles through the resonant layer. The damping is related to, but distinct from spatial Landau damping in which a linear wave–particle resonance produces exponential damping. This paper reports an analytic theory that captures the main signatures reported for this novel damping, namely, that the damping begins when the halo particles reach the resonant radius and that the damping is algebraic in time with nearly constant damping rate. The model also predicts a nonlinear frequency shift. The model provides two ways to think about the damping. It results from a transfer of canonical angular momentum from the mode to halo particles being swept by the mode field through the nonlinear cat's eye orbits of the resonant region. More mechanistically, the electric field produced by the perturbed charge density of the resonant particles acts back on the plasma core causing E × B drift that gives rise to the damping and nonlinear frequency shift.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Nonlinear programming optimization of a single-stack electrodialysis desalination system for cost efficiency

Electrodialysis (ED) presents a competitive method for desalinating brackish waters. In this work, we perform cost optimization of a single-stack ED system across a range of feed salinities and water recoveries while optimizing operating voltage, number of cell pairs, and cell length. The results of our optimization show that the levelized cost of water (LCOW) increases with an increase in feed salinity. The outcomes of our optimization show that cost-optimal design generally increases cell length while decreasing cell pair number and operating voltage with an increase in salinity. These trends are nonlinear, with the number of cell pairs and applied voltage exhibiting local maxima when operating at low salinity and high recovery. We discuss the underlying mechanism for cell length becoming a leveraging design parameter by inspecting the length-dependent profiles of key electrochemical properties of the ED cell. Finally, we present how increasing performance metrics and decreasing costs impact LCOW, demonstrating that innovations that decrease counter-current diffusion have resulted in the highest decrease of LCOW.

42 ENGINEERING↗
DOE OSTI2022

Spin-piston problem for a ferromagnetic thin film: Shock waves and solitons

In this report the unsteady, nonlinear magnetization dynamics induced by spin injection in an easy-plane ferromagnetic channel subject to an external magnetic field are studied analytically. Leveraging a dispersive hydrodynamic description, the Landau-Lifshitz equation is recast in terms of hydrodynamic-like variables for the magnetization's perpendicular component (spin density) and azimuthal phase gradient (fluid velocity). Spin injection acts as a moving piston that generates nonlinear, dynamical spin textures in the ferromagnetic channel with downstream quiescent spin density set by the external field. In contrast to the classical problem of a piston accelerating a compressible gas, here, variable spin injection and field lead to a rich variety of nonlinear wave phenomena from oscillatory spin shocks to solitons and rarefaction (expansion) waves. A full classification of solutions is provided using nonlinear wave modulation theory by identifying two key aspects of the fluid-like dynamics: subsonic/supersonic conditions and convex/nonconvex hydrodynamic flux. Familiar waveforms from the classical piston problem such as rarefaction waves and shocks manifest in their spin-based counterparts as smooth and highly oscillatory transitions, respectively, between two distinct magnetic states. The spin shock is an example of a dispersive shock wave, which arises in many physical systems. New features without a gas dynamics counterpart include composite wave complexes with “contact” spin shocks and rarefactions. Magnetic supersonic conditions lead to two pronounced piston edge behaviors including a stationary soliton and an oscillatory wave train. These coherent wave structures have physical implications for the generation of high frequency spin waves from pulsed injection and persistent, stable stationary and/or propagating solitons in the presence of magnetic damping. The analytical results are favorably compared with numerical simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2024

S-OPT: A Points Selection Algorithm for Hyper-Reduction in Reduced Order Models

While projection-based reduced order models can reduce the dimension of full order solutions, the resulting reduced models may still contain terms that scale with the full order dimension. Hyper-reduction techniques are sampling-based methods that further reduce this computational complexity by approximating such terms with a much smaller dimension. The goal of this work is to introduce the points selection algorithm developed by Shin and Xiu as a hyper-reduction method. The selection algorithm was originally proposed as a stochastic collocation method for uncertainty quantification. Since the algorithm aims at maximizing a quantity $\mathcal{S}$ that measures both the column orthogonality and the determinant, we refer to the algorithm as S-OPT. Numerical examples are provided to demonstrate the performance of S-OPT and to compare its performance with a gappy proper orthogonal decomposition (POD) algorithm. Here, we found that using the S-OPT algorithm is shown to predict the full order solutions with higher accuracy than gappy POD especially when the number of sampling points is small, although we note that S-OPT shows slow asymptotic convergence with respect to the number of samples for some applications, e.g., Lagrangian hydrodynamics.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Operator inference for non-intrusive model reduction of systems with non-polynomial nonlinear terms

Here in this work we present a non-intrusive model reduction method to learn low-dimensional models of dynamical systems with non-polynomial nonlinear terms that are spatially local and that are given in analytic form. In contrast to state-of-the-art model reduction methods that are intrusive and thus require full knowledge of the governing equations and the operators of a full model of the discretized dynamical system, the proposed approach requires only the non-polynomial terms in analytic form and learns the rest of the dynamics from snapshots computed with a potentially black-box full-model solver. The proposed method learns operators for the linear and polynomially nonlinear dynamics via a least-squares problem, where the given non-polynomial terms are incorporated on the right-hand side. The least-squares problem is linear and thus can be solved efficiently in practice. The proposed method is demonstrated on three problems governed by partial differential equations, namely the diffusion–reaction Chafee–Infante model, a tubular reactor model for reactive flows, and a batch-chromatography model that describes a chemical separation process. The numerical results provide evidence that the proposed approach learns reduced models that achieve comparable accuracy as models constructed with state-of-the-art intrusive model reduction methods that require full knowledge of the governing equations.

42 ENGINEERING↗
DOE OSTI2022

Hybrid large-eddy simulation with adaptive downstream anisotropic eddy viscosity model

Turbulence pervades most flows of engineering interest, and its prediction remains a challenge on both accuracy and cost. One promising predictive approach that reduces cost combines large-eddy simulation (LES) with simulation based on Reynolds averaged Navier–Stokes equations (RANS). This study presents a method to overcome stability and accuracy issues associated with these hybrid LES-RANS methods. The method developed involves extracting the Reynolds stresses from the upstream LES solution and then using that information to convert the downstream RANS model from a scalar eddy-viscosity model to an anisotropic nonlinear eddy-viscosity model. The method developed differs from the downstream tensorial eddy-viscosity model by being independent of the coordinate system. The method developed was evaluated by computing film cooling of a flat plate with the coolant injected through one row of circular holes. Results obtained show instabilities at the LES-to-RANS interface to be eliminated. Also, the method developed yielded solutions that compare reasonably well with those from LES, even though a significant portion of the flow is computed by the adapted anisotropic RANS model instead of LES, which significantly reduced the number of grid points and computational time needed. Since modification of the downstream RANS model is based on information extracted from the upstream LES solution, the method developed is adaptive to the problem being studied.

42 ENGINEERING↗
DOE OSTI2024

Even-point multi-loop unitarity and its applications: exponentiation, anomalies and evanescence

Abstract We identify novel structure in newly computed multi-loop amplitudes and quantum actions for even-point effective field theories, including both the nonlinear sigma model (NLSM) and double-copy gauge theories such as Born-Infeld and its supersymmetric generalizations. We exploit special properties of all even-point theories towards establishing an efficient unitarity based amplitude construction. In doing so, we find evidence that the leading IR divergence of NLSM amplitudes exponentiates when the target space is$$ \mathbbm{CP} $$ CP 1 ≅ SU(2)/U(1). We then systematically compute the two-loop anomalous behavior of Born-Infeld, and find that the counterterms needed to restore U(1) invariant behavior at loop-level can be constructed via a symmetric-structure double-copy. We also demonstrate that the divergent part of the one-minus (−+++) two-loop anomaly vanishes upon introducing an evanescent operator. In addition to these purely photonic counterterms, we verify through explicit calculation that the anomalous matrix elements that violate U(1) duality invariance can be alternatively cancelled by summing over internal$$ \mathcal{N} $$ N = 4 DBIVA superfields. Finally we find that$$ \mathcal{N} $$ N = 4 Dirac-Born-Infeld-Volkov-Akulov (DBIVA) amplitudes admit double-copy construction through two-loop order by reproducing our unitarity based result with a double copy between color-dual$$ \mathcal{N} $$ N = 4 super-Yang-Mills and our two-loop NLSM amplitudes. This result supports the possibility of color-dual representations for NLSM beyond one-loop. We conclude with an overview of howD-dimensional four-photon counterterms can be constructed in generality with the symmetric-structure double-copy, and outline a convenient way of counting evanescent operators using Hilbert series as generating functions.

Physics↗
DOE OSTI2024

Revealing the landscape of globally color-dual multi-loop integrands

Abstract We report on progress in understanding how to construct color-dual multi-loop amplitudes. First we identify a cubic theory, semi-abelian Yang-Mills, that unifies many of the color-dual theories studied in the literature, and provides a prescriptive approach for constructingD-dimensional color-dual numerators through one-loop directly from Feynman rules. By a simple weight counting argument, this approach does not further generalize to two-loops. As a first step in understanding the two-loop challenge, we use aD-dimensional color-dual bootstrap to successfully construct globally color-duallocaltwo-loop four-point nonlinear sigma model (NLSM) numerators. The double-copy of these NLSM numerators with themselves, pure Yang-Mills, and$$\mathcal{N}$$= 4 super-Yang-Mills correctly reproduce the known unitarity constructed integrands of special Galileons, Born-Infeld theory, and Dirac-Born-Infeld-Volkov-Akulov theory, respectively. Applying our bootstrap to two-loop four-point pure Yang-Mills, we exhaustively search the space of local numerators and find that it fails to satisfy global color-kinematics duality, completing a search previously initiated in the literature. We pinpoint the failure to the bowtie unitarity cut, and discuss a path forward towards non-local construction of color-dual integrands at generic loop order.

Physics↗
DOE OSTI2022

On-shell correlators and color-kinematics duality in curved symmetric spacetimes

We define a perturbatively calculable quantity — the on-shell correlator — which furnishes a unified description of particle dynamics in curved spacetime. Specializing to the case of flat and anti-de Sitter space, on-shell correlators coincide precisely with on-shell scattering amplitudes and boundary correlators, respectively. Remarkably, we find that symmetric manifolds admit a generalization of on-shell kinematics in which the corresponding momenta are literally the isometry generators of the spacetime acting on the external kinematic data. These isometric momenta are intrinsically non-commutative but exhibit on-shell conditions that are identical to those of flat space, thus providing a common language for computing and representing on-shell correlators which is agnostic about the underlying geometry. Afterwards, we compute tree-level on-shell correlators for biadjoint scalar (BAS) theory and the nonlinear sigma model (NLSM) and learn that color-kinematics duality is manifested at the level of fields under a mapping of the color algebra to the algebra of gauged isometries on the spacetime manifold. Last but not least, we present a field theoretic derivation of the fundamental BCJ relations for on-shell correlators following from the existence of certain conserved currents in BAS theory and the NLSM.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

NS5-brane backgrounds and coset CFT partition functions

Worldsheet string theory is solvable for a variety of backgrounds involving Neveu-Schwarz fivebranes, in terms of gauged nonlinear sigma models on group manifolds. We compute the worldsheet torus partition function of these models, and propose gauging of null isometries as a unifying principle and conceptual framework for this large family of string backgrounds. In the process, we explain how partition functions of asymmetrically gauged Wess-Zumino-Witten models can be computed from the path integral, and organize and systematize various results scattered throughout the literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

BCJ amplitude relations for Anti-de Sitter boundary correlators in embedding space

We generalize the color/kinematics duality of flat-space scattering amplitudes to the embedding space formulation of AdS boundary correlators. Kinematic numerators and propagators are replaced with differential operators acting on a scalar contact diagram that is the AdS generalization of the momentum conserving delta function of flat space scattering amplitudes. We show that color/kinematics duality implies differential relations among AdS boundary correlators that naturally generalize the flat space BCJ amplitude relations and verify them for the correlators of Yang-Mills theory and of the Nonlinear Sigma Model through four- and six-points, respectively. For the latter we also find representations of the four- and six-point correlator that manifest the duality. Possible double-copy procedures in AdS space are also discussed.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗