Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “beam propagation method”

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 199 records · Page 11

A Novel Design of Guiding Stress Wave Propagation

Impact loads can induce a series of undesirable physical phenomena including vibration, acoustical shock, perforation, fracture and fragmentation, etc. The energy associated with the impact loads can lead to severe structure damage and human injuries. A design approach which effectively reduces these negative impacts through shock/stress wave diversion is highly needed. In this paper, a computational model which predicts stress wave propagation by considering different beam geometries and configurations is developed. Additionally, a novel concept of wave guide design which modifies the stress wave propagation path without disturbance is also presented. Furthermore, this design approach is not only useful for material property characterization particularly at intermediate or high strain rates, but also allows stress wave propagation in a desired direction as the shock/impact energy can be redistributed in controllable paths. The numerical results are experimentally verified through a Drop-Hopkinson bar apparatus at Sandia National Laboratories.

36 MATERIALS SCIENCE↗

Simulation of decelerating streamers in inhomogeneous atmosphere with implications for runaway electron generation

The dynamics of positive and negative streamers is numerically simulated in atmospheric pressure air in the range of parameters corresponding to the streamer deceleration and termination in the middle of the discharge gap. A detailed comparison with experiments in air at constant and variable density demonstrates good agreement between the 2D simulation results and the observations. It is shown that positive and negative streamers behave in radically different ways when decelerating and stopping. When the head potential drops, the negative streamer transits to the mode in which the propagation is due to the forward electron drift. In this case, the radius of the ionization wave front increases, whereas the electric field at the streamer head decreases further and the streamer stops. Its head diameter continues to increase due to the slow drift of free electrons in the residual under-breakdown field. On the contrary, the only advancement mechanism for a positive streamer with a decreasing head potential is a decrease in the effective radius of the ionization wave, leading to a local increase in the electric field. This mechanism makes it possible to compensate for the reduction in the efficiency of gas photoionization at small head diameters. A qualitative 1D model is suggested to describe streamer deceleration and stopping for different discharge polarities. Estimates show that, during positive streamer stopping, the local electric field at the streamer head can exceed the threshold corresponding to the transition of electrons to the runaway mode when the head potential (relative to the surrounding space) decreases to~1.2 kV in atmospheric pressure air. In this case, pulsed generation of a beam of runaway electrons directed into the channel of a stopping positive streamer can occur. The energy of the formed pulsed electron beam depends on the intensity of photoionization in front of the streamer head. This energy can vary from 700 V (when increasing the photoionization rate by a factor of 10 with respect to the value in atmospheric pressure air) to 2.6 kV (when decreasing the photoionization rate by a factor of 1000). It is possible that this behavior of decelerating positive streamers can explain the observed bursts of x-ray radiation during the streamer propagation in long air gaps.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Analysis of Nd3+:glass, solar-pumped, high-powr laser systems

The operating characteristics of Nd(3+):glass lasers energized by a solar concentrator were analyzed for the hosts YAG, silicate glass, and phosphate glass. The modeling is based on the slab zigzag laser geometry and assumes that chemical hardening methods for glass are successful in increasing glass hardness by a factor of 4. On this basis, it was found that a realistic 1-MW solar-pumped laser might be constructed from phosphate glass 4 sq m in area and 2 mm thick. If YAG were the host medium, a 1-MW solar-pumped laser need only be 0.5 sq m in area and 0.5 cm thick, which is already possible. In addition, Nd(3+) doped glass fibers were found to be excellent solar-pumped laser candidates. The small diameter of fibers eliminates thermal stress problems, and if their diameter is kept small (10 microns), they propagate a Gaussian single mode which can be expanded and transmitted long distances in space. Fiber lasers could then be used for communications in space or could be bundled and the individual beams summed or phase-matched for high-power operation.

Zapata, L. E.↗

Precise control of flexible manipulators

