SEARCH · Search NASA

Results for “Spatial Behavior”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 1,513 records · Page 84

All sources NASA NTRS DOE OSTI

NASA NTRS1979

Spectral evolution of multiply impulsive solar bursts

Some results from the analysis of a set of multiply impulsive hard X-ray and microwave solar bursts are presented, showing that some bursts can exhibit widely different magnetic-field strengths at different times. Two categories of microwave spectral behavior are identified: those events during which the microwave turnover frequency (MTF) and spectral shape (SS) remain the same from peak to peak, and those during which the MTF and SS change significantly. These categories correspond to two classes of multiply impulsive bursts: those for which the emission can be characterized by a constant magnetic field and therefore a single source region, and those in which groups of component spikes appear to originate in regions of different magnetic-field strengths, corresponding to separate source regions which flare sequentially. With regard to the latter type, examples are presented, the discrete flaring regions are examined, and their spatial separations are estimated.

Karpen, J. T.↗

DOE OSTI2018

Chirped-probe-pulse femtosecond CARS thermometry in turbulent spray flames

This paper presents temperature measurements in turbulent dilute and dense spray flames using single-laser-shot chirped-probe-pulse femtosecond coherent anti-Stokes Raman spectroscopy (CPP-fs-CARS). This ultrafast technique, with a repetition rate of 5 kHz, is applied to the piloted Sydney Needle Spray Burner (SYNSBURN TM ). The burner system features air-blast atomization of liquid injected from a needle that can be translated within a co-flowing air stream. The pilot-stabilized spray flames can range between the two extremes of dense and dilute by physically translating the needle tip relative to the burner's exit plane. The CPP-fs-CARS set-up has achieved integration times of 3 picoseconds (ps) as well as spatial resolution of approximately 800 µm along beam propagation and 60 µm in the transverse dimension. Brief details of the technique, calibration, correction of interferences, and spectral fitting processes are presented along with estimates of the associated error. The measurements are compared against well-established, line Raman–Rayleigh data for temperature collected in a turbulent CH 4 /air jet diffusion flame, which is largely non-sooting. At peak gaseous flame temperatures of up to 2512 K, the relative accuracy and precision were 2.8% and ±3.4%, respectively. Measurements in turbulent spray flames are shown after applying the relevant corrections based on non-resonant background (NRB) behavior and camera saturation effects on the shape of the CARS signal spectrum. Preliminary mapping of the temperature fields demonstrates the wealth of information available in this dataset which will provide insights into the spatio-temporal structure of spray flames once relevant statistical analysis is applied.

42 ENGINEERING↗

DOE OSTI2022

Observation of an intermediate state during lithium intercalation of twisted bilayer MoS2

Abstract Lithium intercalation of MoS 2 is generally believed to introduce a phase transition from H phase (semiconducting) to T phase (metallic). However, during the intercalation process, a spatially sharp boundary is usually formed between the fully intercalated T phase MoS 2 and non-intercalated H phase MoS 2 . The intermediate state, i.e ., lightly intercalated H phase MoS 2 without a phase transition, is difficult to investigate by optical-microscope-based spectroscopy due to the narrow size. Here, we report the stabilization of the intermediate state across the whole flake of twisted bilayer MoS 2 . The twisted bilayer system allows the lithium to intercalate from the top surface and enables fast Li-ion diffusion by the reduced interlayer interaction. The E 2g Raman mode of the intermediate state shows a peak splitting behavior. Our simulation results indicate that the intermediate state is stabilized by lithium-induced symmetry breaking of the H phase MoS 2 . Our results provide an insight into the non-uniform intercalation during battery charging and discharging, and also open a new opportunity to modulate the properties of twisted 2D systems with guest species doping in the Moiré structures.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Quantum fragmentation in the extended quantum breakdown model

