Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Control variates”

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

Modeling Offshore Wind Farm Performance in Coastal Low-Level Jets Using Coupled Mesoscale-Microscale Large Eddy Simulations

Accurately predicting wind farm reliability under complex offshore atmospheric conditions remains a key challenge, particularly during noncanonical meteorological events such as coastal low-level jets (LLJs). LLJs, characterized by strong nonmonotonic vertical shear and directional veer, depart significantly from the simplified inflow assumptions embedded in conventional design standards, low-fidelity engineering models, and microscale large eddy simulations of the atmospheric boundary layer. In this work, we use the virtual wind farm framework—an exascale, graphics processing unit–accelerated large eddy simulation platform coupled with high-fidelity aeroservoelastic turbine models and advanced mesoscale-microscale coupling via the ExaWind software stack—to investigate turbine responses under realistic LLJ forcing. Simulations are performed over the U.S. North Atlantic offshore domain with the use of meteorological inputs from New York State Energy Research and Development Authority buoy data, focusing on a representative LLJ case impacting the International Energy Agency 15 MW reference turbine. Our results show that LLJs can cause up to 50% power deficits in downstream turbine rows and significantly amplify low-speed shaft and tower loads through nonlinear coupling between complex inflow characteristics and turbine structural dynamics. Two primary mechanisms drive these load amplifications: (1) unique LLJ inflow features—including veer and vertical/lateral shear—and (2) the downstream evolution of the flow under stable thermal stratification, which suppresses turbulence mixing and alters wake recovery. These mechanisms produce streamwise variations in turbine loading not captured by standard hub height–based metrics or existing design load case (DLC) definitions. This study highlights the critical role of rotor-scale flow gradients in driving fatigue and system-level aeroelastic responses, challenging current DLC and control strategies. We advocate the integration of full-flow field, environment-aware wind inputs into load modeling and control algorithms. By leveraging exascale computing to resolve mesoscale-microscale coupling, this work lays the groundwork for next-generation offshore wind turbine design and operation in meteorologically complex marine environments.

17 WIND ENERGY↗
DOE OSTI2025

Host-specific adaptation in Fusarium oxysporum correlates with distinct accessory chromosome content in human and plant pathogenic strains

ABSTRACT Fusarium oxysporumis a cross-kingdom pathogen. While some strains cause disseminated fusariosis and blinding corneal infections in humans, others are responsible for devastating vascular wilt diseases in plants. To better understand the distinct adaptations ofF. oxysporumto animal or plant hosts, we conducted a comparative phenotypic and genetic analysis of two strains: MRL8996 (isolated from a keratitis patient) and Fol4287 (isolated from a wilted tomato [Solanum lycopersicum]). Infection of mouse corneas and tomato plants revealed that, while both strains cause symptoms in both hosts, MRL8996 caused more severe corneal disease in mice, whereas Fol4287 induced more pronounced wilting symptoms in tomato plants.In vitroassays using abiotic stress treatments revealed that the human pathogen MRL8996 was better adapted to elevated temperatures, whereas the plant pathogen Fol4287 was more tolerant to osmotic and cell wall stresses. Both strains displayed broad resistance to antifungal treatment, with MRL8996 exhibiting the paradoxical effect of increased tolerance to higher concentrations of the antifungal caspofungin. We identified a set of accessory chromosomes (ACs) that encode genes with different functions and have distinct transposon profiles between MRL8996 and Fol4287. Interestingly, ACs from both genomes also encode proteins with shared functions, such as chromatin remodeling and post-translational protein modifications. Our phenotypic assays and comparative genomics analyses lay the foundation for future studies correlating genotypes with phenotype and for developing targeted antifungals for agricultural and clinical uses. IMPORTANCE Fusarium oxysporumis a cross-kingdom fungal pathogen that infects both plants and animals. In addition to causing many devastating wilt diseases, this group of organisms was recently recognized by the World Health Organization as a high-priority threat to human health. Climate change has increased the risk ofFusariuminfections, asFusariumstrains are highly adaptable to changing environments. Deciphering fungal adaptation mechanisms is crucial to developing appropriate control strategies. We performed a comparative analysis ofFusariumstrains using an animal (mouse) and plant (tomato) host andin vitroconditions that mimic abiotic stress. We also performed comparative genomics analyses to highlight the genetic differences between human and plant pathogens and correlate their phenotypic and genotypic variations. We uncovered important functional hubs shared by plant and human pathogens, such as chromatin modification, transcriptional regulation, and signal transduction, which could be used to identify novel antifungal targets.

