Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamic efficiency”

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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Scale bridging damage model for quasi-brittle metals informed with crack evolution statistics

Computationally efficient methods for bridging length scales, from highly resolved micro/meso-scale models that can explicitly model crack growth, to macro-scale continuum models that are more suitable for modeling large parts, have been of interest to researchers for decades. In this work, an improved brittle damage model is presented for the simulation of dynamic fracture in continuum scale quasi-brittle metal components. Crack evolution statistics, including the number, length, and orientation of individual cracks, are extracted from high-fidelity, finite discrete element method (FDEM) simulations and used to generate effective material moduli that reflect the material’s damaged state over time. This strategy allows for the retention of small-scale physical behaviors such as crack growth and coalescence in continuum scale hydrodynamic simulations. However, the high-fidelity simulations required to generate the crack statistics are computationally expensive. Thus, steps were taken to produce a flexible constitutive model to reduce the number of costly high-fidelity simulations needed to produce accurate results. A new stress based degradation criterion is introduced for the degradation of individual material zones. This allows for the development of a heterogeneous damage distribution within the bulk material. Then a flow stress model is added to the hydrodynamic simulation to account for plasticity in quasi-brittle materials. As a result, the effective moduli model can be applied to a larger range of materials. The effective moduli constitutive model is used to simulate beryllium flyer plate experiments. The results from the continuum scale simulations using statistics from a single high-fidelity simulation are found to be in excellent agreement with numerical and experimental velocity interferometer data. The same set of crack statistics are used to extrapolate the results of a higher rate flyer plate case using the effective moduli model. In conclusion, the extension of this model to higher rate cases shows promise for further reducing the number of costly high-fidelity simulations needed to generate crack statistics.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Boiling Transitions During Droplet Contact on Superheated Nano/Micro-Structured Surfaces

Manipulating surface topography is one of the most promising strategies for increasing the efficiency of numerous industrial processes involving droplet contact with superheated surfaces. In such scenarios, the droplets may immediately boil upon contact, splash and boil, or could levitate on their own vapor in the Leidenfrost state. In this work we report the outcomes of water droplets coming in gentle contact with designed nano/micro-textured surfaces at a wide range of temperatures as observed using high-speed optical and X-ray imaging. Further, we report a paradoxical increase in the Leidenfrost temperature ($T_{\text{LFP}}$) as the texture spacing is reduced below a critical value (~10 μm) that represents a minima in $T_{\text{LFP}}$. Although droplets on such textured solids appear to boil upon contact, our studies suggest that their behavior is dominated by hydrodynamic instabilities implying that the increase in $T_{\text{LFP}}$ may not necessarily lead to enhanced heat transfer. On such surfaces, the droplets display a new regime characterized by splashing accompanied by a vapor jet penetrating through the droplets before they transition to the Leidenfrost state. We provide a comprehensive map of boiling behavior of droplets over a wide range of texture spacings that may have significant implications towards applications such as electronics cooling, spray cooling, nuclear reactor safety and containment of fire calamities.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2016

Predictive Evaluations of Oxygen-Rich Hydrocarbon Combustion Gas-Centered Swirl Coaxial Injectors using a Flamelet-Based 3-D CFD Simulation Approach

NASA Marshall Space Flight Center (MSFC) has embarked upon a joint project with the Air Force to improve the state-of-the-art of space application combustion device design and operational understanding. One goal of the project is to design, build and hot-fire test a 40,000 pound-thrust Oxygen/Rocket Propellant-2 (RP-2) Oxygen-Rich staged engine at MSFC. The overall project goals afford the opportunity to test multiple different injector designs and experimentally evaluate the any effect on the engine performance and combustion dynamics. To maximize the available test resources and benefits, pre-test, combusting flow, Computational Fluid Dynamics (CFD) analysis was performed on the individual injectors to guide the design. The results of the CFD analysis were used to design the injectors for specific, targeted fluid dynamic features and the analysis results also provided some predictive input for acoustic and thermal analysis of the main Thrust Chamber Assembly (TCA). MSFC has developed and demonstrated the ability to utilize a computationally efficient, flamelet-based combustion model to guide the pre-test design of single-element Gas Centered Swirl Coaxial (GCSC) injectors. Previous, Oxygen/RP-2 simulation models utilizing the Loci-STREAM flow solver, were validated using single injector test data from the EC-1 Air Force test facility. The simulation effort herein is an extension of the validated, CFD driven, single-injector design approach applied to single injectors which will be part of a larger engine array. Time-accurate, Three-Dimensional, CFD simulations were performed for five different classes of injector geometries. Simulations were performed to guide the design of the injector to achieve a variety of intended performance goals. For example, two GCSC injectors were designed to achieve stable hydrodynamic behavior of the propellant circuits while providing the largest thermal margin possible within the design envelope. While another injector was designed to purposefully create a hydrodynamic instability in the fuel supply circuit as predicted by the CFD analysis. Future multi-injector analysis and testing will indicate what if any changes occur in the predicted behavior for the single-element injector when the same injector geometry is placed in a multi-element array.

