Search NASA⌕ Search

SEARCH · Search NASA

Results for “COMPUTER MODELS”

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 127 records · Page 7

Computational Model and Characterization of Stacking Faults in ZIF-8 Polymorphs

Degradation of metal–organic frameworks (MOFs) in aqueous, humid, and acid gas environments likely begins at defect sites. Until now, however, theoretical studies of MOFs have widely assumed an ideal defect-free structure. Here we present a computational model for low-energy extended defects in bulk zeolitic imidazolate frameworks (ZIFs) that are analogous to stacking faults in zeolites. We demonstrate the thermodynamic accessibility of stacking faults in ZIFs and examine the impact of these defects on pore diffusion and accessible surface area. We identify strong correlations between the defect density of a structure and its X-ray diffraction spectra. Finally, by examining a topologically isomorphic ZIF that has been reported experimentally we find characteristic defect-induced peak broadening and splitting in the reported powder patterns, giving strong evidence for the existence for stacking faults in this material.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Computational Modeling of Realistic Cell Membranes

Cell membranes contain a large variety of lipid types and are crowded with proteins, endowing them with the plasticity needed to fulfill their key roles in cell functioning. The compositional complexity of cellular membranes gives rise to a heterogeneous lateral organization, which is still poorly understood. Computational models, in particular molecular dynamics simulations and related techniques, have provided important insight into the organizational principles of cell membranes over the past decades. Now, we are witnessing a transition from simulations of simpler membrane models to multicomponent systems, culminating in realistic models of an increasing variety of cell types and organelles. Here, we review the state of the art in the field of realistic membrane simulations and discuss the current limitations and challenges ahead.

59 BASIC BIOLOGICAL SCIENCES↗

COMPUTATIONAL MODELING OF THE CLASS I LOW-MASS PROTOSTAR ELIAS 29 APPLYING OPTICAL CONSTANTS OF ICES PROCESSED BY HIGH ENERGY COSMIC RAY ANALOGS

We present a study of the effects of high energy cosmic rays (CRs) over the astrophysical ices, observed toward the embedded class I protostar Elias 29, by using computational modeling and laboratory data. Its spectrum was observed with the Infrared Space Observatory (ISO) covering 2.3–190 μm. The modeling employed the three-dimensional Monte Carlo radiative transfer code RADMC-3D and laboratory data of bombarded ice grains by CR analogs and unprocessed ices (not bombarded). We are assuming that Elias 29 has a self-irradiated disk with inclination i = 60.°0, surrounded by an envelope with a bipolar cavity. The results show that absorption features toward Elias 29 are better reproduced by assuming a combination between unprocessed astrophysical ices at low temperature (H{sub 2}O, CO, CO{sub 2}) and bombarded ices (H{sub 2}O:CO{sub 2}) by high energy CRs. Evidences of the ice processing around Elias 29 can be observed by the good fitting around 5.5–8.0 μm, by polar and apolar ice segregation in 15.15–15.25 μm, and by the presence of the CH{sub 4} and HCOOH ices. Given that non-nitrogen compounds were employed in this work, we assume that absorption around 5.5–8.0 μm should not be associated with the NH{sub 4}{sup +} ion (see the 2003 work of Shutte and Khanna ), but more probably with aliphatic ethers (e.g., R1-OCH{sub 2}-R2), CH{sub 3}CHO, and related species. The results obtained in this paper are important because they show that the environment around protostars is better modeled considering processed samples and, consequently, demonstrate the chemical evolution of the astrophysical ices.

79 ASTRONOMY AND ASTROPHYSICS↗

Use of computational modeling combined with advanced visualization to develop strategies for the design of crop ideotypes to address food security

Sustainable crop production is a contributing factor to current and future food security. Innovative technologies are needed to design strategies that will achieve higher crop yields on less land and with fewer resources. Computational modeling coupled with advanced scientific visualization enables researchers to explore and interact with complex agriculture, nutrition, and climate data to predict how crops will respond to untested environments. These virtual observations and predictions can direct the development of crop ideotypes designed to meet future yield and nutritional demands. This review surveys modeling strategies for the development of crop ideotypes and scientific visualization technologies that have led to discoveries in “big data” analysis. Combined modeling and visualization approaches have been used to realistically simulate crops and to guide selection that immediately enhances crop quantity and quality under challenging environmental conditions. Lastly, this survey of current and developing technologies indicates that integrative modeling and advanced scientific visualization may help overcome challenges in agriculture and nutrition data as large-scale and multidimensional data become available in these fields.

