SEARCH · Search NASA

Results for “numerical validation”

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 271 records · Page 15

All sources NASA NTRS DOE OSTI

NASA NTRS2016

Exhibition of Stochastic Resonance in Vestibular Perception

Astronauts experience sensorimotor changes during spaceflight, particularly during G-transitions. Post flight sensorimotor changes include spatial disorientation, along with postural and gait instability that may degrade operational capabilities of the astronauts and endanger the crew. A sensorimotor countermeasure that mitigates these effects would improve crewmember safety and decrease risk. The goal of this research is to investigate the potential use of stochastic vestibular stimulation (SVS) as a technology to improve sensorimotor function. We hypothesize that low levels of SVS will improve sensorimotor perception through the phenomenon of stochastic resonance (SR), when the response of a nonlinear system to a weak input signal is enhanced by the application of a particular nonzero level of noise. This study aims to advance the development of SVS as a potential countermeasure by 1) demonstrating the exhibition of stochastic resonance in vestibular perception, a vital component of sensorimotor function, 2) investigating the repeatability of SR exhibition, and 3) determining the relative contribution of the semicircular canals (SCC) and otolith (OTO) organs to vestibular perceptual SR. A constant current stimulator was used to deliver bilateral bipolar SVS via electrodes placed on each of the mastoid processes, as previously done. Vestibular perceptual motion recognition thresholds were measured using a 6-degree of freedom MOOG platform and a 150 trial 3-down/1-up staircase procedure. In the first test session, we measured vestibular perceptual thresholds in upright roll-tilt at 0.2 Hz (SCC+OTO) with SVS ranging from 0-700 μA. In a second test session a week later, we re-measured roll-tilt thresholds with 0, optimal (from test session 1), and 1500 μA SVS levels. A subset of these subjects, plus naive subjects, participated in two additional test sessions in which we measured thresholds in supine roll-rotation at 0.2 Hz (SCC) and upright y-translation at 1 Hz (OTO) with SVS up to 700 μA. A sinusoidal galvanic vestibular stimulation (GVS) perceptual threshold was also measured on each test day and used to normalize the SVS levels across subjects. In roll-tilt thresholds with SVS, the characteristic SR curve was qualitatively exhibited in 10 of 12 subjects, and the improvement in motion threshold was significant in 6 subjects, indicating that optimal SVS improved passive body motion perception in a way that is consistent with classical SR theory. A probabilistic comparison to numeric simulations further validated these experimental results. On the second test session, 4 out of the 10 SR exhibitors showed repeated improvement with SVS compared to the no SVS condition. Data collection is ongoing for the last two test sessions in which SCC and OTO only perceptual motion recognition thresholds are being measured with SVS. The final results of these test sessions will give insight into whether vestibular perceptual SR can occur when only one type of vestibular sensor is sensing motion or if it is more evident when sensory integration between the SCC and OTO is occurring during the motion. The overall purpose of this research is to further quantify the effects of SVS on various sensorimotor tasks and to gain a more fundamental understanding of how SVS causes SR in the vestibular system. In the context of human space flight, results from this research will help in understanding how SVS may be practically implemented in the future as a component of a comprehensive countermeasure plan for G-transition adaptation.

Galvan-Garza, R. C.↗

NASA NTRS2015

Structured Overlapping Grid Simulations of Contra-rotating Open Rotor Noise

Computational simulations using structured overlapping grids with the Launch Ascent and Vehicle Aerodynamics (LAVA) solver framework are presented for predicting tonal noise generated by a contra-rotating open rotor (CROR) propulsion system. A coupled Computational Fluid Dynamics (CFD) and Computational AeroAcoustics (CAA) numerical approach is applied. Three-dimensional time-accurate hybrid Reynolds Averaged Navier-Stokes/Large Eddy Simulation (RANS/LES) CFD simulations are performed in the inertial frame, including dynamic moving grids, using a higher-order accurate finite difference discretization on structured overlapping grids. A higher-order accurate free-stream preserving metric discretization with discrete enforcement of the Geometric Conservation Law (GCL) on moving curvilinear grids is used to create an accurate, efficient, and stable numerical scheme. The aeroacoustic analysis is based on a permeable surface Ffowcs Williams-Hawkings (FW-H) approach, evaluated in the frequency domain. A time-step sensitivity study was performed using only the forward row of blades to determine an adequate time-step. The numerical approach is validated against existing wind tunnel measurements.

