Search NASA⌕ Search

SEARCH · Search NASA

Results for “resolution matrix”

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 235 records · Page 13
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Single-Cell RNA-Seq Reveals Transcriptomic Heterogeneity and Post-Traumatic Osteoarthritis-Associated Early Molecular Changes in Mouse Articular Chondrocytes

Articular cartilage is a connective tissue lining the surfaces of synovial joints. When the cartilage severely wears down, it leads to osteoarthritis (OA), a debilitating disease that affects millions of people globally. The articular cartilage is composed of a dense extracellular matrix (ECM) with a sparse distribution of chondrocytes with varying morphology and potentially different functions. Elucidating the molecular and functional profiles of various chondrocyte subtypes and understanding the interplay between these chondrocyte subtypes and other cell types in the joint will greatly expand our understanding of joint biology and OA pathology. Although recent advances in high-throughput OMICS technologies have enabled molecular-level characterization of tissues and organs at an unprecedented resolution, thorough molecular profiling of articular chondrocytes has not yet been undertaken, which may be in part due to the technical difficulties in isolating chondrocytes from dense cartilage ECM. In this study, we profiled articular cartilage from healthy and injured mouse knee joints at a single-cell resolution and identified nine chondrocyte subtypes with distinct molecular profiles and injury-induced early molecular changes in these chondrocytes. We also compared mouse chondrocyte subpopulations to human chondrocytes and evaluated the extent of molecular similarity between mice and humans. This work expands our view of chondrocyte heterogeneity and rapid molecular changes in chondrocyte populations in response to joint trauma and highlights potential mechanisms that trigger cartilage degeneration.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Structural characterization of HIV-1 matrix mutants implicated in envelope incorporation

During the late phase of HIV-1 infection, viral Gag polyproteins are targeted to the plasma membrane (PM) for assembly. Gag localization at the PM is a prerequisite for the incorporation of the envelope protein (Env) into budding particles. Gag assembly and Env incorporation are mediated by the N-terminal myristoylated matrix (MA) domain of Gag. Nonconservative mutations in the trimer interface of MA (A45E, T70R, and L75G) were found to impair Env incorporation and infectivity, leading to the hypothesis that MA trimerization is an obligatory step for Env incorporation. Conversely, Env incorporation can be rescued by a compensatory mutation in the MA trimer interface (Q63R). The impact of these MA mutations on the structure and trimerization properties of MA is not known. In this study, we employed NMR spectroscopy, X-ray crystallography, and sedimentation techniques to characterize the structure and trimerization properties of HIV-1 MA A45E, Q63R, T70R, and L75G mutant proteins. NMR data revealed that these point mutations did not alter the overall structure and folding of MA but caused minor structural perturbations in the trimer interface. Analytical ultracentrifugation data indicated that mutations had a minimal effect on the MA monomer–trimer equilibrium. The high-resolution X-ray structure of the unmyristoylated MA Q63R protein revealed hydrogen bonding between the side chains of adjacent Arg-63 and Ser-67 on neighboring MA molecules, providing the first structural evidence for an additional intermolecular interaction in the trimer interface. These findings advance our knowledge of the interplay of MA trimerization and Env incorporation into HIV-1 particles.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2020

Gamma ray Imaging Based Methods to Improve the Accuracy of Uranium Holdup Quantification

Deposits in commercial uranium processing facilities of interest to nuclear safeguards are predominantly low enriched uranium (LEU). The shapes and sizes of holdup deposits can vary significantly, thus making it difficult to model the deposits accurately for calibration purposes. This presents a challenge to currently employed methods such as the Generalized Geometry Holdup (GGH) that rely on quantifying 235U mass by simplifying deposit shapes as a point, a line, or an area. In this work, gamma ray imaging using high energy resolution Germanium Gamma Imagers (GeGI™) is employed to determine the distribution of uranium inside the source containment. Two types of imaging methods are employed; coded-aperture imaging and Compton imaging. Gamma-ray emissions seen in different regions of the image are quantified using an inverse gamma-ray transport solver being developed at ORNL. Using measured data, the inverse method will be used to solve for unknown source parameters such as source matrix thickness, density, and attenuation due to container wall and shielding. The intrinsic efficiency of the detector is determined based on a library of built-in response functions and key dimensions of the detector (e.g., thickness and radius of the detector crystal). Measurements are performed in both imaging modalities using uranium sources of well-known masses and enrichments configured inside mocked-up holdup fixtures. Results from the two imaging modalities are intercompared. The paper presents preliminary results of the imaging measurements and the progress made in the computational methods to quantify uranium mass.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗
DOE OSTI2022

