Search NASA⌕ Search

SEARCH · Search NASA

Results for “RANS CFD”

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 37 records · Page 2

Behavior and mechanisms of Doppler wind lidar error in complex terrain: stable flow case study at Perdigão

A numerical experiment is carried out investigating the magnitude of biases in ground-based lidar measurements in complex flow conditions. Biases assessed include those arising from flow curvature and from the interaction of turbulence with the wind field reconstruction (WFR) algorithms used by a WindCube lidars and anemometers. RANS-CFD and WRF-LES simulations were performed for the Perdig˜ao Field Experiment site for a range of atmospheric conditions. Virtual anemometer and lidar data were generated for four locations: two near exposed ridge tops and two in low-speed regions in the valley. The LES data at these four locations show that the scalar inflation terms (the relation between scalar and vector averaged wind speed) for virtual lidar and virtual cups agree very well with predictions using perturbation theory. While the lidar errors vary greatly with location and height, the contribution from the flow curvature tends to be larger than the differences arising from scalar inflation. For one lidar/mast pair near the ridge top, comparisons between simulations and measurements are carried out for a resonant mountain wave event on June 14th, 2017, and for the whole duration of the Perdigão campaign for winds perpendicular to the ridges. The lidar error during the mountain wave, a period of strong stability and low inversion height, is significantly larger than the campaign average. The sensitivity of the lidar error to atmospheric stability is confirmed by the RANS simulations, which suggests strong sensitivity of flow curvature error to stability conditions and to the shape of the wind speed profile near the top of the boundary layer.

17 WIND ENERGY↗

Silent and Efficient Supersonic Bi-Directional Flying Wing

The supersonic bi-directional (SBiDir) flying wing (FW) concept has a great potential to achieve low sonic boom with high supersonic aerodynamic performance due to removal of performance conflict between high speed and low speed by rotating goo in flight. This NIAC Phase 1 research has achieved three objectives: 2) prove the concept based on simulation that it can achieve very low boom with smooth Sine wave ground over-pressure signature and excellent aerodynamic efficiency; 3) conduct trade study to correlate the geometric parameters with sonic boom and aerodynamic performance for further automated design optimization in Phase II. The design methodology developed in Phase I includes three parts: 1) an advanced geometry model, which can vary airfoil meanline angle distribution to control the expansion and shock waves on the airplane surface to mitigate sonic boom and improve aerodynamic efficiency. 2) a validated CFD procedure to resolve near field flow with accurate shock strength. The sonic boom propagation from near field to far field ground is simulated by NASA NF Boom code. The surface friction drag prediction is based on fiat plate correlation adopted by Seebass and supported by the experimental study of Winter and Smith, which is on the conservative side and is more reliable than CFD RANS simulation. 3) a mission analysis tool based on Corke's model that provides design requirements and constraints of supersonic airplanes for range, payload, volume, size, weight, etc. The design mission target is a supersonic transport with cruise Mach number 1.6, 100 passengers, and 4000nm range. The trade study has several very important findings: 1) The far field ground sonic boom signature is directly related to the smoothness of the flow on the airplane surface. The meanline angle distribution is a very effective control methodology to mitigate surface shock and expansion wave strength, and mitigating compression wave coalescing by achieving smooth loading distribution chord-wise. Compared with a linear meanline angle distribution, a design using nonlinear and non-monotonic meanline angle distribution is able to reduce the sonic boom ground loudness by over 20dBP1. The design achieves sonic boom ground loudness less than 70dBP1 and aerodynamic dynamic efficiency 1/D of 8.4. 2) Decreasing sweep angle within the Mach cone will increase 1/D as well as sonic boom. A design with variable sweep from 84° at the very leading edge to 68° at the tip achieves an extraordinarily high 1/D of 10.4 at Mach number 1.6 due to the low wave drag. If no sonic boom constraint is attached, SBiDir-FW concept still has a lot of room to increase the 1/D. 3) The round leading edge and trailing edge under high sweep angle are beneficial to improve aerodynamic performance, sonic boom, and to increase volume of the airplane. 4) Subsonic performance is benefited greatly from the high slenderness of supersonic configuration after rotating goo. A design with excellent supersonic aspect ratio of 0.44, 1/D of 8.g, gives an extraordinary subsonic aspect ration of 10 and 1/D of 1g.7. Two configurations are designed in details to install internal seats, landing gears, and engine installation to demonstrate the feasibility of SBiDir-FW configuration to accommodate all the required volume for realistic airplane. Here we emphasize that the qualitative findings in Phase I are very encouraging, more important than the quantitative results. Qualitative findings give the understanding of physics and provide the path to achieve the ultimate high performance design. The promising quantitative results achieved in Phase I need to be confirmed by wind tunnel testing in Phase II and ultimately proved by flight test. The other important step forward will be made to study the rotation transition from both CFD unsteady simulation and wind tunnel testing.