3D crop modeling↗

Quantifying Impacts of Biomass Pelletization on Fast Pyrolysis Using a Single-Particle Reactor, X-ray Computed Tomography, and Computational Modeling

The pore structure and density of lignocellulosic feedstocks dictate intraparticle transport phenomena and thereby play an important role in thermochemical conversion processes such as fast pyrolysis for biofuel and biochemical production. Variations in microstructure are inherent from different biomass species and can be introduced by preprocessing techniques such as cutting and pelletization. Morphological changes also occur during conversion and lead to vastly different pore structures and behavior during pyrolysis, which impact required conversion times and product distributions. The current work presents a comprehensive comparison of fast pyrolysis of neat and pelletized pine feedstocks, which includes single-particle experiments, modeling, and 3D imaging by X-ray computed tomography (XCT). The particle-scale model included anisotropic heat and mass transport in a shrinking particle with pyrolysis reactions based on the CRECK mechanism with boundary conditions informed by reactor-scale simulations of the single-particle reactor. The models were validated by measurements of the temperature and mass loss from single-particle pyrolysis experiments of neat and pelletized pine. Quantitative analysis of XCT geometries revealed that pyrolytic conversion yielded chars with increased porosity and permeability compared to the unpyrolyzed materials, along with decreased tortuosity and anisotropy. Pelletization of the pine feedstock resulted in a much denser, less permeable material, which converted slower and produced more residual char after pyrolysis compared to neat pine. The results from particle modeling revealed that accounting for the dynamic and anisotropic heat and mass transport caused by differences in pore structure is critical to achieving agreement with experimental results. Overall, this study highlights the dramatic differences in conversion behavior imparted by pelletization and the importance of capturing microstructural attributes in computational models to guide the design and optimization of pyrolysis processes for specific biomass feedstocks.

09 BIOMASS FUELS↗

In Situ X-ray Radiography and Computational Modeling to Predict Grain Morphology in $\beta$-Titanium during Simulated Additive Manufacturing

The continued development of metal additive manufacturing (AM) has expanded the engineering metallic alloys for which these processes may be applied, including beta-titanium alloys with desirable strength-to-density ratios. To understand the response of beta-titanium alloys to AM processing, solidification and microstructure evolution needs to be investigated. In particular, thermal gradients (Gs) and solidification velocities (Vs) experienced during AM are needed to link processing to microstructure development, including the columnar-to-equiaxed transition (CET). In this work, in situ synchrotron X-ray radiography of the beta-titanium alloy Ti-10V-2Fe-3Al (wt.%) (Ti-1023) during simulated laser-powder bed fusion (L-PBF) was performed at the Advanced Photon Source at Argonne National Laboratory, allowing for direct determination of Vs. Two different computational modeling tools, SYSWELD and FLOW-3D, were utilized to investigate the solidification conditions of spot and raster melt scenarios. The predicted Vs obtained from both pieces of computational software exhibited good agreement with those obtained from in situ synchrotron X-ray radiography measurements. The model that accounted for fluid flow also showed the ability to predict trends unobservable in the in situ synchrotron X-ray radiography, but are known to occur during rapid solidification. A CET model for Ti-1023 was also developed using the Kurz–Giovanola–Trivedi model, which allowed modeled Gs and Vs to be compared in the context of predicted grain morphologies. Both pieces of software were in agreement for morphology predictions of spot-melts, but drastically differed for raster predictions. The discrepancy is attributable to the difference in accounting for fluid flow, resulting in magnitude-different values of Gs for similar Vs.

36 MATERIALS SCIENCE↗

Understanding climate-induced migration through computational modeling: A critical overview with guidance for future efforts

