SEARCH · Search NASA

Results for “Lagrangian Particles”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 433 records · Page 24

All sources NASA NTRS DOE OSTI

NASA NTRS1992

A two-phase restricted equilibrium model for combustion of metalized solid propellants

An Eulerian-Lagrangian two-phase approach was adopted to model the multi-phase reacting internal flow in a solid rocket with a metalized propellant. An Eulerian description was used to analyze the motion of the continuous phase which includes the gas as well as the small (micron-sized) particulates, while a Lagrangian description is used for the analysis of the discrete phase which consists of the larger particulates in the motor chamber. The particulates consist of Al and Al2O3 such that the particulate composition is 100 percent Al at injection from the propellant surface with Al2O3 fraction increasing due to combustion along the particle trajectory. An empirical model is used to compute the combustion rate for agglomerates while the continuous phase chemistry is treated using chemical equilibrium. The computer code was used to simulate the reacting flow in a solid rocket motor with an AP/HTPB/Al propellant. The computed results show the existence of an extended combustion zone in the chamber rather than a thin reaction region. The presence of the extended combustion zone results in the chamber flow field and chemical being far from isothermal (as would be predicted by a surface combustion assumption). The temperature in the chamber increases from about 2600 K at the propellant surface to about 3350 K in the core. Similarly the chemical composition and the density of the propellant gas also show spatially non-uniform distribution in the chamber. The analysis developed under the present effort provides a more sophisticated tool for solid rocket internal flow predictions than is presently available, and can be useful in studying apparent anomalies and improving the simple correlations currently in use. The code can be used in the analysis of combustion efficiency, thermal load in the internal insulation, plume radiation, etc.

Sabnis, J. S.↗

DOE OSTI2020

The geometric theory of charge conservation in particle-in-cell simulations

In recent years, several gauge-symmetric particle-in-cell (PIC) methods have been developed whose simulations of particles and electromagnetic fields exactly conserve charge. While it is rightly observed that these methods’ gauge symmetry gives rise to their charge conservation, this causal relationship has generally been asserted via ad hoc derivations of the associated conservation laws. In this work, we develop a comprehensive theoretical grounding for charge conservation in gauge-symmetric Lagrangian and Hamiltonian PIC algorithms. For Lagrangian variational PIC methods, we apply Noether’s second theorem to demonstrate that gauge symmetry gives rise to a local charge conservation law as an off-shell identity. For Hamiltonian splitting methods, we show that the momentum map establishes their charge conservation laws. We define a new class of algorithms – gauge-compatible splitting methods – that exactly preserve the momentum map associated with a Hamiltonian system’s gauge symmetry – even after time discretization. This class of algorithms affords splitting schemes a decided advantage over alternative Hamiltonian integrators. We apply this general technique to design a novel, explicit, symplectic, gauge-compatible splitting PIC method, whose momentum map yields an exact local charge conservation law. Finally, our study clarifies the appropriate initial conditions for such schemes and examines their symplectic reduction.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Analysis of Σ* via isospin selective reaction 𝐾 𝐿 ⁢𝑝 →𝜋 + ⁢Σ 0

The isospin-selective reaction 𝐾 𝐿 ⁢𝑝 →𝜋 + ⁢Σ 0 provides a clean probe for investigating 𝐼 =1 Σ* resonances. In this work, we perform an analysis of this reaction using an effective Lagrangian approach for the first time, incorporating the well-established Σ⁡(1189)⁢1/2 + , Σ⁡(1385)⁢3/2 + , Σ⁡(1670)⁢3/2 − , and Σ⁡(1775)⁢5/2 − states, while also exploring contributions from other unestablished states. By fitting the available differential cross-section and recoil polarization data, adhering to partial-wave phase conventions same as the Particle Data Group, we find that besides the established resonances, contributions from Σ⁡(1660)⁢1/2 + , Σ⁡(1580)⁢3/2 − , and a Σ*⁢(1/2 − ) improve the description. Notably, a Σ*⁢(1/2 − ) resonance with mass around 1.54 GeV, consistent with Σ⁡(1620)⁢1/2 − , is found to be essential for describing the data in this channel, a stronger indication than found in previous analyses focusing on 𝜋⁢Λ final states. While providing complementary support for Σ⁡(1660)⁢1/2 + and Σ⁡(1580)⁢3/2 − , our results highlight the importance of the Σ⁡(1620)⁢1/2 − region in 𝐾 𝐿 ⁢𝑝 →𝜋 + ⁢Σ 0 . Future high-precision measurements are needed to solidify these findings and further constrain the Σ* spectrum.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

