Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “data reduction methods”

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 1,477 records · Page 82

Multigrid-Reduction-in-Time for the Rotating Shallow Water Equations

We consider multilevel time-parallel methods for the numerical solution of the rotating shallow water equations. In particular, the multigrid-reduction-in-time (MGRIT) algorithm is used for the parallel time integration. An asymptotic model is used at the coarse levels while the full model is employed at the finer levels. The asymptotic model is well-suited for highly oscillatory partial differential equations like the rotating shallow water equations because it can accurately and stably take the required large time-steps on coarse levels. Our work exploits the flexibility of the MGRIT algorithm in terms of the number of levels and relaxation schemes to show some computational benefits, especially with respect to FCF-relaxation and data reuse.

97 MATHEMATICS AND COMPUTING↗

Calibration of strain-gage installations in aircraft structures for the measurement of flight loads

A general method has been developed for calibrating strain-gage installations in aircraft structures, which permits the measurement in flight of the shear or lift, the bending moment, and the torque or pitching moment on the principal lifting or control surfaces. Although the stress in structural members may not be a simple function of the three loads of interest, a straightforward procedure is given for numerically combining the outputs of several bridges in such a way that the loads may be obtained. Extensions of the basic procedure by means of electrical combination of the strain-gage bridges are described which permit compromises between strain-gage installation time, availability of recording instruments, and data reduction time. The basic principles of strain-gage calibration procedures are illustrated by reference to the data for two aircraft structures of typical construction, one a straight and the other a swept horizontal stabilizer.

Skopinski, T H↗

Hydrogen Embrittlement

Hydrogen embrittlement (HE) is a process resulting in a decrease in the fracture toughness or ductility of a metal due to the presence of atomic hydrogen. In addition to pure hydrogen gas as a direct source for the absorption of atomic hydrogen, the damaging effect can manifest itself from other hydrogen-containing gas species such as hydrogen sulfide (H2S), hydrogen chloride (HCl), and hydrogen bromide (HBr) environments. It has been known that H2S environment may result in a much more severe condition of embrittlement than pure hydrogen gas (H2) for certain types of alloys at similar conditions of stress and gas pressure. The reduction of fracture loads can occur at levels well below the yield strength of the material. Hydrogen embrittlement is usually manifest in terms of singular sharp cracks, in contrast to the extensive branching observed for stress corrosion cracking. The initial crack openings and the local deformation associated with crack propagation may be so small that they are difficult to detect except in special nondestructive examinations. Cracks due to HE can grow rapidly with little macroscopic evidence of mechanical deformation in materials that are normally quite ductile. This Technical Memorandum presents a comprehensive review of experimental data for the effects of gaseous Hydrogen Environment Embrittlement (HEE) for several types of metallic materials. Common material screening methods are used to rate the hydrogen degradation of mechanical properties that occur while the material is under an applied stress and exposed to gaseous hydrogen as compared to air or helium, under slow strain rates (SSR) testing. Due to the simplicity and accelerated nature of these tests, the results expressed in terms of HEE index are not intended to necessarily represent true hydrogen service environment for long-term exposure, but rather to provide a practical approach for material screening, which is a useful concept to qualitatively evaluate the severity of hydrogen embrittlement. The effects of hydrogen gas on mechanical properties such as tensile strength, ductility, fracture, low and high cycle fatigue, crack growth rate, and creep rupture are analyzed with respect to the general trends established from the HEE index values. It is observed that the severity of HE effects is also influenced by environmental factors such as pressure, temperature, and hydrogen gas purity. The severity of HE effects is also influenced by material factors such as surface finish, heat treatment, and product forms, compositions, grain direction, and crystal orientations.

Jonathan A. Lee↗

Computational Reduction of Specimen Noise to Enable Improved Thermography Characterization of Flaws in Graphite Polymer Composites

Flaw detection and characterization with thermographic techniques in graphite polymer composites are often limited by localized variations in the thermographic response. Variations in properties such as acceptable porosity, fiber volume content and surface polymer thickness result in variations in the thermal response that in general cause significant variations in the initial thermal response. These result in a "noise" floor that increases the difficulty of detecting and characterizing deeper flaws. A method is presented for computationally removing a significant amount of the "noise" from near surface porosity by diffusing the early time response, then subtracting it from subsequent responses. Simulations of the thermal response of a composite are utilized in defining the limitations of the technique. This method for reducing the data is shown to give considerable improvement characterizing both the size and depth of damage. Examples are shown for data acquired on specimens with fabricated delaminations and impact damage.