We introduce a one-dimensional (1D) extended quantum breakdown model comprising a fermionic and a spin degree of freedom per site, and featuring a spatially asymmetric breakdown-type interaction between the fermions and spins. Furthermore, our model resembles the breakdown process of particles incident into a cloud chamber with nonzero quantum amplitudes of both exciting and not exciting the local vapor atoms. We analytically show that, in the absence of any magnetic field for the spins, the model exhibits Hilbert space fragmentation within each symmetry sector into exponentially many Krylov subspaces and hence displays nonthermal dynamics. Here, we demonstrate that the fragmentation naturally occurs in an entangled basis and thus provides an example of “quantum fragmentation.” Besides establishing the nature of fragmentation analytically, we also study the long-time behavior of the entanglement entropy and its deviation from the expected Page value as a probe of ergodicity in the system. Upon introducing a magnetic field for the spins, most of the Krylov subspaces merge and the model becomes chaotic. Finally, we study the effects of strong randomness on the system and observe behavior similar to that of many-body localized systems.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

DOE OSTI2021

Deployment of the Finite Volume Method in Pronghorn for Gas and Salt cooled Pebble Bed Reactors

This report summarizes the activities related to ”Complete FHR and HTGR pebble bed simulator, including initial validation” funded by the NEAMS thermal-hydraulics focus area. The activity revolves around the coarse-mesh thermal-hydraulics code Pronghorn and its application to gas and salt cooled Pebble bed reactor (PBR). The main difference between gas and salt cooled PBR from a thermal-hydraulics perspective is the fluid. To address the difference in fluid behavior, two separate approaches are implemented in the MOOSE Navier Stokes module: 1) a Boussinesq approximation and 2) a fully compressible formulation. The developed finite volume method capabilities are used for improving pre-existing gas-cooled and salt-cooled pebble-bed reactor models. A steady-state, multiphysics gas-cooled pebble-bed reactor model is created that couples the equilibrium core depletion capability developed in previous work, and the finite volume method capability developed for this report. The salt-cooled pebble-bed reactor model is upgraded to use the incompressible finite volume method capability and then extended to three spatial dimensions. Finally, several verification-and-validation exercises performed with Pronghorn are documented using the verification-and-validation report of the MooseDoc system. The goal of this effort to document the verification-and-validation level of Pronghorn and improve stakeholder confidence in the results obtained with Pronghorn.

97 MATHEMATICS AND COMPUTING↗

NASA NTRS1995

Combining Speed Information Across Space

We used speed discrimination tasks to measure the ability of observers to combine speed information from multiple stimuli distributed across space. We compared speed discrimination thresholds in a classical discrimination paradigm to those in an uncertainty/search paradigm. Thresholds were measured using a temporal two-interval forced-choice design. In the discrimination paradigm, the n gratings in each interval all moved at the same speed and observers were asked to choose the interval with the faster gratings. Discrimination thresholds for this paradigm decreased as the number of gratings increased. This decrease was not due to increasing the effective stimulus area as a control experiment that increased the area of a single grating did not show a similar improvement in thresholds. Adding independent speed noise to each of the n gratings caused thresholds to decrease at a rate similar to the original no-noise case, consistent with observers combining an independent sample of speed from each grating in both the added- and no-noise cases. In the search paradigm, observers were asked to choose the interval in which one of the n gratings moved faster. Thresholds in this case increased with the number of gratings, behavior traditionally attributed to an input bottleneck. However, results from the discrimination paradigm showed that the increase was not due to observers' inability to process these gratings. We have also shown that the opposite trends of the data in the two paradigms can be predicted by a decision theory model that combines independent samples of speed information across space. This demonstrates that models typically used in classical detection and discrimination paradigms are also applicable to search paradigms. As our model does not distinguish between samples in space and time, it predicts that discrimination performance should be the same regardless of whether the gratings are presented in two spatial intervals or two temporal intervals. Our last experiment largely confirmed this prediction.

Verghese, Preeti↗

NASA NTRS1995

Combining Speed Information Across Space

Verghese, Preeti↗

NASA NTRS2010

Orbiter BLT Flight Experiment Wind Tunnel Simulations: Nearfield Flowfield Imaging and Surface Thermography

