Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “constraint energy minimization”

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

Multi-Period Active Distribution Network Planning Using Multi-Stage Stochastic Programming and Nested Decomposition by SDDIP

This paper presents a multi-period active distribution network planning (ADNP) with distributed generation (DG). The objective of the proposed ADNP is to minimize the total planning cost, subject to both investment and operation constraints. The paper proposes a multi-stage stochastic optimization model to address DG uncertainties over several periods, in which the decisions are made sequentially by only using the present-stage information. A nested decomposition method is proposed which applies the stochastic dual dynamic integer programming (SDDIP) method to address computational intractabilities of the proposed ADNP approach. The presented numerical results and discussions on a 33-bus distribution system and a large-scale 906-bus system verify the effectiveness of the proposed ADNP method and its solution method.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

NWTC Site 4.0 - NREL ASSIST (SN10) / Thermodynamic retrievals TROPoe

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe v0.12 (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). The profiles are retrieved every 10 minutes from instantaneous observations from the NREL ASSIST-II (SN 10) infrared spectrometer. Observations are noise-filtered but not averaged in time to minimize errors due to non-uniform clouds. Additional input data in TROPoe are cloud base height from a Vaisala CL51 ceilometer. The full pipeline for running the retrieval is available at https://github.com/StefanoWind/TROPoe_processor. Met data were not ingested. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Denver, CO.

17 WIND ENERGY↗
DOE OSTI2025

NWTC Site 3.2 - NREL ASSIST (SN12) / Thermodynamic retrievals TROPoe

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe v0.12 (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). The profiles are retrieved every 10 minutes from instantaneous observations from the NREL ASSIST-II (SN 12) infrared spectrometer. Observations are noise-filtered but not averaged in time to minimize errors due to non-uniform clouds. Additional input data in TROPoe are cloud base height from a Vaisala CL51 ceilometer. The full pipeline for running the retrieval is available at https://github.com/StefanoWind/TROPoe_processor. Met data were not ingested. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Denver, CO.

17 WIND ENERGY↗
DOE OSTI2025

Title NWTC Site 3.2 - NREL ASSIST (SN11) / Thermodynamic retrievals TROPoe

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe v0.12 (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). The profiles are retrieved every 10 minutes from instantaneous observations from the NREL ASSIST-II (SN 11) infrared spectrometer. Observations are noise-filtered but not averaged in time to minimize errors due to non-uniform clouds. Additional input data in TROPoe are cloud base height from a Vaisala CL51 ceilometer. The full pipeline for running the retrieval is available at https://github.com/StefanoWind/TROPoe_processor. Met data were not ingested. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Denver, CO.

17 WIND ENERGY↗
DOE OSTI2023

Optimizing Hydrogen Fueling Infrastructure Plans on Freight Corridors for Heavy-Duty Fuel Cell Electric Vehicles

The development of a future hydrogen energy economy will require the development of several hydrogen market and industry segments including a hydrogen-based commercial freight transportation ecosystem. For a sustainable freight transportation ecosystem, the supporting fueling infrastructure and the associated vehicle powertrains making use of hydrogen fuel will need to be co-established. This article introduces the OR-AGENT (Optimal Regional Architecture Generation for Electrified National Transportation) tool developed at the Oak Ridge National Laboratory, which has been used to optimize the hydrogen refueling infrastructure requirements on the I-75 corridor for heavy-duty (HD) fuel cell electric commercial vehicles (FCEV). This constraint-based optimization model considers existing fueling locations, regional-specific vehicle fuel economy and weight, vehicle origin and destination (O-D), and vehicle volume by class and infrastructure costs to characterize in-mission refueling requirements for a given freight corridor. The authors applied this framework to determine the ideal public access locations for hydrogen refueling (constrained by existing fueling stations), the minimal viable cost to deploy sufficient hydrogen fuel dispensers, and associated equipment, to accommodate a growing population of hydrogen fuel cell trucks. The framework discussed in this article can be expanded and applied to a larger interstate system, expanded regional corridor, or other transportation network. This article is the third in a series of papers that defined the model development to optimize a national hydrogen refueling infrastructure ecosystem for HD commercial vehicles.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2020

Optimization Framework to Assess the Demand Response Capacity of a Water Distribution System