Climate change has the potential to displace large populations in many parts of the developed and developing world. Understanding why, how, and when environmental migrants decide to move is critical to successful strategic planning within organizations tasked with helping the affected groups, and mitigating their systemic impacts. One way to support planning is through the employment of computational modeling techniques. Models can provide a window into possible futures, allowing planners and decision makers to test different scenarios in order to understand what might happen. While modeling is a powerful tool, it presents both opportunities and challenges. This paper builds a foundation for the broader community of model consumers and developers by: providing an overview of pertinent climate-induced migration research, describing some different types of models and how to select the most relevant one(s), highlighting three perspectives on obtaining data to use in said model(s), and the consequences associated with each. It concludes with two case studies based on recent research that illustrate what can happen when ambitious modeling efforts are undertaken without sufficient planning, oversight, and interdisciplinary collaboration. Lastly, we hope that the broader community can learn from our experiences and apply this knowledge to their own modeling research efforts.

54 ENVIRONMENTAL SCIENCES↗

Computer modeling of electrical and thermal performance during bipolar pulsed radiofrequency for pain relief

Purpose: Pulsed RF (PRF) is a nonablative technique for treating neuropathic pain. Bipolar PRF application is currently aimed at creating a “strip lesion” to connect the electrode tips; however, the electrical and thermal performance during bipolar PRF is currently unknown. The objective of this paper was to study the temperature and electric field distributions during bipolar PRF. Methods: The authors developed computer models to study temperature and electric field distributions during bipolar PRF and to assess the possible ablative thermal effect caused by the accumulated temperature spikes, along with any possible electroporation effects caused by the electrical field. The authors also modeled the bipolar ablative mode, known as bipolar Continuous Radiofrequency (CRF), in order to compare both techniques. Results: There were important differences between CRF and PRF in terms of electrical and thermal performance. In bipolar CRF: (1) the initial temperature of the tissue impacts on temperature progress and hence on the thermal lesion dimension; and (2) at 37 °C, 6-min of bipolar CRF creates a strip thermal lesion between the electrodes when these are separated by a distance of up to 20 mm. In bipolar PRF: (1) an interelectrode distance shorter than 5 mm produces thermal damage (i.e., ablative effect) in the intervening tissue after 6 min of bipolar RF; and (2) the possible electroporation effect (electric fields higher than 150 kV m{sup −1}) would be exclusively circumscribed to a very small zone of tissue around the electrode tip. Conclusions: The results suggest that (1) the clinical parameters considered to be suitable for bipolar CRF should not necessarily be considered valid for bipolar PRF, and vice versa; and (2) the ablative effect of the CRF mode is mainly due to its much greater level of delivered energy than is the case in PRF, and therefore at same applied energy levels, CRF, and PRF are expected to result in same outcomes in terms of thermal damage zone dimension.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Computational model development and validation of fuel dispersal phenomena

As the US nuclear industry is renewing efforts to extend the rod average burnup limits in existing pressurized water reactors, it becomes imperative to investigate the safety implications associated with high burnup fuel fragmentation and the potential dispersal of fragmented fuel into the reactor coolant system. Additionally, as there is a growing interest in transitioning to new nuclear fuel designs to accommodate higher fuel burnup levels, it becomes essential to conduct studies on fuel behavior during design basis accident scenarios, particularly during a postulated loss-of-coolant accident, in order to ensure the safe operation of light-water reactors. The present study develops a computational model to simulate the complex three-phase flow of the fuel dispersal phenomena that occur following a breach in the cladding, when the stream of solid fuel particles and carrier fission gas stream may interact with the surrounding liquid medium. Furthermore, the proposed framework considers the solid phase as discrete Lagrangian particles while the interpenetrating gas–liquid continuum is modeled using a Eulerian framework. The developed simulation framework is validated with experimental results performed in a separate-effect test facility by comparing predicted particle settlement locations with the experimental observations reported in the literature. Key highlights of this study include the interpretation of the high-pressure boundary conditions and the three-phase flow coupling strategy. Also discussed are the gas–liquid hydrodynamics and the behavior of particles as they are transported through the gas and liquid phases.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

A computational model for doctoring fluid films in gravure printing