NASA NTRS2002

Transient Mixing Driven by Buoyancy Flows

Mixing driven by buoyancy-induced flows is of particular interest to microgravity processes, as the body force that governs the intensity of flow fields can be directly controlled. We consider a model experimental system to explore the dynamics of mixing which employs two miscible liquids inside a cavity separated initially by a divider. The two liquids are oriented vertically inside a rectangular cavity with constant width and height, and varying depths to span the range of a Hele-Shaw cell to a 3-[) configuration. The two miscible liquids can be sufficiently diluted and died, for example water and deuterium oxide, such that a distinct interface exists across the divider. The transient mixing characteristic of the two fluids is addressed by following the Lagrangian history of the interface for various aspect ratios in the z-plane (depth variation) as well as a range of pulling velocities of the divider. The mixing characteristics of the two fluids are quantified from measurement of the length stretch of the interface and its flow field using respectively image processing techniques and Particle Imaging Velocimetry. Scaling analysis shows that the length stretch depends on four governing parameters, namely the Grashof number (Gr), Schmidt number (Sc), aspect ratio (Ar), and Reynolds number (Re). Variation of the Schmidt number is taken into account through thermophysical property variation. Thus our problem reduces to a codimension three bifurcation in parametric space for Gr, Ar, and Re. Our experimental results show that for Gr on the order of 106 and a nominal cavity aspect ratio Ar = 0.2, the net effect of removal of the divider and the overwhelming buoyancy force causes an overturning motion which stretches and folds the interface to produce ml internal breakwave. The structure of the breakwave is similar to the ubiquitous Rayleigh-Taylor instability morphology. The breakwave is dissipated either through internal or wall collision depending on the impulsive velocity of the divider as prescribed by the Reynolds number. The decay of the collision event occurs through sloshing oscillations over a short time scale. The two fluids then become stably stratified with a diffusive band at the interface indicating local mass transport. The local bifurcation of the internal breakwave is investigated as a function of aspect ratio. Results show that for narrow cavities on the order of 2mm (Ar = 0.04) folding does not occur, the interface only stretches. As the cavity size increases folding occurs through a supercritical bifurcation. Insight into the mechanism of folding is obtained from measurement of the flow field which shows that in the neighborhood of the folding event, there exists hyperbolic points caused by multiple vortex interactions. The global length stretch of the interface as a function of time is nearly Gaussian; calculations of finite-time Liapunov exponents as well as construction of horseshoe maps indicate the likelihood of a chaotic transient.

W.M.B. Duval↗

NASA NTRS2023

Assessment of Oceanographic Conditions during the North Atlantic EXport Processes in the Ocean from RemoTe Sensing (EXPORTS) Field Campaign

This manuscript presents an overview of NASA’s EXport Processes in the Ocean from Remote Sensing 2021 Field Campaign in the North Atlantic (EXPORTS NA) and provides quantitative and dynamical descriptions of the physical processes modulating water transformations during the study. A major programmatic goal was to conduct the sampling in a Lagrangian mode so that ocean ecological and biogeochemical changes can be observed independent from physical advective processes. To accomplish this goal, EXPORTS NA conducted a multi-ship, multi-asset field sampling program within a retentive, anticyclonic mode water eddy. Beneath depths of ∼100 m, Lagrangian sampling assets remained within the eddy core waters (ECWs) throughout the experiment, demonstrating that the ECWs within the mode water eddy were retentive. However, strong westerly winds from four storm events deepened the mixed layer (ML) of the surface core waters (SCWs) above the eddy’s mode water core by 25-40 m and exchanged some of the SCWs with surface waters outside of the eddy via Ekman transport. Estimates of flushing times ranged from 5-8 days, with surface exchange fractions ranging from 20-75% and were consistent with particle tracking advected by combined geostrophic and Ekman velocities. The relative contributions of horizontal and vertical advection on changes in ECW tracers depended on the horizontal and vertical gradients of that tracer. For example, horizontal advection played a large role in ECW salinity fluxes, while vertical entrainment played a larger role in the fluxes of nutrients into SCW ML. Each storm injected nutrients and low oxygen waters into the ML, after which the surface ocean ecosystem responded by reducing nutrient concentrations and increasing %O2 saturation levels. Overall, ECW values of chlorophyll and POC were the largest at the onset of the field program and decreased throughout the campaign. The analysis presented provides a physical oceanographic context for the many measurements made during the EXPORTS NA field campaign while illustrating the many challenges of conducting a production-flux experiment, even in a Lagrangian frame, and the inherent uncertainties of interpreting biological carbon pump observations that were collected in a Eulerian frame of reference.