Grid↗

NASA NTRS2016

Preparation of Morpheus Vehicle for Vacuum Environment Testing

The main objective for this summer 2016 tour was to prepare the Morpheus vehicle for its upcoming test inside Plum Brook's vacuum chamber at NASA John H. Glenn Research Center. My contributions towards this project were mostly analytical in nature, providing numerical models to validate test data, generating computer aided analyses for the structure support of the vehicle's engine, and designing a vacuum can that is to protect the high speed camera used during testing. Furthermore, I was also tasked with designing a tank toroidal spray bar system.

Sandoval, Armando↗

NASA NTRS2019

Solid Fuel Ignition and Extinction (SoFIE) Project on ISS

The Solid Fuel Ignition and Extinction (SoFIE) project studies ignition and flammability of solid spacecraft materials (fuels) in practical geometries and realistic atmospheric conditions. It is an experiment insert designed for use within the existing Combustion Integrated Rack (CIR). The CIR chamber provides a level of containment and permits testing at variable oxygen concentrations and pressures representative of current and planned NASA Space Exploration Atmospheres. The applications of SoFIE include: (1) Determining safer selection of cabin materials and validating NASA materials flammability selection using 1-g test protocols for low-gravity fires, (2) Improving understanding of early fire growth behavior, (3) Validating material flammability numerical models, (4) Determining optimal suppression techniques for burning materials by diluents, flow reduction, and venting, (5) Obtaining burning behavior of actual engineering materials planned for spacecraft, (6) Developing corresponding models of microgravity flame spread, flammability, and extinction, and use the results to improve normal gravity combustion models for terrestrial applications. The hardware permits a wide range of solid-material combustion and fire suppression studies. It supports multiple investigations using common infrastructure including sample holders, flow control, test sections, external radiant heaters, igniters, and diagnostics. SoFIE has been developed to meet the requirements of five unique investigations. It is currently being built and slated to begin operations on the ISS in July 2021. Given the general capabilities of the hardware insert, it is intended to be used as a facility for future researchers who can propose to NASA for related solid combustion studies.

Ferkul, Paul↗

NASA NTRS2000

Transitional mixing in a binary species supercritical, temporal mixing layer

The turbulent mixing of fluids at high pressure is a topic of much interest as it is relavant both to the natural phenomena and to technical applications. In the realm of technical applications, liquid rocket combustion presents a particular challenge as the operating conditions are supercritical with respect to both fuel and oxidizer, making mandatory the understanding of supercritical fluid behavior for a potentially explosive mixture. In this situation, numerical simulations with validated models can contribute information that would be otherwise impossible to obtain.

supercritical↗

NASA NTRS2021

The Scientific Legacy of NASA’s Operation IceBridge

The National Aeronautics and Space Administration (NASA)’s Operation IceBridge (OIB) was a 13-year (2009–2021) airborne mission to survey land and sea ice across the Arctic, Antarctic, and Alaska. Here, we review OIB’s goals, instruments, campaigns, key scientific results, and implications for future investigations of the cryosphere. OIB’s primary goal was to use airborne laser altimetry to bridge the gap in fine-resolution elevation measurements of ice from space between the conclusion of NASA’s Ice, Cloud, and land Elevation Satellite (ICESat; 2003–2009) and its follow-on, ICESat-2 (launched 2018). Additional scientific requirements were intended to contextualize observed elevation changes using a multisensor suite of radar sounders, gravimeters, magnetometers, and cameras. Using 15 different aircraft, OIB conducted 968 science flights, of which 42% were repeat surveys of land ice, 42% were surveys of previously unmapped terrain across the Greenland and Antarctic ice sheets, Arctic ice caps, and Alaskan glaciers, and 16% were surveys of sea ice. The combination of an expansive instrument suite and breadth of surveys enabled numerous fundamental advances in our understanding of the Earth’s cryosphere. For land ice, OIB dramatically improved knowledge of interannual outlet-glacier variability, ice-sheet, and outlet-glacier thicknesses, snowfall rates on ice sheets, fjord and sub-ice-shelf bathymetry, and ice-sheet hydrology. Unanticipated discoveries included a reliable method for constraining the thickness within difficult-to-sound incised troughs beneath ice sheets, the extent of the firn aquifer within the Greenland Ice Sheet, the vulnerability of many Greenland and Antarctic outlet glaciers to ocean-driven melting at their grounding zones, and the dominance of surface-melt-driven mass loss of Alaskan glaciers. For sea ice, OIB significantly advanced our understanding of spatiotemporal variability in sea ice freeboard and its snow cover, especially through combined analysis of fine-resolution altimetry, visible imagery, and snow radar measurements of the overlying snow thickness. Such analyses led to the unanticipated discovery of an interdecadal decrease in snow thickness on Arctic sea ice and numerous opportunities to validate sea ice freeboards from satellite radar altimetry. While many of its data sets have yet to be fully explored, OIB’s scientific legacy has already demonstrated the value of sustained investment in reliable airborne platforms, airborne instrument development, interagency and international collaboration, and open and rapid data access to advance our understanding of Earth’s remote polar regions and their role in the Earth system.