This paper reports a series of wind tunnel tests simulating the near-field behavior of the Space Shuttle Orbiter Boundary Layer Transition Detailed Test Objective (BLT DTO) flight experiment. Hypersonic flow over a flat plate with an attached BLT DTO-shaped trip was tested in a Mach 10 wind tunnel. The sharp-leading-edge flat plate was oriented at an angle of 20 degrees with respect to the freestream flow, resulting in post-shock edge Mach number of approximately 4. The flowfield was visualized using nitric oxide (NO) planar laser-induced fluorescence (PLIF). Flow visualizations were performed at 10 Hz using a wide-field of view and high-resolution NO PLIF system. A lower spatial resolution and smaller field of view NO PLIF system visualized the flow at 500 kHz, which was fast enough to resolve unsteady flow features. At the lowest Reynolds number studied, the flow was observed to be laminar and mostly steady. At the highest Reynolds number, flow visualizations showed streak instabilities generated immediately downstream of the trip. These instabilities transitioned to unsteady periodic and spatially irregular structures downstream. Quantitative surface heating imagery was obtained using the Temperature Sensitive Paint (TSP) technique. Comparisons between the PLIF flow visualizations and TSP heating measurements show a strong correlation between flow patterns and surface heating trends.

Danehy, Paul M.↗

NASA NTRSDate not supplied

Exploring rock-regolith interfaces in Jezero crater with Mars 2020 SHERLOC

The Perseverance rover successfully landed in Jezero crater, Mars in February 2021 at the Octavia E. Butler landing site and began its mission to explore and sample an ancient crater lake basin. Principal goals of the Mars 2020 mission include characterizing the geology of Mars and seeking signs of ancient microbial life via the spacecraft cameras and spectroscopic instruments onboard. The Scanning Habitable Environments with Raman and Luminescence for Organics and Chemicals (SHERLOC) instrument is a deep UV Raman spectrometer that utilizes a 248.6nm pulsed laser. Part of SHERLOC is a color camera known as the Wide Angle Topographic Sensor for Operations and eNgineering (WATSON). The SHERLOC suite provides coordinated, spectroscopic and imaging capabilities at high spatial resolution, to detect minerals and organic molecules in microtextural context. By pairing high spatial resolution (~100 μm) resonance Raman and native fluorescence spectroscopy with microscopic imaging in a novel spacecraft capability, SHERLOC enables texture-specific molecular composition measurements of rock and regolith targets on Mars. Coordinated rock-regolith observations illuminate unique insights into weathering processes and thereby to primary properties of rocks in Jezero crater. This work describes the potential of rock-regolith interfaces to preserve unique records of geological processes in Jezero crater and can powerfully supplement observations of the more general rock record on Mars. Linking observations of local rock texture with associated regolith reveals important lithologic information based on the interrelationship between differential weathering behavior and mineralogy, grain size, and cement chemistry. Preliminary observations indicate that the polygonally fractured lithotype common near the Octavia E. Butler landing site may weather by granular disintegration and/or surface creep, a relation that can be uniquely observed at the rock-regolith interface. SHERLOC -specific observations of microtextural and elemental composition transitions presented here trace rock-regolith boundaries at multiple indurated surfaces adjacent to regolith. At these locales, grain-scale based examinations suggest chemical weathering could be related to a variably distributed coating or rind on dark rock targets that may be mafic in composition. Granule deposits overlying widely distributed fine-grained material are also observable. Mineral identifications of each phase are presented, with cross-scale comparisons to the remote insights gained by the SuperCam instrument.

Emily L. Cardarelli↗

DOE OSTI2020

Local measurements of the spatial magnetic field distribution in a z-pinch plasma during and near stagnation using polarization spectroscopy

We present here the detailed measurements of radial distribution of the magnetic field in a gas-puff z-pinch plasma at the final stages of the implosion phase and at stagnation. While the measurements are chordal, the radial distribution of different charge states was utilized to measure the magnetic field locally for certain radii, so that unlike chordal measurements in general, the magnetic field radial distribution was obtained with no need for the Abel inversion of the data. The distribution was measured using the Zeeman effect via a novel spectroscopic technique, at several axial locations, and demonstrates striking features such as the peak field remaining at a radius much larger than the stagnation radius at all times. Furthermore, while the distribution observed is sometimes monotonic with respect to the radius, it is often not, a behavior that can be linked to 2D features in the plasma column resulting from the Rayleigh–Taylor instability. The current flowing through the stagnating plasma was found to be a small fraction of the total current, resulting in clearly insufficient magnetic pressure to balance the plasma pressure at stagnation. The magnetic field data, taken over several axial positions, are used to obtain the true inductance in the imploding plasma for the first time; it is found that the data cannot explain the current turnover at stagnation. A simulation with the MACH2-Tabular Collisional-Radiative Equilibrium magnetohydrodynamics code in the r–z plane shows that the peak of the magnetic field pinches to a much smaller radius than is observed in the spectroscopic data. Furthermore, the depth of the computed current turnover at stagnation is smaller than the measured one. The two observed features of a radially extended magnetic field at stagnation together with a deep current turnover are a challenge to match in simulations. Various calculations and estimates of the inductive and resistive load voltages are examined to ascertain if they are responsible for the observed current notch. Lastly, the results demonstrate that the knowledge of the true inductance in the driven load requires such magnetic-field-distribution measurements and that imaging data or electrical measurements are insufficient.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