Winfree, William P.↗

Parallel Algorithms for Computing the Tensor-Train Decomposition

The tensor-train (TT) decomposition expresses a tensor in a data-sparse format used in molecular simulations, high-order correlation functions, and optimization. In this paper, we propose four parallelizable algorithms that compute the TT format from various tensor inputs: (1) Parallel-TTSVD for traditional format, (2) PSTT and its variants for streaming data, (3) Tucker2TT for Tucker format, and (4) TT-fADI for solutions of Sylvester tensor equations. We provide theoretical guarantees of accuracy, parallelization methods, scaling analysis, and numerical results. For example, for a d-dimension tensor in $\mathbb{R}$ $n\times∙∙∙$$\times$$n$ a two-sided sketching algorithm PSTT2 is shown to have a memory complexity of $O(n^{[d/2]})$, improving upon $O(n^{d—1})$ from previous algorithms.

97 MATHEMATICS AND COMPUTING↗

Determination of the inflight OTF of orbital earth resources sensors

There are no author-identified significant results in this report. A technique has been investigated that is particularly applicable to evaluation of earth-orbiting multispectral sensors. In all programs employing such sensors, simultaneous underflight photography from aircraft is made during passes of the spacecraft. To evaluate the spacecraft imagery, the two sets of photographs are scanned and digitized with a microdensitometer. The data are then Fourier analyzed, and the spatial frequency spectra is calculated. The spatial frequencies in the underflight trace are scaled to those in the spacecraft trace, and the ratio of the spectra gives the in-flight optical transfer function (OTF) for the orbiting sensor. The method has been applied to an Apollo 9 frame and the sensor OTF has been calculated. The results are good enough to encourage use of the technique and to indicate the accuracies required of the various measurements involved in determining in-flight sensor OTF by this method. A technique for determining correct scan registration and scaling between the two sets of imagery and a promising approach to noise reduction, in the form of weighted averaging of OTF's, are discussed.

Slater, P. N.↗

Tracking Blade Tip Vortices for Numerical Flow Simulations of Hovering Rotorcraft

Blade tip vortices generated by a helicopter rotor blade are a major source of rotor noise and airframe vibration. This occurs when a vortex passes closely by, and interacts with, a rotor blade. The accurate prediction of Blade Vortex Interaction (BVI) continues to be a challenge for Computational Fluid Dynamics (CFD). Though considerable research has been devoted to BVI noise reduction and experimental techniques for measuring the blade tip vortices in a wind tunnel, there are only a handful of post-processing tools available for extracting vortex core lines from CFD simulation data. In order to calculate the vortex core radius, most of these tools require the user to manually select a vortex core to perform the calculation. Furthermore, none of them provide the capability to track the growth of a vortex core, which is a measure of how quickly the vortex diffuses over time. This paper introduces an automated approach for tracking the core growth of a blade tip vortex from CFD simulations of rotorcraft in hover. The proposed approach offers an effective method for the quantification and visualization of blade tip vortices in helicopter rotor wakes. Keywords: vortex core, feature extraction, CFD, numerical flow visualization

Numerical↗

Pyrolysis Vapor and Bio-Oil Preconditioning via Ex Situ Hydrodeoxygenation and Alkylation Using a Heteropolyacid Catalyst

Ex situ catalytic preconditioning of biomass pyrolysis vapors (partial deoxygenation) has the potential for both reducing the oxygen content and beneficially influencing the oxygen speciation of condense bio-oils for enhancement of downstream condensed phase upgrading processes. Here, condensed phase upgrading processes may refer to both coprocessing and hydroprocessing with vacuum gas-oil (VGO). Reduction of oxygen content to produce bio-oils with greater hydrophobicity can enhance VGO-solubility while tailoring of oxygen speciation by removal of reactive carbonyls (e.g., acetic acid) can lessen catalytic deactivation via coking reactions. 1,2 Furthermore, enhanced VGO-solubility and reduced coking leads to greater biogenic carbon incorporation in end-products from condensed phase upgrading processes. By leveraging the hydrodeoxygenation and alkylation activity of a titania-supported molybdenum-heteropolyacid (Mo-HPA/TiO2) catalyst, pyrolysis vapors and resulting condensed bio-oils were effectively preconditioned for downstream upgrading. 3 Data regarding the use of Mo-HPA/TiO2 for pyrolysis vapor preconditioning will be presented on, including method of employment, catalyst characterization, and impacts on bio-oil in relation to downstream condensed phase upgrading processes. Our approach is to: Compare Mo-HPA/TiO2 product distribution to HZSM-5; Investigate Mo-HPA/TiO2 regeneration and characterize changes in acid site character; Assess product differences in terms of enhanced VGO-solubility.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Capacity Design Method for Integrated Convective/Radiant Terminals to Guarantee Overall and Local Environment