The wiping, or doctoring, process in gravure printing presents a fundamental barrier to resolving the micron-sized features desired in printed electronics applications. This barrier starts with the residual fluid film left behind after wiping, and its importance grows as feature sizes are reduced, especially as the feature size approaches the thickness of the residual fluid film. In this work, various mechanical complexities are considered in a computational model developed to predict the residual fluid film thickness. Lubrication models alone are inadequate, and deformation of the doctor blade body together with elastohydrodynamic lubrication must be considered to make the model predictive of experimental trends. Moreover, model results demonstrate that the particular form of the wetted region of the blade has a significant impact on the model's ability to reproduce experimental measurements.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

High-Performance Computational Modeling of Strong Ground Motion for Seismic Response Analysis with Implications for High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSSNLV-015-19, Year 2

Presentation to be given at the FY 2020 SDRD Annual Meeting, to be held Sept. 23-24, on WebEx. There is one video also attached. This presentation is a reduced slide deck of the recently STIPped version, DOE/NV/03624--0875, presented 9-15-20 at the work in progress seminar. High-Performance Computational Modeling of Strong Ground Motion for Seismic Response Analysis with Implications for High-Hazard and/or Nuclear Facilities and Critical Infrastructure at the NNSSNLV-015-19, Year 2

58 GEOSCIENCES↗

Optimal size of the block in block GMRES on GPUs: computational model and experiments

The block version of GMRES (BGMRES) is most advantageous over the single right hand side (RHS) counterpart when the cost of communication is high while the cost of floating point operations is not. This is the particular case on modern graphics processing units (GPUs), while it is generally not the case on traditional central processing units (CPUs). Here, in this paper, experiments on both GPUs and CPUs are shown that compare the performance of BGMRES against GMRES as the number of RHS increases, with a particular focus on GPU performance. The experiments indicate that there are many cases in which BGMRES is slower than GMRES on CPUs, but faster on GPUs. Furthermore, when varying the number of RHS on the GPU, there is an optimal number of RHS where BGMRES is clearly most advantageous over GMRES. A computational model for the GPU is developed using hardware specific parameters, providing insight towards how the qualitative behavior of BGMRES changes as the number of RHS increase, and this model also helps explain the phenomena observed in the experiments.

97 MATHEMATICS AND COMPUTING↗

Computational modeling of grain boundary segregation: A review

Nearly all metals, alloys, ceramics, and their associated composites are polycrystalline in nature, with grain boundaries that separate well-defined crystalline regions that influence materials properties. In all but the most pure elemental systems, intentional solutes or impurities are present and can segregate to, or less commonly away from, the grain boundaries, in turn influencing boundary behavior, their stability, and associated materials properties. In some cases, grain-boundary segregation can also trigger “phase-like” structural transitions that dramatically alter the essential nature of the boundary. With the development of advanced electron microscopy techniques, researchers can directly observe grain-boundary structures and segregation with atomic precision. Despite such spatial resolution, the underlying mechanisms governing grain-boundary segregation remain difficult to characterize. As a result, computational modeling techniques such as density functional theory, molecular dynamics, mesoscale phase-field, continuum defect theory, and others are important complementary tools to experimental observations for studying grain-boundary segregation behavior. In conclusion, these computational methods offer the ability to explore the underlying formation mechanisms of grain-boundary segregation, elucidate complex segregation behavior, and provide insights into solutions to effectively controlling microstructure.

36 MATERIALS SCIENCE↗

FY22 Progress on Computational Modeling of the Water-Based NSTF