Joseph A. MacGregor↗

NASA NTRSDate not supplied

Digitizing Named Entities Found Within Letters of Agreement

Letters of Agreement (LOAs) are text-based air traffic control documents that contain procedures and actions agreed upon by the different parties, typically two or more FAA facilities, that are subject to an agreement. The documents contain among other things generic constraints, which are explicit and implicit combinations of procedures that limit a flight’s trajectory and affects pilot actions. For example, a controller may be required, to assign a specific altitude to an aircraft crossing the boundary between two airspaces. Although LOA generic constraints directly impact the trajectory of an aircraft, they are not currently available in a digital form that can be used for (or directly ingested into automated) flight planning. Instead, the constraints are manually input into an onboard or ground based system. LOA documents are primarily stored at a controlling facility and the generic constraints are implemented by experienced air traffic controllers and pilots primarily using voice instructions. This increases the workload of the controllers, likelihood of error (e.g., due to noisy communication) and makes it impractical for implementation with unmanned aircraft. Therefore, steps must be taken to make existing constraints machine interpretable to enable e.g., automated handoffs which in turn would reduce controller workload. With recent advances in natural language processing, especially the rise in digitization of text documents (e.g., medical documents) and automated extraction of information therein, it is now possible to extract flight specific constraints from LOAs. The goal of this work is to digitize named entities through a combination of natural language processing tasks: named entity disambiguation, toponym resolution, and numeric parsing to extract general constraint components contained within LOAs, herein referred to as Entity Enhancement (EE). Starting with a small list of named entities (e.g., ARTCC, Tower, Altitude and Speed), EE can extract the named entities while simultaneously converting the string-based output into a digital format using an ensemble of processes like rule-based gazetteers and syntactic-lexical patterns. The digital format contains a diverse set of information based on the entity label in question, ranging from standardized facility names to units of measure (e.g., feet) and other numeric information. Upon validating our approach using a truth dataset, we show an overall F1-Score of 0.71 for the extraction process. Looking beyond entity enhancement, we are also working towards the goal of completely digitizing the general constraints by performing EE and fitting them into a standardized exchange model (XM) such as the Aeronautical Information Exchange Model (AIXM). This will allow for easy distribution and dissemination of LOA constraints to air users, better searchability within documents, and enable ingestion into automated flight planning. Finally, we show a preliminary version of the proposed XM architecture and demonstrate how the model can be populated from the EE output.

Stephen S. B. Clarke↗

NASA NTRSDate not supplied

Artemis I Flight Instrumentation Data Quality Assessment and Processing

This paper is in support of the SciTech 2024 Space Launch System Aerosciences Special Sessions being organized by Brent Pomeroy and Jeremy Pinier. On November 16th, 2022, NASA launched the inaugural test flight of the Space Launch System (SLS) carrying the Orion capsule into a high orbit far beyond the Moon. The launch vehicle was instrumented with over three thousand flight instrumentation sensors, which monitored aerodynamic, acoustic, structural, and thermal environments. These data are intended to validate experimental and numerical tools used to predict the design environments which the vehicle experiences during launch and ascent. Prior to launch, tests were performed at the SLS Systems Integration Laboratory (SIL) using flight-like avionics and on the integrated flight hardware of Artemis I at the Vehicle Assembly Building (VAB). The purpose of these tests was to characterize the data acquisition units (DAUs) used to record and telemeter flight data to ground stations in order to assure that flight test objectives can be achieved and to quantify the expected quality of the flight data. In addition, pre-flight assessment and development of tools and methods used to process and disseminate flight data at the Huntsville Operations Support Center (HOSC) were conducted and adjustments made with respect to DAU time-synchronization prior to and after the flight. This paper summarizes these tests and some aspects of the post processing of data are discussed.