To achieve energy efficiency in the operation of heating devices and ensure thermal comfort in indoor environments, the integrated convective/radiant terminals have become an important development direction for heating systems. This study selected a 14 m 2 bedroom hot summer and cold winter region in China as a case study. Computational Fluid Dynamics (CFD) methodology was used to investigate the design of integrated convective/radiant terminals to ensure thermal comfort for both the entire room and partial space. The results indicate that focusing on partial space during the steady-state stage can achieve energy savings of 31.1% compared to guaranteeing the entire room. Additionally, during the start-up stage, there is a significant reduction of convective unit capacity by 17.6% when the start-up time is 15 minutes. These findings provide data support for the design and engineering applications of products of integrated convective/radiant terminals. Besides, the excess convection ratios for the entire room and partial space are specified as 0.69 to 0.80 and 0.73 to 0.83, respectively.

Chi, Junjie↗

Two-phase flow numerical analysis of electrode geometry for alkaline water electrolyzers

Hydrogen is a promising component of a future energy-secure and efficient economy, but its competitiveness depends on reducing production costs. One strategy is to operate alkaline water electrolyzers at higher current densities to increase output. However, this intensifies performance losses due to gas bubble accumulation, which blocks transport pathways and deactivates electrochemically active surfaces. Enhancing bubble evacuation through electrode design is therefore essential. Previous studies have explored various approaches — such as modifying surface morphology, applying sonication or pressure modulation, and introducing surfactants — but these efforts have addressed a limited range of conditions due to the complexity of two-phase flow and electrode geometries. Experiments have also largely been focused on either cell level improvements, which lack the information necessary to isolate each contributing factor, or on modified geometries that are not relevant to practical cell operation. From a modeling perspective, conventional Eulerian multiphase models do not track the complex gas–liquid interfacial dynamics and often neglect surface tension and contact angle effects, reducing their predictive accuracy. To provide insights on the effects of different electrode geometries on the performance of alklaine water electrolyzers this work employs an immersed boundary volume-of-fluid method to simulate bubble behavior in 3D porous electrodes. Multiple base electrode geometries, typically used in practice, with varying porosity are evaluated under a constant surface gas generation rate. Simulation data is analyzed to quantify electrode gas coverage, bubble size dynamics and other relevant metrics. Results show that porosity strongly influences bubble accumulation on electrode surfaces, with higher porosity reducing gas coverage, and its not strictly dependent on the electrode geometry. However, the electrode’s base geometry significantly affects gas accumulation at the separator gap, independent of porosity. A foam electrode geometry resulted in the lowest gas coverage of all electrodes with a median volumetric gas coverage of 11%, but at the cost of a 70% reduction in active area compared with the largest surface area electrode, while gyroid electrodes showed the best trade-off between gas coverage, particularly at the separator surface, and electrochemically active area. In conclusion, the results highlight the need for holistic electrode design strategies.

Alkaline↗

Budget constrained machine learning for early prediction of adverse outcomes for COVID-19 patients

Abstract The combination of machine learning (ML) and electronic health records (EHR) data may be able to improve outcomes of hospitalized COVID-19 patients through improved risk stratification and patient outcome prediction. However, in resource constrained environments the clinical utility of such data-driven predictive tools may be limited by the cost or unavailability of certain laboratory tests. We leveraged EHR data to develop an ML-based tool for predicting adverse outcomes that optimizes clinical utility under a given cost structure. We further gained insights into the decision-making process of the ML models through an explainable AI tool. This cohort study was performed using deidentified EHR data from COVID-19 patients from ProMedica Health System in northwest Ohio and southeastern Michigan. We tested the performance of various ML approaches for predicting either increasing ventilatory support or mortality. We performed post hoc analysis to obtain optimal feature sets under various budget constraints. We demonstrate that it is possible to achieve a significant reduction in cost at the expense of a small reduction in predictive performance. For example, when predicting ventilation, it is possible to achieve a 43% reduction in cost with only a 3% reduction in performance. Similarly, when predicting mortality, it is possible to achieve a 50% reduction in cost with only a 1% reduction in performance. This study presents a quick, accurate, and cost-effective method to evaluate risk of deterioration for patients with SARS-CoV-2 infection at the time of clinical evaluation.