Zha, Gecheng↗

Aerodynamic Conceptual Design of Boundary Layer Ingestion Propulsor Systems: A Quasi-2D Through Flow Analysis Method and Multi-Fidelity Propulsor Design Framework

A propulsor design framework for maximizing the benefits of boundary layer ingestion is presented. The performance of BLI is strongly affected by all propulsor components, including the boundary layer characteristics (displacement thickness and form factor) of the ingested boundary layer at the inlet, the radical loading characteristics if the fan and exit guide vane (EGV), the area contraction from inlet to nozzle, and flow expansion at the exhaust cone. A strategy and it's associated multi-fidelity design framework are proposed for an efficient conceptual design of the BLI proplusor which inherently differs from the conventional engine. In the frame work, a quasi-2D through flow model served as the underlying fidelity model is introduced to incorporate the radial effect of the boundary layer entering the propulsor. Multi-fidelity design work is conducted to maximize the predefined performance metrics. On top of this efficient quasi-2D model. computational fluid dynamics based 3-D propulsor models are implemented to refine and validate the design. A BLI propulsion system integrated with fuselage is designed to showcase the framework. The performance of the resulting BLI propulsor system is evaluated via body-force model in unstructured Reynolds-Averaged Navier-Stoke (RANS) CFD and improvement is presented.

Lee, B. J.↗

Computational Fluid Dynamics Analysis of the Stall Characteristics of a Wing Design Based on Prandtl's Minimum Induced Drag

Stall characteristics of a wing whose design was based on Prandtl’s minimum induced drag analysis is presented. Flow field is resolved using RANS CFD (Computational Fluid Dynamics) solver OVERFLOW-2. Both in freestream and in ground effect are analyzed. In addition, effect of low-Mach preconditioner on the stall characteristic is presented. Results show that simulations that lack preconditioner predicts higher stall angle as well as much more benign behavior near the stall angle. Stall analysis in freestream show that flow begins to separate at the inboard region. The flow at the tip remains attached until approximately 19.0 degrees angle of attack.

PRANDTL↗

A Correction Method for Unsteady Transonic Aerodynamics

This paper presents a new method for transonic pitching airfoils based on a RANS CFD study and the Theodorsen model of an oscillating pitching flat plate. This study quantifies the deviation of the lift coefficient predictions using CFD from that obtained using the Theodorsen model, which is based on the incompressible potential flow assumption. The present method corrects this theoretical model by modulating the Theodorsen functions by coefficient functions that depend on the reduced frequency and the Mach number. It is demonstrated that the modified theoretical model predicts lift coefficient in good agreement with the CFD results in the Mach number range from incompressible (M =0.2) to transonic (M =0.755) flow for a range of reduced frequencies typical of transonic flutter. The simulations are first validated by comparing pitching NACA0012 airfoil results with experimental results at transonic flight conditions, which establishes the requirements for a grid converged unsteady transonic solution. The hysteresis loop, Cl versus α, attains a grid independent solution that compares well with experiment. The present correction method will guide the development of a new state space model for the Variable Camber Continuous Trailing Edge Flap (VCCTEF) system and eventually a new transfer function that will be incorporated in a new aeroelastic framework leading to an appropriate transonic flutter model for use in the future aircraft systems in development under the NASA Advanced Air Transportation Technologies (AATT) project.

Kaul, Upender K.↗

Rapid Prediction of Installed Jet Noise from RANS

