Search NASA⌕ Search

SEARCH · Search NASA

Results for “STRESS MEASUREMENT, X-RAY”

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 91 records · Page 5

Uncovering hidden bias in neutron diffraction residual strain measurements

When calculating residual strain via neutron or X-ray diffraction, uncertainties propagated from the peak fit are often inadequate to describe the true scatter of measurements about a singular strain state, such as one that should describe a macroscopic continuum. Because diffraction is inherently a selective process, orientation-dependent scatter arises from the sub-sampling of strong microstructure and strain gradients. This paper investigates the appropriateness of propagated uncertainties with reference to their original intention, i.e. noise about a mean value. Thirty-six unique orientations of strain measurements are taken at multiple locations within an additive friction-stir deposition component with fine-scale gradients (∼200 µm) of plastic strain, texture and residual elastic strain. Multiple strain and stress calculation pathways are compared: direct substitution of three measurements into Hooke's law, direct inversion of any six unique orientations into the strain state tensor and thirty-six measurement least-squares estimation. For the last two cases, the appropriateness of the uncertainty interval is statistically evaluated on the basis of a physical constraint: common agreement under the strain transformation law. For this sample, the direct inversion of six measurements retains a conservative estimate of the uncertainty. However, propagated uncertainties in the least-squares solution greatly underestimate the true experimental scatter. A simple pathway to estimate appropriate uncertainty intervals is suggested. These results demonstrate that the interpretation of uncertainty in residual strain is strongly dependent on intrinsic sample-dependent effects, and that oversampling orientations and statistical analysis can give more accurate results with realistic uncertainties.

36 MATERIALS SCIENCE↗

Thermal Inspection for the Assessment of Adhesively Bonded Metal Adherends

Certification of adhesively bonded structures is currently a challenge for aircraft manufacturers. The ability to certify bonds in primary structure can reduce dependence on fasteners and this will enable more efficient manufacturing. For example, drilling holes and installing fasteners in bonded joints can be a potentially significant bottleneck in airframe manufacturing. In addition, fasteners increase the airframe weight and can add stress concentration areas. Reducing fastener count can accelerate manufacturing and improve aircraft performance. NASA is currently investigating nondestructive evaluation (NDE) techniques to assess bond integrity. For example, bond thickness influences bond strength and therefore NDE techniques are being investigated to determine bond thickness. For bonded metal structures, there is a large difference in the thermal diffusivity between the 7075 aluminum alloy metal adherends and aircraft grade adhesive. Multi-layered thermal models show a large variation in the thermal response for bondline adhesive thicknesses that vary from 100 to 300 microns for adherend thicknesses of 0.163 cm. Experimental through-transmission thermal measurements reveal promise to quantitatively characterize the bondline thickness. Results were validated with X-ray computed tomography and optical microscopy measurements, and influence of porosity on the thermal model and measurements are investigated.

X-ray computed tomography↗

Thermal Inspection for the Assessment of Adhesively Bonded Metal Adherents

Certification of adhesively bonded structures is currently a challenge for aircraft manufacturers. The ability to certify bonds in primary structure can reduce dependence on fasteners and this will enable more efficient manufacturing. For example, drilling holes and installing fasteners in bonded joints can be a potentially significant bottleneck in airframe manufacturing. In addition, fasteners increase the airframe weight and can add stress concentration areas. Reducing fastener count can accelerate manufacturing and improve aircraft performance. NASA is currently investigating nondestructive evaluation (NDE) techniques to assess bond integrity. For example, bond thickness influences bond strength and therefore NDE techniques are being investigated to determine bond thickness. For bonded metal structures, there is a large difference in the thermal diffusivity between the 7075 aluminum alloy metal adherends and aircraft grade adhesive. Multi-layered thermal models show a large variation in the thermal response for bondline adhesive thicknesses that vary from 100 to 300 microns for adherend thicknesses of 0.163 cm. Experimental through-transmission thermal measurements reveal promise to quantitatively characterize the bondline thickness. Results were validated with X-ray computed tomography and optical microscopy measurements, and influence of porosity on the thermal model and measurements are investigated.

bondline thickness↗

Prediction of residual stresses in additively manufactured parts using lumped capacitance and classical lamination theory