Developmental Flight Instrumentation↗

NASA NTRSDate not supplied

Artemis I Flight Instrumentation Data Quality Assessment and Processing

Developmental Flight Instrumentation↗

NASA NTRSDate not supplied

Validation of Computational Fluid Dynamics Model for Pressurized Expulsion of Liquid Methane using Non-Condensable Pressurant

An important aspect of cryogenic fluid management is the design of the pressurization systems used to maintain and regulate the tank pressure as propellant is expelled from the tanks. To design these systems efficiently and effectively it is important to be able to predict the amount pressurant gas needed and how the pressurant gas affects the propellant as it is added in the tank. Computational fluid dynamics (CFD) modeling is a common tool that is used to help with these predictions, so it is crucial to have CFD models that are validated by comparing numerical results to experimental data. In the 1970’s extensive testing was done at NASA Plum Brook K-Site Facility exploring the pressurant gas requirements needed to maintain tank pressure while liquid methane is discharged from a spherical tank. These tests explored using different pressurant gases (methane, helium, hydrogen, and nitrogen) under static and sloshing conditions. The data gathered from the experiments includes the gas and wall temperatures at the end of the expulsion, the gas concentrations, and the total amount of pressurant added. This study presents a CFD model that accurately predicts similar results to the experimental data specifically for the tests in which helium, a non-condensable gas, was used. This comparison anchors the commercial CFD code Flow3D model and gives confidence for future use in design and analysis of pressurization systems for cryogenic propellant tanks

computational fluid dynamics↗

NASA NTRSDate not supplied

Numerical Investigation of Heat Transfer and Fluid Flow within Electrochemical Hydrogen Peroxide Generation Unit

Long-term manned space missions require the onboard production of disinfectants essential for maintaining crew health and supporting life systems. Currently, disinfection aboard the International Space Station (ISS) relies on disposable wetted wipes, which are regularly resupplied from Earth. This approach imposes a significant burden on resupply logistics, storage, and waste management. To address these challenges and support future missions, efforts are underway to develop an in-situ solution that electrochemically generates hydrogen peroxide disinfectant using onboard resources. In collaboration with NASA, Faraday Technology, Inc. has advanced this concept through a series of Small Business Innovation Research (SBIR) projects, resulting in the development of a Peroxide Generation Unit (PGU). The PGU can produce up to 3 wt.% hydrogen peroxide on-demand at a rate of 1 liter per day, providing a sustainable alternative to Earth-dependent supplies. The resulting aqueous hydrogen peroxide (H₂O₂) is an effective disinfectant, safe for crew use, compatible with spacecraft systems, and free from volatiles, off-gassing, or residues. This innovation offers a reliable, efficient solution for onboard disinfection, reducing dependence on Earth-based resupply while ensuring the health and safety of space crews. Generating hydrogen peroxide at the required rate needs high voltages and currents, exceeding 20V and 2A respectively, which leads to significant heat generation from Joule heating. This temperature rise poses a risk to sensitive system components, especially critical and expensive membranes that can degrade under thermal stress. To mitigate this risk, the thermal, fluid, and electrical flows within the system are modeled computationally using the commercial software COMSOL. The numerical simulations are validated against experimental data from both sub-scale and alpha-scale systems. Once verified, the model is employed to identify thermal hotspots, investigate their underlying causes, and explore solutions to prevent them.

Life Support Systems (LSS)↗

NASA NTRSDate not supplied

Numerical Predictions of the Flow and Heat Transfer Characteristics in the Film Boiling Regime During Tube Quenching