This report summarizes the system-level modeling effort by Argonne National Laboratory (Argonne) of the Natural convection Shutdown heat removal Test Facility (NSTF) in FY22. As an extension of the effort from FY21, this year’s work focuses primarily on the two-phase modeling of the NSTF using RELAP5-3D, particularly with the inclusion of the cavity model. The results from simulations were used to compare against experimental data for benchmarking purposes of the RELAP5 deck. Additionally, RELAP5 was used as a predictive tool to guide planned test operations and identify expected system behaviors. In the first part of this report, details are provided of the cavity omitted model where heat flux is applied directly as a boundary condition to the risers. The general trend predicted by the RELAP5 model matches that from the experimental data when a single-phase natural circulation flow is first established, followed by an oscillatory two-phase period and finally a stable two-phase flow. However, the onset of oscillations is predicted early by the model due to the smaller thermal mixing region in the tank. However, by expanding the simulated thermal mixing region in the tank, the onset of oscillations predicted by the model is able to match that from the experiment. These oscillations where studied in depth and are deduced to be flashing-induced instability. The model was then modified to simulate an accident scenario case where a representative heat load based on the full-scale Framatome’s 625 MW t SC-HTGR was applied directly to the riser channels. The simulated initial and boundary conditions were identical to those performed experimentally, facilitating direct comparisons between the predicted and experimental results. It was determined that the results showed some discrepancies remain, likely due to the overprediction of vapor generation rate by the computer model. In the second part of this report, the cavity model is re-introduced where it is observed that the RELAP5 prediction is now able to capture the major trends of the observed flow commonly observed during two-phase conditions. However, the onset of oscillations is once again predicted early by the model, possibly caused by the underprediction of heat loss from the heater and cavity. This is likely due to the omission of support structures in the cavity that can act as additional pathways for heat to escape to the environment. To overcome the underprediction of heat loss, part of the insulation surrounding the cavity side panels and the back of the heaters are removed to allow heat to escape directly to the environment, which then improves the RELAP5 prediction. Parametric studies are also performed to investigate the effects of heater power, tank inventory level, and tank gas space pressure on flow behaviors, also in direct comparison to conditions tested experimentally. User option-61 in the RELAP5-3D input deck, which changes the heat transfer coefficient correlations used for calculating the vapor generation, is also investigated where it is found that by enabling the option, the overall duration of oscillations is increased and matches that from the experiment better. The RELAP5 model is further benchmarked with a header inlet- throttling case where it is observed that the prediction from the model fails to capture some major features observed in the experiment. By using a modified loss coefficient curve for the valve, the accuracy of the prediction is improved where most of the major features observed in the experiment are predicted by the model. Lastly, the model is benchmarked with an inventory depletion scenario where it is observed that despite the modeling limitation of RELAP5, the prediction shows good agreement with the experimental data where major trends and features are captured by the model. Future work will see continued development of the current RELAP5-3D input deck of the NSTF to both improve the accuracy of the model’s predictive capability and continuing serving the experimental program. The mutually beneficial relationship between analysis and experimental efforts has become integral to the parent NSTF program, and the greater objective to fully understand and accurately predict the heat removal performance of a full scale RCCS concept.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Computational modeling study of the radial line slot antenna microwave plasma source with comparisons to experiments

The radial line slot antenna plasma source is a high-density microwave plasma source comprising a high electron temperature source region within the plasma skin depth from a coupling window and low electron temperature diffusion region far from the window. The plasma is typically comprised of inert gases like argon and mixtures of halogen or fluorocarbon gases for etching. Following the experimental study of Tian et al.[J. Vac. Sci. Technol. A 24, 1421 (2006)], a two-dimensional computational model is used to describe the essential features of the source. A high density argon plasma is described using the quasi-neutral approximation and coupled to a frequency-domain electromagnetic wave solver to describe the plasma-microwave interactions in the source. The plasma is described using a multispecies plasma chemistry mechanism developed specifically for microwave excitation conditions. The plasma is nonlocal by nature with locations of peak power deposition and peak plasma density being very different. The spatial distribution of microwave power coupling depends on whether the plasma is under- or over-dense and is described well by the model. The model predicts the experimentally observed low-order diffusion mode radial plasma profiles. The trends of spatial profiles of electron density and electron temperature over a wide range of power and pressure conditions compare well with experimental results.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NANO.PTML model for read-across prediction of nanosystems in neurosciences. computational model and experimental case of study