NASA NTRS1991

Phase-averaged measurements of perturbations introduced into boundary layers

Large-scale structures in turbulent and transitional wall-bounded flows make a significant contribution to the Reynolds stress and turbulent energy. The behavior of these structures is examined. Small perturbations are introduced into a laminar and a turbulent boundary layer to trigger the formation of large-scale features. Both flows use the same inlet unit Reynolds number, and they experience the same pressure gradient history, i.e. a favorable pressure gradient (FPG) followed by an adverse pressure gradient (APG). The perturbation consists of a small short duration flow repetitively introduced through a hole in the wall located at the C(sub p) minimum. Hot-wire data are averaged on the basis of the phase of the disturbance, and automation of the experiment was used to obtain measurements on large spatially dense grids. In the turbulent boundary, the perturbation evolves into a vortex loop which retains its identity for a considerable streamwise distance. In the laminar layer, the perturbation decays to a very small magnitude before growing rapidly and triggering the transition process in the APG. The 'time-like' animations of the phase-averaged data are used to gain insight into the naturally occurring physical mechanisms in each flow.

Watmuff, Jonathan H.↗

NASA NTRS1993

Stationary eddies in the Mars general circulation as simulated by the NASA-Ames GCM

Quasistationary eddies are prominent in a large set of simulations of the Mars general circulation performed with the NASA-Ames GCM. Various spacecraft observations have at least hinted at the existence of such eddies in the Mars atmosphere. The GCM stationary eddies appear to be forced primarily by the large Mars topography, and (to a much lesser degree) by spatial variations in the surface albedo and thermal inertia. The stationary eddy circulations exhibit largest amplitudes at high altitudes (above 30-40 km) in the winter extratropical regions. In these regions they are of planetary scale, characterized largely by zonal wavenumbers 1 and 2. Southern Hemisphere winter appears to be dominated by a very strong wave 1 pattern, with both waves 1 and 2 being prominent in the Northern Hemisphere winter regime. This difference seems to be basically understandable in terms of differences in the topography in the two hemispheres. The stationary eddies in the northern winter extratropics are found to increase in amplitude with dust loading. This behavior appears to be at least partly associated with changes in the structure of the zonal-mean flow that favor a greater response to wave 1 topographic forcing. There are also strong stationary eddy circulations in the tropics and in the summer hemisphere. The eddies in the summer subtropics and extratropics arc substantially stronger in southern summer than in northern summer. The summer hemisphere stationary circulations are relatively shallow and are characterized by smaller zonal scales than those in the winter extratropics.

Barnes, J. R.↗

NASA NTRS2002

Measurements of Turbulence Attenuation by a Dilute Dispersion of Solid Particles in Homogeneous Isotropic Turbulence

