Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “hydrodynamic control”

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 73 records · Page 4
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Optimal Control of an Oscillating Surge Wave Energy Converter

During this project, we experimentally investigated the hydrodynamics and performance of a laboratory-scale oscillating surge wave energy converter (OSWEC).We looked at how flap buoyancy and driveline losses (primarily in the form of stiction) affected the dynamics and performance of the device. In addition, we assessed the influence of flap profile (rounded vs. square edges) on OSWEC hydrodynamics. Through this, we were able to develop a deeper understanding of OSWEC performance and provide guidance on strategies to counteract artifacts that may be present in laboratory models, but are absent in field-scale devices. To do this, we tested a laboratory-scale OSWEC in the Sea Wave Environmental Lab (SWEL) wave tank at the National Renewable Energy Laboratory (NREL). We ran several types of experiments to investigate the hydrodynamics and performance of the device. Overall, we achieved the overall goal of experimentally investigating the hydrodynamics and performance of this device. We discovered important and unexpected trends in performance, and collected time-resolved data to help us further investigate the underlying hydrodynamics responsible for these trends. In addition, we are currently using the time-resolved data from these experiments to build data-driven models of the dynamics, which can in turn be used to inform data-driven model predictive control of this device and address this objective in the future.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

River Organic Carbon Fluxes Modulated by Hydrodynamic Sorting of Particulate Organic Matter

Rivers regulate the global carbon cycle by transferring particulate organic carbon (POC) from terrestrial landscapes to marine sedimentary basins, but the processes controlling the amount and composition of fluvially exported POC are poorly understood. We propose that hydrodynamic sorting processes modify POC fluxes during fluvial transit. We test this hypothesis by studying POC transported along a ~1,200 km reach of the Rio Bermejo, Argentina. Nanoscale secondary ion mass spectrometry revealed that POC was either fine, mineral-associated organic matter, or coarse discrete organic particles. Mineral-associated POC is more resistant to oxidation and has a lower particle settling velocity than discrete POC. Consequently, hydraulic sorting and downstream fining amplify the proportion of fine, mineral-associated POC from ~55% to ~78% over 1,220 km of downstream transit. This suggests that mineral-associated POC has a greater probability of export and preservation in marine basins than plant detritus, which may be oxidized to CO 2 during transit.

58 GEOSCIENCES↗
DOE OSTI2023

The influence of submarine canyons-related processes on recent benthic foraminiferal distribution, Espírito Santo Basin, Southeastern Brazil

The complex topography of submarine canyons may result in different composition of benthic foraminifera assemblages. To understand how trophic, hydrological and sedimentological conditions in submarine canyons can influence the distribution of benthic foraminifera, and to use this information to corroborate paleoenvironmental interpretations for the Holocene, we investigated recent benthic foraminiferal assemblages (total fauna >63 μm) and sedimentological data in two canyons (CANWN and CAND) in the Espírito Santo Basin (ESB) between 18°20' and 21°20' S. Surface sediment samples (0–2 cm) were collected inside the canyons (150 to 1300 m water depth) and in shelf-slope adjacent transects (50 to 1300 m water depth). The density, taxonomic diversity, and assemblage composition of benthic foraminifera change with depth and location. The distinct ecological preferences of the most abundant taxa allowed us to recognize five benthic foraminiferal groups. Three groups (V, III, and I) are present in different bathymetric sectors; Group V: outer shelf (50 m), Group III: upper, and middle – lower slope (150 to 400 m), and Group I: middle – lower slope (1000 to 1300 m). Groups II and IV show no characteristic bathymetric distribution and are present only in CAND and in CANWN, respectively. Group V consists of Hanzawaia boueana, Peneroplis planatus, and Quinqueloculina lamarckiana; Group III is dominated by Globocassidulina rossensis and Trifarina spp.; Group I consists of Globocassidulina crassa, Bolivina lowmani, Gavelinopsis versiformis, Alabaminella weddellensis, and Epistominella exigua. The main species in Group II (CAND, 150, 1000 to 1300 m) are Trifarina angulosa, Globocassidulina subglobosa, and Discorbis vilardeboanus. Group IV (middle – lower CANWN, 1000 to 1300 m), consists mainly of agglutinated species Glomospira charoides, Rhabdammina abyssorum, and Psammosphaera fusca. Further, our data suggest that the quantity (and quality) of food supply, hydrodynamic conditions and sediment properties are the main drivers controlling the bathymetric distribution of benthic foraminiferal assemblages in both canyons. The middle – lower CANNW revealed unstable trophic conditions, related to terrigenous sediment input due to turbidity currents. In CAND, the foraminiferal assemblages ecology indicated sufficient organic matter supply that favors species establishment and diversity, indicating a more productive and less unstable environment than in CANWN.