oceanography↗

DOE OSTI2024

Hadamard products and BPS networks

We study examples of fourth-order Picard-Fuchs operators that are Hadamard products of two second-order Picard-Fuchs operators. Each second-order Picard-Fuchs operator is associated with a family of elliptic curves, and the Hadamard product computes period integrals on the fibred product of the two elliptic surfaces. We construct 3-cycles on this geometry as the union of 2-cycles in the fibre over contours on the base. We then use the special Lagrangian condition to constrain the contours on the base. This leads to a construction that is reminiscent of spectral networks and exponential networks that have previously appeared in string theory literature.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1994

Mixing by barotropic instability in a nonlinear model

A global, nonlinear, equivalent barotropic model is used to study the isentropic mixing of passive tracers by barotropic instability. Basic states are analytical zonal-mean jets representative of the zonal-mean flow in the upper stratosphere, where the observed 4-day wave is thought to be a result of barotropic, and possibly baroclinic, instability. As is known from previous studies, the phase speed and growth rate of the unstable waves is fairly sensitive to the shape of the zonal-mean jet; and the dominant wave mode at saturation is not necessarily the fastest growing mode; but the unstable modes share many features of the observed 4-day wave. Lagrangian trajectories computed from model winds are used to characterize the mixing by the flow. For profiles with both midlatitude and polar modes, mixing is stronger in midlatitude than inside the vortex; but there is little exchange of air across the vortex boundary. There is a minimum in the Lyapunov exponents of the flow and the particle dispersion at the jet maximum. For profiles with only polar unstable modes, there is weak mixing inside the vortex, no mixing outside the vortex, and no exchange of air across the vortex boundary. These results support the theoretical arguments that, whether wave disturbances are generated by local instability or propagate from other regions, the mixing properties of the total flow are determined by the locations of the wave critical lines and that strong gradients of potential vorticity are very resistant to mixing.

Bowman, Kenneth P.↗

DOE OSTI2021

Fractons and Exotic Symmetries from Branes

The emerging study of fractons, a new type of quasi-particle with restricted mobility, has motivated the construction of several classes of interesting continuum quantum field theories with novel properties. One such class consists of foliated field theories which, roughly, are built by coupling together fields supported on the leaves of foliations of spacetime. Another approach, which we refer to as exotic field theory, focuses on constructing Lagrangians consistent with special symmetries (like subsystem symmetries) that are adjacent to fracton physics. A third framework is that of infinite-component Chern-Simons theories, which attempts to generalize the role of conventional Chern-Simons theory in describing (2+1)D Abelian topological order to fractonic order in (3+1)D. The study of these theories is ongoing, and many of their properties remain to be understood. Historically, it has been fruitful to study QFTs by embedding them into string theory. One way this can be done is via D-branes, extended objects whose dynamics can, at low energies, be described in terms of conventional quantum field theory. Further, QFTs that can be realized in this way can then be analyzed using the rich mathematical and physical structure of string theory. In this paper, we show that foliated field theories, exotic field theories, and infinite-component Chern-Simons theories can all be realized on the world-volumes of branes. We hope that these constructions will ultimately yield valuable insights into the physics of these interesting field theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

NASA NTRS1988

Annual Research Briefs, 1987