The design and experimental testing of end point position controllers for a very flexible one link lightweight manipulator are summarized. The latest upgraded version of the experimental set up, and the basic differences between conventional joint angle feedback and end point position feedback are described. A general procedure for application of modern control methods to the problem is outlined. The relationship between weighting parameters and the bandwidth and control stiffness of the resulting end point position closed loop system is shown. It is found that joint rate angle feedback in addition to the primary end point position sensor is essential for adequate disturbance rejection capability of the closed loop system. The use of a low order multivariable compensator design computer code; called Sandy is documented. A solution to the problem of control mode switching between position sensor sets is outlined. The proof of concept for endpoint position feedback for a one link flexible manipulator was demonstrated. The bandwidth obtained with the experimental end point position controller is about twice as fast as the beam's first natural cantilevered frequency, and comes within a factor of four of the absolute physical speed limit imposed by the wave propagation time of the beam.

Cannon, R. H., Jr.↗

Effects of various cavity designs on the performance of a CO2 TEA laser with an unstable resonator

Unstable resonator modeling has been carried out for an injection-seeded CO2 transversely excited atmosphere (TEA) laser in the NOAA/ERL/Wave Propagation Laboratory (WPL) Doppler lidar to examine the effects of various cavity designs on the quality of the output beam. The results show the effects of an injection pinhole, electrode spacing, mirror tilt, and radial reflectivity function of the output coupler. The electrode spacing in this laser has negligible effect. The injection pinhole, however, produces complicated structures in the output patterns. If the pinhole is removed, the output pattern is much smoother, and the frequency jitter is smaller. Misalignment sensitivity is very closely related to the radial reflectivity function. The superparabolic function provides the highest coupling efficiency, largest beam size, and good collimation, but produces a slightly higher misalignment sensitivity compared with a parabolic function. The Gaussian function provides the lowest misalignment sensitivity, but it produces the smallest beam size and the largest beam divergence. Also, the coupling coefficient is 50 percent lower than the optimum value. Methods for using a flat diffraction grating in unstable resonators are also investigated. The best way is to use a flat grating/positive lens combination to replace the back concave mirror.

Zhao, Yanzeng↗

Observation of a Giant Goos-Hänchen Shift for Matter Waves

The Goos-Hänchen (GH) shift describes a phenomenon in which a specularly reflected beam is translated along the reflecting surface such that the incident and reflected rays no longer intersect at the surface. Using a neutron spin-echo technique and a specially designed magnetic multilayer mirror, we have measured the relative phase between the reflected up and down neutron spin states in total reflection. The relative GH shift calculated from this phase shows a strong resonant enhancement at a particular incident neutron wave vector, which is due to a waveguiding effect in one of the magnetic layers. Calculations based on the observed phase difference between the neutron states indicate a propagation distance along the waveguide layer of 0.65 mm for the spin-down state, which we identify with the magnitude of the giant GH shift. The existence of a physical GH shift is confirmed by the observation of neutron absorption in the waveguide layer. Here, we propose ways in which our experimental method may be exploited for neutron quantum-enhanced sensing of thin magnetic layers.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Development of a Combined Cohesive and Virtual Crack-Closure Technique Approach to Represent R-Curves

Resistance curve (R-curve) effects due to fiber bridging, crack migration, and other blunting mechanisms are common in composite laminates. These mechanisms can dramatically increase the resistance to crack propagation but can be challenging to simulate. Delamination with R-curves can be analyzed using cohesive zone modeling (CZM) or the virtual crack closure technique (VCCT). Large fracture process zones can be simulated with CZM, but they require highly refined meshes. Coarser meshes can be used with VCCT, but this method is only applicable to small fracture process zones with R-curves defined as functions of position. Therefore, a technique with the computational efficiency of VCCT and the natural ability of cohesive elements to represent large fracture process zones is desirable. An approach is proposed that starts as a CZM in which the cohesive traction separation law (TSL) is separated into high-strength (HS) and low-strength (LS) components. The HS part, which is responsible for the mesh requirements of the CZM analysis, is replaced by VCCT. The combined method is evaluated by analyzing the response of double cantilevered beam specimens. Two specimen layup configurations, [0/90/90/0] 3s and [0/90/90/0] 9s , are evaluated. The R-curve response from each configuration is determined and applied to the constitutive properties of models built with CZM and with VCCT. The results and computational efficiency of the CZM, VCCT, and combined approaches are compared. The results are indicative that a combined cohesive/VCCT approach can enable progressive failure analyses to retain the computational efficiency of VCCT with the ability of the cohesive elements to capture R-curve effects.