58 GEOSCIENCES↗
DOE OSTI2022

An Open-Source Frequency-Domain Model for Floating Wind Turbine Design Optimization

A new frequency-domain dynamics model has been developed that uses open-source components to efficiently represent a complete floating wind turbine system. The model, called RAFT (Response Amplitudes of Floating Turbines), incorporates quasi-static mooring reactions, strip-theory and potential-flow hydrodynamics, blade-element-momentum aerodynamics, and linear turbine control. The formulation is compatible with a wide variety of support structure configurations and no manual or time-domain preprocessing steps are required, making RAFT very practical in design and optimization workflows. The model is applied to three reference floating wind turbine designs and its predictions are compared with results from time-domain OpenFAST simulations. There is good agreement in mean offsets as well the statistics and spectra of the dynamic response, verifying RAFT’s general suitability for floating wind analysis. Follow-on work will include verification of potential-flow and turbine-control features and application to optimization problems.

17 WIND ENERGY↗
DOE OSTI2021

The JOREK non-linear extended MHD code and applications to large-scale instabilities and their control in magnetically confined fusion plasmas

JOREK is a massively parallel fully implicit non-linear extended magneto-hydrodynamic (MHD) code for realistic tokamak X-point plasmas. It has become a widely used versatile simulation code for studying large-scale plasma instabilities and their control and is continuously developed in an international community with strong involvements in the European fusion research programme and ITER organization. This article gives a comprehensive overview of the physics models implemented, numerical methods applied for solving the equations and physics studies performed with the code. A dedicated section highlights some of the verification work done for the code. A hierarchy of different physics models is available including a free boundary and resistive wall extension and hybrid kinetic-fluid models. The code allows for flux-surface aligned iso-parametric finite element grids in single and double X-point plasmas which can be extended to the true physical walls and uses a robust fully implicit time stepping. Particular focus is laid on plasma edge and scrape-off layer (SOL) physics as well as disruption related phenomena. Among the key results obtained with JOREK regarding plasma edge and SOL, are deep insights into the dynamics of edge localized modes (ELMs), ELM cycles, and ELM control by resonant magnetic perturbations, pellet injection, as well as by vertical magnetic kicks. Also ELM free regimes, detachment physics, the generation and transport of impurities during an ELM, and electrostatic turbulence in the pedestal region are investigated. Regarding disruptions, the focus is on the dynamics of the thermal quench (TQ) and current quench triggered by massive gas injection and shattered pellet injection, runaway electron (RE) dynamics as well as the RE interaction with MHD modes, and vertical displacement events. Also the seeding and suppression of tearing modes (TMs), the dynamics of naturally occurring TQs triggered by locked modes, and radiative collapses are being studied.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Resistive wall mode stability and resonant field amplification in MAST high beta plasma