Lagrangian techniques have found widespread application to the prediction and understanding of turbulent transport phenomena and have yielded satisfactory results for different cases of shear flow problems. However, it must be kept in mind that in most experiments what is really available are Eulerian statistics, and it is far from obvious how to extract from them the information relevant to the Lagrangian behavior of the flow; in consequence, Lagrangian models still include some hypothesis for which no adequate supporting evidence was until now available. Direct numerical simulation of turbulence offers a new way to obtain Lagrangian statistics and so verify the validity of the current predictive models and the accuracy of their results. After the pioneering work of Riley (Riley and Patterson, 1974) in the 70's, some such results have just appeared in the literature (Lee et al, Yeung and Pope). The present contribution follows in part similar lines, but focuses on two particle statistics and comparison with existing models.

Moin, Parviz↗

NASA NTRS2024

Modeling Soft X‐Ray Emissions at the Dayside Magnetopause

In this study, we simulate the Solar Wind Charge Exchange (SWCX) soft X-ray emissions at dayside magnetosheath and cusps by using magnetohydrodynamic (MHD) and LAtmos TEst Particle (LaTeP) models. MHD models are unable to resolve the particle kinetic effects, such as the different behaviors of ions with different q/m, or distinguish the magnetospheric plasma from the solar wind plasma. We investigate these effects with the LaTeP model. As the LaTeP model does not self-compute magnetic and electric field, the magnetic and electric field data obtained from Open Geospace General Circulation Model (OpenGGCM) and Lagrangian version of the piecewise parabolic method (PPMLR) MHD model are used as the input to LaTeP model. The soft X-ray emissivity maps simulated from pure OpenGGCM and PPMLR MHD approaches and from LaTeP-OpenGGCM and LaTeP-PPMLR approaches are presented and compared. The results indicate that the LaTeP model can well resolve the kinetic effects and can be used to investigate the individual spectral characteristics. Therefore, the LaTeP model is a complementary approach for simulating the X-ray emissions near the dayside magnetopause. We also calculate the ratio of integrated OVII/OVIII line intensities, produced by charge exchange of O7+ ions and O8+ ions, respectively. We find a relatively higher ratio at the bow shock compared to the surrounding areas, suggesting that this ratio can be an effective parameter to identify the bow shock location.

Qiuyu Xu↗

DOE OSTI2021

Energy Exascale Earth System Model v2.0

First release of version 2 of the Energy Exascale Earth System Model. The atmosphere component remains EAM. Major changes since version 1 include: all column-physics parameterizations are computed on a separate grid that has approximately half the number of points of the dynamics grid, a new nonhydrostatic dynamical core (running in hydrostatic mode) with semi-Lagrangian tracer transport, CLUBB updated from v1 to v2, a new convective trigger (dCAPE/ULL) based on the dynamic Convective Available Potential Energy (CAPE) (dCAPE) and the Unrestricted Launch Level (ULL) concepts is used in ZM. minimum cloud droplet number changed, gravity wave drag energy conservation fixed and new tunings used, dust emission size distribution changed to emit more coarse dust particles The land component is still ELM. Major changes since version 1 include: using SNICAR-AD for radiation in snow to match the sea-ice model and fixing bugs in snow compaction and water state calculation. The ocean component remains MPAS-ocean. Major change since version 1 include: Redi isopycnal mixing has been updated, tested, and tuned in combination with the Gent-McWilliams parameterization, a sign error was fixed in the 3rd-order flux routines, the mesh used in low-resolution coupled cases was modified and the time steps adjusted, new regionally refined meshes were created, one focused on North America and another on the Southern Ocean, including ice shelf cavities. The sea-ice component remains MPAS-seaice. Major changes since version 1 include: A new heat- and freshwater-conserving coupling of frazil ice, turning off of SSH filtering, addition of SNICAR-AD and snow grain aging. The land-ice component remains MPAS-Albany-landIce (MALI) and is a static ice sheet. There is more out-of-the-box support for cryosphere configurations including the new regionally refined configuration around Antarctica with ice shelf cavities. The river model is MOSART. The half-degree river mesh was redone so it no longer treats Black and Caspian seas as ocean. The coupler remains cpl7/MCT. Major changes since version 1 include handling of ice shelf melt fluxes (heat / freshwater exchange with the ocean), and data icebergs. All components allow regional refinement of their meshes and two separate example of refinement, one in and around North America and one around Antarctica and the Southern Ocean, are provided.

E3SM Project, DOE↗

