Search NASA⌕ Search

SEARCH · Search NASA

Results for “compact model”

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 487 records · Page 27
All sourcesNASA NTRSDOE OSTI
NASA NTRS1994

The radio emission from the ultraluminous far-infrared galaxy NGC 6240

We present new radio observations of the 'prototypical' ultraluminous far-infrared galaxy NGC 6240, obtained using the Very Large Array (VLA) at lambda = 20 cm in B-configuration and at lambda = 3.6 cm in A-configuration. These data, along with those from four previous VLA observations, are used to perform a comprehensive study of the radio emission from NGC 6240. Approximately 70% (approximately 3 x 10(exp 23) W/Hz) of the total radio power at 20 cm originates from the nuclear region (approximately less than 1.5 kpc), of which half is emitted by two unresolved (R approximately less than 36 pc) cores and half by a diffuse component. The radio spectrum of the nuclear emission is relatively flat (alpha approximately equals 0.6; S(sub nu) proportional to nu(exp -alpha). The supernova rate required to power the diffuse component is consistent with that predicted by the stellar evolution models of Rieke et al. (1985). If the radio emission from the two compact cores is powered by supernova remnants, then either the remnants overlap and form hot bubbles in the cores, or they are very young (approximately less than 100 yr.) Nearly all of the remaining 30% of the total radio power comes from an 'armlike' region extending westward from the nuclear region. The western arm emission has a steep spectrum (alpha approximately equals 1.0), suggestive of aging effects from synchrotron or inverse-Compton losses, and is not correlated with starlight; we suggest that it is synchrotron emission from a shell of material driven by a galactic superwind. Inverse Compton scattering of far-infrared photons in the radio sources is expected to produce an X-ray flux of approximately 2 - 6 x 10(exp -14) ergs/s/sq cm in the 2 - 10 keV band. No significant radio emission is detected from or near the possible ultramassive 'dark core'.

Colbert, Edward J. M.↗
NASA NTRS2009

The Morphology of Passively Evolving Galaxies at Z-2 from HST/WFC3 in the Hubble Ultra Deep Field

We discuss near-IR images of six passive galaxies (SSFR< 10(exp -2)/Gyr) at redshift 1.3 < z < 2.4 with stellar mass M approx 10(exp 11) solar mass, selected from the Great Observatories Origins Deep Survey (GOODS), obtained with WFC3/IR and the Hubble Space Telescope (HST). These WFC3 images provide the deepest and highest angular resolution view of the optical rest-frame morphology of such systems to date. We find that the light profile of these; galaxies is generally regular and well described by a Sersic model with index typical of today's spheroids. We confirm the existence of compact and massive early-type galaxies at z approx. 2: four out of six galaxies have T(sub e) approx. 1 kpc or less. The WFC3 images achieve limiting surface brightness mu approx. 26.5 mag/sq arcsec in the F160W bandpass; yet there is no evidence of a faint halo in the five compact galaxies of our sample, nor is a halo observed in their stacked image. We also find very weak "morphological k-correction" in the galaxies between the rest-frame UV (from the ACS z band), and the rest-frame optical (WFC3 H band): the visual classification, Sersic indices and physical sizes of these galaxies are independent or only mildly dependent on the wavelength, within the errors.

Cassata, P.↗
NASA NTRS2010

The Morphology of Passively Evolving Galaxies at Z approximately 2 from HST/WFC3 Deep Imaging in the Hubble Ultradeep Field

We present near-IR images of six passive galaxies (SSFR< 10(exp -2)/ Gyr) at redshift 1.3 < z < 2.4 with stellar mass M approximately 10(exp 11) solar M, selected from the Great Observatories Origins Deep Survey (GOODS), obtained with the Hubble Space Telescope (HST) and the WFC3/IR camera. These images provide the deepest and highest angular resolution view of the optical rest-frame morphology of such systems to date. We find that the light profile of these galaxies is generally regular and well described by a Sersic model with index typical of today's spheroids. We confirm the existence of compact and massive early-type galaxies at z approximately 2: four out of six galaxies have r(sub e) approximately 1 kpc or less. The images reach limiting surface brightness mu approximates 26.5 mag/square arcsec in the F160W bandpass; yet there is no evidence of a faint halo in the galaxies of our sample, even in their stacked image. We also find very weak "morphological k-correction" in the galaxies between the rest-frame UV (from the ACS z-band), and the rest-frame optical (WFC3 H-band): the visual classification, Sersic indices and physical sizes of these galaxies are independent or only mildly dependent on the wavelength, within the errors. The presence of an active nucleus is suspected in two out of six galaxies (33%), opening the intriguing possibility that a large, presently unaccounted population of AGN is hosted in these galaxies, possibly responsible for the cessation of star formation.

Cassata, P.↗
NASA NTRS2011

Thermal Properties and Evolution of a Coronal Cavity as Observed by the X-Ray Telescope on Hinode

Coronal cavities are voids in coronal emission often observed above high latitude filament channels. sometimes, these cavities have areas of bright x-ray emission in their centers (i.e. Hudson et al 1999). In this study, we use data from the x-ray Telescope (XRT) on Hinode to examine the thermal emission properties of a cavity observed during July 2008 that contains bright x-ray emission in its center. using ratios of XRT filters, we find evidence for elevated temperatures in the cavity center. The area of elevated temperature evolves from a ring-shaped structure at the be~inning of the observation, to an elongated structure two days later, finally appearing as a compact round source four days after the initial observation. We use a morphological model to fit the cavity emission, and find that a uniform structure running through the cavity does not fit the observations well. Instead, the observations are reproduced by modeling several short cavity "cores" with different parameters on different days. These changing core parameters may be due to some observed activity heating different parts of the cavity core at different times. We also find that core temperatures of 1.75 MK, 1.7 MK and 2.25 MK (for 19 July, 21 July and 23 July, respectively) in the model lead to structures that are consistent with the data, and that line-of-sight effects serve to lower the effective temperature derived from the filter ratio. page

Reeves, K. K.↗
NASA NTRS2011

Thermal Properties of A Solar Coronal Cavity Observed with the X-Ray Telescope on Hinode

Coronal cavities are voids in coronal emission often observed above high latitude filament channels. Sometimes, these cavities have areas of bright X-ray emission in their centers. In this study, we use data from the X-ray Telescope (XRT) on the Hinode satellite to examine the thermal emission properties of a cavity observed during July 2008 that contains bright X-ray emission in its center. Using ratios of XRT filters, we find evidence for elevated temperatures in the cavity center. The area of elevated temperature evolves from a ring-shaped structure at the beginning of the observation, to an elongated structure two days later, finally appearing as a compact round source four days after the initial observation. We use a morphological model to fit the cavity emission, and find that a uniform structure running through the cavity does not fit the observations well. Instead, the observations are reproduced by modeling several short cylindrical cavity "cores" with different parameters on different days. These changing core parameters may be due to some observed activity heating different parts of the cavity core at different times. We find that core temperatures of 1.75 MK, 1.7 MK and 2.0 MK (for July 19, July 21 and July 23, respectively) in the model lead to structures that are consistent with the data, and that line-of-sight effects serve to lower the effective temperature derived from the filter ratio.

Reeves, Katherine K.↗
NASA NTRS1999

Seasonal-to-Interannual Variability in Antarctic Sea-Ice Dynamics, and Its Impact on Surface Fluxes and Water Mass Production

Strong seasonal and interannual signals in Antarctic bottom-water outflow remain unexplained yet are highly correlated with anomalies in net sea-ice growth in coastal polynyas. The mechanisms responsible for driving salination and replenishment and rejuvenation of the dense shelf "source" waters likely also generate pulses of bottom water outflow. The objective of this research is to investigate time-scales of variability in the dynamics of sea-ice in the Southern Ocean in order to determine the primary sites for production of dense shelf waters. We are using a merged satellite/buoy sea-ice motion data set for the period 1978-present day to compute the dynamics of opening and closing of coastal polynyas over the continental shelf. The Ocean Circulation and Climate Advanced Model (OCCAM) ocean general circulation model with coupled sea-ice dynamics is presently forced using National Center for Environmental Prediction (NCEP) data to simulate fluxes and the salination impact of the ocean shelf regions. This work is relevant in the context of measuring the influence of polar sea-ice dynamics upon polar ocean characteristics, and thereby upon global thermohaline ocean circulation. Interannual variability in simulated net freezing rate in the Southern Weddell Sea is shown for the period 1986-1993. There is a pronounced maximum of ice production in 1988 and minimum in 1991 in response to anomalies in equatorward meridional wind velocity. This follows a similar approximate 8-year interannual cycle in Sea Surface Temperature (SST) and satellite-derived ice-edge anomalies reported elsewhere as the "Antarctic Circumpolar Wave." The amplitude of interannual fluctuations in annual net ice production are about 40% of the mean value, implying significant interannual variance in brine rejection and upper ocean heat loss. Southward anomalies in wind stress induce negative anomalies in open water production, which are observed in passive microwave satellite images. Thus, cycles of enhanced poleward wind stress reduce ice growth by compacting the ice along the coastline and closing open water in leads and polynyas. Model simulations confirm that years of low ice production, such as 1991, coincide with years of lower than normal bottom water outflow. Future plans include the assimilation of satellite ice concentrations and ice drift dynamics to more accurately constrain boundary conditions in the model.

Drinkwater, Mark R.↗
DOE OSTI2025

Supporting New Advanced Nuclear Technologies for Commercial-Maritime Applications

The ANS summary doesn't require an abstract, but I will produce one for the purpose of LRS: The large demand for maritime nuclear power underscores the need for experimental campaigns and modeling and simulation of new advanced reactors, which offer numerous advantages in terms of safety, efficiency, and compactness. INL, through the work conducted by NRIC and ABS, has addressed some of the technical, regulatory, and economic aspects of potential nuclear commercial maritime applications. However, on the technical side, there remain important physical phenomena, particularly for advanced reactors, that are not yet fully understood. Addressing these knowledge gaps requires a combination of experiments and advanced modeling and simulation techniques. INL possesses significant expertise in Multiphysics modeling and simulation. By collaborating with INL, the maritime nuclear sector can leverage this expertise to advance the development and deployment of innovative nuclear technologies.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS2002

Magnetohydrodynamic MACH Code Used to Simulate Magnetoplasmadynamic Thrusters

The On-Board Propulsion program at the NASA Glenn Research Center is utilizing a state of-the-art numerical simulation to model the performance of high-power electromagnetic plasma thrusters. Such thrusters are envisioned for use in lunar and Mars cargo transport, piloted interplanetary expeditions, and deep-space robotic exploration of the solar system. The experimental portion of this program is described in reference 1. This article describes the numerical modeling program used to guide the experimental research. The synergistic use of numerical simulations and experimental research has spurred the rapid advancement of high-power thruster technologies for a variety of bold new NASA missions. From its inception as a U.S. Department of Defense code in the mid-1980's, the Multiblock Arbitrary Coordinate Hydromagnetic (MACH) simulation tool has been used by the plasma physics community to model a diverse range of plasma problems--including plasma opening switches, inertial confinement fusion concepts, compact toroid formation and acceleration, z-pinch implosion physics, laser-target interactions, and a variety of plasma thrusters. The MACH2 code used at Glenn is a time-dependent, two-dimensional, axisymmetric, multimaterial code with a multiblock structure. MACH3, a more recent three-dimensional version of the code, is currently undergoing beta tests. The MACH computational mesh moves in an arbitrary Lagrangian-Eulerian (ALE) fashion that allows the simulation of diffusive-dominated and dispersive-dominated problems, and the mesh can be refined via a variety of adaptive schemes to capture regions of varying characteristic scale. The mass continuity and momentum equations model a compressible viscous fluid, and three energy equations are used to simulate nonthermal equilibrium between electrons, ions, and the radiation field. Magnetic fields are modeled by an induction equation that includes resistive diffusion, the Hall effect, and a thermal source for magnetic fields. Various models of plasma resistivity are included, along with ablation models and multiport circuit solvers. The set of equations is closed using either an ideal gas or real equation of state.

Mikellides, Pavlos G.↗
NASA NTRS1984

Spectral formation in compact X-ray sources

The gas flows responsible for X-ray spectral emission in compact X-ray sources are complex and likely to be optically thick. In galactic sources these flows include: stellar wind, accretion disk and corona, magnetosphere, and shock at the compact object surface. Active galactic nuclei are surrounded by many dense clouds responsible for the broad optical and UV emission lines. The nature of these flows can be discovered through observations of their X-ray spectra and time variability. Theoretical models for the formation of line and continuum spectra in these sources are described, emphasizing atomic photoabsorption, line emission, and Comptonization. Current observations of X-ray and UV spectra of compact galactic sources and active galactic nuclei are interpreted in terms of these models. Great advances in the understanding of these systems will result from future missions capable of observing their spectra at higher resolution. Observations of the complex of iron K-alpha lines at 6.4-6.7 keV will be particularly valuable.

Mccray, R.↗
NASA NTRS1986

The evolution of methods for noise prediction of high speed rotors and propellers in the time domain

Linear wave equation models which have been used over the years at NASA Langley for describing noise emissions from high speed rotating blades are summarized. The noise sources are assumed to lie on a moving surface, and analysis of the situation has been based on the Ffowcs Williams-Hawkings (FW-H) equation. Although the equation accounts for two surface and one volume source, the NASA analyses have considered only the surface terms. Several variations on the FW-H model are delineated for various types of applications, noting the computational benefits of removing the frequency dependence of the calculations. Formulations are also provided for compact and noncompact sources, and features of Long's subsonic integral equation and Farassat's high speed integral equation are discussed. The selection of subsonic or high speed models is dependent on the Mach number of the blade surface where the source is located.

Farassat, F.↗
NASA NTRS2007

The Chaotic Light Curves of Accreting Black Holes

We present model light curves for accreting Black Hole Candidates (BHC) based on a recently developed model of these sources. According to this model, the observed light curves and aperiodic variability of BHC are due to a series of soft photon injections at random (Poisson) intervals and the stochastic nature of the Comptonization process in converting these soft photons to the observed high energy radiation. The additional assumption of our model is that the Comptonization process takes place in an extended but non-uniform hot plasma corona surrounding the compact object. We compute the corresponding Power Spectral Densities (PSD), autocorrelation functions, time skewness of the light curves and time lags between the light curves of the sources at different photon energies and compare our results to observation. Our model reproduces the observed light curves well, in that it provides good fits to their overall morphology (as manifest by the autocorrelation and time skewness) and also to their PSDs and time lags, by producing most of the variability power at time scales 2 a few seconds, while at the same time allowing for shots of a few msec in duration, in accordance with observation. We suggest that refinement of this type of model along with spectral and phase lag information can be used to probe the structure of this class of high energy sources.

Kazanas, Demosthenes↗
DOE OSTI2025

Multiphysics Modeling of Microreactors with NEAMS codes, and Validation Based on KRUSTY Reactivity Insertion

The NEAMS Multiphysics Applications team continues to assess code usability and functionality for microreactor design and safety analyses, while demonstrating that NEAMS tools capture both steady-state and transient behavior across distinct microreactor concepts. In FY2025, the team advanced full-core, high-fidelity, multiphysics models that solve more complex problems and strengthen verification/validation for several microreactor systems: heat-pipe microreactor (HPMR), gas-cooled microreactor (GCMR), and the KRUSTY experiment. These models employ the MOOSE MultiApp/Transfers architecture with Griffin for neutronics, BISON for heat conduction/thermomechanics, Sockeye for heat pipes, SAM/THM for coolant channels and loops, and SWIFT for hydride behavior, with meshes generated via the MOOSE Reactor Module. The graphite models available in the Grizzly code were also investigated for future analyses. For the HPMR, a Na-HPMR variant was constructed to align with recently validated heat-pipe experiments and Sockeye’s LCVF capability, enabling mechanistic heat-pipe transients and startup modeling. The Na-HPMR will serve as the primary model for HPMR investigations in upcoming tasks. The load-following and single heat-pipe failure scenarios (Griffin/BISON/Sockeye), which were previously modeled for the K-HPMR, were replicated for the Na-HPMR, showing strong negative temperature feedback and highly localized thermal effects, respectively, while the startup case captured vapor-front progression and heat-removal activation. Solid mechanics was added to the previously built K-HPMR full-core model in BISON, showing minimal impact on steady-state reactivity yet enabling stress-field predictions that prepare the path for full-core TRISO performance analyses. For the GCMR, automated steady-state and four transient scenarios were executed using Griffin/BISON/SAM/SWIFT. Results confirm robust inherent safety: power collapses promptly in loss-of-cooling events, the inlet-temperature drop settles to a new equilibrium, and a single-channel blockage yields only a ~30 K local fuel-temperature rise with <0.4% power decrease. SWIFT-predicted hydrogen redistribution affects reactivity during both steady-state and transient conditions, underscoring its importance. A Brayton-cycle balance of plant (BOP) model in SAM/THM demonstrated stable startup behavior, and xenon-driven reactivity during load following was analyzed. To improve TRISO-compact temperature fidelity, a fast multiscale Heat Source Decomposition (HSD) treatment was implemented. Against heterogeneous benchmarks, HSD reduces underprediction of kernel temperatures and lowers predicted peak powers in reactivity-insertion transients compared to previous homogenized models. KRUSTY warm-critical validation progressed from FY2024 baselines: the 15Ȼ insertion shows excellent agreement in peak power (~2% high) and temperature trends, and the 30Ȼ case was automated via a feedback controller that maintained power near 3 kW for ~150 s with close agreement to data. The successful modeling of the warm critical tests has laid a strong foundation for simulating more complex nuclear system tests in the years ahead. Throughout FY2025, developer feedback was provided (e.g., MOOSE batch mesh generation, distributed pre-split meshes, Griffin sweeper on displaced meshes), several new models were contributed to the Virtual Test Bed, and an OECD-NEA WPRS multiphysics benchmark based on the HPMR was initiated to enable broader cross-comparison and best-practice development with the nuclear community at large.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
NASA NTRS1988

Compact imaging spectrometer for induced emissions

On the basis of spectral measurements made from the Space Shuttle and on models of the possible Space Station external environment, it appears likely that, even at the planned altitudes of Space Station, photon emissions will be induced. These emissions will occur to some degree throughout the ultraviolet-visible-infrared spectrum. The emissions arise from a combination of processes including gas phase collisions between relatively energetic ambient and surface emitted or re-emitted atoms or molecules, where the surface raises some species to excited energy states. At the present time it is not possible to model these processes or the anticipated intensity levels with any accuracy, as a number of fundamental parameters needed for such calculations are still poorly known or unknown. However, it is possible that certain spectral line and band features will exceed the desired goal that concomitant emissions not exceed the natural zodiacal background. Also, in the near infrared and infrared, it appears that this level will be exceeded to a significant degree. Therefore it will be necessary to monitor emission levels in the vicinity of Space Station, both in order to establish the levels and to better model the environment. A small spectrometer is briefly described which is suitable for monitoring the spectrum from 1200 A to less than or equal to 12,000 A. The instrument uses focal plane array detectors to image this full spectral range simultaneously. The spectral resolution is 4 to 12 A, depending on the portion of the wavelength range.

Torr, Douglas G.↗
DOE OSTI2024

Inverse Design of Photonic Surfaces via High throughput Femtosecond Laser Processing and Tandem Neural Networks

Abstract This work demonstrates a method to design photonic surfaces by combining femtosecond laser processing with the inverse design capabilities of tandem neural networks that directly link laser fabrication parameters to their resulting textured substrate optical properties. High throughput fabrication and characterization platforms are developed that generate a dataset comprising 35280 unique microtextured surfaces on stainless steel with corresponding measured spectral emissivities. The trained model utilizes the nonlinear one‐to‐many mapping between spectral emissivity and laser parameters. Consequently, it generates predominantly novel designs, which reproduce the full range of spectral emissivities (average root‐mean‐squared‐error < 2.5%) using only a compact region of laser parameter space 25 times smaller than what is represented in the training data. Finally, the inverse design model is experimentally validated on a thermophotovoltaic emitter design application. By synergizing laser‐matter interactions with neural network capabilities, the approach offers insights into accelerating the discovery of photonic surfaces, advancing energy harvesting technologies.

36 MATERIALS SCIENCE↗
DOE OSTI2026

A fractional calculus framework for open quantum dynamics: From Liouville to Lindblad to memory kernels

Open quantum systems exhibit dynamics ranging from unitary evolution to irreversible dissipation. While the Gorini–Kossakowski–Sudarshan–Lindblad equation uniquely characterizes Markovian completely positive and trace-preserving (CPTP) evolution, many physical platforms display non-Markovian features such as algebraic relaxation and coherence backflow. Fractional calculus provides a natural way to model such long-memory behavior through power-law temporal kernels introduced by fractional time derivatives. Here, we develop a unified framework that embeds fractional master equations within the broader hierarchy of open-system formalisms. The fractional equation forms a structured subclass of memory-kernel models, reduces to the Lindblad form at unit order, and, through Bochner–Phillips subordination, admits a CPTP representation as an average over Lindblad semigroups. Its resolvent structure further connects fractional dynamics to established non-Markovian approaches, including Nakajima–Zwanzig kernels and hierarchical equations of motion, providing a compact surrogate for long-memory effects. This formulation positions fractional calculus as a rigorous and practical language for modeling non-Markovian quantum dynamics in chemical physics and physical chemistry, providing a CPTP-preserving, computationally efficient surrogate for structured condensed-phase environments where long-time memory and dissipation play a central role.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Waveform modelling for the Laser Interferometer Space Antenna

LISA, the Laser Interferometer Space Antenna, will usher in a new era in gravitational-wave astronomy. As the first anticipated space-based gravitational-wave detector, it will expand our view to the millihertz gravitational-wave sky, where a spectacular variety of interesting new sources abound: from millions of ultra-compact binaries in our Galaxy, to mergers of massive black holes at cosmological distances; from the early inspirals of stellar-mass black holes that will ultimately venture into the ground-based detectors’ view to the death spiral of compact objects into massive black holes, and many sources in between. Central to realising LISA’s discovery potential are waveform models, the theoretical and phenomenological predictions of the pattern of gravitational waves that these sources emit. This White Paper is presented on behalf of the Waveform Working Group for the LISA Consortium. It provides a review of the current state of waveform models for LISA sources, and describes the significant challenges that must yet be overcome.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
NASA NTRS1998

Performance of Low Dissipative High Order Shock-Capturing Schemes for Shock-Turbulence Interactions

Accurate and efficient direct numerical simulation of turbulence in the presence of shock waves represents a significant challenge for numerical methods. The objective of this paper is to evaluate the performance of high order compact and non-compact central spatial differencing employing total variation diminishing (TVD) shock-capturing dissipations as characteristic based filters for two model problems combining shock wave and shear layer phenomena. A vortex pairing model evaluates the ability of the schemes to cope with shear layer instability and eddy shock waves, while a shock wave impingement on a spatially-evolving mixing layer model studies the accuracy of computation of vortices passing through a sequence of shock and expansion waves. A drastic increase in accuracy is observed if a suitable artificial compression formulation is applied to the TVD dissipations. With this modification to the filter step the fourth-order non-compact scheme shows improved results in comparison to second-order methods, while retaining the good shock resolution of the basic TVD scheme. For this characteristic based filter approach, however, the benefits of compact schemes or schemes with higher than fourth order are not sufficient to justify the higher complexity near the boundary and/or the additional computational cost.

Sandham, N. D.↗
NASA NTRS1998

Performance of Low Dissipative High Order Shock-Capturing Schemes for Shock-Turbulence Interactions

Abstract Accurate and efficient direct numerical simulation of turbulence in the presence of shock waves represents a significant challenge for numerical methods. The objective of this paper is to evaluate the performance of high order compact and non-compact central spatial differencing employing total variation diminishing (TVD) shock-capturing dissipations as characteristic based filters for two model problems combining shock wave and shear layer phenomena. A vortex pairing model evaluates the ability of the schemes to cope with shear layer instability and eddy shock waves, while a shock wave impingement on a spatially-evolving mixing layer model studies the accuracy of computation of vortices passing through a sequence of shock and expansion waves. A drastic increase in accuracy is observed if a suitable artificial compression formulation is applied to the TVD dissipations. With this modification to the filter step the fourth-order non-compact scheme shows improved results in comparison to second-order methods, while retaining the good shock resolution of the basic TVD scheme. For this characteristic based filter approach, however, the benefits of compact schemes or schemes with higher than fourth order are not sufficient to justify the higher complexity near the boundary and/or the additional computational cost.

Sandham, N. D.↗