The n=1 (n is the toroidal mode number) resistive wall mode (RWM) instability and the resonant field amplification (RFA) due to a stable RWM response are numerically investigated for a MAST high-pressure plasma scenario, utilizing the MARS-F (Liu Y.Q. et al 2000 Phys. Plasmas 7 3681) and MARS-K (Liu Y.Q. et al 2008 Phys. Plasmas 15 112503) codes, for the purpose of understanding the potential RWM behavior in MAST-U plasmas. Ideal kink stability analysis identifies a target plasma, with parameters similar to that of the reference equilibrium reconstructed from a MAST high-pressure discharge, that accesses the RWM regime. The unstable n=1 RWM for the target equilibrium is subject to strong damping by the plasma toroidal flow and/or the drift kinetic effects from thermal particles. As a result, the mode is found to be stable under the experimental flow conditions, or even without flow stabilization if drift kinetic stabilization is included. The stability prediction is robust against variation of the assumed resistive wall minor radius. Active magneto-hydrodynamic (MHD) spectroscopy modeling, using the magnetic coils designed for controlling the edge localized modes (ELMs) in MAST-U as the antenna, shows strong resonant field amplification due to a stable RWM response in the target plasma. Maximal amplification, of 6-7 times larger than the applied vacuum field within the plasma, is obtained assuming a coil phasing of 120-180 degrees between the upper and lower rows of the ELM control coils. The MHD-kinetic hybrid model predicts generally higher RFA than the fluid model, in particular near the low-field and high-field sides of the plasma boundary surface. Furthermore, the MARS-F/K modeling thus shows that the RWM in MAST-U is subject to strong flow and kinetic stabilization, but can nevertheless be detected by active MHD spectroscopy.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

RAFT (Response Amplitudes of Floating Turbines) [SWR-22-40]

RAFT is a Python code for frequency-domain analysis of floating wind turbines. RAFT works in the frequency domain to provide efficient computations of a floating wind system's response spectra accounting for platform hydrodynamics, quasi-static mooring reactions, rotor aerodynamics, and linearized control. It can be used to calculate response amplitude operators (RAO's), power spectral densities, mean/static properties, and response metrics based on mean and root-mean-square values of each output. RAFT serves as the Level 1 model in the WEIS (Wind Energy with Integrated Servo-control) toolkit, and features an OpenMDAO wrapper.

Hall, Matthew↗
DOE OSTI2026

Sensitivity Analysis of Numerical Modeling Input Parameters on Wind Turbine Loads in Deterministic Transient Load Cases

Aero-hydro-elastic-servo numerical models used to design and analyze wind turbines are based on thousands of variable input parameters that dictate the inflow, aerodynamic, structural, and control characteristics of the system as well as sea state, hydrodynamic, and mooring characteristics for fixed-bottom and floating offshore wind turbines. Each of these parameters has some level of uncertainty, which can significantly impact the predicted loads. Understanding the uncertainty in the inputs is critical to understanding the uncertainty in the outputs. This work demonstrates a screening technique to identify which parameters ultimate loads are most sensitive to so that more focus can be given to quantifying the possible range of those parameters. This technique has been demonstrated previously for different turbine and load case types and is extended here for a floating offshore wind turbine in design load cases with transient events both in the inflow and operations. Each load case features a deterministic gust, including variations in wind speed, direction, and shear. Load cases are considered with an operating turbine as well as with prescribed fault, startup, and shutdown procedures. The study found that key input parameters with a large impact on loads include the length of the gust, the magnitude of direction change and speed in the gust, the initial wind speed, and the shape of the gust profile.

17 WIND ENERGY↗
DOE OSTI2025

Hydrodynamic modeling of plasma channel systems for laser plasma accelerators

Structured plasma channels are an essential technology for driving high-gradient, plasma-based acceleration and control of electron and positron beams for advanced concepts accelerators. Laser and gas technologies can permit the generation of long plasma columns known as hydrodynamic, optically-field-ionized (HOFI) channels, which feature low on-axis densities and steep walls. By carefully selecting the background gas and laser properties, one can generate narrow, tunable plasma channels for guiding high intensity laser pulses. Here, we present on the development of simulations of HOFI channels using the FLASH code, a publicly available radiation hydrodynamics code. We explore sensitivities of the channel evolution to laser profile, intensity, and background gas conditions, and identify relevant scalings with laser intensity through a range of practical channel delays.

Accelerators↗
DOE OSTI2021

Relativistic spin hydrodynamics with torsion and linear response theory for spin relaxation

Using the second law of local thermodynamics and the first-order Palatini formalism, we formulate relativistic spin hydrodynamics for quantum field theories with Dirac fermions, such as QED and QCD, in a torsionful curved background. We work in a regime where spin density, which is assumed to relax much slower than other non-hydrodynamic modes, is treated as an independent degree of freedom in an extended hydrodynamic description. Spin hydrodynamics in our approach contains only three non-hydrodynamic modes corresponding to a spin vector, whose relaxation time is controlled by a new transport coefficient: the rotational viscosity. We study linear response theory and observe an interesting mode mixing phenomenon between the transverse shear and the spin density modes. We propose several field-theoretical ways to compute the spin relaxation time and the rotational viscosity, via the Green-Kubo formula based on retarded correlation functions.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Control Design for a Marine Hydrokinetic Cycloturbine Vehicle