Cryogenic fluid management plays a major role in refueling of spacecrafts while in space for NASA’s future human space exploration missions. Due to the low boiling points of cryogens, storage, transport and handling of these fluids becomes difficult and may result in inefficient operation of the space propulsion systems. For refueling applications in space, the cryogenic fluids have to be transported across different locations and hence, the transfer of cryogenic fluids through pipes become critical. The cryogenic chill-down process is characterized by different regimes of flow boiling, viz., film boiling, transition boiling and nucleate boiling. The prediction of these regimes in a single CFD framework available in the literature is challenging and the present work attempts to address this challenge by initially modeling the film boiling regime accurately and to incorporate an user-defined function for transition and nucleate boiling at a later stage. Hence, the aim of the present work is to numerically model and validate the film boiling regime of the chilldown curve for liquid nitrogen experiments available in the literature. The validations are carried out at different inlet mass fluxes to have a robust simulation methodology. A dispersed mixture model is used to predict the vapor-liquid interface dynamics with the phase change phenomena modeled using the Lee model.

line chilldown↗

NASA NTRSDate not supplied

Validation Experiments of Incipient Turbulent Separation over an Axisymmetric Afterbody

Historically, the flow physics involved with most turbulent separated flows have presented fundamental challenges in validation between experimental and numerical approaches. As recognized by the CFD Vision 2030 study commissioned by NASA, validation of Reynolds-averaged Navier-Stokes (RANS) models and other scale-resolving methods for turbulent separated flow requires the support of advanced, high-fidelity experiments designed specifically for CFD implementation. In accordance with this effort, a new test platform, referred to as the NASA Axisymmetric Afterbody, was designed to obtain detailed information of the flow field over a smooth body, undergoing adverse pressure gradient induced separation for a fixed Reynolds number, Re = 180,000. The parametric body offers a range of flow states from fully attached,incipient separation, to fully separated flow based on variable afterbody geometries. In an initial effort to evaluate RANS turbulence model capabilities, the present configuration of the axisymmetric model features a mild adverse pressure gradient over the contoured boattail section, inducing incipient turbulent separation. Experiments include steady pressure measurements to serve as preliminary comparisons to simulation studies, which examine the effect of variable grid domains and RANS turbulence models, in an effort to understand and evaluate the critical variability between solutions for the present configuration. Results indicate potential discrepancies may be due to the effect of the square tunnel test section walls, relatively large blockage ratio, and inconsistent reference parameters. Ongoing work will focus on the experimental campaigns to obtain surface flow visualizations and high-resolution, off-body flow field measurements using Stereoscopic Particle Image Velocimetry (SPIV) and Laser Doppler Velocimetry (LDV) techniques to provide a detailed benchmark dataset to aid turbulence modelers.

axisymmetric afterbody↗

NASA NTRS1987

Numerical simulation of shock wave diffraction by TVD schemes

An upwind total variation diminishing (TVD) scheme and a predictor-corrector symmetric TVD scheme were used to numerically simulate the blast wave diffraction on a stationary object. The objective is to help design an optimum configuration so that lateral motion is minimized and at the same time vortex shedding and flow separation are reduced during a blast wave encounter. Results are presented for a generic configuration for both a coarse grid and a fine grid to illustrate the global and local diffraction flow fields. Numerical experiments for the shock wave reflection on a wedge are also included to validate the current approach. Numerical study indicated that these TVD schemes are more stable and produced higher shock resolution than classical shock capturing methods such as the explicit MacCormack scheme.

Young, Victor Y. C.↗

NASA NTRS2020

Computational Simulations of a Mach 0.745 Transonic Truss-Braced Wing Design

A joint effort between the NASA Ames and Langley Research Centers was undertaken to analyze the Mach 0.745 variant of the Boeing Transonic Truss-Braced Wing (TTBW) Design. Two different flow solvers, LAVA and USM3D, were used to predict the TTBW flight performance. Sensitivity studies related to mesh resolution and numerical schemes were conducted to define best practices for this type of geometry and flow regime. Validation efforts compared the numerical simulation results of various modeling methods against experimental data taken from the NASA Ames 11-foot Unitary Wind Tunnel experimental data. The fidelity of the computational representation of the wind tunnel experiment, such as utilizing a porous wall boundary condition to model the ventilated test section, was varied to examine how different tunnel effects influence CFD predictions. LAVA and USM3D results both show an approximate 0.5 angle of attack shift from experimental lift curve data. This drove an investigation that revealed that the trailing edge of the experimental model was rounded in comparison to the CAD model, due to manufacturing tolerances, which had not been accounted for in the initial simulations of the experiment. Simulating the TTBW with an approximation of this rounded trailing-edge reduces error by approximately 60%. An accurate representation of the tested TTBW geometry, ideally including any wing twists and deflections experienced during the test under various loading conditions, will be necessary for proper validation of the CFD.