The Timepix4 analog front-end design: Lessons learnt on fundamental limits to noise and time resolution in highly segmented hybrid pixel detectors

This manuscript describes the optimization of the front-end readout electronics for high granularity hybrid pixel detectors. The theoretical study aims at minimizing the noise and jitter. The model presented here is validated with both circuit post layout simulations and measurements on the Timepix4 Application Specific Integrated Circuit (ASIC). The analog front-end circuit and the procedure to optimize the dimensions of the main transistors are described with detail. The Timepix4 is the most recent ASIC designed in the framework of the Medipix4 Collaboration. It was manufactured in 65nm CMOS process, and consists of a four side buttable matrix of 448 × 512 pixels with 55µm pitch. The analog front-end has a gain of ~36mV/ke - when configured in High Gain Mode, and ~20mV/ke - when configured in Low Gain Mode. The Equivalent Noise Charge (ENC) is ~68e - rms and ~80e - rms in High Gain Mode and in Low Gain Mode respectively. In event driven mode the incoming hits can be time stamped within a ~ 200ps time bin and the chip can deal with a maximum flux of ~ 3.6MHzmm –2 s –1 . In photon counting mode, the chip can deal with up to ~ 5GHzmm –2 s –1 . The routine designed to optimize the Timepix4 front-end is then used to analyze the performance limits in terms of jitter and noise for Charge Sensitive Amplifiers in pixel detectors.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2023

A mesoscale crystal plasticity model to predict room-temperature deformation and martensitic transformation of high-strength Quenching and Partitioning (Q&P) Steels and validation with synchrotron X-ray diffraction

Renowned for the superior mechanical properties and adeptness at cold-forming, Quenching and Partitioning (QP) steels have gained prominence as a promising candidate material in fabricating safety-critical components in various industries. The pertinent research on QP steels focus on the martensitic transformation of the Retained Austenite (RA) phase during cold-forming, a crucial mechanism that substantially influences the overall strength and ductility of QP steels. The austenite stability and transformation rate heavily rely on the local strain path and the initial microstructure, which is challenging for analytical prediction. In this paper, a mesoscale model is developed to capture the deformation and transformation kinetics of QP steels inside the microstructure. The model integrates the detailed explicit microstructure, acquired from characterization experiments, into a high-resolution finite element (FE) mesh. It distinctly model the deformation and interaction between the various phases and the effect on the transformation of RA. The model is validated with high energy X-ray diffraction (HEXRD) data, and shows excellent capability in predicting the asymmetric stress-strain behavior under uniaxial tension and compression, as well as the martensitic transformation rate. The model is used to investigate the strain and load partitioning effect of surrounding matrix to the transformation of RA, offering insights into the complex behavior of QP980 and facilitates further material development.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Toward Revealing Microcystin Distribution in Mouse Liver Tissue Using MALDI-MS Imaging