Abstract Marine Hydrokinetic (MHK) cycloturbines generate sustainable power by exploiting tidal currents. By powering the turbines and using pitching foils for control, a vehicle comprised of MHK cycloturbines also has the ability to station keep and maneuver. The vehicle consists of four counter-rotating cycloturbines, with hydrofoils oriented perpendicular to the flow in a paddlewheel configuration. Lift and drag generated from these foils sum together to produce thrust. An experimentally tuned simulation model that solves the six-degrees-of-freedom rigid body equations of motion for the MHK vehicle subject to hydrodynamic, hydrostatic, and propulsive forces is used to aid the design of vehicle controllers. Global feedback controllers are initially designed by applying classical control methods to an approximate linear model of the system dynamics. A higher performing nonlinear controller is designed using the nonlinear dynamic inversion (NDI) method. NDI accounts for the nonlinearities of the MHK system and therefore is suitable for a wide range of operating conditions. The response of the classical and NDI controllers to speed, depth, roll, pitch, and yaw commands are evaluated and compared in simulation. The classical controller outperforms the NDI controller for small amplitude maneuvers, although the degradation with NDI is minor. However, in the nonlinear operating regime the NDI controller outperforms the classical controller and the classical controller exhibits instability.

Automation & Control Systems↗
DOE OSTI2023

Influence on Structural Loading of a Wave Energy Converter by Controlling Variable-Geometry Components and the Power Take-Off

Oceans are harsh environments and can impose significant loads on deployed structures. A wave energy converter (WEC) should be designed to maximize the energy absorbed while ensuring the operating wave condition does not exceed the failure limits of the device itself. Therefore, the loads endured by the support structure are a design constraint for the system. Furthermore, the WEC should be adaptable to different sea states. Herethis work uses a WEC-Sim model of a variable-geometry oscillating wave energy converter (VGOSWEC) mounted on a support structure simulated under different wave scenarios. A VGOSWEC resembles a paddle pitching about a fixed hinge perpendicular to the incoming wave fronts. The geometry of the VGOSWEC is varied by opening a series of controllable flaps on the pitching paddle when the structure experiences threshold loads. It is hypothesized that opening the flaps should result in load shedding at the base of the support structure by reducing the moments about the hinge axis. This work compares the hydrodynamic coefficients, natural periods, and response amplitude operators from completely closed to completely open configurations of the controllable flaps. This work shows that the completely open configuration can reduce the pitch and surge loads on the base of the support structure by as much as 80%. Increased loads at the structure’s natural period can be mitigated by an axial power take-off damping acting as an additional design parameter to control the loads at the WEC’s support structure.

16 TIDAL AND WAVE POWER↗
DOE OSTI2022

Influence on Structural Loading of a Wave Energy Converter by Controlling Variable-Geometry Components and the Power Take-Off

Oceans are harsh environments and can impose significant loads on deployed structures. The deployment of wave energy converters (WECs) faces a design challenge with apparently contradictory goals. A WEC should be designed to maximize the energy absorbed while ensuring the operating wave condition does not exceed the failure limits of the device itself. Therefore, the loads endured by the support structure are a design constraint for the system. Adaptability to different sea states is, therefore, highly desirable. This work uses a WEC-Sim model of a variable-geometry oscillating wave energy converter (VGOSWEC) mounted on a support structure simulated under different wave scenarios. A VGOSWEC resembles a paddle pitching about a fixed hinge perpendicular to the incoming wave fronts. Therefore, the hinge experiences loads perpendicular to its axis as it maintains its position. The geometry of the VGOSWEC is varied by opening a series of controllable flaps on the pitching paddle when the structure experiences threshold loads. Because opening the flaps lets the waves transmit through the paddle, it is hypothesized that opening the flaps should result in load shedding at the base of the support structure. The load shedding is achieved by reducing the moments about the hinge axis. This work compares the hydrodynamic coefficients, natural periods, and response amplitude operators from completely closed to completely open configurations of the controllable flaps. The comparisons quantify the effects of letting the waves transmit through the VGOSWEC. This work shows that the completely open configuration can reduce the pitch and surge loads on the base of the support structure by as much as 80%. It was observed that at the paddle's resonance frequency, the loads on the structure increased substantially. This increase in loads can be mitigated by a rotational power take-off damping about the hinge axis. Changing the rotational power take-off damping was identified as an additional design parameter that can be used to control the loads experienced by the WEC's support structure.