DOE OSTI2023

Mixing in Low Reynolds Number Reacting Impinging Jets in Crossflow

Previous efforts to model uranyl fluoride formation in an impinging jet gas reactor underpredicted spatial mixing and overpredicted chemical conversion into particulates. The previous fluid dynamics model was based on the solution of the Reynolds Averaged Navier Stokes equations. After simulating fluid dynamics, aerosol dynamics were superimposed onto CFD-simulated gas reactant species concentrations. The current work explores the influence of complex unsteady flow features on the overall flow physics and chemistry for a low Reynolds number, opposed flow, impinging jet gas reactor where there is a low Reynolds number cross flow. The objective of this study was to assess the impact of model formulation on scalar mixing and transport. Here, transient flow simulations were performed using Scale Resolving Simulations. Large-Eddy Simulations with the dynamic Smagorinsky turbulence model were performed along with simulations which directly resolved the flow. Average and root-mean-square (RMS) velocities and species concentrations were computed along with modeled and resolved turbulence kinetic energy (TKE), modeled turbulence dissipation, and modeled turbulent viscosity. Lagrangian flow tracers were also used to quantify species concentrations along path lines emanating from the jet tips. Transient simulation data were compared to results from RANS simulations using the k-ω shear stress transport (SST) model and Reynolds Stress Model (RSM). Transient simulations showed spatial mixing patterns which were more consistent with experimental data and helped elucidate the process of particle formation observed in experiments.

42 ENGINEERING↗

DOE OSTI2022

Pore-Scale Modeling of Mineral Growth and Nucleation in Reactive Flow

A fundamental understanding of mineral precipitation kinetics relies largely on microscopic observations of the dynamics of mineral surfaces exposed to supersaturated solutions. Deconvolution of tightly bound transport, surface reaction, and crystal nucleation phenomena still remains one of the main challenges. Particularly, the influence of these processes on texture and morphology of mineral precipitate remains unclear. This study presents a coupling of pore-scale reactive transport modeling with the Arbitrary Lagrangian-Eulerian approach for tracking evolution of explicit solid interface during mineral precipitation. It incorporates a heterogeneous nucleation mechanism according to Classical Nucleation Theory which can be turned “on” or “off.” This approach allows us to demonstrate the role of nucleation on precipitate texture with a focus at micrometer scale. In this work precipitate formation is modeled on a 10 micrometer radius particle in reactive flow. The evolution of explicit interface accounts for the surface curvature which is crucial at this scale in the regime of emerging instabilities. The results illustrate how the surface reaction and reactive fluid flow affect the shape of precipitate on a solid particle. It is shown that nucleation promotes the formation of irregularly shaped precipitate and diminishes the effect of the flow on the asymmetry of precipitation around the particle. The observed differences in precipitate structure are expected to be an important benchmark for reaction-driven precipitation in natural environments.

58 GEOSCIENCES↗

DOE OSTI2020

Markovian transport processes in a heterogeneous, variably saturated watershed: A multi-domain spatial Markov model

Transport processes in watersheds remain challenging to describe accurately with upscaled models because of their intricate complexities at multiple scales that can lead to multi-modal, non-Fickian breakthrough curves. One of the recent advances in solute transport modeling has been the Lagrangian spatial Markov model (SMM), which describes transport using velocity increments and models the correlated transitions between them. The method has been applied previously to saturated porous media and fractured media, but these successful applications suggest it may be useful for describing watershed scale processes as well. An existing, 3-D, heterogeneous variably saturated model of a headwaters catchment and classical random walk particle tracking were used to generate a forward simulation and the resulting transport was analyzed with a spatial Markov framework. Existing SMM methods were used to analyze the transport paths through the complex, unconfined domain and we found that the vadose zone plays an important role in the correlation structure. As such, a new kind of SMM is proposed, termed a multi-domain SMM (MD-SMM), that uses separate correlation models for auto-transitions in the vadose zone and the saturated zone, and for the cross-transitions between them. The MD-SMM improves representation of transport relative to the single domain model and can naturally delineate times when a particle is in the vadose zone or the saturated zone, adding an additional degree of realism to the upscaled model. Lastly, we show that the transition probability matrix for transport in the saturated portion of this unconfined system is approximately Markovian, which is an important validation of the proposed framework.