Microbiology↗
DOE OSTI2026

Hardware Aware Mitigation of Timing Side-Channel Vulnerabilities in Critical Infrastructure Software

Program runtime/timing attacks exploit variations in a program’s execution times to extract sensitive information from the program (e.g. encryption keys, sensitive variable data, intellectual property). State-of-the-art solutions to runtime sidechannel attacks attempt to balance the execution time of the sensitive code for different control flow paths to eliminate the timing leakage. However, during the mitigation process, most techniques do not consider the underlying hardware/device on which the target program is supposed to run on. This can lead to over-fixing (unnecessary extra operations), under-fixing (not solving the imbalance properly), and even failures. We propose DISARM, a joint hardware-software methodology (unlike any existing solution) for mitigating runtime side-channel vulnerabilities that utilizes timing values from real embedded devices to generate targeted software fixes. We implement DISARM to support C/C++/Java source codes and validate it across 22 standard benchmarks. DISARM outperforms state-of-the-art solutions such as PENDULUM and DifFuzzAR in terms of execution time overhead (up to −46%), code size overhead (up to −10%), and correctness (no failures) on five different embedded/edge devices.

Suha, Tasneem [University of Maine]↗
DOE OSTI2022

Root angle is controlled by EGT1 in cereal crops employing an antigravitropic mechanism

Root angle in crops represents a key trait for efficient capture of soil resources. Root angle is determined by competing gravitropic versus antigravitropic offset (AGO) mechanisms. Here we report a root angle regulatory gene termed ENHANCED GRAVITROPISM1 ( EGT1 ) that encodes a putative AGO component, whose loss-of-function enhances root gravitropism. Mutations in barley and wheat EGT1 genes confer a striking root phenotype, where every root class adopts a steeper growth angle. EGT1 encodes an F-box and Tubby domain-containing protein that is highly conserved across plant species. Haplotype analysis found that natural allelic variation at the barley EGT1 locus impacts root angle. Gravitropic assays indicated that Hvegt1 roots bend more rapidly than wild-type. Transcript profiling revealed Hvegt1 roots deregulate reactive oxygen species (ROS) homeostasis and cell wall-loosening enzymes and cofactors. ROS imaging shows that Hvegt1 root basal meristem and elongation zone tissues have reduced levels. Atomic force microscopy measurements detected elongating Hvegt1 root cortical cell walls are significantly less stiff than wild-type. In situ analysis identified HvEGT1 is expressed in elongating cortical and stele tissues, which are distinct from known root gravitropic perception and response tissues in the columella and epidermis, respectively. We propose that EGT1 controls root angle by regulating cell wall stiffness in elongating root cortical tissue, counteracting the gravitropic machinery’s known ability to bend the root via its outermost tissues. We conclude that root angle is controlled by EGT1 in cereal crops employing an antigravitropic mechanism.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2020

C{sub 2}N{sub 2} Vertical Profile in Titan’s Stratosphere

In this paper, we present the first measurements of the vertical distribution of cyanogen (C{sub 2}N{sub 2}) in Titan’s lower atmosphere at different latitudes and seasons, using Cassini’s Composite Infrared Spectrometer far-infrared data. We also study the vertical distribution of three other minor species detected in our data: methylacetylene (C{sub 3}H{sub 4}), diacetylene (C{sub 4}H{sub 2}), and H{sub 2}O, in order to compare them to C{sub 2}N{sub 2}, but also to get an overview of their seasonal and meridional variations in Titan’s lower stratosphere from 85 km to 225 km. We measured an average volume mixing ratio of C{sub 2}N{sub 2} of 6.2±0.8×10{sup −11} at 125 km at the equator, but poles exhibit a strong enrichment in C{sub 2}N{sub 2} (up to a factor 100 compared to the equator), greater than what was measured for C{sub 3}H{sub 4} or C{sub 4}H{sub 2}. Measuring C{sub 2}N{sub 2} profiles provides constraints on the processes controlling its distribution, such as bombardment by Galactic Cosmic Rays which seem to have a smaller influence on C{sub 2}N{sub 2} than predicted by photochemical models.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Enabling BWR fuel rod analysis in the BISON fuel performance code