Cyanotoxins can be found in water and air during cyanobacterial harmful algal blooms (cHABs) in lakes and rivers. Therefore, it is very important to monitor their potential uptake by animals and humans as well as their health effects and distribution in affected organs. Herein, the distribution of hepatotoxic peptide microcystin-LR (MC-LR) is investigated in liver tissues of mice gavaged with this most common MC congener. Preliminary matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS) imaging experiments performed using a non-automated MALDI matrix deposition device and a MALDI-time-of-flight (TOF) mass spectrometer yielded ambiguous results in terms of MC-LR distribution in liver samples obtained from MC-LR-gavaged mice. The tissue preparation for MALDI-MS imaging was improved by using an automated sprayer for matrix deposition, and liver sections were imaged using an Nd:YAG MALDI laser coupled to a 15 Tesla Fourier-transform ion cyclotron resonance (FT-ICR)-mass spectrometer. MALDI-FT-ICR-MS imaging provided unambiguous detection of protonated MC-LR (calculated m/z 995.5560, z = +1) and the sodium adduct of MC-LR (m/z 1017.5380, z = +1) in liver sections from gavaged mice with great mass accuracy and ultra-high mass resolution. Since both covalently bound and free MC-LR can be found in liver of mice exposed to this toxin, the present results indicate that the distribution of free microcystins in tissue sections from affected organs, such as liver, can be monitored with high-resolution MALDI-MS imaging.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Synergic grain boundary segregation and precipitation in W- and W-Mo-containing high-entropy borides

The structures of W- and W-Mo-containing high-entropy borides (HEBs) are systematically studied by combining atomic-resolution transmission electron microscopy imaging, electron diffraction, and chemical analysis. Here, we reveal that W or W-Mo addition in HEBs leads to segregation of these elements to the grain boundaries (GBs). In the meantime, W- or W-Mo-rich precipitates also form along the GBs. Crystallographic analysis and atomic-scale imaging show that the GB precipitates in both W- and W-Mo-containing HEBs have a cube-on-cube orientation relationship with the matrix. With further strain analysis, the coherency of the precipitate/matrix interface is validated. Nanoindentation tests show that the simultaneous GB segregation and coherent precipitation, as a supplement to the grain hardening, provide additional hardening of the HEBs. Our work provides an in-depth understanding of the GB segregation and precipitation behaviors of HEBs. It suggests that GB engineering could be potentially used for optimizing the performance of high-entropy ceramics.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2021

Inversion of dynamical Bragg intensities to complex structure factors by iterated projections. For Ultramic. 2020. ("Pico" Festschrift, May 2021)

We discuss a method for recovering complex structure factors from many simultaneously excited Bragg beam intensities is described. The method is applied to simulated transmission electron diffraction data over a wide range of crystal thickness and beam energies. The method is based on iterated projections between structure and scattering matrices, which are related by a matrix unitary transformation, exponential, which we invert. The algorithm removes multiple-scattering perturbations from diffraction data and might be extended to other fields, including X-ray and neutron diffraction and cryo-electron microscopy. Because coherent multiple scattering involves interference between Bragg beams, the method also solves the phase problem. Unlike dynamical inversion from electron microscope images or ptychography data, the method, which starts with Bragg beam intensities, provides complex structure factors unaffected by focusing errors or resolution limitations imposed by lenses. We provide inversions from simulated data with 441 simultaneously excited Bragg beams over a range of thickness and beam energy. We discuss the retrieval of chirality information from enantiomorphs, the efficient incorporation of symmetry information using the irreducible representation of the group of structure matrices, and the effect of HOLZ lines to provide three-dimensional information.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A General Method for Estimating Zonal Transmission Interface Limits from Nodal Network Data: Preprint

Capacity expansion models for the electric power system often employ zonal (rather than nodal) resolution, necessitating estimates of aggregate power transfer limits across the interfaces between model zones. Interface limits between planning areas are sometimes published, but they are not generalizable to arbitrary zone shapes. There is thus a need for a reproducible method for estimating interface transfer limits (ITLs) between user-defined zones directly from nodal transmission system data. Here, we present a simple method for estimating ITLs using a DC power flow approximation via the power transfer distribution factor (PTDF) matrix. Linear optimization is performed to identify the distribution of power flows that maximizes the total flow on interface-crossing lines, subject to individual line ratings, limits on bus injection/withdrawal, and the relationships among flows, injections, and withdrawals imposed by the PTDF matrix. We demonstrate the application of the method on a 134-zone ~65000-bus system, and we explore the influence of flow direction, contingency level, and zone size on the estimated ITLs. There is significant heterogeneity in the ratio of the ITL to the sum of interface-crossing line ratings, which highlights the importance of accounting for the physical constraints on power flows imposed by Kirchhoff's laws when estimating zonal ITLs.

