Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Lagrangian approach”

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 307 records · Page 17
All sourcesNASA NTRSDOE OSTI
NASA NTRSDate not supplied

Modeling of Cavity Residence Time via Modeled Lagrangian Particle Tracking

One of the critical components in a scramjet is the cavity flame holder. In evaluating the sizing of this component, the residence time of the fuel in the cavity is a common parameter used. However, to obtain the residence time via CFD is an expensive undertaking. One method is to perform an LES simulation of the cavity and use Lagrangian particle tracking to obtain both the mean and probability density function of the time that the particles remain in the cavity. While this method affords a significant amount of information, it is very costly in CPU run time. A simpler method involves performing a time accurate RANS simulation of the cavity, and when converged, overlay a traceable scalar in the cavity region, and then record the decay of this scalar’s flux at a downstream location. The time scale of the decay rate can then be estimated and a residence time inferred. This method is less computationally expensive than the LES approach, but still requires a time-accurate solution while only providing a mean cavity residence time evaluation. In this paper a method of determining the mean cavity residence time as well as the Probability Density Function (PDF) of the residence time is proposed and demonstrated. This method is based on Lagrangian particle tracking modeled by the Generalized Langevin Equation and is implemented as a post-processing step for a steady RANS solution. Computational cost is on the order of minutes and the results are in agreement with values obtained by the LES and scalar decay methods.

Langevin↗
DOE OSTI2020

Photon scattering by a 4π -spherically-focused ultrastrong electromagnetic wave

The scattering of high-power probe laser pulses by tightly focused ultrastrong laser pulses has been investigated through a semiclassical approach using a perturbative method based on the Born approximation. Under a 4π-spherically-focused ultrastrong light field, the electric permittivity and magnetic permeability tensors for vacuum are calculated from the Euler-Heisenberg Lagrangian to show the nonlinear birefringent property of vacuum. And, from permittivity and permeability tensors, the scattering potential is derived for the Born approximation. The first-order solution of the Born approximation is taken as an electric field scattered from the 4π-spherically-focused laser pulse when a probe laser pulse propagates through the focused laser field. The differential cross section of the nonlinear birefringent vacuum is derived, and the number of photons scattered from the nonlinear birefringent vacuum is analyzed in the laser power range of 10-1000 PW.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2024

An asymptotic-preserving semi-Lagrangian algorithm for the anisotropic heat transport equation with arbitrary magnetic fields