Maldonado, Daniel↗

NASA NTRS2020

Computational Simulations of a Mach 0.745 Transonic Truss-Braced Wing Design

A joint effort between the NASA Ames and Langley Research Centers was undertaken to analyze the Mach 0.745 variant of the Boeing Transonic Truss-Braced Wing (TTBW) Design. Two different flow solvers, LAVA and USM3D, were used to predict the TTBW flight performance. Sensitivity studies related to mesh resolution and numerical schemes were conducted to define best practices for this type of geometry and flow regime. Validation efforts compared the numerical simulation results of various modeling methods against experimental data taken from the NASA Ames 11-foot Unitary Wind Tunnel experimental data. The fidelity of the computational representation of the wind tunnel experiment, such as utilizing a porous wall boundary condition to model the ventilated test section, was varied to examine how different tunnel effects influence CFD predictions. LAVA and USM3D results both show an approximate 0.5o angle of attack shift from experimental lift curve data. This drove an investigation that revealed that the trailing edge of the experimental model was rounded in comparison to the CAD model, due to manufacturing tolerances, which had not been accounted for in the initial simulations of the experiment. Simulating the TTBW with an approximation of this rounded trailing-edge reduces error by approximately 60%. An accurate representation of the tested TTBW geometry, ideally including any wing twists and deflections experienced during the test under various loading conditions, will be necessary for more thorough validation of the CFD.

Maldonado, Daniel↗

NASA NTRS1993

A systems approach to a DSMC calculation of a control jet interaction experiment

This study deals with the development of a methodology for numerically simulating the interaction of a reaction control system (RCS) jet with a low-density external flow. A European Space Agency (ESA) experiment was chosen as a test case since it provided experimental data that could validate some of the numerical results. The initial approach was to focus on several subproblems having direct relevance to the full interaction problem. This enabled different numerical methods to be separately investigated and validated for each part of the interaction problem. In this manner, the best methodology for solving the full interaction problem was developed. The subproblems considered in this study included typical RCS nozzle and plume flows, a flat plate at zero incidence, and the flow past the experimental test model without the control jet firing. Once these calculations were completed, a simulation was performed with a control jet operating at a relatively low density. The results from this latter simulation provided qualitative insight into the complex interaction process.

Tartabini, Paul V.↗

NASA NTRS1992

Validation of three-dimensional incompressible spatial direct numerical simulation code: A comparison with linear stability and parabolic stability equation theories for boundary-layer transition on a flat plate

Spatially evolving instabilities in a boundary layer on a flat plate are computed by direct numerical simulation (DNS) of the incompressible Navier-Stokes equations. In a truncated physical domain, a nonstaggered mesh is used for the grid. A Chebyshev-collocation method is used normal to the wall; finite difference and compact difference methods are used in the streamwise direction; and a Fourier series is used in the spanwise direction. For time stepping, implicit Crank-Nicolson and explicit Runge-Kutta schemes are used to the time-splitting method. The influence-matrix technique is used to solve the pressure equation. At the outflow boundary, the buffer-domain technique is used to prevent convective wave reflection or upstream propagation of information from the boundary. Results of the DNS are compared with those from both linear stability theory (LST) and parabolized stability equation (PSE) theory. Computed disturbance amplitudes and phases are in very good agreement with those of LST (for small inflow disturbance amplitudes). A measure of the sensitivity of the inflow condition is demonstrated with both LST and PSE theory used to approximate inflows. Although the DNS numerics are very different than those of PSE theory, the results are in good agreement. A small discrepancy in the results that does occur is likely a result of the variation in PSE boundary condition treatment in the far field. Finally, a small-amplitude wave triad is forced at the inflow, and simulation results are compared with those of LST. Again, very good agreement is found between DNS and LST results for the 3-D simulations, the implication being that the disturbance amplitudes are sufficiently small that nonlinear interactions are negligible.

Joslin, Ronald D.↗