As large electricity consumers, water distribution system (WDS) pumping stations have the potential to become meaningful participants in demand response (DR) programs. The authors propose an optimization framework for assessing the DR capacity of a WDS and identifying the optimal bidding strategy for maximizing WDS revenue in the DR spot market. The proposed mixed integer linear programming (MILP) model overcomes computational constraints of previous DR optimization models by adopting a preprocessing procedure to minimize the number of binary variables and implementing a convex relaxation technique to linearize the hydraulic equations. The proposed MILP model also explicitly accounts for varying levels of risk tolerance of WDS operators by varying the recovery period over which pumping returns to business-as-usual operation. The optimization framework is implemented on a skeletonized 48-node WDS model that includes 7 pumps, 6 tanks, and 39 pipes. Using a simulated DR event and water consumption profile, the authors derive the optimal DR supply curves (i.e., compensation price versus load curtailment quantity) and revenue potential of the WDS under six scenarios for DR participation.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2025

BLOC Site - ASSIST Thermodynamic retrievals TROPoe / Derived Data

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). The profiles are retrieved every 10 minutes from instantaneous radiances observed with an Atmospheric Sounder Spectrometer by Infrared Spectral Technology (ASSIST, Rochette et al. 2009) operated by NOAA Physical Sciences Laboratory on Block Island for WFIP3. The spectral bands used in the retrieval are in the wavenumber range from 612 - 905.4 cm-1 and are specified in Turner and Löhnert (2021). Radiances are noise-filtered but not averaged in time to minimize errors due to non-uniform clouds. Additional input data in TROPoe are cloud base height from a collocated ceilometer operated by NOAA GML and temperature, water vapor mixing ratio, and pressure from a sensor attached to the ASSIST housing. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Upton, NY.

17 WIND ENERGY↗
NASA NTRS1991

Convection in three dimensions with surface plates - Generation of toroidal flow

This work presents numerical calculations of mantle convection that incorporate some of the basic observational constraints imposed by plate tectonics. The model is three-dimensional and includes surface plates; it allows plate velocity to change dynamically according to the forces which result from convection. It is shown that plates are an effective means of introducing a toroidal component into the flow field. After initial transients the plate motion is nearly parallel to transform faults and in the direction that tends to minimize the toroidal flow field. The toroidal field decays with depth from its value at the surface; the poloidal field is relatively constant throughout the layer but falls off slightly at the top and bottom boundaries. Layered viscosity increasing with depth causes the toroidal field to decay more rapidly, effectively confining it to the upper, low-viscosity layer. The effect of viscosity layering on the poloidal field is relatively small, which is attributed to its generation by temperature variations distributed throughout the system. The generation of toroidal flow by surface plates would seem to account for the observed nearly equal energy of toroidal and poloidal fields of plate motions on the earth. A low-viscosity region in the upper mantle will cause the toroidal flow to decay significantly before reaching the lower mantle. The resulting concentration of toroidal flow in the upper mantle may result in more thorough mixing there and account for some of the geochemical and isotopic differences proposed to exist between the upper and lower mantles.

Gable, Carl W.↗
DOE OSTI2026

FC Site 4.0 - NLR Thermodynamic profiler (ASSIST II-11) Thermodynamic Retrievals TROPoe

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe v0.19 (Turner and Löhnert 2014; Turner and Blumberg 2019; Turner and Löhnert 2021). The profiles are retrieved every 10 minutes from instantaneous observations from the NLR ASSIST II infrared spectrometer. Observations are noise-filtered but not averaged in time to minimize errors due to non-uniform clouds. Additional input data in TROPoe are cloud base height (CBH) from co-located scanning lidar. The full pipeline for running the retrieval is available at https://github.com/StefanoWind/TROPoe_processor. Met data was not ingested. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al. 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches in Denver, CO.

17 WIND ENERGY↗
NASA NTRS2019

Space Launch System and Missions to Near Rectilinear Halo Orbits