Richardson, Brian R.↗
DOE OSTI2018

Creation of quark–gluon plasma droplets with three distinct geometries

Experimental studies of the collisions of heavy nuclei at relativistic energies have established the properties of the quark–gluon plasma (QGP), a state of hot, dense nuclear matter in which quarks and gluons are not bound into hadrons. In this state, matter behaves as a nearly inviscid fluid that efficiently translates initial spatial anisotropies into correlated momentum anisotropies among the particles produced, creating a common velocity field pattern known as collective flow. In recent years, comparable momentum anisotropies have been measured in small-system proton–proton (p+p) and proton–nucleus (p+A) collisions, despite expectations that the volume and lifetime of the medium produced would be too small to form a QGP. Here we report on the observation of elliptic and triangular flow patterns of charged particles produced in proton–gold (p+Au), deuteron–gold (d+Au) and helium–gold ( 3 He+Au) collisions at a nucleon–nucleon centre-of-mass energy s N N = 200 GeV. The unique combination of three distinct initial geometries and two flow patterns provides unprecedented model discrimination. Hydrodynamical models, which include the formation of a short-lived QGP droplet, provide the best simultaneous description of these measurements.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS2000

Fast Particle Methods for Multiscale Phenomena Simulations

We are developing particle methods oriented at improving computational modeling capabilities of multiscale physical phenomena in : (i) high Reynolds number unsteady vortical flows, (ii) particle laden and interfacial flows, (iii)molecular dynamics studies of nanoscale droplets and studies of the structure, functions, and evolution of the earliest living cell. The unifying computational approach involves particle methods implemented in parallel computer architectures. The inherent adaptivity, robustness and efficiency of particle methods makes them a multidisciplinary computational tool capable of bridging the gap of micro-scale and continuum flow simulations. Using efficient tree data structures, multipole expansion algorithms, and improved particle-grid interpolation, particle methods allow for simulations using millions of computational elements, making possible the resolution of a wide range of length and time scales of these important physical phenomena.The current challenges in these simulations are in : [i] the proper formulation of particle methods in the molecular and continuous level for the discretization of the governing equations [ii] the resolution of the wide range of time and length scales governing the phenomena under investigation. [iii] the minimization of numerical artifacts that may interfere with the physics of the systems under consideration. [iv] the parallelization of processes such as tree traversal and grid-particle interpolations We are conducting simulations using vortex methods, molecular dynamics and smooth particle hydrodynamics, exploiting their unifying concepts such as : the solution of the N-body problem in parallel computers, highly accurate particle-particle and grid-particle interpolations, parallel FFT's and the formulation of processes such as diffusion in the context of particle methods. This approach enables us to transcend among seemingly unrelated areas of research.

Koumoutsakos, P.↗
DOE OSTI2026

Modeling and Experimental Validation of a Direct-Contact Counter-Flow Fluidized Bed Heat Exchanger for Thermal Energy Storage (TES) Applications