Nuclear fuel vendors around the world are pursuing approaches to sustain the existing nuclear reactor fleet consisting primarily of pressurized-water reactors (PWRs) and boiling-water reactors (BWRs). To support the industry's efforts, advanced modeling and simulation tools need to be capable of analyzing both legacy reactor concepts. BWR fuel rods are significantly different than those used in PWRs, which can affect fuel performance analysis. BWR fuel rods include: (1) an extensive use of gadolinia dopant as a burnable absorber, (2) an axial variation in fuel enrichment and gadolinia content, (3) the inclusion of a liner on the inner cladding surface to mitigate the impact of pellet-clad mechanical interaction (which impacts hydrogen and hydride distribution), (4) a lower initial fill gas pressure, (5) bottom-entry control rods, and (6) a lower coolant pressure that results in the two-phase flow boiling phenomenon. Although the primary focus of BISON has been in the area of PWR and advanced reactor fuel analyses, this paper presents the developments in BISON to support BWR fuel performance analysis. An overview of the models that account for the effects of gadolinia is highlighted. Internal mesh generation capabilities to include a liner is presented. Normal operating and transient (reactivity insertion accident) demonstration problems are presented to illustrate the impact of gadolinia, the hydrogen and hydride evolution due to the presence of the liner, and BISON's ability to simulate axially varying enrichments and dopant concentration. Bottom-entry control effects are captured by the axial power peaking factors present in the demonstration cases. Comparisons to integral experiments from the Halden IFA-681 experiments are discussed as initial validation. Reasonable comparisons are obtained for fuel centerline temperature and rod internal pressure as a function of time. In conclusion, simulations of additional experiments containing Gd 2 O 3 -bearing fuel are necessary to completely validate the code for BWR applications.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2024

Measurements of Penning-Malmberg trap patch potentials and associated performance degradation

Antiprotons created by laser ionization of antihydrogen are observed to rapidly escape the ALPHA trap. Further, positron plasmas heat more quickly after the trap is illuminated by laser light for several hours. These phenomena can be caused by patch potentials—variations in the electrical potential along metal surfaces. A simple model of the effects of patch potentials explains the particle loss, and an experimental technique using trapped electrons is developed for measuring the electric field produced by the patch potentials. The model is validated by controlled experiments and simulations. Published by the American Physical Society 2024

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

Magnetoelectric properties at the Co/AlN(0001) interface

AlN nitride is a large-band-gap polar material that has gained interest due to its ability to become ferroelectric and its compatibility with Si- and GaN-based technologies. This compound could, therefore, be a good candidate to design new artificial multiferroics if combined with a ferromagnetic electrode. For this work, we performed first-principles calculations to investigate the set up of the magnetoelectric coupling at the Co/AlN(0001) interface. Our results describe a complex interfacial atomic structure with variable local magnetic properties as a function of the atom alignment. We predict an average variation of interface spin magnetization, when the polarization is reversed, of 3.83 𝜇B nm −2 , close to the values already reported in the literature at metal/oxide interfaces. This confirms the potential of this AlN-based ferroelectric compound to be used in future voltage-controlled spintronic devices.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Kensor: Coordinated Intelligence from Co-Located Sensors