97 MATHEMATICS AND COMPUTING↗

Progressive Failure And Life Prediction of Ceramic and Textile Composites

An engineering approach to predict the fatigue life and progressive failure of multilayered composite and textile laminates is presented. Analytical models which account for matrix cracking, statistical fiber failures and nonlinear stress-strain behavior have been developed for both composites and textiles. The analysis method is based on a combined micromechanics, fracture mechanics and failure statistics analysis. Experimentally derived empirical coefficients are used to account for the interface of fiber and matrix, fiber strength, and fiber-matrix stiffness reductions. Similar approaches were applied to textiles using Repeating Unit Cells. In composite fatigue analysis, Walker's equation is applied for matrix fatigue cracking and Heywood's formulation is used for fiber strength fatigue degradation. The analysis has been compared with experiment with good agreement. Comparisons were made with Graphite-Epoxy, C/SiC and Nicalon/CAS composite materials. For textile materials, comparisons were made with triaxial braided and plain weave materials under biaxial or uniaxial tension. Fatigue predictions were compared with test data obtained from plain weave C/SiC materials tested at AS&M. Computer codes were developed to perform the analysis. Composite Progressive Failure Analysis for Laminates is contained in the code CPFail. Micromechanics Analysis for Textile Composites is contained in the code MicroTex. Both codes were adapted to run as subroutines for the finite element code ABAQUS and CPFail-ABAQUS and MicroTex-ABAQUS. Graphic user interface (GUI) was developed to connect CPFail and MicroTex with ABAQUS.

Xue, David Y.↗

Seismic Study of the Dynamics of the Solar Subsurface from SoHO Observations

In collaboration with Dr. Baudin, we have developed and refined the new observational methodology for local helioseismology known as time-distance analysis. Global helioseismology study the solar oscillations as a superposition of resonant modes, whose properties (mode frequencies) reflect the global structure of the sun (sound speed stratification, rotation rate, etc). In contrast, local helioseismology look at the solar oscillations as wave packets whose propagation will be affected by perturbations of the media sampled. These local perturbations (sound speed or velocity flows) will modify the propagation time, that in turn can be used as a diagnostic tool for a given region. From a data reduction perspective, the processing of solar dopplergrams that result in time-distance maps, i.e. propagation times as a function of distance between bounces at the surface, is radically different from the methodology used for global mode analysis. We have, in a first step, further develop the programs needed to carry out such analysis. We have then applied them to NMI data set, and explore the trade-off between various averaging and filtering approaches - steps required to improve the signal-to-noise ratio of correlation maps - and the resulting stability and precision of the fitted propagation times. While excessive averaging (whether over space, propagation distance, or time) will reduce the diagnostic potential of the method, insufficient averaging lead to unstable fits, or uncertainties so large as to hide the information we seek. In a second phase, we have developed the analysis methodology required to infer local properties from perturbation in time propagation. Namely, we have developed time-distance inversion techniques, with an emphasis on inferences of velocity flows from time anomalies. Note also that during the period covered by this grant, all the investigators on this proposal (i.e., Drs. Baudin, Eff-Darwich, Korzennik, and Noyes) took part in the organization of the SOHO 6 /GONG 99 Workshop: Structure and Dynamics of the Interior of the Sun and Sunlike Stars, held on June 1-4 1999 at the Boston Park Plaza Hotel in Boston, Massachusetts, USA. it was very well attended by more than 160 participants from 26 countries from all over the world. The proceedings were published in two volumes as ESA SP-418, with Sessions I-III in Volume 1, and Sessions IV-VI in Volume 2 (1,000 pages in total). The complete contents are also included in digital form on a CD-ROM included with Volume 1. This CD-ROM also contains additional multi-media material that complements some of the contributions.

Korzennik, Sylvain G.↗

Simulating the Effects of Irrigation over the U.S. in a Land Surface Model Based on Satellite Derived Agricultural Data

