Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “orientational order”

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 1,117 records · Page 62

Nanolayered Features of Collagen-like Peptides

We have been investigating collagen-like model oligopeptides as molecular bases for complex ordered biomimetic materials. The collagen-like molecules incorporate aspects of native collagen sequence and secondary structure. Designed modifications to native primary and secondary structure have been incorporated to control the nanostructure and microstructure of the collagen-like materials produced. We find that the collagen-like molecules form a number of lyotropic rod liquid crystalline phases, which because of their strong temperature dependence in the liquid state can also be viewed as solvent intercalated thermotropic liquid crystals. The liquid crystalline phases formed by the molecules can be captured in the solid state by drying off solvent, resulting in solid nanopatterned (chemically and physically) thermally stable (to greater than 100 C) materials. Designed sequences which stabilize smectic phases have allowed a variety of nanoscale multilayered biopolymeric materials to be developed. Preliminary investigations suggest that chemical patterns running perpendicular to the smectic layer plane can be functionalized and used to localize a variety of organic, inorganic, and organometallic moieties in very simple multilayered nanocomposites. The phase behavior of collagen-like oligopeptide materials is described, emphasizing the correlation between mesophase, molecular orientation, and chemical patterning at the microscale and nanoscale. In many cases, the textures observed for smectic and hexatic phase collagens are remarkably similar to the complex (and not fully understood) helicoids observed in biological collagen-based tissues. Comparisons between biological morphologies and collagen model liquid crystalline (and solidified materials) textures may help us understand the molecular features which impart order and function to the extracellular matrix and to collagen-based mineralized tissues. Initial studies have utilized synthetic collagen-like peptides while future work will also focus on similar sequences generated via genetic engineering methods.

Valluzzi, Regina↗

Scan‐Path‐ and Initial‐State‐Dependent Superdomain Switching in (111)‐Oriented PZT

Polarization switching in ferroelectric materials arises from the collective evolution of complex domain hierarchies, yet deterministic control over these processes remains challenging. Here, we investigate scan-path- and initial-state-dependent switching in epitaxial (111)-oriented PbZr 0.2 Ti 0.8 O 3 thin films using automated AFM-based writing combined with quantitative 3D piezoresponse force microscopy. We show that the scan trajectory acts as an experimentally accessible control parameter for superdomain formation. Box-in-box raster scans reproducibly stabilize ordered stripe superdomains with a reduced subset of symmetry-allowed variants, whereas spiral trajectories generate frustrated mixed-variant states with a broader distribution of final microstructures. Automated pulsing experiments further show that the local superdomain configuration at the nucleation site strongly influences the final written morphology. Phase-field modeling qualitatively reproduces the contrast between representative initial-state geometries and supports the role of compatibility constraints among competing ferroelastic pathways. These findings establish scan-path and initial-state engineering as practical handles to program ferroic order in hierarchical ferroelectric domain structures.