Internet of Things (IoT) is becoming more pervasive in many installations, including homes, manufacturing plants, and industrial facilities of all kinds. The data that IoT produces is a reflection of usual behavior such as daily routines and scheduled tasks, but also from unexpected behavior due to unintentional or undesirable abnormalities. Here, we focus on achieving coordinated intelligence about normal and abnormal phenomena from multiple sensors that are geographically colocated in close proximity, monitoring and controlling a set of co-located devices. Given a set of co-located sensors, we seek an intelligent approach that would automatically determine the “normal” patterns of behaviors among the correlated sensors. After normal behavior is extracted, later monitoring should detect any deviant variations over time. An example application is an entry monitoring and alert system for facilities such as nuclear reactors, where badge readers, door locks, lights, weight trackers and other co-located sensors at the entry point are collectively tracked. To address this problem, we identify the possible solution approach that can be used to solve its different variants. The implemented model is developed as a combination of rules and Markov Chain methods.

Kotevska, Olivera↗
DOE OSTI2025

Demonstration of a Novel Technology to Manage Electricity Demand in Grid-Independent Military Microgrids

This research was conducted by the National Renewable Energy Laboratory (NREL) in collaboration with the S&C Electric Inc. through funding provided by the ESTCP. The project demonstrates use of cybersecure Automated Demand Response (ADR) technology to effectively manage microgrid loads during grid-independent, also known as "islanded," operation. When military microgrids become isolated from the main electrical grid, they are required to balance electricity supply and demand locally. Given that local generation may be constrained, the prevailing strategy involves shedding all but the most critical loads by tripping smart circuit breakers, which then necessitate manual resetting. This approach is generally implemented at the building level, which means that the buildings with mission-critical activities are exempt from load management and remain fully powered, whereas those deemed non-critical can experience a complete loss of service. In this research we developed a method that allows building automation systems to selectively control their assets in response to load shedding request from a microgrid controller, avoiding total loss of service in contrast to the conventional control approach. A commercial OpenADR client server by GridFabric is used for communication between the microgrid controller and the building management system (BMS). The microgrid controller monitors both generation capacity and various assets within the microgrid and issues a demand reduction request when necessary. This request is communicated to the OpenADR server via Modbus. Upon receiving the request, the OpenADR server forwards it to the BMS utilizing the OpenADR protocol. The BMS is pre-configured with various levels of load reduction strategies based on the controllable assets available, allowing for a nuanced approach to demand reduction. Both lab and field tests were performed that considered load shedding needed to achieve closed transition into island mode and to accommodate changing loads and power source availability while islanded. A commercial microgrid controller was used for these tests with normal programming within the expected constraints of the system capabilities. That is, the solution did not require any specialized modification to the code base of the controller. Given the latency of the round-trip communication path between the microgrid controller and the various devices involved with the load shed processes, there are certain scenarios for which the demonstrated solution are appropriate and some which are not. The methods described in this report can be used for load shedding/restoration during transitions between islanded and grid-tied modes of operation, as well as accommodating normal variations in load and the need to remove a power source from operation for maintenance. These methods should not be used for scenarios that require load shedding within a second or two such as sudden and unanticipated significant load increases or loss of power sources through equipment faults.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2023

Climatic factors and human population changes in Eurasia between the Last Glacial Maximum and the early Holocene

Archaeological records document a significant expansion of populations from the Last Glacial Maximum (LGM, ~23–19 ka) to the early Holocene (EH, ~9 ka) in Eurasia, which is often attributed to the influence of orbital-scale climate changes. Yet, information remains limited concerning the climatic factor(s) which were responsible for conditioning demographic patterns. Here, in this work, we present results from an improved Minimalist Terrestrial Resource Model (MTRM), forced by a transient climate simulation from the LGM to the EH. Simulated potential hunter-gatherer population densities and spatial distributions across Eurasia are supported by observed archaeological sites in Europe and China. In the low latitudes, potential population size change was predominantly controlled by precipitation and its strong influence on plant and animal resources. In the middle-high latitudes, temperature was the dominant driver in influencing potential population size change and animal resource availability. Different regional responses of potential populations to climate change across Eurasia - owing to variations in available food resources between the LGM and EH - provide a better understanding of human dispersal during the Late Pleistocene.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

Understanding the effect of minor alloying elements on helium bubble formation in ferritic-martensitic steels