A novel method is introduced for integrating satellite derived irrigation data and high-resolution crop type information into a land surface model (LSM). The objective is to improve the simulation of land surface states and fluxes through better representation of agricultural land use. Ultimately, this scheme could enable numerical weather prediction (NWP) models to capture land-atmosphere feedbacks in managed lands more accurately and thus improve forecast skill. Here we show that application of the new irrigation scheme over the continental US significantly influences the surface water and energy balances by modulating the partitioning of water between the surface and the atmosphere. In our experiment, irrigation caused a 12% increase in evapotranspiration (QLE) and an equivalent reduction in the sensible heat flux (QH) averaged over all irrigated areas in the continental US during the 2003 growing season. Local effects were more extreme: irrigation shifted more than 100 W/m from QH to QLE in many locations in California, eastern Idaho, southern Washington, and southern Colorado during peak crop growth. In these cases, the changes in ground heat flux (QG), net radiation (RNET), evapotranspiration (ET), runoff (R), and soil moisture (SM) were more than 3 W/m(sup 2), 20 W/m(sup 2), 5 mm/day, 0.3 mm/day, and 100 mm, respectively. These results are highly relevant to continental- to global-scale water and energy cycle studies that, to date, have struggled to quantify the effects of agricultural management practices such as irrigation. Based on the results presented here, we expect that better representation of managed lands will lead to improved weather and climate forecasting skill when the new irrigation scheme is incorporated into NWP models such as NOAA's Global Forecast System (GFS).

Ozdogan, Mutlu↗

Thermal Performance Testing of Single Channel GRCop-84 SLM Components

The surface finish found on components manufactured by sinter laser manufacturing (SLM) is rougher (0.013 - 0.0006 inches) than parts made using traditional fabrication methods. Internal features and passages built into SLM components do not readily allow for roughness reduction processes. Alternatively, engineering literature suggests that the roughness of a surface can enhance thermal performance within a pressure drop regime. To further investigate the thermal performance of SLM fabricated pieces, several GRCop‐84 SLM single channel components were tested using a thermal conduction rig at MSFC. A 20 kW power source running at 25% duty cycle and 25% power level applied heat to each component while varying water flow rates between 2.1 - 6.2 gallons/min (GPM) at a supply pressure of 500 to 750 psi. Each test was allowed to reach quasi‐steady state conditions where pressure, temperature, and thermal imaging data were recorded. Presented in this work are the heat transfer responses compared to a traditional machined OHFC Copper test section. An analytical thermal model was constructed to anchor theoretical models with the empirical data.

Garcia, Chance P.↗

Deciphering Reaction Mechanisms of Molecular Proton Reduction Catalysts with Cyclic Voltammetry: Kinetic vs Thermodynamic Control

The kinetics and thermodynamics of elementary reaction steps involved in the catalytic reduction of protons to hydrogen define the reaction landscape for catalysis. The mechanisms can differ in the order of the elementary proton transfer, electron transfer, and bond-forming steps and can be further differentiated by the sites at which protons and electrons localize. Access to fully elucidated mechanistic, kinetic, and thermochemical details of molecular catalysts is crucial to facilitate the development of new catalysts that operate with optimal efficiency, selectivity, and durability. The mechanism by which a catalyst operates, as well as the kinetics and thermodynamics associated with the individual steps, can often be accessed through electroanalytical studies. Here, this Account details the application of cyclic voltammetry to interrogate reaction mechanisms and quantify the kinetics and thermodynamics of elementary reaction steps for a series of molecular catalysts that mediate electrochemical proton reduction. I distinguish the limiting scenarios wherein a catalyst operates under kinetic control vs thermodynamic control, with a focus on detecting how cyclic voltammetry features shift with proton source strength and concentration, as well as scan rate. For systems that operate under kinetic control, catalytic currents are observed at, or slightly positive toward, the formal potential for the redox process that triggers catalysis. Under thermodynamic control, catalytic responses shift as a function of the proton source pKa and effective pH of the solution. After drawing this distinction, we introduce the appropriate voltammetry experiments and accompanying analytical expressions for extracting key metrics from the data. To illustrate analytical strategies to quantify elementary reaction steps of catalysts operating under kinetic control, I describe our studies of proton reduction catalysts Co(dmgBF 2 ) 2 (CH 3 CN) 2 (dmgBF 2 = difluoroboryl-dimethylglyoxime) and [Ni(P 2 Ph N 2 Ph ) 2 ] 2+ (P 2 Ph N 2 Ph = 1,5-phenyl-3,7-phenyl-1,5-diaza-3,7-diphosphacyclooctane). Here, peak shift analysis, foot-of-the-wave analysis, and plateau current analysis are applied to data sets wherein voltammetric response are recorded as a function of catalyst concentration, proton source concentration, proton source strength, and scan rate to quantify rate constants for elementary proton transfer and bond-forming steps in a catalytic cycle. Further, the case study of [Ni(P 2 Ph N 2 Ph ) 2 ] 2+ illustrates how complementary spectroscopic methods can bolster the mechanistic assignment. Collectively, these two studies showcase how detailed mechanistic studies inform on rate-limiting elementary steps in catalysis and other key processes underpinning catalysis. Second, I present analytical strategies to interrogate catalysts operating under thermodynamic control, centered on the case study of [Ni II (P 2 Ph N 2 Bn ) 2 ] 2+ (P 2 Ph N 2 Bn = 1,5-dibenzyl-3,7-diphenyl-1,5-diaza-3,7-diphosphacyclooctane). Here, the application of nonaqueous Pourbaix theory to extract thermodynamic information is introduced, and the construction of a coupled Pourbaix diagram is detailed. This study identifies ligand-based protonation as the key process that places catalysis under thermodynamic control and influences the reaction mechanism. Together, the work detailed in this Account showcases the utility of electroanalytical methods to disentangle complex reaction mechanisms and extract key thermochemical and kinetic parameters for elementary steps of catalysis. Through detailed presentation of the key analytical expressions that underpin these analyses, this Account seeks to facilitate the adoption of cyclic voltammetry by the community to fully extract kinetic, thermochemical, and mechanistic information on electrochemical small-molecule activation.