17 WIND ENERGY↗
DOE OSTI2020

Measurements of the Flow in the Vicinity of an Additively Manufactured Turbine Leading-Edge Using X-Ray Particle Tracking Velocimetry

X-ray particle tracking velocimetry (PTV) is performed, for the first time, to measure the velocity field inside a leading-edge of a turbine blade made by laser-additive-manufacturing (LAM) process. The traditional showerhead holes were replaced by a porous matrix in the leading-edge. The flow through such a leading-edge piece cannot be faithfully recreated by traditional prototype testing methods due to the surface roughness and imperfections caused by LAM process. Hence, direct measurement is the only option. However, it is difficult to measure flow inside such pieces with traditional velocimetry measurements due to the existence of metallic walls. Moreover, small internal size and high flow speeds call for a measurement technique with high spatial and temporal resolutions. To address these issues, we performed time-resolved X-ray PTV using the Advanced Photon Source (APS) synchrotron facility at the Argonne National Laboratory (ANL). A hydraulic system was constructed to run water, mixed with seeding particles, through the leading-edge piece. A high-speed camera captured the images of the seeding particles, which were later processed to create particle tracks. The time-averaged velocity field showed distinct pairs of vortices located in front of the porous outlet inside the leading-edge piece. The inlet channel showed reversed flow due to partial obstruction by the porous inlet of the test piece. Such knowledge of the flow field inside a leading-edge of a turbine blade will help us to design better cooling paths leading to higher cooling efficiency and increased life-span of a turbine blade.

42 ENGINEERING↗
DOE OSTI2024

A method for crystallographic mapping of an alpha-beta titanium alloy with nanometre resolution using scanning precession electron diffraction and open-source software libraries

An approach for the crystallographic mapping of two-phase alloys on the nanoscale using a combination of scanned precession electron diffraction and open-source python libraries is introduced in this paper. This method is demonstrated using the example of a two-phase α/β titanium alloy. The data were recorded using a direct electron detector to collect the patterns, and recently developed algorithms to perform automated indexing and analyse the crystallography from the results. Very high-quality mapping is achieved at a 3 nm step size. The results show the expected Burgers orientation relationships between the α laths and β matrix, as well as the expected misorientations between α laths. A minor issue was found that one area was affected by 180° ambiguities in indexing occur due to this area being aligned too close to a zone axis of the α with twofold projection symmetry (not present in 3D) in the zero-order Laue Zone, and this should be avoided in data acquisition in the future. Nevertheless, this study demonstrates a good workflow for the analysis of nanocrystalline two- or multi-phase materials, which will be of widespread use in analysing two-phase titanium and other systems and how they evolve as a function of thermomechanical treatments.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Mass Spectrometer Transient Analysis

This software implements a complete preprocessing pipeline for transient mass spectrometry (MS) data collected during TAP (Temporal Analysis of Products) experiments. It is designed to extract chemically meaningful fluxes from overlapping ion signals by applying a calibrated defragmentation matrix and solving the resulting linear system using non-negative least squares (NNLS) regression. The core script, preprocess_mass_spec.py, performs the following operations: Gain correction: Applies amplifier gain scalars derived from inert-packed calibration pulses to normalize signal intensities across AMUs and acquisition settings. Background subtraction: Removes experiment baselines using user-defined time windows, ensuring compatibility with slow-diffusing species and preventing negative values that would interfere with NNLS. Options to subtract before and after defragmentation. Defragmentation: Constructs a fragmentation matrix A from zeroth moments of calibration pulses (equal molar gas:inert mixtures) and solves Ax=b at each time point, where b is the raw MS signal and x is the estimated species flux. The matrix is normalized to inert signals and accounts for instrument-specific fragmentation behavior. Pulse-mode handling: Supports both averaged and individual pulse modes, enabling statistical treatment of fluxes and calculation of standard deviations. Integration and output: Computes zeroth moments (integrated fluxes) and exports time-resolved and integrated data in CSV format, suitable for downstream kinetic modeling. The software is validated using both virtual TAP simulations (VTAP) and experimental data from propane dehydrogenation (PDH) on CrOx/Al2O3 catalysts. It preserves temporal resolution by applying NNLS point-by-point across the pulse duration (typically 6,000+ time slices per pulse), leveraging the linear superposition principle to reconstruct full flux profiles. The defragmented outputs are compatible with kinetic extraction methods such as the G and Y procedures, which are used to derive rate–concentration relationships from TAP data. The details of these validations are discussed in detail in the supporting manuscript and supporting information. Example data and output files are also included. The methodology is robust to experimental noise and drift, with calibration protocols that account for pulse size effects, MS aging, and inert gas normalization. The software is modular, reproducible, and tailored for high-throughput TAP-MS workflows in catalysis research.