A new method of computing jet noise, called mSrc, was developed on the general principles of acoustic analogies. In the method, the problem of translating turbulent flow energy into acoustic energy at a far-field observer is broken into two parts, a calculation of acoustic source strengths and then their propagation. The acoustic sources are related to turbulent quantities in the jet plume in a robust manner. The propagation, which is more properly computed using a Green's function accounting for nonuniform speed of sound and solid surfaces, is instead modeled using commonly observed features of jet acoustic far field directivity, and by diffraction barrier theory for surfaces. The mSrc method does not require gradients of the predicted flow field, allowing it to make use of robust unstructured RANS CFD methods, including embedded boundary codes. Such codes do not require specification of surface meshes, and auto-refine their grid to resolve flow gradients, putting resolution where it is required without a priori user input. The ability of mSrc to use such radically unstructured flow input results in an efficient method of estimating noise from jet flows from complex nozzles installed on aircraft. Many validation cases are presented to demonstrate the accuracy and range of applicability of the mSrc method for representative jet noise applications.

Bridges, James↗

Rapid Prediction of Installed Jet Noise From RANS

A new method of computing jet noise, called mSrc, was developed on the general principles of acoustic analogies. In the method, the problem of translating turbulent flow energy into acoustic energy at a far-field observer is broken into two parts, a calculation of acoustic source strengths and then their propagation. The acoustic sources are related to turbulent quantities in the jet plume in a robust manner. The propagation, which is more properly computed using a Green's function accounting for nonuniform speed of sound and solid surfaces, is instead modeled using commonly observed features of jet acoustic far field directivity, and by diffraction barrier theory for surfaces. The mSrc method does not require gradients of the predicted flow field, allowing it to make use of robust unstructured RANS CFD methods, including embedded boundary codes. Such codes do not require specification of surface meshes, and auto-refine their grid to resolve flow gradients, putting resolution where it is required without a priori user input. The ability of mSrc to use such radically unstructured flow input results in an efficient method of estimating noise from jet flows from complex nozzles installed on aircraft. Many validation cases are presented to demonstrate the accuracy and range of applicability of the mSrc method for representative jet noise applications.

Bridges, James↗

Extending A Correction Method for Unsteady Transonic Aerodynamics to Variable Camber Continuous Trailing Edge Flap

In this study, we have extended our correction method for transonic aerodynamics to baseline and variable camber continuous trailing edge (VCCTEF) geometries of super-critical airfoils at cruise conditions. This correction method modifies Theodorsen theory, which is based on incompressible potential flow over a flat plate. The method is partly based on CFD RANS simulations of oscillating NASA generic transport model (GTM) baseline cruise and VCCTEF airfoil geometries in pitch and shows promise for addressing transonic flutter problems of aeroelasticity. In the earlier study1, the Theodorsen functions were modulated with non-linear functions of Mach number, M, and reduced frequency, k, and thus the correction was demonstrated for the symmetrical NACA0012 airfoil. In the present study, the modulation functions are also allowed to account for the effects of camber and thickness of the airfoil. This is demonstrated by applying the correction method to various symmetric NACA00xx airfoil geometries (NACA001, NACA0002, NACA0012, NACA0015 and NACA0018), the NASA GTM super-critical baseline and VCCTEF airfoil geometries, e.g., the circular airfoil, VCCTEF222. First results indicate that the correction method accounts for a large variation of thickness of airfoils and camber. The present correction method will guide the development of a new state space model for the VCCTEF system and eventually a new transfer function that will be incorporated in a new aeroelastic framework leading to an appropriate transonic flutter model for use in the future aircraft systems in development under the NASA Advanced Air Transportation Technologies (AATT) project.

Kaul, Upender K.↗

Performance Analysis of Optimized STARC-ABL Designs Across the Entire Mission Profile