Virtual Crack-Closure Technique↗

Compensation of heat load deformations using adaptive optics for the ALS upgrade: a wave optics study

In this paper, a realistic wave optics simulation method has been developed to study how wavefront distortions originating from heat load deformations can be corrected using adaptive X-ray optics. Several planned soft X-ray and tender X-ray insertion-device beamlines in the Advanced Light Source upgrade rely on a common design principle. A flat, first mirror intercepts the white beam; vertical focusing is provided by a variable-line-space monochromator; and horizontal focusing comes from a single, pre-figured, adaptive mirror. A variety of scenarios to cope with thermal distortion in the first mirror are studied by finite-element analysis. The degradation of the intensity distribution at the focal plane is analyzed and the adaptive optics that correct it is modeled. The range of correctable wavefront errors across the operating range of the beamlines is reported in terms of mirror curvature and spatial frequencies. The software developed is a one-dimensional wavefront propagation package made available in the OASYS suite, an adaptable, customizable and efficient beamline modeling platform.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Efficient end-to-end simulation of time-dependent coherent X-ray scattering experiments

Physical optics simulations for beamlines and experiments allow users to test experiment feasibility and optimize beamline settings ahead of beam time in order to optimize valuable beam time at synchrotron light sources like NSLS-II. Further, such simulations also help to develop and test experimental data processing methods and software in advance. The Synchrotron Radiation Workshop ( SRW ) software package supports such complex simulations. We demonstrate how recent developments in SRW significantly improve the efficiency of physical optics simulations, such as end-to-end simulations of time-dependent X-ray photon correlation spectroscopy experiments with partially coherent undulator radiation (UR). The molecular dynamics simulation code LAMMPS was chosen to model the sample: a solution of silica nanoparticles in water at room temperature. Real-space distributions of nanoparticles produced by LAMMPS were imported into SRW and used to simulate scattering patterns of partially coherent hard X-ray UR from such a sample at the detector. The partially coherent UR illuminating the sample can be represented by a set of orthogonal coherent modes obtained by simulation of emission and propagation of this radiation through the coherent hard X-ray (CHX) scattering beamline followed by a coherent-mode decomposition. GPU acceleration is added for several key functions of SRW used in propagation from sample to detector, further improving the speed of the calculations. The accuracy of this simulation is benchmarked by comparison with experimental data.

36 MATERIALS SCIENCE↗

Grain boundary zirconia-modified garnet solid-state electrolyte

Here, we report a method for promoting electrochemical stability in garnet Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 solid-state electrolyte based on a composite two-phase oxide–oxide microstructure. Grain boundary precipitation of the controlled distribution of amorphous zirconium oxide microparticles is achieved through the addition of reactive tantalum carbide. During ambient-atmosphere sintering, the carbide decomposes through an in situ reaction, the ‘extra’ Ta substituting for Zr within the Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 lattice. Density functional theory (DFT) calculations identify a thermodynamically favourable reaction path and show how substituting Ta 5+ at Zr 4+ sites affects the crystal structure as well as bulk ionic and electronic conductivities. Quantitative stereology highlights that zirconia also acts as a sintering aid, reducing compact porosity. Cryogenic focused-ion-beam scanning electron microscopy and fractography analysis of cycled solid-state electrolytes illustrates that near-universally observed intergranular Li-metal dendrite propagation is suppressed by the two-phase microstructure, favouring transgranular dendrites instead. Importantly, DFT demonstrates that compared with the Li 6.4 La 3 Zr 1.4 Ta 0.6 O 12 surface, the zirconium oxide surface per se is less electronically conductive and does not trap excess electrons to reduce Li ions. This is a key reason for the substantial improvement in the electrochemical properties over the single-phase baseline.

36 MATERIALS SCIENCE↗

Radar Differential Phase Signatures of Ice Orientation for the Prediction of Lightning Initiation and Cessation