Vasudevan, Rama K. [Oak Ridge National Laboratory ↗

Dose conversion factors for absorbed dose in a mobile phone to absorbed dose in critical organs in an anthropomorphic phantom for emergency dosimetry applications: OSL and TL experimental results, and Monte Carlo simulations

Here, an anthropomorphic phantom was equipped with three smartphones and irradiated in five separate exposure geometries in air using a 137 Cs source. The phantom was fitted internally with thermoluminescence dosimeters (TLD-700 chips) to measure the absorbed dose to the individual organs, while the absorbed doses in the phones were evaluated using optically stimulated luminescence (OSL) from the phones’ surface mount resistors (SMRs), along with reference dosimeters, namely Luxel OSL films and TLD-700 chips. The results indicate a strong dependence on the phone location with respect to the source due to shielding and scattering from the body. Dose conversion factors were calculated for several plausible exposure scenarios in order to correlate the absorbed doses measured in the phones to the average body dose received by the phantom. The measured conversion factors were compared with calculated factors using MCNP6.2 Monte Carlo simulations and an ADAM voxel phantom. Experiments were conducted to assess how the phone orientation and angle impacts the dose received by the phone components. The overall results indicate that the doses measured can vary significantly depending upon exposure geometry (including phone position and angle). However, in rotational irradiation symmetry, the phone doses lead to an over-estimate of the average whole-body dose by an average of ~10%, independent of the phone locations tested in this study. This result, if confirmed in future studies, may be acceptable for triage.

61 RADIATION PROTECTION AND DOSIMETRY↗

Controlling the Magneto-Optical Response in Ultrathin Films of EuO 1– x via Interface Engineering with Ferroelectric BaTi 2 O 5

Utilizing pulsed laser deposition, a film of EuO 1–x was deposited onto a Si(001) substrate with MgO buffer and compared to the same heterostructure with an additional BaTi 2 O 5 thin film on top of the EuO 1–x surface. X-ray diffraction (XRD) indicates the films crystallize into a preferred EuO(111) orientation; it also reveals the clear presence of EuSi 2 , which suggests Si or Eu diffuses across the MgO buffer layer. EuO 1–x films exhibit a ferromagnetic (FM) signature and temperature-dependent exchange bias, indicated by MOKE measurements, suggesting the presence of a magnetic order well above the EuO Curie temperature with possible origins in charge carrier density near the interface. In comparison, an antiferromagnetic character persists well above the EuO Curie temperature of 69 K and the enhanced Curie temperature of 150 K for BaTi 2 O 5 films grown on the EuO 1–x films. Furthermore, the antiferromagnetic behavior is not seen in thicker EuO 1–x thin films when integrated with other ferroelectric (FE) phases of the BaO–TiO 2 system, suggesting an origin in the perturbed charge population at the BaTi 2 O 5 /EuO 1–x interface.

36 MATERIALS SCIENCE↗

Phase-field framework with constraints and its applications to ductile fracture in polycrystals and fatigue

Abstract Modeling of ductile fracture in polycrystalline structures is challenging, since it requires integrated modeling of cracks, crystal plasticity, and grains. Here we extend the typical phase-field framework to the situations with constraints on the order parameters, and formulate two types of phase-field models on ductile fracture. The Type-I model incorporates three sets of order parameters, which describe the distributions of cracks, plastic strain, and grains, respectively. Crystal plasticity is employed within grain interiors accommodated by J 2 plasticity at grain boundaries. The applications of the Type-I model to single crystals and bicrystals demonstrate the influences of grain orientations and grain boundaries on crack growth. In the Type-II model, J 2 plasticity is assumed for the whole system and grain structures are neglected. Taking advantage of the efficiency of the fast Fourier transform, our Type-II model is employed to study low cycle fatigue. Crack closure and striation-like patterning of plastic strain are observed in the simulations. Crack growth rate is analyzed as a function of the J -integral, and the simulated fatigue life as a function of plastic strain agrees with the Coffin–Manson relation without a priori assumption.

36 MATERIALS SCIENCE↗

Magnetic Weyl Semimetallic Phase in Thin Films of Eu 2 Ir 2 O 7

The interplay between electronic interactions and strong spin-orbit coupling is expected to create a plethora of fascinating correlated topological states of quantum matter. Of particular interest are magnetic Weyl semimetals originally proposed in the pyrochlore iridates, which are only expected to reveal their topological nature in thin film form. To date, however, direct experimental demonstrations of these exotic phases remain elusive, due to the lack of usable single crystals and the insufficient quality of available films. Here, in this study, we report on the discovery of signatures for the long-sought magnetic Weyl semimetallic phase in (111)-oriented Eu 2 Ir 2 O 7 high-quality epitaxial thin films. We observed an intrinsic anomalous Hall effect with colossal coercivity but vanishing net magnetization, which emerges right below the onset of a peculiar magnetic phase with all-in-all-out (AIAO) antiferromagnetic ordering. The anomalous Hall conductivity obtained experimentally is consistent with the theoretical prediction, likely arising from the nonzero Berry curvature emanated by Weyl node pairs near the Fermi level that act as sources and sinks of Berry flux, activated by broken cubic crystal symmetry at the top and bottom terminations of the thin film.

36 MATERIALS SCIENCE↗

Iron-doped lithium niobate as a read-write holographic storage medium

The response of iron-doped lithium niobate under conditions corresponding to hologram storage and retrieval is described, and the material characteristics are discussed. The optical sensitivity can be improved by heavy chemical reduction of lightly doped crystals such that most of the iron is in the divalent state, the remaining part being trivalent. The best reduction process found to be reproducible so far is the anneal of the doped crystal in the presence of a salt such as lithium carbonate. It is shown by analysis and simulation that a page-oriented read-write holographic memory with 1000 bits per page would have a cycle time of about 60 msec and a signal-to-noise ratio of 27 dB. This cycle time, although still too long for a practical memory, represents an improvement of two orders of magnitude over that of previous laboratory prototypes using a thermoplastic storage medium

Alphonse, G. A.↗

The distribution of earthquakes with depth and stress in subducting slabs

The global variation of Benioff zone seismicity with depth and the orientation of stress axes of deep and intermediate earthquakes is explained using numerical models of subducting slabs. Models that match the seismicity and stress require a barrier to flow at the 670 km seismic discontinuity. The barrier may be a viscosity increase of at least an order of magnitude or a chemical discontinuity. Instantaneous flow is subparallel to the slabs for models with a viscosity increase but contorted for models with a chemical barrier. Log N (number of earthquakes) decreases linearly to 250-300 km depth and increases thereafter. Stress magnitude in the models shows the same pattern, in accord with experiments showing N proportional to e(k-sigma), with k a constant and sigma the stress magnitude. The models predict downdip compression in the slabs at depths below 300-400 km, as observed for earthquake stress axes.

Vassiliou, M. S.↗

The dynamic aspects of thermo-elasto-viscoplastic snap-through and creep buckling phenomena

Use of a mathematical model and solution methodology, to examine dynamic buckling and dynamic postbuckling behavior of shallow arches and spherical caps made of a realistic material and undergoing non-isothermal, elasto-viscoplastic deformation was examined. Thus, geometric as well as material type nonlinearities of higher order are included in this analysis. The dynamic stability problem is studied under impulsive loading and suddenly applied loading with loads of constant magnitude and infinite duration. A finite element model was derived directly from the incrementally formulated nonlinear shell equations, by using a tensor-oriented procedure. As an example of the results, the time history of the midspan displacement of a damped shallow circular arch is presented.

Riff, R.↗

A different view of plasma flow inside P/Halley

The Giotto spacecraft carried two different instruments - the JPA and the IMS - for the observation of hot ions in the coma of P/Halley. Although there are many similarities in the time and distance profiles of the plasma flow parameters (bulk velocity, number density, and temperature) computed from the two data sets, there are also some significant differences, especially at cometocentric distances less than 500,000 km. The principal discrepancies between the JPA results presented by Formisano et al. (1990) and the IMS observations are: (1) the IMS did not detect the levelling off of the speed and temperature profiles that Formisano et al. interpreted as flow stabilization; (2) the IMS detected differential north-south flow between the solar wind and cometary ions for only a brief interval when the magnetic field was oriented nearly southward, whereas Formisano et al. reported more extensive differential north-south flow that was independent of the direction of the field; (3) the JPA ion densities were factors of 2 to 4 higher than the IMS ion densities which, in turn, were an order of magnitude greater than theoretical values.

Neugebauer, M.↗

OSSE observations of the ultraluminous infrared galaxies ARP 220 and MRK 273

The results of oriented scintillation spectrometer experiment (OSSE) observations of the ultraluminous infrared galaxies Arp 220 and Mrk 273 are reported. The pointings of Arp 220 and Mrk 273 concentrated on their upper limits. The gamma ray luminosities from these sources were found to be between one and two orders of magnitude smaller than the infrared luminosities. Multiwavelength luminosity spectra are produced from the radio to the gamma ray regime, and are compared with the typical multiwavelength spectra of active galactic nuclei. The lack of measured gamma ray emission provides no evidence for the existence of buried active galactic nuclei in these ultraluminous infrared galaxies, but is consistent with an origin of the infrared luminosity from starburst activity.

Dermer, C. D.↗

The role of the epidermis and cortex in gravitropic curvature of maize roots

In order to determine the role of the epidermis and cortex in gravitropic curvature of seedling roots of maize (Zea mays L. cv. Merit), the cortex on the two opposite flanks was removed from the meristem through the growing zone; gravitropic curvature was measured with the roots oriented horizontally with the cut flanks either on the upper and lower side, or on the lateral sides as a wound control. Curvature was slower in both these treatments (53 degrees in 5 h) than in intact roots (82 degrees), but there was no difference between the two orientations in extent and rate of curvature, nor in the latent time, showing that epidermis and cortex were not the site of action of the growth-regulating signal. The amount of cortex removed made no difference in the extent of curvature. Curvature was eliminated when the endodermis was damaged, raising the possibility that the endodermis or the stele-cortex interface controls gravitropic curvature in roots. The elongation rate of roots from which just the epidermis had been peeled was reduced by 0.01 mM auxin (indole-3-acetic acid) from 0.42 to 0.27 mm h-1, contradicting the hypothesis that only the epidermis responds to changes in auxin activity during gravistimulation. These observations indicate that gravitropic curvature in maize roots is not driven by differential cortical cell enlargement, and that movement of growth regulator(s) from the tip to the elongating zone is unlikely to occur in the cortex.

NASA Program Space Biology↗

Three-dimensional analysis of vestibular efferent neurons innervating semicircular canals of the gerbil

Anterograde labeling techniques were used to examine peripheral innervation patterns of vestibular efferent neurons in the crista ampullares of the gerbil. Vestibular efferent neurons were labeled by extracellular injections of biocytin or biotinylated dextran amine into the contralateral or ipsilateral dorsal subgroup of efferent cell bodies (group e) located dorsolateral to the facial nerve genu. Anterogradely labeled efferent terminal field varicosities consist mainly of boutons en passant with fewer of the terminal type. The bouton swellings are located predominately in apposition to the basolateral borders of the afferent calyces and type II hair cells, but several boutons were identified close to the hair cell apical border on both types. Three-dimensional reconstruction and morphological analysis of the terminal fields from these cells located in the sensory neuroepithelium of the anterior, horizontal, and posterior cristae were performed. We show that efferent neurons densely innervate each end organ in widespread terminal fields. Subepithelial bifurcations of parent axons were minimal, with extensive collateralization occurring after the axons penetrated the basement membrane of the neuroepithelium. Axonal branching ranged between the 6th and 27th orders and terminal field collecting area far exceeds that of the peripheral terminals of primary afferent neurons. The terminal fields of the efferent neurons display three morphologically heterogeneous types: central, peripheral, and planum. All cell types possess terminal fields displaying a high degree of anisotropy with orientations typically parallel to or within +/-45 degrees of the longitudinal axis if the crista. Terminal fields of the central and planum zones predominately project medially toward the transverse axis from the more laterally located penetration of the basement membrane by the parent axon. Peripheral zone terminal fields extend predominately toward the planum semilunatum. The innervation areas of efferent terminal fields display a trend from smallest to largest for the central, peripheral, and planum types, respectively. Neurons that innervate the central zone of the crista do not extend into the peripheral or planum regions. Conversely, those neurons with terminal fields in the peripheral or planum regions do not innervate the central zone of the sensory neuroepithelium. The central zone of the crista is innervated preferentially by efferent neurons with cell bodies located in the ipsilateral group e. The peripheral and planum zones of the crista are innervated preferentially by efferent neurons with cell bodies located in the contralateral group e. A model incorporating our anatomic observations is presented describing an ipsilateral closed-loop feedback between ipsilateral efferent neurons and the periphery and an open-loop feed-forward innervation from contralateral efferent neurons. A possible role for the vestibular efferent neurons in the modulation of semicircular canal afferent response dynamics is proposed.

Non-NASA Center↗

Magnetospheric Multiscale (MMS) Mission Attitude Ground System Design

This paper describes the attitude ground system (AGS) design to be used for support of the Magnetospheric MultiScale (MMS) mission. The AGS exists as one component of the mission operations control center. It has responsibility for validating the onboard attitude and accelerometer bias estimates, calibrating the attitude sensors and the spacecraft inertia tensor, and generating a definitive attitude history for use by the science teams. NASA's Goddard Space Flight Center (GSFC) in Greenbelt, Maryland is responsible for developing the MMS spacecraft, for the overall management of the MMS mission, and for mission operations. MMS is scheduled for launch in 2014 for a planned two-year mission. The MMS mission consists of four identical spacecraft flying in a tetrahedral formation in an eccentric Earth orbit. The relatively tight formation, ranging from 10 to 400 km, will provide coordinated observations giving insight into small-scale magnetic field reconnection processes. By varying the size of the tetrahedron and the orbital semi-major axis and eccentricity, and making use of the changing solar phase, this geometry allows for the study of both bow shock and magnetotail plasma physics, including acceleration, reconnection, and turbulence. The mission divides into two phases for science; these phases will have orbit dimensions of 1.2 x 12 Earth radii in the first phase and 1.2x25 Earth radii in the second in order to study the dayside magnetopause and the nightside magnetotail, respectively. The orbital periods are roughly one day and three days for the two mission phases. Each of the four MMS spacecraft will be spin stabilized at 3 revolutions per minute (rpm), with the spin axis oriented near the ecliptic north pole but tipped approximately 2.5 deg towards the Sun line. The main body of each spacecraft will be an eight-sided platform with diameter of 3.4 m and height of 1.2 m. Several booms are attached to this central core: two axial booms of 14.9 m length, two radial magnetometer booms of 5 m length, and four radial wire booms of 60 m length. Attitude and orbit control will use a set of axial and radial thrusters. A four-head star tracker and a slit-type digital Sun sensor (DSS) provide input for attitude determination. In addition, an accelerometer will be used for closed-loop orbit maneuver control. The primary AGS product will be a daily definitive attitude history. Due to power limitations, the star tracker and accelerometer data will not be available at all times. However, tracker data from at least 10 percent of each orbit and continuous DSS data will be provided. An extended Kalman filter (EKF) will be used to estimate the three-axis attitude (i.e., spin axis orientation and spin phase) and rotation rate for all times when the tracker data is valid. For other times, the attitude is generated by assuming a constant angular momentum vector in the inertial frame. The DSS sun pulse will provide a timing signal to maintain an accurate spin phase. There will be times when the Sun is occulted and DSS data is not available. If this occurs at the start or end of a definitive attitude product, then the spin phase will be extrapolated using the mean rate determined by the EKF.

Sedlak, Joseph E.↗

Characterization of Delaminations and Transverse Matrix Cracks in Composite Laminates Using Multiple-Angle Ultrasonic Inspection

Delaminations and transverse matrix cracks often appear concurrently in composite laminates. Normal-incidence ultrasound is excellent at detecting delaminations, but is not optimum for matrix cracks. Non-normal incidence, or polar backscattering, has been shown to optimally detect matrix cracks oriented perpendicular to the ultrasonic plane of incidence. In this work, a series of six composite laminates containing slots were loaded in tension to achieve various levels of delamination and ply cracking. Ultrasonic backscattering was measured over a range of incident polar and azimuthal angles, in order to characterize the relative degree of damage of the two types. Sweptpolar- angle measurements were taken with a curved phased array, as a step toward an array-based approach to simultaneous measurement of combined flaws.

Johnston, Patrick H.↗

Electric Sail Propulsion for Exploring Nearby Interstellar Space

An Electric Sail is a revolutionary propellant-less propulsion system that is ideal for deep space missions to the outer planets, the Heliopause, and beyond. It is revolutionary in that it uses momentum exchange with the hypersonic solar wind to propel a spacecraft within the heliosphere. The momentum exchange is affected by the deflection of charged solar wind particles by an array of electrically biased wires that extend outward up to 30 km from a slowly rotating spacecraft. A high-voltage, positive bias on the wires, which are oriented normal to the solar wind flow, deflects the streaming protons, resulting in a reaction force on the wires that is also directed radially away from the sun. Over a period of months, this small force can accelerate the spacecraft to enormous speeds-on the order of 100-150 km/s (approximately 20 to 30 AU/yr). Unlike solar sails, Electric Sails do not rely on a fixed area to produce thrust. In fact, as they move away from the Sun and solar wind pressure decreases, the area for solar proton momentum transfer becomes larger, increasing system efficiency. As a result, thrust decreases at ≈1/r**(7/6) instead of the ≈1/r**2 rate typical for solar sails. The net effect is that an increased radial range of operation, together with increased thrust, both contribute to higher velocities and shorter total trip times to distant destinations. The MSFC Advanced Concepts Office (ACO) was awarded a Phase II NASA Innovative Advanced Concepts (NIAC) study to mature the technology for possible future demonstration and implementation. Preliminary results indicate that the physics of the system is viable and that a spacecraft propelled by an Electric Sail could reach the Heliopause in less than 15 years - and could be developed within a decade.

Johnson, Les↗

Magnetospheric Multiscale (MMS) Mission Attitude Ground System Design

This paper describes the attitude ground system (AGS) design to be used for support of the Magnetospheric MultiScale (MMS) mission. The AGS exists as one component of the mission operations control center. It has responsibility for validating the onboard attitude and accelerometer bias estimates, calibrating the attitude sensors and the spacecraft inertia tensor, and generating a definitive attitude history for use by the science teams. NASA's Goddard Space Flight Center (GSFC) in Greenbelt, Maryland is responsible for developing the MMS spacecraft, for the overall management of the MMS mission, and for mission operations. MMS is scheduled for launch in 2014 for a planned two-year mission. The MMS mission consists of four identical spacecraft flying in a tetrahedral formation in an eccentric Earth orbit. The relatively tight formation, ranging from 10 to 400 km, will provide coordinated observations giving insight into small-scale magnetic field reconnection processes. By varying the size of the tetrahedron and the orbital semi-major axis and eccentricity, and making use of the changing solar phase, this geometry allows for the study of both bow shock and magnetotail plasma physics, including acceleration, reconnection, and turbulence. The mission divides into two phases for science; these phases will have orbit dimensions of l.2xl2 Earth radii in the first phase and l.2x25 Earth radii in the second in order to study the dayside magnetopause and the nightside magnetotail, respectively. The orbital periods are roughly one day and three days for the two mission phases. Each of the four MMS spacecraft will be spin stabilized at 3 revolutions per minute (rpm), with the spin axis oriented near the ecliptic north pole but tipped approximately 2.5 deg towards the Sun line. The main body of each spacecraft will be an eight-sided platform with diameter of 3.4 m and height of 1.2 m. Several booms are attached to this central core: two axial booms of 14.9 m length, two radial magnetometer booms of 5 m length, and four radial -wire booms of 60 m length. Attitude and orbit control will use a set of axial and radial thrusters. A four-head star tracker and a slit-type digital Sun sensor (DSS) provide input for attitude determination. In addition, an accelerometer will be· used for closed-loop orbit maneuver control. The primary AGS product will be a daily definitive attitude history. Due to power limitations; the star tracker and accelerometer data will not be available at all times. However, tracker data from at least 10 percent of each orbit and continuous DSS data will be provided. An extended Kalman filter (EKF) will be used to estimate the three-axis attitude (i.e., spin axis orientation and spin phase) and rotation rate for all times when the tracker data is valid. For other times, the attitude is generated by assuming a constant angular momentum vector in the inertial frame. The DSS sun pulse will provide a timing signal to maintain an accurate spin phase. There will be times when the Sun is occulted and DSS data is not available. If this occurs at the start or end of a definitive attitude product, then the spin phase will be extrapolated using the mean rate determined by the EKF.

Sedlak, Joseph E.↗

An arbitrarily high-order three-dimensional Cartesian-grid method for reconstructing interfaces from volume fraction fields

Here Tthis work describes a newly developed, arbitrarily high-order Cartesian-grid method for reconstructing material interfaces from a volume fraction field. The method begins by identifying all of the grid cells in the volume fraction field that are intersected by the interface and need to be approximated by the reconstruction scheme. Finite-differences are used to calculate the gradient of the volume fraction field and provide an estimate of the surface normal in all of the interfacial grid cells. Groups of connected grid cells are then identified which all have the same dominant component of the normal vector. This grouping by orientation determines the proper dependent variable to use in the surface reconstruction (e.g. for a 2D curve, this step determines if the surface will be approximated by a function of x or y). A cumulative integral over the surface is constructed and fit using b-splines for two-dimensional problems or tensor-product b-splines for three-dimensional problems. This construction allows for the interface to be recovered through application of the second fundamental theorem of calculus. Fitting the cumulative integral with $\mathscr{N}$ th-order b-splines (or tensor-product b-splines) yields an ($\mathscr{N}$-1) th-order convergence rate of the interface shape. Differentiation of the b-spline interface function(s) allows for the high-order approximation of the normal vector and curvature to be obtained directly anywhere along b-spline. Together, the proposed reconstruction technique can achieve arbitrarily high mesh convergence rates. Validation tests are presented with mesh convergence rates ranging from fourth- to tenth-order.

97 MATHEMATICS AND COMPUTING↗