The first major evolution of NASA’s Space Launch System (SLS) will begin its flights starting in the mid-2020’s. This new configuration, called Block 1B, replaces the Interim Cryogenic Propulsion Stage with a larger Exploration Upper Stage (EUS). The additional capability provided by the new upper stage will allow SLS to send heavier payloads into deep space. One destination of interest to the SLS program is called a Near Rectilinear Halo Orbit (NRHO). This is a type of lunar orbit with multiple advantages for deep space exploration. These benefits include Earth/lunar access, low station-keeping requirements, and high communication potential with Earth. Therefore, it is a leading candidate for the proposed Lunar Orbital Platform-Gateway (LOP-G). This paper will provide a detailed assessment of the SLS Block 1B requirements and capabilities for sending payloads to an NRHO. Analysts at Marshall Space Flight Center are producing a multi-year mission availability scan for the SLS Block 1B configuration to a predefined NRHO orbit. The analysis produces an optimized trajectory for each day of the scan window. A maximized payload and minimized propellant requirement are determined for each day. All maneuvers from launch to the end of the Trans-Lunar-Injection (TLI) are being modeled as finite burns. Injections into an NRHO are being modeled as impulsive maneuvers. The payload element is arbitrary, but includes sufficient mass to represent a large habitat or propulsion module. The resulting parameters of payload capability, delta-v requirements, and launch windows length vary over the course of the scan. Many launch days in the scan are eliminated in post-processing as they violate mission constraints such as payload mass and propellant usage to insert into an NRHO. Based off previous one-year scan results for the SLS Block 1B Design Analysis Cycle 2 (DAC-2), it is expected that there will be one to three days per week where the payload is able to insert into an NRHO within SLS constraints objectives. This scan provides results for longer than one-year, allowing analysts to better understand the launch availability and energy requirements of SLS Block 1B over time. The in-space mission design and scans utilize Copernicus, an n-body trajectory optimization tool originally developed out of the University of Texas at Austin with further development at Johnson Space Center in Houston, TX. To seed the in-space trajectory, Copernicus uses a plugin to call a database of SLS ascent trajectories optimized in the Program to Simulate Optimized Trajectories II (POST2). The ascent trajectories are developed using a framework that parametrizes payload mass and LEO inclination, and targets a 100 nmi (nautical mile) altitude circular parking orbit. Though this analysis is specific to the SLS program, it will provide a summary of mission design benefits and constraints associated with generic NRHO access, and may be applied to other programs or concepts that will utilize this orbit.

Statham, T. L.↗
DOE OSTI2023

New Theories for Neutrino Masses and Dark Matter

The main goal of the project is to investigate new theories for neutrino masses and dark matter, and understand their testability in different experiments. The PI proposed several mechanisms for neutrino masses in the context of B-L theories, left-right symmetric theories, supersymmetric theories and grand unified theories. The PI would like to investigate new ways to understand the origin of neutrino masses where the seesaw scale is in the multi-TeV region. One can have a multi-TeV seesaw scale in scenarios where the same symmetry relevant for neutrino masses defines the dark matter relic density, or in the minimal supersymmetric theory based on local B-L. In these theories one can have new ways to understand the testability of the origin of neutrino masses at colliders and low energy experiments. In these theories one predicts the possibility to observe lepton number violating signatures at the LHC, one can predict large contributions to lepton number violating processes such as mu to e conversion, neutrinoless double beta decay experiments and others. The nature of the dark matter in the Universe is one of the most important problems in cosmology. Two of the most popular candidates are the Axions and the Weakly Interacting Massive Particles (WIMPs). The PI would like to investigate the implications of a new electroweak theta term similar to the QCD vacuum angle but in the SU(2) gauge sector of the SM which is physical if the baryon and lepton numbers are broken symmetries. We will investigate the implications of having a new electroweak theta term dark matter portal, the implications for baryogenesis and understand the constraints coming from different experiments such as the searches for electric dipole moments and axion experiments. The implications of having the electroweak theta term in theories with different sources of baryon and lepton number violation will be investigated. WIMPs are perhaps the most appealing candidates for many reasons. The PI would like to investigate the properties of WIMPs in different gauge theories, study the implications of the cosmological bounds on the WIMPs relic density for the symmetry breaking scale in models for new physics. The PI would like to investigate in great detail the predictions for gamma lines, neutrino lines and other striking signatures which can be present in these theories. These studies can have a profound impact in the testability of different theories for physics beyond the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2021

PHOENIX Electrostatic Design

PHOENIX (Portable, High-efficiency, Optimal ENergy Imaging X-rays) is a quasi-DC, electrostatic, vacuum-diode designed as a portable x-ray source with national defense and commercial applications. The patent-pending PHOENIX concept combines a megavoltage, Cockroft-Walton voltage multiplier with a Van de-Graaff electrostatic charge-storage dome to create a vacuum-diode suitable for x-ray production. Naturally this structure must minimize internal electric fields to reduce electrical breakdown while simultaneously reducing size and weight to enhance portability. In this paper we describe the optimization process and model results obtained using the COMSOL multi-physics code. We describe three models: a prototype model built to demonstrate the PHOENIX concept as part of Laboratory Directed Research and Development (LDRD) Mission Foundation Research (MFR) Phase-I , a “back-of-the-envelope” design used as a starting point for further COMSOL optimization, and finally, the optimized geometry implemented in the MFR Phase-II. In all cases compromises resulting from cost, schedule, and manufacturing constraints were taken into account as the design progressed.