Particle-based thermal energy storage (TES) systems are an emerging energy storage technology. The technological advances have reduced costs, making TES more competitive and reliable in the marketplace but an efficient and reliable operation is heavily dependent on coherent heat transfer between air to particles or vice versa. The particle-based TES technologies provide an intermediate system that can store energy for short (0-10 h), long (10-200 h) and seasonal (> 200 h) timescales. The TES systems store energy by converting electricity to thermal energy; electricity can be directly sourced intermittent generation technologies and/or the grid, helping manage peak loads and other mismatches in supply and demand. The overall efficiency of the TES system depends on the performance of system components (particle storage silos and particle transfer mechanism etc.). The particle heat exchanger is one of the key system components that affects the system efficiency. The pressurized fluidized bed heat exchanger (PFB HX) performance is challenging to predict due to the chaotic behavior of particle and fluid interaction. This research presents a computational study of a novel direct-contact, counter-flow and air-to-particles PFB HX, that contributes in advancing the particle-based long-duration TES technologies. For the current analysis an unsteady Eulerian-Eulerian CFD model was developed and validated against experiments performed at the National Laboratory of the Rockies for two particle sizes (600 ..mu..m and 825 ..mu..m ). Following validation, parametric simulations were conducted to evaluate the effects of interphase drag models (Syamlal-O'Brien and Gidaspow), particle size, bed height and the influence of a frictional-viscosity term on hydrodynamics and heat transfer between the air & particles. The key findings from the analysis are: (1) for the studied operating window Syamlal-O'Brien provides superior agreement with measured gas temperatures (errors generally < 10%) while Gidaspow shows large deviations for the coarse particle case; (2) model predictions are most sensitive in the lower 0.2 m above the air distributor where bubble initiation and local mixing dominate interphase heat transfer; (3) representation of the distributor (number of inlet ports) materially affects predicted local mixing and temperature stratification; and (4) the Eulerian-Eulerian framework reproduces bulk thermal trends but shows regime dependent limitations for coarse particles, motivating mesoscale informed closures for scale-up analysis for future studies. These results provide validated guidance for drag selection and distributor design in particle-based thermal energy storage applications. Collectively, the validated model and parametric results quantify key drivers of PHB-HX performance and provide practical guidance for design and optimization. The results provide confidence in the model predictability and provide a step forward to improve on heat exchange performance. The demonstrated performance and modeling approach support the deployment and further development of this novel PHB-HX concept for robust, particle-based long-duration thermal energy storage systems.

25 ENERGY STORAGE↗
DOE OSTI2014

PPPL-5081

This paper formulates the Tokamak Magneto-Hydrodynamics (TMHD), initially outlined by X. Li and L.E. Zakharov [Plasma Science and Technology, accepted, ID:2013-257 (2013)] for proper simulations of macroscopic plasma dynamics. The simplest set of magneto-hydrodynamics equations, sufficient for disruption modeling and extendable to more refined physics, is explained in detail. First, the TMHD introduces to 3-D simulations the Reference Magnetic Coordinates (RMC), which are aligned with the magnetic field in the best possible way. The numerical implementation of RMC is adaptive grids. Being consistent with the high anisotropy of the tokamak plasma, RMC allow simulations at realistic, very high plasma electric conductivity. Second, the TMHD splits the equation of motion into an equilibrium equation and the plasma advancing equation. This resolves the 4 decade old problem of Courant limitations of the time step in existing, plasma inertia driven numerical codes. The splitting allows disruption simulations on a relatively slow time scale in comparison with the fast time of ideal MHD instabilities. A new, efficient numerical scheme is proposed for TMHD.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2017

Stochastic simulation of reaction-diffusion systems: A fluctuating-hydrodynamics approach

Here, we develop numerical methods for stochastic reaction-diffusion systems based on approaches used for fluctuating hydrodynamics (FHD). For hydrodynamic systems, the FHD formulation is formally described by stochastic partial differential equations (SPDEs). In the reaction-diffusion systems we consider, our model becomes similar to the reaction-diffusion master equation (RDME) description when our SPDEs are spatially discretized and reactions are modeled as a source term having Poisson fluctuations. However, unlike the RDME, which becomes prohibitively expensive for an increasing number of molecules, our FHD-based description naturally extends from the regime where fluctuations are strong, i.e., each mesoscopic cell has few (reactive) molecules, to regimes with moderate or weak fluctuations, and ultimately to the deterministic limit. By treating diffusion implicitly, we avoid the severe restriction on time step size that limits all methods based on explicit treatments of diffusion and construct numerical methods that are more efficient than RDME methods, without compromising accuracy. Guided by an analysis of the accuracy of the distribution of steady-state fluctuations for the linearized reaction-diffusion model, we construct several two-stage (predictor-corrector) schemes, where diffusion is treated using a stochastic Crank-Nicolson method, and reactions are handled by the stochastic simulation algorithm of Gillespie or a weakly second-order tau leaping method. We find that an implicit midpoint tau leaping scheme attains second-order weak accuracy in the linearized setting and gives an accurate and stable structure factor for a time step size of an order of magnitude larger than the hopping time scale of diffusing molecules. We study the numerical accuracy of our methods for the Schlögl reaction-diffusion model both in and out of thermodynamic equilibrium. We demonstrate and quantify the importance of thermodynamic fluctuations to the formation of a two-dimensional Turing-like pattern and examine the effect of fluctuations on three-dimensional chemical front propagation. Furthermore, by comparing stochastic simulations to deterministic reaction-diffusion simulations, we show that fluctuations accelerate pattern formation in spatially homogeneous systems and lead to a qualitatively different disordered pattern behind a traveling wave.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

PIAFS: A 2D nonlinear hydrodynamics code to model gaseous optics

The survivability of final optics is expected to be a major challenge for all future inertial fusion energy concepts. Due to their higher damage threshold, gaseous optics have been identified as a promising solution to this problem. Gaseous optics can be created through the photoabsorption of spatially modulated UV light, which induces various chemical processes that heat the gas. This heating leads to a pressure perturbation, which in turn launches a density perturbation that can imprint a refractive index modulation such as a grating. In this article, we introduce a parallel C/C++ code to simulate gaseous optics. PIAFS2D is a high-order conservative finite-difference code to solve the compressible Navier–Stokes equations along with the photochemical heating sources on Cartesian grids. The simulations are validated by the linear theory derived in a previous paper [Michel et al., Phys. Rev. Appl. 22, 024014 (2024)]. For larger perturbations, the behavior of the system—particularly the evolution of the generated acoustic wave—demonstrates strong nonlinearity. PIAFS2D allows the study of nonlinear behaviors and can be used for the design of high-efficiency gaseous optics elements in realistic experimental conditions.

Oudin, A. [Lawrence Livermore National Laboratory ↗
DOE OSTI2022

Pathways towards break even for low convergence ratio direct-drive inertial confinement fusion

Following indirect-drive experiments which demonstrated promising performance for low convergence ratios (below 17), previous direct-drive simulations identified a fusion-relevant regime which is expected to be robust to hydrodynamic instability growth. This paper expands these results with simulated implosions at lower energies of 100 and 270 kJ, and ‘hydrodynamic equivalent’ capsules which demonstrate comparable convergence ratio, implosion velocity and in-flight aspect ratio without the need for cryogenic cooling, which would allow the assumptions of one-dimensional-like performance to be tested on current facilities. A range of techniques to improve performance within this regime are then investigated, including the use of two-colour and deep ultraviolet laser pulses. Finally, further simulations demonstrate that the deposition of electron energy into the hotspot of a low convergence ratio implosion through auxiliary heating also leads to significant increases in yield. Results include break even for 1.1 MJ of total energy input (including an estimated 370 kJ of short-pulse laser energy to produce electron beams for the auxiliary heating), but are found to be highly dependent upon the efficiency with which electron beams can be created and transported to the hotspot to drive the heating mechanism.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2023

Gelation and Re-entrance in Mixtures of Soft Colloids and Linear Polymers of Equal Size

Liquid mixtures composed of colloidal particles and much smaller non-adsorbing linear homopolymers can undergo a gelation transition due to polymer-mediated depletion forces. We now show that the addition of linear polymers to suspensions of soft colloids having the same hydrodynamic size yields a liquid-togel-to-re-entrant liquid transition. In particular, the dynamic state diagram of 1,4-polybutadiene star-linear polymer mixtures was determined with the help of linear viscoelastic and small-angle Xray scattering experiments. While keeping the star polymers below their nominal overlap concentration, a gel was formed upon increasing the linear polymer content. Further addition of linear chains yielded a re-entrant liquid. This unexpected behavior was rationalized by the interplay of three possible phenomena: (i) depletion interactions, driven by the size disparity between the stars and the polymer length scale which is the mesh size of its entanglement network; (ii) colloidal deswelling due to the increased osmotic pressure exerted onto the stars; and (iii) a concomitant progressive suppression of the depletion efficiency on increasing the polymer concentration due to reduced mesh size, hence a smaller range of attraction. Our results unveil an exciting new way to tailor the flow of soft colloids and highlight a largely unexplored path to engineer soft colloidal mixtures.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2026

Mixed-anion electrolytes: from bulk speciation to interfacial dynamics in divalent metal electrodeposition

Mixing anions is emerging as a promising strategy for multivalent electrolyte design, allowing for adjustment of the solvation structure of bulk cations and enhancing the efficiency of electrochemical processes (e.g. metal deposition for batteries and catalysis). Further progress in electrolyte development requires a fundamental understanding of how tailored electrolyte speciation in mixed anion systems can modify the dynamic electrochemical interface during metal cycling. In this study, we present an anode-focused mechanistic study of exemplar Mg electrolytes containing three different secondary anions, correlating electrochemical behavior with bulk speciation and operando interfacial dynamics. Electrospray Ionization-Mass Spectrometry (ESI-MS) results reveal a general trend of forming mixed anion contact ion pairs (CIPs) across various anions, with the extent of ion pairing influenced by the association strength of the secondary anion. Operando multiharmonic electrochemical quartz crystal microbalance with dissipation (EQCM-D) reveals how these bulk species influence interfacial mass uptake, viscoelasticity, and solvent-coupled hydrodynamic behavior during deposition and stripping. The results indicate that Mg-containing ion pairs and solvated complexes shape adsorption, nucleation, and deposit growth, leading to distinct anion-dependent interphases ranging from more permeable, solvent-coupled layers to relatively compact and rigid deposits. This work establishes a quantitative link between bulk speciation and interfacial dynamics in divalent metal electrodeposition and provides mechanistic guidance for electrolyte design.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Effects of stratification on overshooting and waves atop the convective core of M ⊙ main-sequence stars

As a massive star evolves along the main sequence, its core contracts, leaving behind a stable stratification in helium. We simulate two-dimensional convection in the core at three different stages of evolution of a $5\,\mathrm{ M}_{\odot }$ star, with three different stratifications in helium atop the core. We study the propagation of internal gravity waves in the stably stratified envelope, along with the overshooting length of convective plumes above the convective boundary. We find that the stratification in helium in evolved stars hinders radial motions and effectively shields the radiative envelope against plume penetration. This prevents convective overshooting from being an efficient mixing process in the radiative envelope. In addition, internal gravity waves are less excited in evolved models compared to the zero-age-main-sequence model, and are also more damped in the stratified region above the core. As a result, the wave power is several orders of magnitude lower in mid- and terminal-main-sequence models compared to zero-age-main-sequence stars.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Simulation of Mesoscale Cellular Convection in Marine Stratocumulus: 2. Nondrizzling Conditions

This is the second of two large-eddy simulation studies on the mechanisms of mesoscale cellular organization in drizzling (open cells) and nondrizzling marine stratocumulus (closed cells). This study uses a hard nudging approach which maintains fixed horizontal-mean temperature and humidity profiles for a well-mixed boundary layer with a constant boundary layer depth. For the case studied, closed cells develop and broaden by 32 hr to an aspect ratio of 25. Simulations show that the closed-cell mesoscale cellular convection is driven by positive feedback from cloud-induced mesoscale perturbations of longwave radiative cooling. A conceptual model for closed-cell stratocumulus as a mesoscale wavelength hydrodynamic instability in which mesoscale moist and dry anomalies spontaneously grow is presented. In simulations in which long-wavelength sinusoidal moisture anomalies are initially imposed, these anomalies evolve into amplifying closed cells. The cell structure is visualized with a compositing approach based on sorting grid columns by their mesoscale-smoothed total water path. A thermally direct mesoscale circulation pattern develops in the interior of the boundary layer with buoyant mesoscale updrafts, thicker cloud, and a slightly higher capping inversion in the moister columns. There is a mesoscale flow of above-inversion air down the slightly sloping capping inversion from the moist to the dry regions, reinforced by cloud top radiative cooling. This strengthens the mesoscale anomalies by preferentially cooling and drying the already dry regions. The sloping inversion flow is not driven as efficiently if the radiative cooling is artificially horizontally homogenized, partly disrupting this positive feedback and the resulting closed-cell development.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2014

A general hybrid radiation transport scheme for star formation simulations on an adaptive grid

Radiation feedback plays a crucial role in the process of star formation. In order to simulate the thermodynamic evolution of disks, filaments, and the molecular gas surrounding clusters of young stars, we require an efficient and accurate method for solving the radiation transfer problem. We describe the implementation of a hybrid radiation transport scheme in the adaptive grid-based FLASH general magnetohydrodyanmics code. The hybrid scheme splits the radiative transport problem into a raytracing step and a diffusion step. The raytracer captures the first absorption event, as stars irradiate their environments, while the evolution of the diffuse component of the radiation field is handled by a flux-limited diffusion solver. We demonstrate the accuracy of our method through a variety of benchmark tests including the irradiation of a static disk, subcritical and supercritical radiative shocks, and thermal energy equilibration. We also demonstrate the capability of our method for casting shadows and calculating gas and dust temperatures in the presence of multiple stellar sources. Our method enables radiation-hydrodynamic studies of young stellar objects, protostellar disks, and clustered star formation in magnetized, filamentary environments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2021

Full-Loop Simulation of the Combustion of Biomass in a Circulating Fluidized Bed Combustor

Interest in circulating fluidized bed (CFB) boilers as a power generation technology has skyrocketed in recent years because of several advantages this technology offers over conventional boilers, such as increased gas-solid mixing, which results in higher combustion efficiency, and the ability to use lower rank fuels. CFB combustors are operated at lower temperatures than conventional thermal power generation combustors, thus reducing NO x emissions, while SO 2 emissions can be conveniently controlled through the addition of Ca-based sulfur sorbents within the combustor. This report summarizes the development of a comprehensive computational fluid dynamics (CFD) model for biomass combustion in a 50 kW th CFB combustor employing the multiphase particle-in-cell (PIC) approach in the open-source Multiphase Flow with Interphase eXchanges (MFiX) Software Suite at the U.S. Department of Energy (DOE) National Energy Technology Laboratory (NETL). The experimental combustor was designed, built, and operated at CanmetENERGY, Natural Resources Canada in Ottawa, Canada. This work is the culmination of a multi-year collaboration between NETL and CanmetENERGY. The hydrodynamics in the full loop are validated against cold-flow experiments conducted at CanmetENERGY. To mimic the dense packing in the standpipe in the absence of a collisional model in the MFiX-PIC approach, the ram valve in the return leg is partially closed to obstruct the flow. A combustion reaction scheme was developed that considers the volatiles released by the fuel as a lumped species that fully converts to CO 2 and H 2 O. The so-called enhanced simplified combustion scheme produces an excellent match with experimental data in terms of the species concentrations of CO 2 and O 2 at the outlet and the freeboard, demonstrating the utility of the comprehensive model for a range of flow conditions and operating temperatures. Furthermore, the model can be readily adapted to simulate combustion of a variety of biomass fuels as well as co-fired systems under air and oxy-fuel conditions.

09 BIOMASS FUELS↗
DOE OSTI2023

Utilizing the maximum likelihood estimator for flow analysis

We explore the possibility of evaluating flow harmonics by employing the maximum likelihood estimator (MLE). For a given finite multiplicity, the MLE simultaneously furnishes estimations for all the parameters of the underlying distribution function while efficiently suppressing the variance of measures. Also, the method provides a means to assess a specific class of mixed harmonics, which is not straightforwardly feasible by the approaches primarily based on particle correlations. The results are analyzed using the Wald, likelihood ratio, and score tests of hypotheses. Besides, the resultant flow harmonics obtained using MLE are compared with those derived using particle correlations and event plane methods. Here, the dependencies of extracted flow harmonics on the multiplicity of individual events and the total number of events are analyzed. It is shown that the proposed approach works efficiently to deal with the deficiency in detector acceptability. Moreover, we elaborate on a fictitious scenario where the event plane is not a well-defined quantity in the distribution function. For the latter case, the MLE is shown to largely perform better than the two-particle correlation estimator. In this regard, one concludes that the MLE furnishes a meaningful alternative to the existing approaches for flow analysis.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