Boundary layer ingestion (BLI) offers the potential for significant fuel burn reduction by exploiting strong aeropropulsive interactions. NASA’s STARC–ABL concept uses an electri- cally powered BLI tail cone thruster on what is otherwise a traditional airframe. Despite the traditional airframe of this configuration, aeropropulsive integration is critical to the perfor- mance of the BLI propulsor. Furthermore, due to being electrically powered, the fan pressure ratio and efficiency of the BLI tail cone thruster vary widely across the flight envelope, and this variation in fan performance must be accounted for with the aeropropulsive integration of the BLI system. Thus, accurate performance prediction for this novel propulsion configu- ration requires the use of a coupled aeropropulsive model across the flight envelope. In this work, we analyze the off-design performance of 18 optimized designs using an aeropropulsive model that is built with the OpenMDAO framework to couple 3-D RANS CFD simulations to 1-D thermodynamic cycle analyses. The designs are created via high-fidelity aeropropulsive design optimizations that span a range of fan pressure ratio and thrust values at the cruise conditions for the STARC-ABL concept, which was chosen as the aerodynamic design point for the propulsor. Performance analyses we present herein are then performed at a range of off-design flight conditions that span the flight envelope, including low-speed and low-altitude flight conditions. This study provides the first set of high-fidelity data for the STARC–ABL configuration at off-design conditions, and the results quantify the power savings through BLI compared to a traditional propulsion system across the entire mission profile.

optimization↗

Transonic Correction to Theodorsen's Theory for Oscillating Airfoil in Pitch and Plunge Toward Flutter

This paper presents a transonic correction method for an oscillating airfoil in pitch and plunge. The proposed method applies correction functions to the Theodorsen’s theory to capture the transonic nonlinear aero- dynamics. These correction functions apply necessary corrections to the amplitudes and the phase angles of the unsteady lift and pitching moment coefficients to account for transonic aerodynamics. The proposed method also postulates a correction for the motion of the aerodynamic center which could be induced by moving shocks. A series of unsteady RANS CFD simulations of the airfoil at the mean aerodynamic chord of the Transonic Truss-Braced Wing aircraft are conducted using FUN3D to provide data to construct these transonic correction functions. The computed responses of the unsteady lift and pitching moment coefficients using these transonic correction functions match the CFD simulation results very well even when the pitching moment coefficient is highly nonlinear. A flutter analysis of an airfoil in pitch and plunge illustrates the potential use of the proposed transonic correction method.

Transonic Flutter↗

High-Fidelity Flight Dynamic Analysis of Transonic Truss-Braced Wing

This paper presents a high-fidelity flight dynamic analysis of the Mach 0.8 Transonic Truss-Braced Wing (TTBW). Unsteady RANS CFD simulations of the Mach 0.8 TTBW in pitch, plunge, roll, and yaw oscillations are conducted in FUN3D. The time-domain data are transformed into the frequency-domain data by Fourier series. Transfer functions of the dynamic stability derivatives are then estimated by a frequency-domain re- gression. The dynamic stability derivatives with respect to the angle of attack are determined by the regression of the unsteady aerodynamic coefficients for the plunge motion. The dynamic stability derivatives with re- spect to the pitch rate are determined by the regression of the differential unsteady aerodynamic coefficients for the pitch motion upon the removal of the angle of attack contribution by the plunge motion. Similarly, dynamic stability derivatives with respect to the angle of sideslip, roll rate, and yaw rates are determined from the frequency domain regression. The longitudinal and lateral-direction flight dynamic models of the Mach 0.8 TTBW are constructed from these dynamic stability derivatives. The eigenvalues of the aircraft modes are analyzed to determine the vehicle stability.

Aircraft Stability and Control↗

Nonlinear Dynamic Control Derivative Analysis for Aircraft with Application to Transonic Truss-Braced Wing

This paper presents the development of a nonlinear dynamic control derivative estimation method. A nonlinear aerodynamic model is developed to account for the effects of large control surface deflections on aircraft aerodynamic forces and moments. A series of unsteady RANS CFD simulations is performed to simulate the control surface oscillations at various reduced frequencies. The time-domain data are transformed into the frequency-domain data by a Fourier series analysis. Transfer functions of the dynamic control derivatives are then estimated by a frequency-domain regression. The method is applied to the Transonic Truss-Braced Wing (TTBW) to estimate the dynamic control derivatives for the elevator, aileron, and rudder.

Aircraft Control Derivatives↗

Computationally Efficient Frequency Domain Method for Dynamic Control Derivative Estimation with Application to Transonic Truss-Braced Wing