Several industries are interested in Laser Powder Bed Fusion (L-PBF) Additively Manufactured (AM) metal parts because their designs can be made arbitrarily complex while retaining bulk-type material properties. However, the residual stresses (RS) and distortions caused by the heat gradients inherent to L-PBF processes are detrimental to the structural integrity of the parts and must be taken into consideration during the part design cycle. Predicting the state of stresses in as-built 3D printed parts is a difficult problem that is typically approached with the use of transient thermomechanical Finite Element Models (FEMs). However, the nonlinearities associated with AM processes are difficult to capture in these FEMs without increasing the computational cost of the simulation, limiting their ability to be incorporated into practical design cycles. This work presents a novel analytical framework that combines lumped capacitance nonlinear heat transfer with time dependent classical lamination theory to efficiently and accurately predict RS in as-built L-PBF parts without the need of FEMs. The simulation was compared to Neutron Diffraction (ND) residual strain measurements taken at Oak Ridge National Laboratories (ORNL) as well as Synchrotron X-ray Diffraction (XRD) strain data published by the National Institute of Standards and Technology (NIST). The simulation predictions and the experimental data showed excellent agreement for the in-plane strain directions, and general agreement for the out of plane strain component, highlighting an area where further development can be implemented.

42 ENGINEERING↗

Artificial Soiling Replication of Field Losses on Commercial Photovoltaic Modules

Here, this study demonstrates the capabilities of an indoor artificial soiling approach developed to closely replicate the natural, cyclic soil accumulation processes in the field-dust suspension, deposition, and sedimentation/cementation-for a subtropical climate. In this work, a near-field environment is replicated in an artificial soiling cubic chamber through controlled regulation of humidity, temperature, dust type, and dust concentration, based on site-specific historical climate data. Two different models (MA and MB) of full-size commercial photovoltaic modules from a single manufacturer, installed side by side in the mid-Atlantic United States, were retrieved and subjected to artificial soiling experiments and various characterization measurements, including short-circuit current, colorimetry, reflectance, X-ray fluorescence, laser diffraction, and optical microscopy. Both in the field and in our improved field-representative artificial soiling tests, the MA model experienced roughly twice the soiling loss as the MB model. To closely replicate the field soiling losses for a site-specific climate, it is critical to include: 1) The use of field-collected dust with identical dust chemistry and particle distribution instead of standardized ISO 12103 Arizona Road dusts, 2) the use of only a small amount of field-collected dust inside the chamber during the deposition process (e.g., 0.15 g), and 3) the preconditioning of the surface coating for the partial/full dose of UV stress as experienced in the field during sunlight exposure and the abrasion as experienced in the field during regular module cleaning activities, if/as needed. The field-representative artificial soiling method developed here could potentially be adopted for rank ordering of various antisoiling coatings developed by researchers and industry stakeholders.

14 SOLAR ENERGY↗

ETP-0492, Measured Residual Stresses in CYL S/N 53 Fretted Area

This test report presents the results of a residual stress survey of the inner clevis leg of lightweight cylinder SIN 053 as described by ETP-0492. The intent of this testing was to evaluate the residual stresses that occur in and around the inner clevis leg at the capture feature contact zone during a normal flight cycle. Lightweight case cylinder segment IU50717, S/N L053 from Flight STS-27 exhibited fretting around the contact zone of the inner clevis leg and the capture feature of the field joint. Post flight inspection revealed several large fitting pits on the inside of the inner clevis leg. This cylinder was assigned for both residual stress and metallurgical evaluation. This report is concerned only with the residual so= evaluations. The effects of glass bead cleaning and fi=ing were evaluated using the x-ray diffraction method.

Webster, Ronald L.↗

Analysis of strain in ion implanted 4H-SiC by fringes observed in synchrotron X-ray topography