This research addresses turbulent gas flows laden with fine solid particles at sufficiently large mass loading that strong two-way coupling occurs. By two-way coupling we mean that the particle motion is governed largely by the flow, while the particles affect the gas-phase mean flow and the turbulence properties. Our main interest is in understanding how the particles affect the turbulence. Computational techniques have been developed which can accurately predict flows carrying particles that are much smaller than the smallest scales of turbulence. Also, advanced computational techniques and burgeoning computer resources make it feasible to fully resolve very large particles moving through turbulent flows. However, flows with particle diameters of the same order as the Kolmogorov scale of the turbulence are notoriously difficult to predict. Some simple flows show strong turbulence attenuation with reductions in the turbulent kinetic energy by up to a factor of five. On the other hand, some seemingly similar flows show almost no modification. No model has been proposed that allows prediction of when the strong attenuation will occur. Unfortunately, many technological and natural two-phase flows fall into this regime, so there is a strong need for new physical understanding and modeling capability. Our objective is to study the simplest possible turbulent particle-laden flow, namely homogeneous, isotropic turbulence with a uniform dispersion of monodisperse particles. We chose such a simple flow for two reasons. First, the simplicity allows us to probe the interaction in more detail and offers analytical simplicity in interpreting the results. Secondly, this flow can be addressed by numerical simulation, and many research groups are already working on calculating the flow. Our detailed data can help guide some of these efforts. By using microgravity, we can further simplify the flow to the case of no mean velocity for either the turbulence or the particles. In fact the addition of gravity as a variable parameter may help us to better understand the physics of turbulence attenuation. The experiments are conducted in a turbulence chamber capable of producing stationary or decaying isotropic turbulence with nearly zero mean flow and Taylor microscale Reynolds numbers up to nearly 500. The chamber is a 410 mm cubic box with the corners cut off to make it approximately spherical. Synthetic jet turbulence generators are mounted in each of the eight corners of the box. Each generator consists of a loudspeaker forcing a plenum and producing a pulsed jet through a 20 mm diameter orifice. These synthetic jets are directed into ejector tubes pointing towards the chamber center. The ejector tubes increase the jet mass flow and decrease the velocity. The jets then pass through a turbulence grid. Each of the eight loudspeakers is forced with a random phase and frequency. The resulting turbulence is highly Isotropic and matches typical behavior of grid turbulence. Measurements of both phases are acquired using particle image velocimetry (PIV). The gas is seeded with approximately 1 micron diameter seeding particles while the solid phase is typically 150 micron diameter spherical glass particles. A double-pulsed YAG laser and a Kodak ES-1.0 10-bit PIV camera provide the PIV images. Custom software is used to separate the images into individual images containing either gas-phase tracers or large particles. Modern high-resolution PIV algorithms are then used to calculate the velocity field. A large set of image pairs are acquired for each case, then the results are averaged both spatially and over the ensemble of acquired images. The entire apparatus is mounted in two racks which are carried aboard NASA's KC-135 Flying Microgravity Laboratory. The rack containing the turbulence chamber, the laser head, and the camera floats freely in the airplane cabin (constrained by competent NASA personnel) to minimize g-jitter.

Eaton, John↗

DOE OSTI2022

Machine-Learning of Nonlocal Kernels for Anomalous Subsurface Transport from Breakthrough Curves

Anomalous behavior is ubiquitous in subsurface solute transport due to the presence of high degrees of heterogeneity at different scales in the media. Although fractional models have been extensively used to describe the anomalous transport in various subsurface applications, their application is hindered by computational challenges. Simpler nonlocal models characterized by integrable kernels and finite interaction length represent a computationally feasible alternative to fractional models; yet, the informed choice of their kernel functions still remains an open problem. We propose a general data-driven framework for the discovery of optimal kernels on the basis of very small and sparse data sets in the context of anomalous subsurface transport. Using spatially sparse breakthrough curves recovered from fine-scale particle-density simulations, we learn the best coarse-scale nonlocal model using a nonlocal operator regression technique. Predictions of the breakthrough curves obtained using the optimal nonlocal model show good agreement with fine-scale simulation results even at locations and time intervals different from the ones used to train the kernel, confirming the excellent generalization properties of the proposed algorithm. A comparison with trained classical models and with black-box deep neural networks confirms the superiority of the predictive capability of the proposed model.

97 MATHEMATICS AND COMPUTING↗

DOE OSTI2024

Object-oriented analysis as a foundation for building climate storylines of compounding short-term drought and crop heat stress