The majority of lightning-related casualties typically occur during thunderstorm initiation (e.g., first flash) or dissipation (e.g., last flash). The physics of electrification and lightning production during thunderstorm initiation is fairly well understood. As such, the literature includes a number of studies presenting various radar techniques (using reflectivity and, if available, other dual-polarimetric parameters) for the anticipation of initial electrification and first lightning flash. These radar techniques have shown considerable skill at forecasting first flash. On the other hand, electrical processes and lightning production during thunderstorm dissipation are not nearly as well understood and few, if any, successful techniques have been developed to anticipate the last flash and subsequent cessation of lightning. One promising approach involves the use of dual-polarimetric radar variables to infer the presence of oriented ice crystals in lightning producing storms. In the absence of strong vertical electric fields, ice crystals fall with their largest (semi-major) axis in the horizontal associated with gravitational and aerodynamic forces. In thunderstorms, strong vertical electric fields (100-200 kV m(sup -1)) have been shown to orient small (less than 2 mm) ice crystals such that their semi-major axis is vertical (or nearly vertical). After a lightning flash, the electric field is typically relaxed and prior radar research suggests that ice crystals rapidly resume their preferred horizontal orientation. In active thunderstorms, the vertical electric field quickly recovers and the ice crystals repeat this cycle of orientation for each nearby flash. This change in ice crystal orientation from primarily horizontal to vertical during the development of strong vertical electric fields prior to a lightning flash forms the physical basis for anticipating lightning initiation and, potentially, cessation. Research has shown that radar reflectivity (Z) and other co-polar back-scattering radar measurements like differential reflectivity (Z(sub dr)) typically measured by operational dual-polarimetric radars are not sensitive to these changes in ice crystal orientation. However, prior research has demonstrated that oriented ice crystals cause significant propagation effects that can be routinely measured by most dual-polarimetric radars from X-band (3 cm) to S-band (10 cm) wavelengths using the differential propagation phase shift (often just called differential phase, phi(sub dp)) or its range derivative, the specific differential phase (K(sub dp)). Advantages of the differential phase include independence from absolute or relative power calibration, attenuation, differential attenuation and relative insensitivity to ground clutter and partial beam occultation effects (as long as the signal remains above noise). In research mode, these sorts of techniques have been used to anticipate initial cloud electrification, lightning initiation, and cessation. In this study, we develop a simplified model of ice crystal size, shape, orientation, dielectric, and associated radar scattering and propagation effects in order to simulate various idealized scenarios of ice crystals responding to a hypothetical electric field and their dual-polarimetric radar signatures leading up to lightning initiation and particularly cessation. The sensitivity of the K(sub dp) ice orientation signature to various ice properties and radar wavelength will be explored. Since K(sub dp) is proportional to frequency in the Rayleigh- Gans scattering regime, the ice orientation signatures should be more obvious at higher (lower) frequencies (wavelengths). As a result, simulations at radar wavelengths from 10 cm down to 1 cm (Ka-band) will be conducted. Resonance effects will be considered using the T-matrix method. Since most K(sub dp) Vbased observations have been shown at S-band, we will present ice orientation signatures from C-band (UAH/NASA ARMOR) and X-bd (UAH MAX) dual-polarimetric radars located in Northern Alabama. Issues related to optimal radar scanning for the detection of oriented ice will be discussed. Preliminary suggestions on how these differential phase signatures of oriented ice could contribute to lightning initiation and cessation algorithms will be presented.

Carey, L.D.↗

3013 Baseline Inner Container Report: Description, Imaging, and Analysis