Abstract Neurodegenerative diseases involve progressive neuronal death. Traditional treatments often struggle due to solubility, bioavailability, and crossing the Blood-Brain Barrier (BBB). Nanoparticles (NPs) in biomedical field are garnering growing attention as neurodegenerative disease drugs (NDDs) carrier to the central nervous system. Here, we introduced computational and experimental analysis. In the computational study, a specific IFPTML technique was used, which combined Information Fusion (IF) + Perturbation Theory (PT) + Machine Learning (ML) to select the most promising Nanoparticle Neuronal Disease Drug Delivery (N2D3) systems. For the application of IFPTML model in the nanoscience, NANO.PTML is used. IF-process was carried out between 4403 NDDs assays and 260 cytotoxicity NP assays conducting a dataset of 500,000 cases. The optimal IFPTML was the Decision Tree (DT) algorithm which shown satisfactory performance with specificity values of 96.4% and 96.2%, and sensitivity values of 79.3% and 75.7% in the training (375k/75%) and validation (125k/25%) set. Moreover, the DT model obtained Area Under Receiver Operating Characteristic (AUROC) scores of 0.97 and 0.96 in the training and validation series, highlighting its effectiveness in classification tasks. In the experimental part, two samples of NPs (Fe 3 O 4 _A and Fe 3 O 4 _B) were synthesized by thermal decomposition of an iron(III) oleate (FeOl) precursor and structurally characterized by different methods. Additionally, in order to make the as-synthesized hydrophobic NPs (Fe 3 O 4 _A and Fe 3 O 4 _B) soluble in water the amphiphilic CTAB (Cetyl Trimethyl Ammonium Bromide) molecule was employed. Therefore, to conduct a study with a wider range of NP system variants, an experimental illustrative simulation experiment was performed using the IFPTML-DT model. For this, a set of 500,000 prediction dataset was created. The outcome of this experiment highlighted certain NANO.PTML systems as promising candidates for further investigation. The NANO.PTML approach holds potential to accelerate experimental investigations and offer initial insights into various NP and NDDs compounds, serving as an efficient alternative to time-consuming trial-and-error procedures.

60 APPLIED LIFE SCIENCES↗

Introduction to Focus Issue: Rhythms and Dynamic Transitions in Neurological Disease: Modeling, Computation, and Experiment

Rhythmic neuronal oscillations across a broad range of frequencies, as well as spatiotemporal phenomena, such as waves and bumps, have been observed in various areas of the brain and proposed as critical to brain function. While there is a long and distinguished history of studying rhythms in nerve cells and neuronal networks in healthy organisms, the association and analysis of rhythms to diseases are more recent developments. Indeed, it is now thought that certain aspects of diseases of the nervous system, such as epilepsy, schizophrenia, Parkinson's, and sleep disorders, are associated with transitions or disruptions of neurological rhythms. This focus issue brings together articles presenting modeling, computational, analytical, and experimental perspectives about rhythms and dynamic transitions between them that are associated to various diseases.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Computational modeling of the N-terminus of the human dopamine transporter and its interaction with PIP 2 -containing membranes: Modeling of DAT N-Terminus/PIP 2 Interactions

The dopamine transporter (DAT) is a transmembrane protein belonging to the family of neurotransmitter:sodium symporters (NSS). Members of the NSS are responsible for the clearance of neurotransmitters from the synaptic cleft, and for their translocation back into the presynaptic nerve terminal. The DAT contains long intracellular N- and C-terminal domains that are strongly implicated in the transporter function. The N-terminus (N-term), in particular, regulates the reverse transport (efflux) of the substrate through DAT. Currently, the molecular mechanisms of the efflux remain elusive in large part due to lack of structural information on the N-terminal segment. Here we report a computational model of the N-term of the human DAT (hDAT), obtained through an ab initio structure prediction, in combination with extensive atomistic molecular dynamics (MD) simulations in the context of a lipid membrane. Our analysis reveals that whereas the N-term is a highly dynamic domain, it contains secondary structure elements that remain stable in the long MD trajectories of interactions with the bilayer (totaling >2.2 μs). Combining MD simulations with continuum mean-field modeling we found that the N-term engages with lipid membranes through electrostatic interactions with the charged lipids PIP 2 (phosphatidylinositol 4,5-Biphosphate) or PS (phosphatidylserine) that are present in these bilayers. We identify specific motifs along the N-term implicated in such interactions and show that differential modes of N-term/membrane association result in differential positioning of the structured segments on the membrane surface. These results will inform future structure-based studies that will elucidate the mechanistic role of the N-term in DAT function.

59 BASIC BIOLOGICAL SCIENCES↗