catalysts↗

In Situ XAFS, XRD, and DFT Characterization of the Sulfur Adsorption Sites on Cu and Ce Exchanged Y Zeolites

Adsorptive desulfurization with Cu and Ce ion-exchanged Y zeolite (CuCeY) has proven to be an effective method for the removal of sulfur compounds from hydrocarbon fuels. In this study, Cu and Ce exchanged Y materials including CuY, CeY, and CuCeY were prepared and examined to investigate the mechanism behind the superior sulfur adsorption and selectivity of CuCeY. In situ conditions were used to study the materials as prepared for optimal desulfurization. X-ray diffraction (XRD) confirmed the absence of large well-ordered crystalline phases from metallic or oxide Cu and Ce after the reduction of the samples. The oxidation states and local environments of Cu and Ce were determined using X-ray adsorption fine structure (XAFS) analysis and correlated to theoretical findings obtained from density functional theory (DFT) calculations. XAFS data indicate the successful reduction of Cu species to Cu + and Cu o , and Ce to Ce 3+ . Analysis of XAFS spectra located Cu and Ce within the Y zeolite framework with Cu cations in the six-member ring sites and as small metallic Cu clusters. Ce cations were found to occupy both six-member ring and hexagonal prism sites. Furthermore, these results reveal the structure of CuCeY as prepared for desulfurization and provide insight into its superior sulfur adsorption performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Optimization of foreground moment deprojection for semi-blind CMB polarization reconstruction

Abstract Upcoming Cosmic Microwave Background (CMB) experiments, aimed at measuring primordial CMB polarization B-modes, require exquisite control of instrumental systematics and Galactic foreground contamination. Blind minimum-variance techniques, like the Needlet Internal Linear Combination (NILC), have proven effective in reconstructing the CMB polarization signal and mitigating foregrounds and systematics across diverse sky models without suffering from foreground mismodelling errors. Still, residual foreground contamination from NILC may bias the recovered CMB polarization at large angular scales when confronted with the most complex foreground scenarios.By adding constraints to NILC to deproject statistical moments of the Galactic emission, the Constrained Moment ILC (cMILC) method has been demonstrated to further enhance foreground subtraction, albeit with an associated increase in overall noise variance. Faced with this trade-off between foreground bias reduction and overall variance minimization, there is still no recipe on which moments to deproject and which are better suited for blind variance minimization. To address this, we introduce the optimized cMILC (ocMILC) pipeline, which performs full automated optimization of the required number and set of foreground moments to deproject, pivot parameter values, and deprojection coefficients across the sky and angular scales, depending on the actual sky complexity, available frequency coverage, and experiment sensitivity. The optimal number of moments for deprojection, before paying significant noise penalty, is determined through a data diagnosis inspired by the Generalized NILC (GNILC) method.Validated on B-mode simulations of thePICOspace mission concept with four challenging foreground models, ocMILC exhibits lower Galactic foreground contamination compared to NILC and cMILC at all angular scales, with limited noise penalty. This multi-layer optimization enables the ocMILC pipeline to achieve unbiased posteriors of the tensor-to-scalar ratio, regardless of foreground complexity.

Astronomy & Astrophysics↗