This paper presents a computationally efficient method for dynamic control derivative estimation technique via control surface oscillation numerical experiments. Unsteady RANS CFD simulations in FUN3D are performed to simulate the control surface oscillations with a prescribed truncated square wave containing sufficient frequencies of interest. The truncated square wave oscillation offers the computational efficiency which reduces the computational cost by almost an order of magnitude compared to a sine wave oscillation. The nonlinear effect of large control surface oscillation amplitudes creates a spillover effect whereby the frequency response at the same input frequencies contains not only the linear aerodynamic response but also nonlinear aerodynamic response. A correction procedure is developed to remove the spillover effect from the linear aerodynamic response. A frequency-domain regression is performed to estimate the dynamic control derivatives after the correction. The results generally agree with the previous results obtained from the sine wave oscillation.

Stability and Control↗

Surrogates for Valve-Controlled Pipe Flow: Accelerating Nuclear Reactor Design

Neural surrogate models are developed to replace expensive steady-state RANS CFD simulations for valve-controlled pipe flow in nuclear reactor design. Using parametric CFD data generated with MOOSE Pronghorn across a range of valve geometry and flow conditions, three approaches are compared: a POD-based reduced-order model, a structured UNet on a cylindrical grid, and unstructured models (DeepONet and BiStride MeshGraphNet) on nondimensionalized point clouds. POD achieves the highest accuracy (99%) with fast inference but requires storing all solution snapshots, while the DeepONet and BSMS-GNN both achieve ~89% accuracy at sub-second inference, with the BSMS-GNN offering superior geometric generalizability. These surrogates enable rapid ranking of candidate valve designs and can warm-start CFD solvers to accelerate convergence, supporting agentic design iteration on the Prometheus platform.

42 - ENGINEERING↗

Establishing Consensus Turbulence Statistics for Hot Subsonic Jets

Many tasks in fluids engineering require knowledge of the turbulence in jets. There is a strong, although fragmented, literature base for low order statistics, such as jet spread and other meanvelocity field characteristics. Some sources, particularly for low speed cold jets, also provide turbulence intensities that are required for validating Reynolds-averaged Navier-Stokes (RANS) Computational Fluid Dynamics (CFD) codes. There are far fewer sources for jet spectra and for space-time correlations of turbulent velocity required for aeroacoustics applications, although there have been many singular publications with various unique statistics, such as Proper Orthogonal Decomposition, designed to uncover an underlying low-order dynamical description of turbulent jet flow. As the complexity of the statistic increases, the number of flows for which the data has been categorized and assembled decreases, making it difficult to systematically validate prediction codes that require high-level statistics over a broad range of jet flow conditions. For several years, researchers at NASA have worked on developing and validating jet noise prediction codes. One such class of codes, loosely called CFD-based or statistical methods, uses RANS CFD to predict jet mean and turbulent intensities in velocity and temperature. These flow quantities serve as the input to the acoustic source models and flow-sound interaction calculations that yield predictions of far-field jet noise. To develop this capability, a catalog of turbulent jet flows has been created with statistics ranging from mean velocity to space-time correlations of Reynolds stresses. The present document aims to document this catalog and to assess the accuracies of the data, e.g. establish uncertainties for the data. This paper covers the following five tasks: Document acquisition and processing procedures used to create the particle image velocimetry (PIV) datasets. Compare PIV data with hotwire and laser Doppler velocimetry (LDV) data published in the open literature. Compare different datasets acquired at roughly the same flow conditions to establish uncertainties. Create a consensus dataset for a range of hot jet flows, including uncertainty bands. Analyze this consensus dataset for self-consistency and compare jet characteristics to those of the open literature. One final objective fulfilled by this work was the demonstration of a universal scaling for the jet flow fields, at least within the region of interest to aeroacoustics. The potential core length and the spread rate of the half-velocity radius were used to collapse of the mean and turbulent velocity fields over the first 20 jet diameters in a highly satisfying manner.

Bridges, James↗

Jet Noise Diagnostics Supporting Statistical Noise Prediction Methods