54 ENVIRONMENTAL SCIENCES↗

DOE OSTI2023

Neutrino Properties and Interactions

In this chapter, we present basic concepts of neutrino physics. We start with a brief introduction to the standard model electroweak sector, followed by calculations of some relevant neutrino interaction cross sections. Here, we obtain the oscillation formalism from the Dirac equation, alongside with a self-contained derivation of matter effects from the standard model electroweak Lagrangian. Some key features of neutrino oscillation phenomenology are discussed as well. We then review the most precise oscillation measurements and finalize with some broad discussion of the current open questions in neutrino physics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

DOE OSTI2024

Gauge invariance of radiative jet functions in the position-space formulation of SCET

In subleading powers of soft-collinear effective theory (SCET), the Lagrangian contains couplings between soft quarks and hard-collinear quarks. Matrix elements of the hard-collinear parts of these couplings are radiative jet functions. In the position-space formulation of SCET, the Lagrangians are constructed from operators that appear to be gauge invariant. Nevertheless, we find violations of gauge invariance arise in the hard-collinear sector because gauge transformations can shift the momentum of a hard-collinear quark field from the hard-collinear sector to the soft sector, where the hard-collinear fields, by definition, have no support. The violations of gauge invariance are manifested in perturbation theory in the hard-collinear sector through the absence of certain Feynman diagrams that would be present in full QCD. A consequence of the absence of these diagrams is that the radiative jet functions that follow directly from the position-space Lagrangians are not gauge invariant, and we demonstrate this through explicit calculations in lower-order perturbation theory. We obtain gauge-invariant Lagrangians by adding to existing position-space Lagrangians terms that are proportional to the soft-quark equation of motion. These gauge-invariant Lagrangians are valid for nonzero, as well as zero, quark masses. We also remark briefly on the gauge invariance of certain Lagrangians that have been constructed in the label-momentum formulation of SCET. Published by the American Physical Society 2024

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

NASA NTRS1993

Streaming vorticity flux from oscillating walls with finite amplitude

How to describe vorticity creation from a moving wall is a long standing problem. This paper discusses relevant issues at the fundamental level. First, it is shown that the concept of 'vorticity flux due to wall acceleration' can be best understood by following fluid particles on the wall rather than observing the flow at fixed spatial points. This is of crucial importance when the time-averaged flux is to be considered. The averaged flux has to be estimated in a wall-fixed frame of reference (in which there is no flux due to wall acceleration at all); or, if an inertial frame of reference is used, the generalized Lagrangian mean (GLM) also gives the same result. Then, for some simple but typical configurations, the time-averaged vorticity flux from a harmonically oscillating wall with finite amplitude is analyzed, without appealing to small perturbation. The main conclusion is that the wall oscillation will produce an additional mean vorticity flux (a fully nonlinear streaming effect), which is partially responsible for the mechanism of vortex flow control by waves. The results provide qualitative explanation for some experimentally and/or computationally observed phenomena.

Wu, J. Z.↗

DOE OSTI2026

LeWRON: Agentic Analysis of Electroweak Phase Transitions

The electroweak phase transition (EWPT) is a central topic in particle physics and cosmology, connecting collider phenomenology, baryogenesis, and gravitational-wave observatories. Its analysis requires a technically demanding, convention-sensitive, and model-dependent pipeline, from constructing the finite-temperature effective potential to tracking thermal histories, computing bubble nucleation rates, and predicting gravitational-wave spectra. We present LeWRON (Learning ElectroWeak phase tRansitiON), an agentic framework that orchestrates this pipeline starting from an input Lagrangian. LeWRON combines audited toolbox construction with an Explorer module that uses the generated model-specific code for further analysis, including scans and plots. Intermediate analytic outputs are checked by auditor agents and stored as structured artifacts, enabling reproducible human inspection and downstream use through both a command-line interface and a public Python API. The framework supports a reproduction mode, which infers conventions from the literature and reproduces published results, and a discovery mode, which guides users through structured checkpoints for new models. We demonstrate LeWRON across representative beyond-the-Standard-Model scenarios and release the code on GitHub.

Wang, Isaac R. [Fermilab] (ORCID:000000030789218X)↗