fixed and floating structures↗
DOE OSTI2023

Wind/Wave Testing of a 1:70-Scale Performance-Matched Model of the IEA Wind 15 MW Reference Wind Turbine with Real-Time ROSCO Control and Floating Feedback

Experimental results from the Floating Offshore-wind and Controls Advanced Laboratory (FOCAL) experimental program, which tested a performance-matched model of the IEA Wind 15 MW Reference Turbine on a 1:70 scale floating semisubmersible platform, are compared with OpenFAST simulations. Four experimental campaigns were performed, and data from the fourth campaign, which focused on wind and wave testing of the scaled floating wind turbine system, are considered. Simulations of wave-only, wind-only, and wind/wave environments are performed in OpenFAST, and results for key metrics are compared with the experiment. Performance of the real-time Reference OpenSource COntroller (ROSCO) in above-rated wind conditions, including the effects of the floating feedback loop, are investigated. Results show good agreement in mean values for key metrics, and hydrodynamic effects are matched well. Differences in the surge resonant behavior of the platform are identified and discussed. The effect of the controller and floating feedback loop is evident in both the experiment and OpenFAST, showing significant reduction in platform pitch response and tower base bending load near the platform pitch natural frequency.

17 WIND ENERGY↗
DOE OSTI2026

Influence of the drive design on ablation front hydrodynamic instability growth in a capsule implosion at the National Ignition Facility

A critical aspect of inertial confinement fusion experiments lies in the control of ablation front instabilities during the implosion process. The growth of Rayleigh–Taylor and Richtmyer–Meshkov instabilities, seeded by target defects, can significantly degrade the performance of the implosion. Here, this study explores the influence of drive design on ablation front stability using the hydro-growth radiography platform at the National Ignition Facility. This platform allows the measurement of the ablation front hydrodynamic instability growth. Because these measurements are restricted to early convergence, their primary role is to constrain the initial instability growth and benchmark the simulations used to infer shell integrity and performance at peak velocity and ignition. Three ignition designs, Hybrid-E (HyE), SQ-n (“S” for scaling and “Q” for quality), and High temperature, High thickness (HiT) were analyzed. The results demonstrated a high dependence of the growth factor on the choice of design, capsule scale as well as the hohlraum conditions, mainly the radiative temperature and the gold M-band emission. HyE and SQ-n display a similar growth factor range, but their evolution is different due to the design differences in hohlraum conditions. HiT has the lowest growth factor of the studied designs, which is a result of a higher radiative temperature. These results highlight the importance of design choices in controlling instability dynamics. These insights inform future drive design strategies to enhance stability and efficiency in fusion ignition experiments.

Physics↗
DOE OSTI2022

Influence on Structural Loading of a Wave Energy Converter by Controlling Variable-Geometry Components and the Power Take-Off: Preprint