Ferritic-martensitic steels are promising structural materials for advanced nuclear reactors. To minimize long-term radioactivity, reduced-activation ferritic-martensitic steels have been developed by substituting high-activation elements like Ni and Mo with low-activation elements such as W. However, the impact of these alloying modifications on helium bubble formation, which plays a key role in material swelling, remains unclear. Here, in this study, we compared helium bubble formation in ferritic-martensitic steel T91 and reduced-activation ferritic-martensitic steel F82H. Both materials were irradiated with sequential 100 keV, 150 keV, and 200 keV helium ions to a dose of 0.5 dpa and a helium concentration of 9,000 appm at 500°C. The helium bubbles in F82H exhibited a larger average size and a lower density than those in T91, suggesting differences in minor alloying elements may influence the bubble growth. Here, to investigate the effects of these alloying elements, we characterized radiation-induced segregation near bubbles and grain boundaries. Prominent Ni-Mn-Si enriched clusters were found near bubbles in T91, while only Mn-Si enriched clusters were found near bubbles in F82H. In addition, the obvious Cr enrichment near grain boundaries was absent around bubbles in both steels. The different segregation trends among elements revealed the variations in element diffusion mechanisms and the different sink biases between bubbles and grain boundaries. Cr enrichment near grain boundaries is mostly driven by interstitial-mediated diffusion. However, since bubble growth relies on net vacancy flux, vacancy-mediated diffusion plays a dominant role in controlling element segregation near bubbles. Therefore, Cr enrichment was not found near bubbles. Because of preferential vacancy-drag diffusion for Ni, Si and Mn, these elements were enriched near bubbles. Due to the strong binding energies of vacancies with these solute atoms, the vacancy diffusivity can be reduced near these solutes. Therefore, the more prominent Ni-Si-Mn clustered near helium bubbles in T91 lead to stronger suppression of helium bubble growth compared to F82H.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Improving Data and Prediction Quality of High-Throughput Perovskite Synthesis with Model Fusion

Combinatorial fusion analysis (CFA) is an approach for combining multiple scoring systems using the rank-score characteristic function and cognitive diversity measure. One example is to combine diverse machine learning models to achieve better prediction quality. In this work, we apply CFA to the synthesis of metal halide perovskites containing organic ammonium cations via inverse temperature crystallization. Using a data set generated by high-throughput experimentation, four individual models (support vector machines, random forests, weighted logistic classifier, and gradient boosted trees) were developed. We characterize each of these scoring systems and explore 66 possible combinations of the models. When measured by the precision on predicting crystal formation, the majority of the combination models improves the individual model results. The best combination models outperform the best individual models by 3.9 percentage points in precision. In addition to improving prediction quality, we demonstrate how the fusion models can be used to identify mislabeled input data and address issues of data quality. In particular, we identify example cases where all single models and all fusion models do not give the correct prediction. Experimental replication of these syntheses reveals that these compositions are sensitive to modest temperature variations across the different locations of the heating element that can hinder or enhance the crystallization process. In summary, we demonstrate that model fusion using CFA can not only identify a previously unconsidered influence on reaction outcome but also be used as a form of quality control for high-throughput experimentation.

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

Understanding and Manipulating Helical Nanofilaments in Binary Systems with Achiral Dopants

Here, in this paper, we report the relationship between helical pitch of the helical nanofilament (HNF) phase formed by bent-core molecule NOBOW and the concentration of achiral dopants 5CB and octane, using linearly polarized resonant soft X-ray scattering (RSoXS). Utilizing theory-based simulation, which fits well with the experiments, the molecular helices in the filament were probed and the superstructure of helical 5CB directed by groove of HNFs was observed. Quantitative pitch determination with RSoXS reveals that helical pitch variation is related to 5CB concentration with no temperature dependence. Doping rodlike immiscible 5CB led to a pitch shortening of up to 30%, which was attributed to a change in interfacial tension. By shedding light not only on phase behavior of binary systems but also enabling control over pitch length, our work may benefit various applications of HNF-containing binary systems, including optical rotation devices, circularly polarized light emitters, and chirality transfer agents.

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

Transcription factor binding divergence drives transcriptional and phenotypic variation in maize