Kristy, Stephen [Idaho National Laboratory (INL), ↗
DOE OSTI2021

Asymptotically Compatible Reproducing Kernel Collocation and Meshfree Integration for Nonlocal Diffusion

Reproducing kernel (RK) approximations are meshfree methods that construct shape functions from sets of scattered data. We present an asymptotically compatible (AC) RK collocation method for nonlocal diffusion models with Dirichlet boundary condition. The numerical scheme is shown to be convergent to both nonlocal diffusion and its corresponding local limit as nonlocal interaction vanishes. The analysis is carried out on a special family of rectilinear Cartesian grids for a linear RK method with designed kernel support. The key idea for the stability of the RK collocation scheme is to compare the collocation scheme with the standard Galerkin scheme, which is stable. In addition, assembling the stiffness matrix of the nonlocal problem requires costly computational resources because high-order Gaussian quadrature is necessary to evaluate the integral. We thus provide a remedy to the problem by introducing a quasi-discrete nonlocal diffusion operator for which no numerical quadrature is further needed after applying the RK collocation scheme. The quasi-discrete nonlocal diffusion operator combined with RK collocation is shown to be convergent to the correct local diffusion problem by taking the limits of nonlocal interaction and spatial resolution simultaneously. The theoretical results are then validated with numerical experiments. We additionally illustrate a connection between the proposed technique and an existing optimization based approach based on generalized moving least squares.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

Direct observation on supersonic microprojectile penetration of carbon fiber composites with ultrafast synchrotron X-ray phase contrast imaging

High-speed penetration into carbon fiber composites is of fundamental importance to materials science and impact engineering, but research along this line suffers considerably from the lack of direct experimental observations. Here we investigate such penetration dynamics of a unidirectional carbon fiber reinforced epoxy (UCFRE) composite, with a combination of in situ, ultrafast, synchrotron phase contrast imaging and finite element (FE) analysis. The experiments yield the first direct observation on projectile trajectories and fiber-scale deformation and damage in the UCFRE composites during supersonic microprojectile penetration, for different fiber orientations (0 degrees - 90 degrees from the impact direction) and projectile velocities (600 - 850 m.s -1 ), at unprecedented temporal (~ 100 ps) and spatial (5 μm) resolutions. The maximum penetration depth decreases with increasing fiber orientation angles, as a result of anisotropic damage evolution in the composite sample. Strain localizations are prone to develop along a direction perpendicular to the fiber orientation, while the damage or cavity region, along the fiber direction. FE modeling with a three-dimensional Hashin criterion yields consistent projectile trajectory and cavity morphology with the experimental results. Finally, with increasing fiber orientation angles, damage analyses show a transition in the damage mode from fiber compression to matrix compression damage, in line with the increasing maximum penetration depth.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Evolution of Efimov states

The Efimov phenomenon manifests itself as an emergent discrete scaling symmetry in the quantum three-body problem. In the unitarity limit, it leads to an infinite tower of three-body bound states with energies forming a geometric sequence. Here in this work, we study the evolution of these so-called Efimov states using relativistic scattering theory. We identify them as poles of the three-particle S matrix and trace their trajectories in the complex energy plane as they evolve from virtual states through bound states to resonances. We dial the scattering parameters toward the unitarity limit and observe the emergence of the universal scaling of energies and couplings - a behavior known from the non-relativistic case. Interestingly, we find that Efimov resonances follow unusual, cyclic trajectories accumulating at the three-body threshold and then disappear at some values of the two-body scattering length. We propose a partial resolution to this ?missing states? problem.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Time-dependent approach to inelastic scattering spectroscopies in and away from equilibrium: Beyond perturbation theory

We propose a nonperturbative numerical approach to calculate the spectrum of a many-body Hamiltonian with time and momentum resolution by exactly recreating a scattering event using the time-dependent Schrödinger equation. Akin to an actual inelastic scattering experiment, we explicitly account for the incident and scattered particles (e.g., photons, neutrons, electrons, etc.) in the Hamiltonian and obtain the spectrum by measuring the energy and momentum lost by the particle after interacting with the sample. We illustrate the method by calculating the spin excitations of a Mott-insulating Hubbard chain after a sudden quench with the aid of the time-dependent density matrix renormalization group method. Our formalism can be applied to different forms of spectroscopies, such as neutron and Compton scattering and electron-energy-loss spectroscopy, as an example.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2024

3D-Reconstruction of Tau Neutrinos in LArTPC Detectors

The Deep Underground Neutrino Experiment (DUNE) is a next-generation neutrino experiment currently under construction. DUNE will consist of two high-resolution neutrino interaction imaging detectors exposed to the world’s most intense neutrino beam, with the Near Detector at Fermilab and the Far Detector 1,300 km away in the Sanford Underground Research Facility in South Dakota, US. The high statistics and excellent resolution capabilities of DUNE's $^{40}$Ar detector will allow us to make precision studies of oscillation parameters capable of searching for CP violation in the lepton sector, testing interaction models, and studying phenomena that have until now, seemed too complex to measure, like $\nu_\tau$ detection and therefore, providing the completion of the 3-flavor neutrino paradigm. Knowledge of the $\nu_\tau$ detection can impact a broad spectrum of open questions. These include searching for non-standard neutrino interactions, constraining the unitarity of the PMNS matrix, searching for sterile neutrinos, and studying neutrino interactions. In the case of LArTPC data, the detector hits can be considered nodes in a graph, and the edges represent the spatial and temporal relationships between them. By using graph neural networks, it is possible to exploit these relationships and improve the accuracy of particle identification and reconstruction. During my presentation and specifically for tau neutrino reconstruction, I will show the effectiveness and reliability of our in-house developed graph neural network (GNN), NuGraph. This GNN classifies detector hits based on the particle type responsible for their production, assuring that the system accurately identifies and categorizes information based on its unique characteristics.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2020

Transport Simulations on Scanning Transmission Electron Microscope Images of Nanoporous Shale

Digital rock physics is an often-mentioned approach to better understand and model transport processes occurring in tight nanoporous media including the organic and inorganic matrix of shale. Workflows integrating nanometer-scale image data and pore-scale simulations are relatively undeveloped, however. In this paper, a workflow is demonstrated progressing from sample acquisition and preparation, to image acquisition by Scanning Transmission Electron Microscopy (STEM) tomography, to volumetric reconstruction to pore-space discretization to numerical simulation of pore-scale transport. Key aspects of the workflow include (i) STEM tomography in high angle annular dark field (HAADF) mode to image three-dimensional pore networks in µm-sized samples with nanometer resolution and (ii) lattice Boltzmann method (LBM) simulations to describe gas flow in slip, transitional, and Knudsen diffusion regimes. It is shown that STEM tomography with nanoscale resolution yields excellent representation of the size and connectivity of organic nanopore networks. In turn, pore-scale simulation on such networks contributes to understanding of transport and storage properties of nanoporous shale. Interestingly, flow occurs primarily along pore networks with pore dimensions on the order of tens of nanometers. Smaller pores do not form percolating pathways in the sample volume imaged. Apparent gas permeability in the range of 10−19 to 10−16 m2 is computed.

42 ENGINEERING↗