Oceans are a harsh environment and can impose significant loads on deployed structures. Deployment of wave energy converters (WECs) faces a design challenge with apparently contradictory goals. A WEC should be designed to maximize the energy absorbed while ensuring the operating wave condition does not exceed the failure limits of the device itself. Therefore, the loads endured by the support structure are a design constraint for the system. Adaptability to different sea states is, therefore, highly desirable. This work uses a WEC-Sim model of a Variable Geometry Oscillating Wave Energy Converter (VGOSWEC) mounted on a support structure simulated under different wave scenarios. A VGOSWEC resembles a paddle pitching about a fixed hinge perpendicular to the incoming waves-fronts. Therefore, the hinge experiences loads perpendicular to its axis as it maintains its position. The geometry of the VGOSWEC was varied by opening a series of controllable flaps on the pitching paddle when the structure experiences threshold loads. Since opening the flaps lets the waves transmit through the paddle, it is hypothesized that opening the flaps should result in load shedding at the base of the support structure. This load shedding is achieved by reducing the moments about the hinge axis. This work compares the hydrodynamic coefficients, natural periods, and Response Amplitude Operators (RAOs) from completely closed to completely open configurations of the controllable flaps. These comparisons quantify the effects of letting the waves transmit through the VGOSWEC. This work shows that the completely open configuration can reduce the pitch and surge loads on the base of the support structure by as much as 80%. It was observed that at the paddle's resonance frequency, the loads on the structure increased substantially. This increase in loads can be mitigated by a rotational Power Take-off (PTO) damping about the hinge axis. Changing the rotational PTO damping was identified as an additional design parameter that can be used to control the loads experienced by the WEC's support structure.

fixed and floating structures↗
DOE OSTI2024

Precision Loading and Delivery of Molecular Cargo by Size-Controlled Coacervation of Gold Nanoparticles Functionalized with Elastin-like Peptides

Thermoresponsive elastin-like peptides (ELPs) have been extensively investigated in biotechnology and medicine, but little attention has been paid to the process by which coacervation causes ELP-decorated particles to aggregate. Here, using gold nanoparticles (AuNPs) functionalized with a cysteine-terminated 96-repeat of the VPGVG sequence (V96-Cys), we show that the size of the clusters that reversibly form above the ELP transition temperature can be finely controlled in the 250 to 930 nm range by specifying the concentration of free V96-Cys in solution and using AuNPs of different sizes. We further find that the localized surface plasmon resonance peak of the embedded AuNPs progressively red-shifts with cluster size, likely due to an increase in particle–particle contacts. We exploit this fine control over size to homogeneously load precise amounts of the dye Nile Red and the antibiotic Tetracycline into clusters of different hydrodynamic diameters and deliver cargos near-quantitatively by deconstructing the aggregates below the ELP transition temperature. Beyond establishing a key role for free ELPs in the agglomeration of ELP-functionalized particles, our results provide a path for the thermally controlled delivery of precise quantities of molecular cargo. This capability might prove useful in combination photothermal therapies and theranostic applications, and to trigger spatially and temporally uniform responses from biological, electronic, or optical systems.

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

Soft matter approach for creating novel protein hydrogels using fractal whey protein assemblies as building blocks

In this study, mesoscopic sized fractal assembly (FA) particles were prepared using whey protein isolate; the coldset gelation properties of FA particles were investigated in-depth. Two types of FA particles (FA -62 and FA -90) with different mean sizes were synthesized through controlled thermal treatment of whey protein solutions at two concentrations (62 g/L and 90 g/L). Particle characteristics e.g. hydrodynamic radius, ζ-potential and surface hydrophobicity were dependent on pH and structure of FA. Transmission electron microscopy (TEM) observation confirmed the fractal morphology and small-angle X-ray scattering (SAXS) analysis suggested an internal fractal structure for the obtained FA particles. Eight cold-set FA protein gels (2 % w/v) were manufactured by controlling two gelling factors at two levels: pH (5.8 and 7.0) and Ca 2+ content (5 mM and 10 mM). Rheological characteristics in the large amplitude oscillatory shear regime revealed that pH 7.0 gels were softer and elastic while pH 5.8 gels were harder and brittle. Rheology Pipkin diagrams demonstrated that the strain softening/stiffening and the shear-thinning/thickening behaviors may be fine-tuned by manipulating the key gelation factors: e.g. FA structure, pH, and Ca 2+ . The entrainment speed-dependent friction coefficient curves showed that at an intermediate velocity regime (6-250 mm s -1 ), FA -90 particles induced hydrogels had superior lubrication effect compared to FA -62 gels. Further, this work demonstrated a food structure design approach regarding tuning texture and lubrication properties of protein gels without changing protein content and protein composition. The optimized protein hydrogels may be used as texturizer for "cleaner label" food formulas and/or as delivery system for carrying micronutrients.

60 APPLIED LIFE SCIENCES↗