Regulatory elements are essential components of plant genomes that have shaped the domestication and improvement of modern crops. However, their identity, function and diversity remain poorly characterized, limiting our ability to harness their full power for agricultural advances using induced or natural variation. Here, in this study, we mapped transcription factor (TF) binding for 200 TFs from 30 families in two distinct maize inbred lines historically used in maize breeding. TF binding comparison revealed widespread differences between inbreds, driven largely by structural variation, that correlated with gene expression changes and explained complex quantitative trait loci such as Vgt1, an important determinant of flowering time, and DICE, an herbivore resistance enhancer. CRISPR–Cas9 editing of TF binding regions validated the function and structure of regulatory regions at various loci controlling plant architecture and biotic resistance. Our maize TF binding catalogue identifies functional regulatory regions and enables collective and comparative analysis, highlighting its value for agricultural improvement.

Galli, Mary [Rutgers Univ., Piscataway, NJ (United↗
DOE OSTI2024

Real-time tracking of structural evolution in 2D MXenes using theory-enhanced machine learning

In situ Electron Energy Loss Spectroscopy (EELS) combined with Transmission Electron Microscopy (TEM) has traditionally been pivotal for understanding how material processing choices affect local structure and composition. However, the ability to monitor and respond to ultrafast transient changes, now achievable with EELS and TEM, necessitates innovative analytical frameworks. Here, we introduce a machine learning (ML) framework tailored for the real-time assessment and characterization of in operando EELS Spectrum Images (EELS-SI). We focus on 2D MXenes as the sample material system, specifically targeting the understanding and control of their atomic-scale structural transformations that critically influence their electronic and optical properties. This approach requires fewer labeled training data points than typical deep learning classification methods. By integrating computationally generated structures of MXenes and experimental datasets into a unified latent space using Variational Autoencoders (VAE) in a unique training method, our framework accurately predicts structural evolutions at latencies pertinent to closed-loop processing within the TEM. This study presents a critical advancement in enabling automated, on-the-fly synthesis and characterization, significantly enhancing capabilities for materials discovery and the precision engineering of functional materials at the atomic scale.

47 OTHER INSTRUMENTATION↗
DOE OSTI2026

DISARM: Target Electronic Device Informed Mitigation of Software Runtime Side-Channel Vulnerabilities

Program runtime/timing attacks exploit variations in a program’s execution times to extract sensitive information from the program (e.g. encryption keys, sensitive variable data, intellectual property). State-of-the-art solutions to runtime side-channel attacks attempt to balance the execution time of the sensitive code for different control flow paths to eliminate the timing leakage. However, during the mitigation process, most techniques do not consider the underlying hardware/device on which the target program is supposed to run on. This can lead to over-fixing (unnecessary extra operations), under-fixing (not solving the imbalance properly), and even failures. Here, we propose DISARM, a joint hardware-software methodology (unlike any existing solution) for mitigating runtime side-channel vulnerabilities that utilizes timing values from real embedded devices to generate targeted software fixes. We implement DISARM to support C/C++/Java source codes and validate it across 22 standard benchmarks. DISARM outperforms state-of-the-art solutions such as PENDULUM and DifFuzzaR in terms of execution time overhead, code size overhead, and correctness on five different embedded/edge devices.

Timing/runtime side-channel↗
DOE OSTI2021

Inter-Subunit Dynamics Controls Tunnel Formation During the Oxygenation Process in Hemocyanin Hexamers

Hemocyanin from horseshoe crab in its active form is a homo-hexameric protein. It exists in open and closed conformations when transitioning between deoxygenated and oxygenated states. Here, we present a detailed dynamic atomistic investigation of the oxygenated and deoxygenated states of the hexameric hemocyanin using explicit solvent molecular dynamics simulations. We focus on the variation in solvent cavities and the formation of tunnels in the two conformational states. By employing principal component analysis and CVAE-based deep learning, we are able to differentiate between the dynamics of the deoxy- and oxygenated states of hemocyanin. Finally, our results identify the deoxygenated open conformation, which adopts a stable, closed conformation after the oxygenation process.

59 BASIC BIOLOGICAL SCIENCES↗