This report presents consolidated findings from three microscopy studies of baseline Bagless Transfer Containers (BTCs), which are the inner containers used in DOE Standard 3013-compliant plutonium storage packages. The studies were conducted by Los Alamos National Laboratory (LANL), Savannah River National Laboratory (SRNL), and DNV, a subcontractor specializing in corrosion science. The baseline BTCs examined had not been exposed to corrosive environments, serving as control samples for comparison with destructively examined (DE) containers. Additionally, the Background section provides an overview of the BTC design, including manufacturing methods, material properties, and container configurations that may influence corrosion behavior. The primary objective of the LANL and SRNL studies was to establish reference characteristics of BTCs to distinguish manufacturing artifacts from corrosion-induced features observed in DE containers. Both laboratories used advanced microscopy techniques, including Wide Area Microscopy (WAMS), Laser Confocal Microscopy (LCM), Scanning Electron Microscopy (SEM), and Focused Ion Beam (FIB), to reveal that most surface anomalies were shallow and attributable to flow-forming processes, not corrosion. Subsurface impurities, identified as non-metallic inclusions, were also observed. DNV’s investigation focused on characterizing the mechanical and chemical properties of the BTC material to support crack growth modeling. Their findings indicate higher hardness and strength but reduced ductility—consistent with flow-formed 304L stainless steel and a higher-than-expected martensite content, in contrast to manufacturer claims. Ongoing and future studies, including Electron Backscatter Diffraction (EBSD), aim to further understand the microstructural factors influencing crack initiation and propagation. The combined studies provide a crucial baseline for distinguishing corrosion effects in DE containers, contributing to the safe and reliable long-term storage of plutonium-bearing materials.

36 MATERIALS SCIENCE↗

Probing drift velocity dispersion in MAPbI 3 photovoltaic cells with nonlinear photocurrent spectroscopy

Conventional time-of-flight (TOF) measurements yield charge carrier mobilities in photovoltaic cells with time resolution limited by the RC time constant of the device, which is on the order of 0.1–1 µs for the systems targeted in the present work. We have recently developed an alternate TOF method, termed nonlinear photocurrent spectroscopy (NLPC), in which carrier drift velocities are determined with picosecond time resolution by applying a pair of laser pulses to a device with an experimentally controlled delay time. In this technique, carriers photoexcited by the first laser pulse are “probed” by way of recombination processes involving carriers associated with the second laser pulse. Here, we report NLPC measurements conducted with a simplified experimental apparatus in which synchronized 40 ps diode lasers enable delay times up to 100 µs at 5 kHz repetition rates. Carrier mobilities of ∼0.025 cm 2 /V/s are determined for MAPbI 3 photovoltaic cells with active layer thicknesses of 240 and 460 nm using this instrument. Here, our experiments and model calculations suggest that the nonlinear response of the photocurrent weakens as the carrier densities photoexcited by the first laser pulse trap and broaden while traversing the active layer of a device. Based on this aspect of the signal generation mechanism, experiments conducted with co-propagating and counter-propagating laser beam geometries are leveraged to determine a 60 nm length scale of drift velocity dispersion in MAPbI 3 films. Contributions from localized states induced by thermal fluctuations are consistent with drift velocity dispersion on this length scale.

14 SOLAR ENERGY↗

GPR Robotic Assisted Non-Destructive Evaluation of Concrete Structures - 20196

Reinforced concrete structures are always under exposure of internal or external deteriorations due to many factors such as environmental, nuclear radiation or even construction deficiencies. The use of nondestructive test is one of the reliable methods to inspect structural members. In this case, GPR has a good capability to detect inside deteriorations of concrete structures. However, in many cases it is hard to inspect structural members by hand due to safety issues or inaccessibility. Thus, an alternative way to reduce such risks is to employ robot to carry nondestructive equipment. In this research, twelve reinforced concrete RC beams (8 x 16 x 96, width, height and length) are cast in the lab. Three RC beams have no defect and nine others various defects such as corrosion, de-bonded bars, void, delamination and honeycombing. The one-point load is applied at the center of the RC beam. The loading level is carried out at three levels: 50%, 75% and 100% of maximum load and total load-displacement diagram are performed. The load is held at each level and the robotic GPR is used to scan the side of RC beams. Then, RADAN software is employed to find out any defects at different loading conditions and visualize the crack propagation. The RADAN results show that the pattern of crack propagation for corrosion and debonding are same. Also, the load-displacement results present that inside defect cause reduction of load-capacity of RC beam. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

Electron Modulational Instability in the Strong Turbulent Regime for an Electron Beam Propagating in a Background Plasma

We study collective processes for an electron beam propagating through a background plasma using simulations and analytical theory. A new regime where the instability of a Langmuir wave packet can grow locally much faster than ion frequency is clearly identified. The key feature of this new regime is an electron modulational instability that rapidly creates a local Langmuir wave packet, which in its turn produces local charge separation and strong ion density perturbations because of the action of the ponderomotive force, such that the beam-plasma wave interaction stops being resonant. Three evolution stages of the process and observed periodic burst features are discussed. In conclusion, different physical regimes in the plasma and beam parameter space are demonstrated for the first time.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Electromagnetic Properties of Filled, Pixellated Arrays