A novel high energy implantation system has been successfully developed to fabricate 4H-SiC superjunction devices for medium and high voltage via implantation of dopant atoms with multi-energies ranging from 13 to 66 MeV. The significantly higher levels of energy used compared to conventional implantation processes, necessitates detailed characterization of the lattice damage caused by implantation. To achieve this by employing the novel high energy system, 4H-SiC wafer with 12 μm epilayers were blanket implanted by 13.8–65.7 MeV Al atoms. The lattice damages induced by the implantation were primarily characterized by Synchrotron X-ray Plane Wave Topography (SXPWT) and Reciprocal Space Mapping (RSM). Topographs reveal fringe contrast akin to multiple asymmetric diffraction peaks with an angular separation of only 2″ (arcseconds) observed on rocking curves, indicating inhomogeneous strain distribution across the implanted layer. The strain profile of the implanted layer was extracted from the fringe contrast by applying Rocking-curve Analysis by Dynamical Simulation (RADS). In conclusion, the maximum strain value is similar to that measured on the RSM.

A1. Characterization↗

Anisotropic structure and optical engineering of strontium titanate and zirconia responding to sequential hydrogen and helium irradiation

Unraveling the correlation between strain engineering with anisotropic optical properties via nanoscale defects gradually evolved into a strategy for fundamental studies and technological applications, yet it remains understudied in functional complex oxides compared to metals and semiconductors. Here, the methodology of strain engineering for the determination of the lattice parameters, parallel/normal to the sample surface, in the individual layers of single-crystalline superlattices is derived, which is based on analysis of high-angle X-ray diffraction measurements in combination with diffraction reciprocal space mapping. With modeling that takes into account the effect of elastic properties, two elastically anisotropic materials, SrTiO 3 and ZrO 2 , have been compared in terms of defect-induced elastic strain caused by individual and sequential H + and He 2+ irradiation. The anisotropic lattice swelling with corresponding refractive index, and obstructive behavior of elastic strain recovery are demonstrated in SrTiO 3 , while approaching strain behaviors accompanied by isotropic refractive index distribution in both orientations are confirmed in ZrO 2 . Under sequential He 2+ /H + irradiation, pre-existing He-vacancy complexes acted as vacancy traps, preferentially capturing H + clusters to induce lattice distortion and enhance absorption. Nanohardness increments (ΔH) calculated by the DBH model matched nanoindentation results, confirming that sequential irradiation generated higher indentation yield stress than individual irradiations.

Anisotropic expansion↗

Statistical fracture behavior of doped UO 2 using a ball-on-ring equibiaxial flexure test method

Metal oxide dopants, such as titanium and chromium oxides, have garnered considerable attention for their potential to increase grain size (≥ 30 µm) in UO 2 fuel, purportedly enhancing fission gas retention during reactor operation. Fuel performance is significantly impacted by fuel fracture behavior, so it is important to understand the effects of enhanced grain size and dopant content on UO 2 fuel fracture. UO 2 pellets were doped with 0.1 wt% TiO 2 and 0.3 wt% Cr 2 O 3 to alter density and grain size. Inductively coupled plasma mass spectroscopy measured dopant levels pre- and post-sintering. X-ray diffraction revealed lattice changes and microstrain via Rietveld refinement. Field emission scanning electron microscopy determined grain sizes of approximately 30 µm for TiO 2 doping and 7 µm for Cr 2 O 3 doping. Transverse rupture strength tests were performed on over 30 samples per dataset to obtain characteristic strength and Weibull modulus. Results indicate no statistical difference in fracture strength between 0.1 wt% TiO 2 doped UO 2 and undoped UO 2 , while 0.3 wt% Cr 2 O 3 doped UO 2 exhibited a 20% decrease in fracture strength. Doped UO 2 samples also showed reduced Weibull modulus compared to undoped UO 2 , suggesting increased scatter in fracture strength. This study's findings suggest that titanium and chromium oxide doping in UO 2 , regardless of grain size, induce residual stresses, decreasing fracture strength and increasing variability in fracture behavior.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Quaternary i-MAX Phases (Mo 2/3 RE 1/3 ) 2 AlC (RE: Dy, Tb, Er): Experimental Characterization and First-Principles Insights into their Fundamental Properties