42 ENGINEERING↗
DOE OSTI2023

Optimizing Long Term Hydrogen Fueling Infrastructure Plans on Freight Corridors for Heavy Duty Fuel Cell Electric Vehicles

The development of a future hydrogen energy economy will require the development of several hydrogen market and industry segments including a hydrogen based commercial freight transportation ecosystem. For a sustainable freight transportation ecosystem, the supporting fueling infrastructure and the associated vehicle powertrains making use of hydrogen fuel will need to be co-established. This paper develops a long-term plan for refueling infrastructure deployment using the OR-AGENT (Optimal Regional Architecture Generation for Electrified National Transportation) tool developed at the Oak Ridge National Laboratory, which has been used to optimize the hydrogen refueling infrastructure requirements on the I-75 corridor for heavy duty (HD) fuel cell electric commercial vehicles (FCEV). This constraint-based optimization model considers existing fueling locations, regional specific vehicle fuel economy and weight, vehicle origin and destination (OD), vehicle volume by class and infrastructure costs to characterize in-mission refueling requirements for a given freight corridor. The authors applied this framework to determine the ideal long term public access locations for hydrogen refueling (constrained by existing fueling stations and dispensing technology), the minimal viable cost to deploy sufficient hydrogen fuel dispensers, and associated equipment, to accommodate a growing population of hydrogen fuel cell trucks. So the framework discussed in this paper can be expanded and applied to additional electrified powertrains as well as a larger interstate system, expanded regional corridor, or other transportation networks.

08 HYDROGEN↗
DOE OSTI2025

Rhode Island Site - NREL ASSIST Thermodynamic Retrievals TROPoe v0.19 / Derived Data

This dataset contains daily files with thermodynamic profiles retrieved with the optimal estimation physical retrieval TROPoe v0.19 (Turner and Löhnert, 2014; Turner and Blumberg, 2019; Turner and Löhnert, 2021). The profiles are retrieved every 10 minutes from instantaneous observations from the NREL ASSIST-II (SN 12) infrared spectrometer. Observations are noise-filtered but not averaged in time to minimize errors due to non-uniform clouds. Additional input data in TROPoe are cloud base height from the collocated NREL upgraded Galion lidar (Newsom et al., 2019), and surface temperature, relative humidity, and pressure from the collocated PNNL met tower. The full pipeline to run the retrieval is available at https://github.com/StefanoWind/TROPoe_processor. In addition to these temporally resolved input data, TROPoe requires an a priori dataset (prior) that provides mean climatological estimates of thermodynamic profiles and specifies how temperature and humidity covary with height as an input (for details see, e.g., Djalalova et al., 2022). The prior is a key component of the retrieval and provides a constraint on the ill-posed inversion problem. A monthly prior was computed from operational radiosonde launches at Upton, NY.

17 WIND ENERGY↗
NASA NTRS2007

Development of a Deployable Nonmetallic Boom for Reconfigurable Systems of Small Modular Spacecraft

Launch vehicle payload capacity and the launch environment represent two of the most operationally limiting constraints on space system mass, volume, and configuration. Large-scale space science and power platforms as well as transit vehicles have been proposed that greatly exceed single-launch capabilities. Reconfigurable systems launched as multiple small modular spacecraft with the ability to rendezvous, approach, mate, and conduct coordinated operations have the potential to make these designs feasible. A key characteristic of these proposed systems is their ability to assemble into desired geometric (spatial) configurations. While flexible and sparse formations may be realized by groups of spacecraft flying in close proximity, flyers physically connected by active structural elements could continuously exchange power, fluids, and heat (via fluids). Configurations of small modular spacecraft temporarily linked together could be sustained as long as needed with minimal propellant use and reconfigured as often as needed over extended missions with changing requirements. For example, these vehicles could operate in extremely compact configurations during boost phases of a mission and then redeploy to generate power or communicate while coasting and upon reaching orbit. In 2005, NASA funded Phase 1 of a program called Modular Reconfigurable High-Energy Technology Demonstrator Assembly Testbed (MRHE) to investigate reconfigurable systems of small spacecraft. The MRHE team was led by NASA's Marshall Space Flight Center and included Lockheed Martin's Advanced Technology Center (ATC) in Palo Alto and its subcontractor, ATK. One of the goals of Phase 1 was to develop an MRHE concept demonstration in a relevant 1-g environment to highlight a number of requisite technologies. In Phase 1 of the MRHE program, Lockheed Martin devised and conducted an automated space system assembly demonstration featuring multipurpose free-floating robots representing Spacecraft in the newly built Controls and Automation Laboratory (CAL) at the ATC. The CAL lab features a 12' x 24' granite air-bearing table and an overhead simulated starfield. Among the technologies needed for the concept demo were mating interfaces allowing the spacecraft to dock and deployable structures allowing for adjustable separation between spacecraft after a rigid connection had been established. The decision to use a nonmetallic deployable boom for this purpose was driven by the MRHE concept demo requirements reproduced in Table 1.