Filled arrays of bolometers are currently being employed for use in astronomy from the far-infrared through millimeter parts of the electromagnetic spectrum. Because of the large range of wavelengths for which this is applicable, the number of modes supported by a bolometric pixel will vary according to specific application. Withington et al. (2003) have built a formalism for treating the electromagnetic properties of such bolometers by propagating the second order statistical properties of the radiation through a canonical optical system. In this work, we use this formalism to construct beam pattern images of square pixels for various ratios of p/lambda where p is the pixel size and lambda is the wavelength. In the low mode limit, the diffraction effects cause the beam pattern to be circular with a quadrupole dependence of Stokes Q and U. High mode cases approach the geometric limit. The polarization in these cases can be seen to trace the pixel edges. The effective size of the beam has a direct impact on the inter-pixel coupling and sets the number of independent detectors in an astronomical focal plane. This technique illustrates and quantifies the relationship between pixel size and angular resolution limits for a given wavelength and telescope. This is especially true in the limit of low p/lambda. In this case, the diffraction due to pixelization is non-negligible for the calculation of the overall angular resolution of the telescope. In addition, for instruments that are polarization sensitive, this method also provides a quantitative method for determining the contribution of the instrument to the measured polarization.

Chuss, David T.↗

Endothelium Preserving Microwave Treatment for Atherosclerosis

Method and apparatus are provided to treat atherosclerosis wherein the artery is partially closed by dilating the artery while preserving the vital and sensitive endothelial layer thereof Microwave energy having a frequency from 3 GHz to 300 GHz is propagated into the arterial wall to produce a desired temperature profile therein at tissue depths sufficient for thermally necrosing connective tissue and softening fatty and waxy plaque while limiting heating of surrounding tissues including the endothelial laser and/or other healthy tissue, organs, and blood. The heating period for raising the temperature a potentially desired amount, about 20 C., within the atherosclerotic lesion may be less than about one second. In one embodiment of the invention, a radically beveled waveguide antenna is used to deliver microwave energy at frequencies from 25 GHz or 30 GHz to about 300 GHz and is focused towards a particular radial sector of the artery. Because the atherosclerotic lesions are often asymmetrically disposed, directable of focussed heating preserves healthy sectors of the artery and applies energy to the asymmetrically positioned lesion faster than a non-directed beam. A computer simulation predicts isothermic temperature profiles for the given conditions and man be used in selecting power, pulse duration, beam width, and frequency of operation to maximize energy deposition and control heat rise within the atherosclerotic lesion without harming healthy tissues or the sensitive endothelium cells.

Carl, James R.↗

Endothelium Preserving Microwave Treatment for Atherosclerois

Method and apparatus are provided to treat atherosclerosis wherein the artery is partially closed by dilating the artery while preserving the vital and sensitive endothelial layer thereof. Microwave energy having a frequency from 3 GHz to 300 GHz is propagated into the arterial wall to produce a desired temperature profile therein at tissue depths sufficient for thermally necrosing connective tissue and softening fatty and waxy plaque while limiting heating of surrounding tissues including the endothelial layer and/or other healthy tissue, organs, and blood. The heating period for raising the temperature a potentially desired amount about 20 C, with the atherosclerotic lesion may be less than about one second. In one embodiment of the invention, a radically beveled waveguide antenna is used to deliver microwave energy at frequencies from 25 GHz or 30 GHz to about 300 GHz and is focused towards a particular radical sector of the artery. Because the atherosclerotic lesions are often asymmetrically disposed, directable or focussed heating preserves healthy sectors of the artery and applies energy to the asymmetrically positioned lesion faster than a non-directed beam. A computer simulation predicts isothermic temperature profiles for the given conditions and may be used in selecting power, pulse duration, beam width, and frequency of operation to maximize energy deposition and control heat rise within the atherosclerotic lesion without harming healthy tissues or the sensitive endothelium cells.

Carl, James R.↗