Here, we extend the recently proposed semi-Lagrangian algorithm for the extremely anisotropic heat transport equation [Chacón et al., J. Comput. Phys ., 272 (2014)] to deal with arbitrary magnetic field topologies. The original scheme (which showed remarkable numerical properties) was valid for the so-called tokamak-ordering regime, in which the magnetic field magnitude was not allowed to vary much along field lines. The proposed extension maintains the attractive features of the original scheme (including the analytical Green's function, which is critical for tractability) with minor modifications, while allowing for completely general magnetic fields. The accuracy and generality of the approach are demonstrated by numerical experiment with an analytical manufactured solution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2020

A Scalable and Distributed Algorithm for Managing Residential Demand Response Programs using Alternating Direction Method of Multipliers (ADMM)

For effective engagement of residential demand-side resources and to ensure efficient operation of distribution networks, we must overcome the challenges of controlling and coordinating residential components and devices at scale. To overcome this challenge, we present a distributed and scalable algorithm with a three-level hierarchical information exchange architecture for managing the residential demand response programs. First, a centralized optimization model is formulated to maximize community social welfare. Then, this centralized model is solved in a distributed manner with alternating direction method of multipliers (ADMM) by decomposing the original problem to utility-level and house-level problems. The information exchange between the different layers is limited to the primary residual (i.e., supply-demand mismatch), Lagrangian multipliers, and the total load of each house to protect each customer’s privacy. Simulation studies are performed on the IEEE 33 bus test system with 605 residential customers. The results demonstrate that the proposed approach can save customers’ electricity bills and reduce the peak load at the utility level without much affecting customers’ comfort and privacy. Finally, a quantitative comparison of the distributed and centralized algorithms shows the scalability advantage of the proposed ADMM-based approach, and it gives benchmarking results with achievable value for future research works.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2024

Impacts of stochastic coalescence variability on warm rain initiation using Lagrangian microphysics in box and large-eddy simulations

Various coalescence methods for Lagrangian microphysics schemes are tested in box and large-eddy simulation (LES) models, including the stochastic all-or-nothing super-droplet method (SDM) and a version of SDM (dSDM) that applies a fractional approach similar to the average impact method. In LES, variabilities driven by microphysics and by flow realizations are separated using the “piggybacking” technique. Rain initiation averaged over many realizations of the box model is delayed and rain variability increases as the number of super-drops per collision volume ( N SD ) is decreased using SDM. In contrast, rain initiation time using SDM in LES is insensitive to N SD for 32 ≤ N SD ≤ 512. This is explained through the interaction between LES grid boxes, each acting as a separate collision volume. Variability across the ensemble of LES collision volumes using SDM results in rain quickly initiating in some of the LES grid cells at low N SD and leading to a similar overall timing of rain initiation from the cloud compared to simulations with high N SD . There is a ~20% decrease in the total rain mass and mean rain flux as N SD is increased from 32 to 256, with little additional change as N SD is increased from 256 to 512. The fractional coalescence approach in dSDM leads to reduced microphysical variability and a 15-18 min delay in rain initiation compared to SDM. Finally, an additional LES ensemble with microphysical variability feeding back to the dynamics shows that flow variability dominates the impact of microphysical variability on rain properties. Thus, flow variability must be constrained to isolate impacts of microphysical variability.

54 ENVIRONMENTAL SCIENCES↗
NASA NTRS1976

Structural dynamics, stability, and control of helicopters

The dynamic synthesis of a helicopter is reported. The method of approach is a variation of the component mode synthesis in the sense that it regards the aircraft as an assemblage of interconnected substructures. The equations of motion are derived in general form by means of the Lagrangian formulation in conjunction with an orderly kinematical procedure that takes into account the superposition of motion of various substructures, thus circumventing constraint problems.

Meirovitch, L.↗
NASA NTRS1980

Solutions of contact problems by the assumed stress hybrid model

A method was developed for contact problems which may be either frictional or frictionless and may involve extensive sliding between deformable bodies. It was based on an assumed stress hybrid approach and on an incremental variational principle for which the Euler's equations of the functional include the equilibrium and compatibility conditions at the contact surface. The tractions at an assumed contact surface were introduced as Lagrangian multipliers in the formulation. It was concluded from the results of several example solutions that the extensive sliding contact between deformable bodies can be solved by the present method.

Kubomura, K.↗
NASA NTRS1998

Current Status on the use of Parallel Computing in Turbulent Reacting Flow Computations Involving Sprays, Monte Carlo PDF and Unstructured Grids

The state of the art in multidimensional combustor modeling as evidenced by the level of sophistication employed in terms of modeling and numerical accuracy considerations, is also dictated by the available computer memory and turnaround times afforded by present-day computers. With the aim of advancing the current multi-dimensional computational tools used in the design of advanced technology combustors, a solution procedure is developed that combines the novelty of the coupled CFD/spray/scalar Monte Carlo PDF (Probability Density Function) computations on unstructured grids with the ability to run on parallel architectures. In this approach, the mean gas-phase velocity and turbulence fields are determined from a standard turbulence model, the joint composition of species and enthalpy from the solution of a modeled PDF transport equation, and a Lagrangian-based dilute spray model is used for the liquid-phase representation. The gas-turbine combustor flows are often characterized by a complex interaction between various physical processes associated with the interaction between the liquid and gas phases, droplet vaporization, turbulent mixing, heat release associated with chemical kinetics, radiative heat transfer associated with highly absorbing and radiating species, among others. The rate controlling processes often interact with each other at various disparate time 1 and length scales. In particular, turbulence plays an important role in determining the rates of mass and heat transfer, chemical reactions, and liquid phase evaporation in many practical combustion devices.

Raju, M. S.↗
NASA NTRS1991

Efficient high-resolution rotor wake calculations using flow field reconstruction

Recent computational studies of helicopter flowfields for acoustic applications has led to the development of exceptionally efficient analysis methods for carrying out Lagrangian simulations of the vortex wake. The capabilities of these methods are demonstrated here through calculations involving a variety of wake/rotor interactions associated with both main rotor and tail rotor flowfields. These methods incorporate a flowfield reconstruction approach efficiently to compute vortex wake interactions within arbitrary regions located near the main rotor. This paper summarizes the development of the reconstruction approach and outlines some of its potential uses, while focusing primarily on its application to rotor noise calculations. The development of the nearfield velocity corrections used in the analysis of close vortex interactions is described, and the overall structure of the current reconstruction approach is discussed.

Quackenbush, Todd R.↗
DOE OSTI2023

Ocean & Geohazard Analysis Tool

The Ocean & Geohazard Analysis (OGA) software tool is designed to summarize insights into key offshore hazards drawing from a diverse set of approaches, including artificial intelligence, machine learning, probabilistic and statistical, and offshore data sources. The offshore hazards that can be analyzed include submarine landslides, extreme wind/wave/current event probabilities, earthquakes, and metocean pathways (CIAM Climatological Isolation and Attraction Model–Climatological Lagrangian Coherent Structures - Submissions - EDX (doe.gov)). Currently, the tool is developed for use in the Gulf of Mexico. The data underlying the offshore hazard analyses can be found here: https://edx.netl.doe.gov/dataset/gulf-of-mexico-risk-analysis-database-gomrad This work was conducted under the Advanced Offshore Research Portfolio, FWP Number 1022409 at National Energy Technology Laboratory, U.S. Dept. of Energy. Disclaimer This project was funded by the United States Department of Energy, National Energy Technology Laboratory, in part, through a site support contract. Neither the United States Government nor any agency thereof, nor any of their employees, nor the support contractor, nor any of their employees, makes any warranty, express or implied, or assumes any legal liability or responsibility for the accuracy, completeness, or usefulness of any information, apparatus, product, or process disclosed, or represents that its use would not infringe privately owned rights. Reference herein to any specific commercial product, process, or service by trade name, trademark, manufacturer, or otherwise does not necessarily constitute or imply its endorsement, recommendation, or favoring by the United States Government or any agency thereof. The views and opinions of authors expressed herein do not necessarily state or reflect those of the United States Government or any agency thereof.

AIML↗
NASA NTRS1988

Predicting the velocity and azimuth of fragments generated by the range destruction or random failure of rocket casings and tankage

A calculational method is described which provides a powerful tool for predicting solid rocket motor (SRM) casing and liquid rocket tankage fragmentation response. The approach properly partitions the available impulse to each major system-mass component. It uses the Pisces code developed by Physics International to couple the forces generated by an Eulerian-modeled gas flow field to a Lagrangian-modeled fuel and casing system. The details of the predictive analytical modeling process and the development of normalized relations for momentum partition as a function of SRM burn time and initial geometry are discussed. Methods for applying similar modeling techniques to liquid-tankage-overpressure failures are also discussed. Good agreement between predictions and observations are obtained for five specific events.

Eck, Marshall↗
DOE OSTI2021

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↗
NASA NTRS1985

Three-dimensional hybrid-stress finite element analysis of composite laminates with cracks and cutouts

A three-dimensional hybrid-stress finite element analysis of composite laminates containing cutouts and cracks is presented. Fully three-dimensional, hexahedral isoparametric elements of the hybrid-stress model are formulated on the basis of the Hellinger-Reissner variational principle. Traction-free edges, cutouts, and crack surfaces are modeled by imposition of exact traction boundary conditions along element surfaces. Special boundary and surface elements are constructed by introducing proper constraints on assumed stress functions. The Lagrangian multiplier technique is used to enforce ply-interface continuity conditions in hybrid bimaterial composite elements for modeling the interface region in a composite laminate. Two examples are given to illustrate the capability of the present method of approach: (1) the well-known delamination problem in an angle-ply laminate, and (2) the important problem of a composite laminate containing a circular hole. Results are presented in detail for each case. Implications of interlaminar and intralaminar crack initiation, growth and fracture in composites containing cracks and cutouts are discussed.

Wang, S. S.↗
NASA NTRS1990

The pdf approach to turbulent flow

This paper provides a detailed discussion of the theory and application of probability density function (pdf) methods, which provide a complete statistical description of turbulent flow fields at a single point or a finite number of points. The basic laws governing the flow of Newtonian fluids are set up in the Eulerian and the Lagrangian frame, and the exact and linear equations for the characteristic functionals in those frames are discussed. Pdf equations in both frames are derived as Fourier transforms of the equations of the characteristic functions. Possible formulations for the nonclosed terms in the pdf equation are discussed, their properties are assessed, and closure modes for the molecular-transport and the fluctuating pressure-gradient terms are reviewed. The application of pdf methods to turbulent combustion flows, supersonic flows, and the interaction of turbulence with shock waves is discussed.

Kollmann, W.↗
DOE OSTI2024

On small black holes in string theory

Abstract We discuss the worldsheet sigma-model whose target space is thed+1 dimensional Euclidean Schwarzschild black hole. We argue that in the limit where the Hawking temperature of the black hole,T, approaches the Hagedorn temperature,T H , it can be described in terms of a generalized version of the Horowitz-Polchinski effective theory. Ford≥ 6, where the Horowitz-Polchinski EFT [1, 2] does not have suitable solutions, the modified effective Lagrangian allows one to study the black hole CFT in an expansion in powers ofd− 6 andT H −T. AtT=T H , the sigma model is non-trivial for alld> 6. It exhibits an enhanced SU(2) symmetry, and is described by a non-abelian Thirring model with a radially dependent coupling. The resulting picture connects naturally to the results of [3–5], that relate Schwarzschild black holes in flat spacetime at largedto the two dimensional black hole. We also discuss an analogous open string system, in which the black hole is replaced by a system of two separated D-branes connected by a throat. In this system, the asymptotic separation of the branes plays the role of the inverse temperature. At the critical separation, the system is described by a Kondo-type model, which again exhibits an enhanced SU(2) symmetry. At larged, the brane system gives rise to the hairpin brane [6].

Physics↗
DOE OSTI2021

Quasinormal modes for non-minimally coupled scalar fields in regular black hole spacetimes: Grey-body factors, area spectrum and shadow radius

Black hole perturbation theory is a useful approach to study interactions between black holes and fundamental fields. A particular class of black hole solutions arising out of modification of Einstein’s general theory of relativity are regular black holes (RBHs) which can be constructed using a nonlinear electrodynamic Lagrangian. Because of their importance, we are interested in studying the behavior of three kinds of such RBHs under perturbations generated by an external field. Indeed, we investigate the quasinormal modes (QNMs) of a massive scalar field propagating near the RBHs which is non-minimally coupled to the Ricci scalar-tensor of background geometry. We will attempt to find the low-lying quasinormal frequencies of the perturbations by using WKB approximation. We shall also study the relationship between the QNMs and some characteristic properties of black holes such as grey-body factors, area quantization, and shadow radius.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2024

Distributed optimization for multi-commodity urban traffic control

A distributed method for concurrent traffic signal and routing control of traffic networks is proposed. The method is based on the multi-commodity store-and-forward model, in which the destinations are the commodities. The system benefits from the communication between vehicles and infrastructure, providing optimal signal timings to intersections and routes to vehicles on a link-by-link basis. Using the augmented Lagrangian to model the constraints into the objective, the baseline centralized problem is decomposed into a set of objective-coupled subproblems, one for each intersection, enabling the solution to be computed by a distributed- gradient projection algorithm. Further, the intersection agents only need to communicate and coordinate with neighboring intersections to ensure convergence to the optimal solution while tolerating suboptimal iterations that offer more flexibility, unlike other distributed approaches. Through microsimulation, we demonstrate the effectiveness of the proposed algorithm in traffic networks with time-varying demand. Computational analysis shows that the distributed problem is suitable for real-time applications. A robustness analysis show that the distributed formulation enables a graceful degradation of the system in case of failure.

Augmented Lagrangian↗
DOE OSTI2024

Pion Nucleon Scattering in BCHPT Combined with 1/Nc Expansion

Pion nucleon scattering has played a crucial role in advancing our understanding of the low-energy regime of strong interactions. This dissertation presents the development of a new theoretical framework that combines Baryon Chiral Perturbation Theory (BChPT) with the 1/Nc Expansion to analyze pion nucleon scattering. The resulting effective Lagrangian incorporates both chiral and spin-flavor symmetries. Scattering amplitudes are calculated up to the one-loop level (NNLO), which includes four tree-level Feynman diagrams and 45 non-zero loop diagrams. The effective Lagrangian is renormalized at the one-loop level. With spin-3/2 baryons naturally included, the framework demonstrates strong convergence at higher energies. The NNLO pion nucleon scattering amplitude shows good agreement with experimental data from the SAID database. This novel approach, BChPT combined with 1/Nc Expansion, demonstrates substantial predictive power at higher energies where other low-energy theories fall short, achieving the goals of this combined framework. The contribution from loop diagrams is essential for achieving good agreement with the data, highlighting the importance of conducting the calculation up to NNLO. Extracted values for physically measurable quantities agree reasonably well with independently extracted values, underscoring the robustness of the theory.

Jayakodige, Dulitha [Hampton Univ., Hampton, VA (U↗