Rare earth (RE)-based materials have unique electronic, magnetic, and optical properties, leading to the recent discovery of atomically layered solids with the chemical formula (M' 2/3 RE 1/3 ) 2 AlC, which have since garnered significant attention in the scientific community. This study aims to synthesize, characterize, and investigate the structural and thermal stability of the RE i-MAX phases. We prepared i-MAX phases using molybdenum (Mo) as M′ and RE elements as Dy, Tb, and Er, namely (Mo 2/3 Dy 1/3 ) 2 AlC, (Mo 2/3 Tb 1/3 ) 2 AlC, and (Mo 2/3 Er 1/3 ) 2 AlC. Structural characterization through x-ray diffraction (XRD) and Raman spectroscopy confirms the formation of the RE-based i-MAX phase, along with the presence of minor impurity phases in the alloys. Thermogravimetric analysis (TGA) conducted up to 1000°C under ambient conditions reveals that the i-MAX phases remain thermally stable up to approximately 450°C, beyond which oxidation leads to a noticeable weight gain in all samples. Differential scanning calorimetry (DSC) measurements during heating and cooling cycles show endothermic and exothermic peaks for (Mo 2/3 Dy 1/3 ) 2 AlC i-MAX in the 410–420°C range, indicating a temperature-induced minor atomic arrangement. In contrast, these peaks are absent in the Tb- and Er-based i-MAX phases. These findings offer valuable insights into the thermal behavior and stability of these i-MAX phases under thermal stress, contributing to a deeper understanding of their unique properties. Furthermore, first-principles density functional theory (DFT) calculations were performed to investigate the electronic and optical properties of the i-MAX phases. The results reveal their metallic nature, with pronounced contributions from Mo and RE elements near the Fermi level and within the conduction band.

Rare earth↗

Compression–tension cell with sample manipulator for in situ X‐ray nanotomography experiments

In situ X-ray nanotomography experiments where tensile or compressive force is applied on the sample require specialized equipment. A compression-tension device with fluid flow-through capability has been designed for X-ray nanotomography beamlines. The compression-tension cell is equipped with a triaxial stage for sample alignment and a high sensitivity loadcell for measurement of applied force. To handle the <100 µm samples used for X-ray nanotomography imaging and for loading samples on the compression-tension cell a sample manipulator has been built. The sample manipulator is capable of selecting a single <100 µm particle for nanotomography scanning while viewing multiple samples under an optical microscope. To test the functionality of these two devices an initial compression experiment involving two glass beads was performed. To demonstrate instrument stability two spherical glass beads were compressed from a no load condition until one of the beads fractured. Nanotomography data were collected at each step of increasing compressive force. The experimentally observed contact area of the spherical glass beads was compared with the theoretical estimate using the Hertz analysis. To demonstrate the fluid flow capability, two calcite grains were compressed against each other under a calcite saturated solution. Surface topological changes were observed for the stressed grain contact area.

X-ray tomography↗

Using X-ray radiography to study oxygen flow in a proton exchange membrane electrolyzer operating under balanced pressure conditions

Of the various water electrolyzer technologies, the proton exchange membrane electrolyzer (PEMWE) is one of the best solutions for producing clean hydrogen without releasing CO 2 . In order to allow for widespread use of clean hydrogen, it is necessary to decrease its cost, which is intrinsically related to system operation. Current PEMWE plants operate in differential mode, directly pressurizing hydrogen and benefiting from thermodynamic compression, which increases overall system efficiency. However, high differential pressure above 30 bar can cause membrane stress, resulting in membrane creeping and failure. Pressurizing the water and operating at balanced pressure allows hydrogen to be produced at higher pressures while preserving the integrity of the membrane and porous layers. Nevertheless, the impact of pressurizing water on PEMWE performance must be better understood to maximize performance under balanced pressure conditions. Here, this study examined the impact of water pressure on two-phase flow. A high-pressure electrolyzer setup was developed to perform operando X-ray radiography and examine oxygen transport with high temporal resolution. The imaging segmentation process, developed to capture bubble properties in the channel, was applied to a specific experiment. The results clearly showed that as pressure increased up to 30 bars, the initial bubbly flow transitioned to slug flow, which led to channel saturation with oxygen. This work demonstrates that two-phase flow in an electrolyzer can be studied using X-ray radiography, which has the advantages of fast measurements and the ability to probe dense materials, such as those required for pressurized electrolyzers.

Balanced high pressure operation↗

Microstructural and rheological training and memory of nanocolloidal soft glasses under cyclic shear

An intrinsic feature of disordered and out-of-equilibrium materials, such as glasses, is the dependence of their properties on their history. An important example is rheological memory, in which disordered solids obtain properties based on their deformation history. Here, in this study, we employ x-ray photon correlation spectroscopy with in situ rheometry to characterize memory formation in a nanocolloidal soft glass due to cyclic shear. During a cycle, particles undergo irreversible displacements composed of a combination of shear-induced diffusion and heterogeneous, residual strain fields. At lower shear amplitudes, the displacements resemble a random walk in which the directions in each cycle are independent of those in preceding cycles, while at high amplitude, the irreversible displacements in consecutive cycles become correlated. The magnitudes of the displacements decrease with each cycle before reaching a steady state where the microstructure has been trained to achieve enhanced reversibility even at shear amplitudes well above yielding and despite the presence of thermal fluctuations. At amplitudes below and near yielding, these decreases are monotonic, while well above yielding, they are nonmonotonic, suggesting evidence of shear banding. Accompanying this microstructural training are corresponding decreases in the dissipation during each cycle and the magnitude of the residual stress toward steady-state values. Memory of the training is revealed by measurements in which the amplitude of the shear is changed after steady state is reached. The magnitude of the particle displacements, as well as the dissipation and the change in residual stress, vary nonmonotonically with the new shear amplitude, having minima near the training amplitude, thereby revealing correlated microscopic and macroscopic signatures of memory.

Chen, Yihao [Johns Hopkins Univ., Baltimore, MD (U↗

Multimodal Defect Imaging of Pure Tungsten Components Fabricated via Electron Beam Powder Bed Fusion

The utilization of additive manufacturing (AM) techniques for refractory materials in high-temperature environments has significantly expanded because of the ability to fabricate geometrically complex components. Electron beam powder bed fusion (EB-PBF), which provides lower residual stress, a cleaner vacuum environment, and better efficiency for high melting point, is one of the best-suited AM methods to produce advanced refractory components. However, the property variation attributed to the heterogeneous microstructure and process-induced defects has hindered the widespread adoption of EB-PBF-produced material like tungsten. While numerous in-situ monitoring and defect detection methods have been demonstrated for EB-PBF, a workflow that compares and evaluates process-induced abnormalities from different imaging perspectives is still limited. This study examines a feature-embedded tungsten component manufactured via the EB-PBF process to demonstrate the defect detection capabilities of a multimodal defect imaging workflow. The predefined and process-induced defects are evaluated by harnessing various imaging techniques, including in-situ electron imaging, layerwise near-infrared (NIR) imaging, post-build high-energy x-ray computed tomography (CT), and conventional destructive metallography. The results highlight the strengths and limitations of distinctive defect imaging techniques concerning specific defect types, sizes, and conditions. It was found that electron imaging can provide more abnormal detection capabilities while maintaining a higher measuring accuracy, against the conventional metallography in this case study, compared with NIR and CT imaging techniques.

36 MATERIALS SCIENCE↗

Residual Stress Measurements After Proof and Flight: ETP-0403

The intent of this testing was to evaluate the residual stresses that occur in and around the attachment details of a case stiffener segment that has been subjected to flight/recovery followed by proof loading. Not measured in this test were stresses relieved at joint disassembly due to out-of-round and interference effects, and those released by cutting the specimens out of the case segment. The test article was lightweight case stiffener segment 1U50715, S/N L023 which was flown in the forward stiffener position on flight SRM 14A and in the aft position on flight SRM24A. Both of these flights were flown with the 3 stiffener ring configuration. Stiffener L023 had a stiffener ring installed only on the aft stub in its first flight, and it had both rings installed on its second flight. No significant post flight damage was found on either flight. Finally, the segment was used on the DM-8 static test motor in the forward position. No stiffener rings were installed. It had only one proof pressurization prior to assignment to its first use, and it was cleaned and proof tested after each flight. Thus, the segment had seen 3 proof tests, two flight pressurizations, and two low intensity water impacts prior to manufacturing for use on DM-8. On DM-8 it received one static firing pressurization in the horizontal configuration. Residual stresses at the surface and in depth were evaluated by both the x-ray diffraction and neutron beam diffraction methods. The x-ray diffraction evaluations were conducted by Technology for Energy Corporation (TEC) at their facilities in Knoxville, TN. The neutron beam evaluations were done by Atomic Energy of Canada Limited (AECL) at the Chalk River Nuclear Laboratories in Ontario. The results showed general agreement with relatively high compressive residual stresses on the surface and moderate to low subsurface tensile residual stresses.

Webster, Ronald L..↗

Apollo 6, Anomaly Report No. 6: Abnormal Structural Performance During Launch Phase

Approximately 2 minutes 13 seconds after lift-off of the Apollo 6 mission, abrupt changes of strain, vibration, and acceleration measurements were indicated in the S-IVB, instrument unit, adapter, lunar module, and command and service modules; photographs showed objects coming from the area of the adapter. The adapter, however, continued to sustain the required loads with no impairment of the mission. The investigation was first focused upon the understanding of the coupled vibration modes and characteristics of the launch vehicle and spacecraft.' Extensive test programs were conducted. It was eventually concluded that the adapter failure was not caused by vibration. Extensive study of the airborne photography and other evidence indicated that a large area of the adapter had lost inner facesheet from the honeycomb sandwich panels. Loads and stresses resulting from vibration were determined to be insufficient to initiate such a failure. The investigation was then directed toward determining the range of pressures that could have been trapped in the Apollo 6 adapter sandwich panels, and toward determining the tolerance of the panels to withstand pressure with various degrees of flaws such as adhesive voids and facesheet dents. The degradation effects of moisture and heat exposure on the adhesive strength were also studied and tested. These tests and analyses led to the conclusion that pressure internal to the sandwich panels could have caused the failure, if a large flaw existed. The pressure buildup would have been caused by aerodynamic heating effects on air and moisture trapped in the panel. The probable cause of the failure was found in the original ultrasonic inspection scan record of the affected adapter panel. In the center of the region where the adapter failed, horizontally along the station 709 panel splice, the record contained two thick anomalous lines extending several feet. Without an X-ray record of this region, the significance of this particular scan record cannot be fully understood. However, since all other evidence had indicated that the adhesive had to be weakened in a rather large area to initiate the failure, the investigation was focused intently upon the station 709 splices of other adapters. Sufficient information was developed to verify that deficient assembly techniques have consistently resulted in abnormalities in the structure at this station. These abnormalities were identified in adapters 12, 13, 14, 15, and 16. Before the splice abnormalities were pinpointed, corrective action was taken to reduce pressure "buildup in the honeycomb panels and to reduce heat degrading effects on the adhesive. This was done by drilling vent holes in the inner facesheet and covering the outer facesheet with cork. The adapters having identified abnormalities in the station 709 splice are being repaired, and the contractor is investigating ways of avoiding these -abnormalities in panels yet to be bonded.

Source record↗

Strain-induced topological transitions in epitaxial films of cadmium arsenide

Strain in epitaxial thin films can have a profound influence on their electronic bands. Here, we investigate the influence of epitaxial strains on the electronic states of (001)-oriented cadmium arsenide (Cd 3 As 2 ) thin films, which are two-dimensional topological insulators. Cd 3 As 2 films were grown coherently on Al 1−x In x Sb buffer layers with varying degrees of lattice mismatch to Cd 3 As 2 , producing in-plane stresses that ranged from compressive to tensile. Using magneto-transport measurements of the Landau levels of the lowest subbands, we show that films under compression remain in the two-dimensional topological insulator state, while a sufficiently large tensile stress induces a transition to a topologically trivial phase. These results qualitatively agree with expectations from theoretical models for these films. At the topological transition, the gap avoids completely closing, and we discuss possible origins.

36 MATERIALS SCIENCE↗

Optimizing ablator thickness for laser shock experiments

In laser shock experiments, a well-defined, flat-top shock wave at the ablator/sample interface is important for accurately probing material response under uniaxial strain compression. However, the relationship between the ablator thickness and the resulting shock wave characteristics is insufficiently understood, limiting the ability to design optimal experiments. To address this need, we conducted a systematic experimental study using a 100 J laser to shock-compress polyimide ablators to peak stresses ranging from 20.4 to 111.6 GPa. Laser interferometry diagnostics measured the transmitted wave profiles at the ablator/sample interface, consistently showing a single jump followed by a constant peak state before the arrival of release waves. Here, by analyzing shock transit time, flat-top duration, and stress, our results establish a framework for selecting ablator thickness to maximize the flat-top duration, improving the precision and reproducibility of laser shock experiments.

design of experiments↗