The primary focus of my presentation is the development of the jet noise prediction code JeNo with most examples coming from the experimental work that drove the theoretical development and validation. JeNo is a statistical jet noise prediction code, based upon the Lilley acoustic analogy. Our approach uses time-average 2-D or 3-D mean and turbulent statistics of the flow as input. The output is source distributions and spectral directivity. NASA has been investing in development of statistical jet noise prediction tools because these seem to fit the middle ground that allows enough flexibility and fidelity for jet noise source diagnostics while having reasonable computational requirements. These tools rely on Reynolds-averaged Navier-Stokes (RANS) computational fluid dynamics (CFD) solutions as input for computing far-field spectral directivity using an acoustic analogy. There are many ways acoustic analogies can be created, each with a series of assumptions and models, many often taken unknowingly. And the resulting prediction can be easily reverse-engineered by altering the models contained within. However, only an approach which is mathematically sound, with assumptions validated and modeled quantities checked against direct measurement will give consistently correct answers. Many quantities are modeled in acoustic analogies precisely because they have been impossible to measure or calculate, making this requirement a difficult task. The NASA team has spent considerable effort identifying all the assumptions and models used to take the Navier-Stokes equations to the point of a statistical calculation via an acoustic analogy very similar to that proposed by Lilley. Assumptions have been identified and experiments have been developed to test these assumptions. In some cases this has resulted in assumptions being changed. Beginning with the CFD used as input to the acoustic analogy, models for turbulence closure used in RANS CFD codes have been explored and compared against measurements of mean and rms velocity statistics over a range of jet speeds and temperatures. Models for flow parameters used in the acoustic analogy, most notably the space-time correlations of velocity, have been compared against direct measurements, and modified to better fit the observed data. These measurements have been extremely challenging for hot, high speed jets, and represent a sizeable investment in instrumentation development. As an intermediate check that the analysis is predicting the physics intended, phased arrays have been employed to measure source distributions for a wide range of jet cases. And finally, careful far-field spectral directivity measurements have been taken for final validation of the prediction code. Examples of each of these experimental efforts will be presented. The main result of these efforts is a noise prediction code, named JeNo, which is in middevelopment. JeNo is able to consistently predict spectral directivity, including aft angle directivity, for subsonic cold jets of most geometries. Current development on JeNo is focused on extending its capability to hot jets, requiring inclusion of a previously neglected second source associated with thermal fluctuations. A secondary result of the intensive experimentation is the archiving of various flow statistics applicable to other acoustic analogies and to development of time-resolved prediction methods. These will be of lasting value as we look ahead at future challenges to the aeroacoustic experimentalist.

Bridges, James E.↗

ECAR-7932 Rev 0 Large Eddy Simulation of MARVEL Reactor Core Subchannel to Evaluate Model Uncertainty of Reynolds-Averaged Navier-Stokes Equation Based Computational Fluid Dynamics Analysis

In the previous work (ECAR-7210), the peak cladding temperature of the MARVEL microreactor has been evaluated by steady-state Reynolds-Averaged Navier-Stokes (RANS) based computational fluid dynamics (CFD) simulations. Although numerical uncertainties of RANS-based CFD simulations has been assessed in ECAR-7210, the model uncertainty of RANS turbulence models must be investigated to resolve the issues related to inaccurate prediction of turbulent heat flux and flow pulsation in a tight lattice rod bundle using the steady-state RANS simulations. Consequently, this ECAR conducted a high-fidelity CFD analysis utilizing Large Eddy Simulation (LES) to generate reference data and investigated the model uncertainty of RANS-based CFD simulations.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Improvements to RANS Modeling for Aeroheating Predictions on Blunt Bodies

Accurate predictions of aeroheating are critical for designing thermal protection systems for planetary entry vehicles. For larger vehicles, turbulence in the boundary layer can substantially increase convective heating. This turbulence must be accurately modeled to ensure the thermal protection system is sufficient. The majority of hypersonic turbulence model development and validation focuses on boundary layers developing over flat-plates or sharp cones; these cases are substantially different than the boundary layer that develops over the heatshield of a blunt body traveling at hypersonic speeds. Planetary missions often use blunt body geometries, such as the 70-degree sphere-cone favored by Mars missions or the 45-degree sphere-cone planned for the upcoming DAVINCI mission. Due to smaller vehicle size and the lower velocities in the stagnation region, planetary entry vehicles have relatively low Reynolds numbers. Surface curvature and high enthalpy gradients create additional challenges. These difficulties must be addressed to obtain high accuracy needed for the ambitious planetary missions in the upcoming decade. This work focuses on both assessing and improving Reynolds-averaged Navier-Stokes (RANS) turbulence models for blunt-body geometries typical of planetary entry vehicles, with a focus on one-equation and two-equation formulations.

Mars2020↗