Rehnmark, Fredrik↗
NASA NTRS2006

Small Body GN&C Research Report: A G-Guidance Enhancement to Increase Mission Functionality

G-Guidance is a robust G&C (guidance and control) algorithm developed under the small-body GN&C task. The G-Guidance scheme utilizes a model predictive control approach, along with a convexification of the governing dynamics, control constraints, and trajectory/state constraints. This report details an enhancement to the FY2005 G-Guidance algorithm; the addition of a fire-second approach helps to nullify velocity errors and hit desired final velocities much more precisely than the original Fire-first scheme developed in FY2005. The enhancement preserves the primary benefit of the algorithm, which is to ensure required thruster silent times during trajectory maneuvers. The fire-second scheme increases the versatility of the algorithm for missions employing G-Guidance. For instance, a landing sequence could employ the fire-second scheme to ensure a null final velocity at landing, whereas an ascent sequence could utilize a fire-first scheme to ensure there is no drift of the spacecraft toward the ground at the beginning of the maneuver. Examples are provided within to demonstrate a fire-first versus fire-second guidance scheme. As in the existing G-Guidance algorithm, the examples and schemes incorporate gravity models and thruster firing times into discrete dynamics that are solved as a optimal control problem to minimize fuel consumption or thruster energy expenditure.

velocity errors↗
DOE OSTI2026

Minimal dark 𝑆⁢𝑈⁡(2) origin of a massless Dirac neutrino

We propose a gauge-symmetry origin of a rank-two Dirac neutrino mass matrix that enforces one exactly massless neutrino, while being consistent with the oscillation data, as well as cosmological constraints. The mechanism relies on a minimal dark 𝑆⁢𝑈⁢(2) 𝐷 gauge symmetry under which one right-handed neutrinolike Weyl fermion is charged, thereby forbidding its Standard Model Yukawa coupling. Quantum consistency then fixes the minimal dark-sector completion: Cancellation of the Witten anomaly requires a second fermionic 𝑆⁢𝑈⁢(2) 𝐷 doublet, while a discrete 𝑍 4 symmetry that forbids Majorana masses allows the two dark doublets to form a vectorlike pair. This anomaly-free completion gives rise to a secluded, confining dark sector that can contain a potential dark matter candidate, linking the protected neutrino texture to dark infrared dynamics.

confinement↗
DOE OSTI2025

Lepton flavor asymmetries: from the early Universe to BBN

Large primordial lepton flavor asymmetries with almost vanishing total baryon-minus-lepton number can evade the usual BBN and CMB constraints if neutrino oscillations lead to perfect flavor equilibration. Solving the momentum averaged quantum kinetic equations (QKEs) describing neutrino oscillations and interactions, we perform the first systematic investigation of this scenario, uncovering a rich flavor structure in stark contradiction to the assumption of simple flavor equilibration. We find (i) a particular direction in flavor space, ∆ne ≃ – 2/3 (– 1)∆n μ for normal (inverted) neutrino mass hierarchy, in which the flavor equilibration is efficient and primordial asymmetries are essentially unconstrained, (ii) a minimal washout factor, ∆$n_{e}^{2}$| BBN ≤ 0.03 (0.016) ∑ α ∆$n_{α}^{2}$| ini yielding a conservative estimate for the allowed primordial asymmetries in a generic flavor direction, and (iii) particularly strong or weak washout if one of the initial flavor asymmetries vanishes due to non-adiabatic muon- or electron-driven MSW transitions. These results open up the possibility of a first-order QCD phase transition facilitated by large lepton asymmetries as well as baryogenesis from large and compensated ∆n e = ∆n μ asymmetries. Our systematic approach of deriving momentum averaged QKEs includes collision terms beyond the damping approximation, energy transfer between the neutrino and electron-photon plasma, and provides a fast and reliable way to investigate the impact of primordial lepton asymmetries at the time of BBN. We publicly release the Mathematica code COFLASY-M on https://github.com/mariofnavarro/COFLASY which solves the QKEs numerically.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