Introduction: Crops are vulnerable to precipitation and heat extremes during late spring through summer. Methods: We analyzed for a north-central U.S. region short-term drought and agricultural heat stress during April-May-June-July. We used the 4-km Parameter Elevation Regression on Independent Slopes Model (PRISM) for observations, aggregated to a 25-km grid, and two 25-km Regional Climate Model version 4 (RegCM4) simulns used either GFDL- or MPI-GCM boundary conditions. We chose 1981-2000 as our contemporary time period, and 2041- 2060 as our scenario time period, which used the Representative Concentration Pathway 8.5 emissions scenario. We used object-oriented analysis to identify events of interest in observations and simulations by identifying objects in a space-time domain that meet specified criteria, such as exceeding a heat-stress temperature threshold. The event diagnosis allowed analysis of compound events, occurring when temperature and drought objects overlap. Results: Identified objects yielded events that can undermine agricultural productivity and which are thus relevant to decision makers, making them building blocks for possible climate storylines. The observations and simulations showed similar spatial distributions of event frequencies across the analysis region. However, the simulations attained this distribution by having fewer events that tend to cover larger areas compared to observed events, suggesting that the effective resolution of the simulations was coarser than their 25-km grids. Short-term drought frequency increased and heat-stress frequency decreased in transitioning to the scenario climate. When compounding occurred heat-stress events generally preceded the short-term drought events. The overlapping, compound events tended to be more extreme compared to non-overlapping events of either type. Discussion: The information yielded projected changes in these agriculturally motivated events. One prominent conditional behavior emerging from the work was that a heat-stress event should be a warning to watch for potential drought, as both could compound each other to more intense levels.

54 ENVIRONMENTAL SCIENCES↗

NASA NTRS1980

A heuristic mathematical model for the dynamics of sensory conflict and motion sickness

The etiology of motion sickness is explained in terms of a qualitatively formulated sensory conflict hypothesis. By consideration of the information processing task faced by the central nervous system in estimating body spatial orientation and in controlling active body movement using an internal model referenced control strategy, a mathematical model for sensory conflict generation is developed. The model postulates a major dynamic functional role for sensory conflict signals in movement control, as well as in sensory-motor adaptation. It accounts for the role of active movement in creating motion sickness symptoms in some experimental circumstances, and in alleviating them in others. The relationship between motion sickness produced by sensory rearrangement and that resulting from external motion disturbances is explicitly defined. A nonlinear conflict averaging model is proposed which describes dynamic aspects of experimentally observed subjective discomfort sensation, and suggests resulting behaviors.

Oman, C. M.↗

NASA NTRS1984

Pattern selection in solidification

Directional solidification of alloys produces a wide variety of cellular or lamellar structures which, depending upon growth conditions, may be reproducibly regular or may behave chaotically. It is not well understood how these patterns are selected and controlled or even whether there ever exist sharp selection mechanisms. A related phenomenon is the spatial propagation of a pattern into a system which has been caused to become unstable against pattern-forming deformations. This phenomenon has some features in common with the propagation of sidebranching modes in dendritic solidification. In a class of one-dimensional models, the nonlinear system can be shown to select the propagating mode in which the leading edge of the pattern is just marginally stable. This stability principle, when applicable, predicts both the speed of propagation and the geometrical characteristics of the pattern which forms behind the moving front. A boundary-layer model for fully two or three dimensional solidification problems appears to exhibit similar mathematical behavior.

Langer, J. S.↗

NASA NTRS1995

Diamond thin film temperature and heat-flux sensors

Diamond film temperature and heat-flux sensors are developed using a technology compatible with silicon integrated circuit processing. The technology involves diamond nucleation, patterning, doping, and metallization. Multi-sensor test chips were designed and fabricated to study the thermistor behavior. The minimum feature size (device width) for 1st and 2nd generation chips are 160 and 5 micron, respectively. The p-type diamond thermistors on the 1st generation test chip show temperature and response time ranges of 80-1270 K and 0.29-25 microseconds, respectively. An array of diamond thermistors, acting as heat flux sensors, was successfully fabricated on an oxidized Si rod with a diameter of 1 cm. Some problems were encountered in the patterning of the Pt/Ti ohmic contacts on the rod, due mainly to the surface roughness of the diamond film. The use of thermistors with a minimum width of 5 micron (to improve the spatial resolution of measurement) resulted in lithographic problems related to surface roughness of diamond films. We improved the mean surface roughness from 124 nm to 30 nm by using an ultra high nucleation density of 10(exp 11)/sq cm. To deposit thermistors with such small dimensions on a curved surface, a new 3-D diamond patterning technique is currently under development. This involves writing a diamond seed pattern directly on the curved surface by a computer-controlled nozzle.

Aslam, M.↗