Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Closure 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 649 records · Page 36
All sourcesNASA NTRSDOE OSTI
NASA NTRS2017

Interaction of Delaminations and Matrix Cracks in a CFRP Plate, Part II: Simulation Using an Enriched Shell Finite Element Model

Numerical simulations are presented of a recently developed test which creates multiple delaminations in a CFRP laminate specimen that grow and interact via transverse matrix cracks. A novel shell element enriched with the Floating Node Method, and a damage algorithm based on the Virtual Crack Closure Technique, were used to successfully simulate the tests. Additionally, a 3D high mesh fidelity model based on cohesive zones and continuum damage mechanics was used to simulate the tests and act as a representative of other similar state-of-the-art high mesh fidelity modeling techniques to compare to the enriched shell element. The enriched shell and high mesh fidelity models had similar levels of accuracy and generally matched the experimental data. With runtimes of 36 minutes for the shell model and 55 hours for the high mesh fidelity model, the shell model is 92 times faster than the high- fidelity simulation.

Mark W McElroy↗
NASA NTRS1987

The role of turbulent convection in the primitive solar nebula. I - Theory. II - Results

A model of convective turbulence which takes radiative dissipation, rotation, and convective motion anisotropy into account, on the basis of a closure for the nonlinear interactions that employs the growth rates of hydrodynamic instabilities, is used to obtain a theoretical framework for modeling the primordial solar nebula. It is assumed that convection is the sole source of turbulence causing the solar nebula to evolve. Vertical structure equations in the thin disk approximation are developed and a detailed comparison with the previous solar nebula convective models of such workers as Lin et al. (1981, 1982) is undertaken. The present values for the turbulent efficiency are much lower and more sensitive to opacity and surface density, resulting in low turbulent speeds, a more massive disk, a lower accretion rate 'best value', and a longer characteristic dispersal time for the disk. It is concluded that convection may not be the dominant source of turbulence needed to evolve young solar/stellar nebulae.

Cabot, W.↗
NASA NTRS1996

[?]Nonlinear Issues in the Aerothermochemistry of Gases and Materials and the Associated Physics and Dynamics of Interfaces

Our research and technology are focused on nonlinear issues in the aerothermochemistry of gases and materials and the associated physics and dynamics of interfaces. Our program is now organized to aggressively support the NASA Aeronautics Enterprise so as to: (a) develop a new generation of environmentally compatible, economic subsonic aircraft; (b) develop the technology base for an economically viable and environmentally compatible high-speed civil transport; (c) develop the technology options for new capabilities in high-performance aircraft; (d) develop hypersonic technologies for air-breathing flight; and (e) develop advanced concepts, understanding of physical phenomena, and theoretical, experimental, and computational tools for advanced aerospace systems. The implications from our research for aeronautical and aerospace technology have been both broad and deep. For example, using advanced computational techniques, we have determined exact solutions for the Schrodinger equation in electron-molecule scattering allowing us to evaluate atmospheric models important to reentry physics. We have also found a new class of exact solutions for the Navier Stokes equations. In experimental fluid dynamics, we have found explicit evidence of turbulence modification of droplet sizes in shock tube flow with condensation. We have developed a new diagnostic tool for the direct estimation of flow velocities at MHz sampling rates in quasi-one dimensional turbulent flow. This procedure suggests an unexpected confirmation of the possibility of 'natural' closure in Reynolds stresses with deep implications for the development of turbulent models. A transient increase is observed in both the spectral energy decay rate and the degree of chaotic complexity at the interface of a shock wave and a turbulent ionized gas. Even though the gas is apparently brought to rest by the shock wave, no evidence is found either of the expected relaminarization. A unique diamond-shaped nozzle has been designed for a detailed investigation of the effect of significant streamwise vorticity on the acoustic and IR characteristics of supersonic jets. Our results provide convincing evidence of the significant effect of vorticity on the far-field noise for the diamond jet as compared to the conventional round jets. We have found that the countercurrent shear layer mixes much more efficiently than conventional coflowing shear layers. We also developed the fluid thrust vectoring procedures which use counter flow to vector a jet. Our materials research has shown that the steep stress gradients at the fiber-matrix interface could be the primary cause of interface cracks after the processing of metallic and intermetallic matrix composites. New techniques have been evolved for: the microcharacterization of materials including microplastic strain and, point by point, the misorientation and plasticity for matrix composites; thermally induced stress measurements and load relaxation; the growth and characterization of metallic matrix composite interfaces; and for the growth of ferrite materials by pulsed laser deposition. The FAMU commitment to the HBCU Research Center also continues to be broad and deep.

Johnson, Joseph A., III↗
NASA NTRS1993

Turbulence modelling of flow fields in thrust chambers

Following the consensus of a workshop in Turbulence Modelling for Liquid Rocket Thrust Chambers, the current effort was undertaken to study the effects of second-order closure on the predictions of thermochemical flow fields. To reduce the instability and computational intensity of the full second-order Reynolds Stress Model, an Algebraic Stress Model (ASM) coupled with a two-layer near wall treatment was developed. Various test problems, including the compressible boundary layer with adiabatic and cooled walls, recirculating flows, swirling flows, and the entire SSME nozzle flow were studied to assess the performance of the current model. Detailed calculations for the SSME exit wall flow around the nozzle manifold were executed. As to the overall flow predictions, the ASM removes another assumption for appropriate comparison with experimental data to account for the non-isotropic turbulence effects.

Chen, C. P.↗
NASA NTRS1993

Development of an algebraic stress/two-layer model for calculating thrust chamber flow fields

Following the consensus of a workshop in Turbulence Modeling for Liquid Rocket Thrust Chambers, the current effort was undertaken to study the effects of second-order closure on the predictions of thermochemical flow fields. To reduce the instability and computational intensity of the full second-order Reynolds Stress Model, an Algebraic Stress Model (ASM) coupled with a two-layer near wall treatment was developed. Various test problems, including the compressible boundary layer with adiabatic and cooled walls, recirculating flows, swirling flows and the entire SSME nozzle flow were studied to assess the performance of the current model. Detailed calculations for the SSME exit wall flow around the nozzle manifold were executed. As to the overall flow predictions, the ASM removes another assumption for appropriate comparison with experimental data, to account for the non-isotropic turbulence effects.

Chen, C. P.↗
NASA NTRS2003

Macroscopic Models of Radiative Transfer as Applied to Computation of the Radiation Field in the Solar Atmosphere

In this paper a new three dimensional half-moment model for radiative transfer is presented for a gray medium. It describes the evolution of the zeroth and first directional half moments of the radiative intensity. The closure is provided, similarly to Dubroca and Klar, by the maximum entropy concept. This work generalizes that model to three dimensions. The model presented here (the derivation being done in Ripoll and Wray, called the M(sup 1/2)(sub 1) model, is a hyperbolic system consisting of a total of eight equations in three dimensions, four equations for each direction. Each half model has the classical form of a macroscopic moment model in which the pressure tensor is constructed from the well-known Eddington tensor with a particular Eddington factor. Moreover, different source and border terms occur. The latter introduce couplings between the macroscopic and microscopic quantities and between the + and - streams, through the intensity in the plane perpendicular to the flux. The main theoretical application of the half moment model, treated in this paper, is its reduction to a full moment model, called M(sup +)(sub 1), for the particular but important case of a hot, opaque source radiating in a cold transparent (or semi-transparent) medium for very specific applications, such as stellar interiors or atmospheres, or combustion problems. The structure of the paper is as follows. In section 2, the model M(sup 1/2)(sub 1) is presented. In section 3, for the particular case of a hot, opaque source radiating into a cold medium, the half moment model is reduced to the M(sup +)(sub 1) model. In section 4, we first solve a simple and academic problem to validate the models, followed by a simplified solar atmosphere.

Ripoll, J.-F.↗
NASA NTRS2003

Modeling and Control of Airport Queueing Dynamics under Severe Flow Restrictions

Based on field observations and interviews with controllers at BOS and EWR, we identify the closure of local departure fixes as the most severe class of airport departure restrictions. A set of simple queueing dynamics and traffic rules are developed to model departure traffic under such restrictions. The validity of the proposed model is tested via Monte Carlo simulation against 10 hours of actual operations data collected during a case-study at EWR on June 29,2000. In general, the model successfully reproduces the aggregate departure congestion. An analysis of the average error over 40 simulation runs indicates that flow-rate restrictions also significantly impact departure traffic; work is underway to capture these effects. Several applications and what-if scenarios are discussed for future evaluation using the calibrated model.

Carr, Francis↗
NASA NTRS1992

Application of fiber bridging models to fatigue crack growth in unidirectional titanium matrix composites

Several fiber bridging models were reviewed and applied to study the matrix fatigue crack growth behavior in center notched (0)(sub 8) SCS-6/Ti-15-3 and (0)(sub 4) SCS-6/Ti-6Al-4V laminates. Observations revealed that fatigue damage consisted primarily of matrix cracks and fiber matrix interfacial failure in the (0)(sub 8) SCS-6/Ti-15-3 laminates. Fiber-matrix interface failure included fracture of the brittle reaction zone and cracking between the two carbon rich fiber coatings. Intact fibers in the wake of the matrix cracks reduce the stress intensity factor range. Thus, an applied stress intensity factor range is inappropriate to characterize matrix crack growth behavior. Fiber bridging models were used to determine the matrix stress intensity factor range in titanium metal matrix composites. In these models, the fibers in the wake of the crack are idealized as a closure pressure. An unknown constant frictional shear stress is assumed to act along the debond or slip length of the bridging fibers. The frictional shear stress was used as a curve fitting parameter to available data (crack growth data, crack opening displacement data, and debond length data). Large variations in the frictional shear stress required to fit the experimental data indicate that the fiber bridging models in their present form lack predictive capabilities. However, these models provide an efficient and relatively simple engineering method for conducting parametric studies of the matrix growth behavior based on constituent properties.

Bakuckas, J. G., Jr.↗
NASA NTRS1994

National facilities study. Volume 4: Space operations facilities task group

The principal objectives of the National Facilities Study (NFS) were to: (1) determine where U.S. facilities do not meet national aerospace needs; (2) define new facilities required to make U.S. capabilities 'world class' where such improvements are in the national interest; (3) define where consolidation and phase-out of existing facilities is appropriate; and (4) develop a long-term national plan for world-class facility acquisition and shared usage. The Space Operations Facilities Task Group defined discrete tasks to accomplish the above objectives within the scope of the study. An assessment of national space operations facilities was conducted to determine the nation's capability to meet the requirements of space operations during the next 30 years. The mission model used in the study to define facility requirements is described in Volume 3. Based on this model, the major focus of the Task Group was to identify any substantive overlap or underutilization of space operations facilities and to identify any facility shortfalls that would necessitate facility upgrades or new facilities. The focus of this initial study was directed toward facility recommendations related to consolidations, closures, enhancements, and upgrades considered necessary to efficiently and effectively support the baseline requirements model. Activities related to identifying facility needs or recommendations for enhancing U.S. international competitiveness and achieving world-class capability, where appropriate, were deferred to a subsequent study phase.

Source record↗
NASA NTRSDate not supplied

Progressive Damage Failure Analysis of Post-Buckled Composite Single-Stringer Panel with Teflon Inserts

The high strength-to-weight and stiffness-to-weight ratio materials, such as laminated composites, are advantageous for modern aircraft. Laminated composites with initial flaws are susceptible to delamination under buckling loads. PDA tools help enhance the industry’s understanding of the mechanisms for damage initiation and growth in composite structures while assisting in the design, analysis, and sustainment methods of these composite structures. The global-local modeling approach for the single-stringer post-buckled panel was evaluated through this effort, using Teflon inserts to simulate the defect of damage during manufacturing. This understanding is essential for designing the post-buckled structure, reducing weight while predicting damage initiation location, and addressing a potential design review for future aircraft repairs. In this work, the initial damage was captured with Teflon inserts as the starting configuration; and any reference to the damage initiation refers to any damage beyond the “initial unbonded region.” The effort aims to develop, evaluate, and enhance methods to predict damage initiation and progression and the failure of post-buckled hat-stiffened panels using multiple Abaqus FEA Virtual Crack Closure Technique (VCCT) definitions. Validation of the PDA using the VCCT material model was performed on a large single-stringer panel subjected to compressive loading. The compressive loading of the panel caused the skin to buckle before any damage began to occur locally. In addition, comparisons are made for critical aspects of the damage morphology, such as a growth pattern that included delamination from the skin-stiffener interface to the skin and ply interfaces. When compared against the experimental data produced through the NASA Advanced Composites Project (ACP), the present model captured damage migration from one surface to another, and model validations were within 4% of the experimental data.

PDA↗
NASA NTRS1996

Aspirating Seal Development: Analytical Modeling and Seal Test Rig

This effort is to develop large diameter (22 - 36 inch) Aspirating Seals for application in aircraft engines. Stein Seal Co. will be fabricating the 36-inch seal(s) for testing. GE's task is to establish a thorough understanding of the operation of Aspirating Seals through analytical modeling and full-scale testing. The two primary objectives of this project are to develop the analytical models of the aspirating seal system, to upgrade using GE's funds, GE's 50-inch seal test rig for testing the Aspirating Seal (back-to-back with a corresponding brush seal), test the aspirating seal(s) for seal closure, tracking and maneuver transients (tilt) at operating pressures and temperatures, and validate the analytical model. The objective of the analytical model development is to evaluate the transient and steady-state dynamic performance characteristics of the seal designed by Stein. The transient dynamic model uses a multi-body system approach: the Stator, Seal face and the rotor are treated as individual bodies with relative degrees of freedom. Initially, the thirty-six springs are represented as a single one trying to keep open the aspirating face. Stops (Contact elements) are provided between the stator and the seal (to compensate the preload in the fully-open position) and between the rotor face and Seal face (to detect rub). The secondary seal is considered as part of the stator. The film's load, damping and stiffness characteristics as functions of pressure and clearance are evaluated using a separate (NASA) code GFACE. Initially, a laminar flow theory is used. Special two-dimensional interpolation routines are written to establish exact film load and damping values at each integration time step. Additionally, other user-routines are written to read-in actual pressure, rpm, stator-growth and rotor growth data and, later, to transfer these as appropriate loads/motions in the system-dynamic model. The transient dynamic model evaluates the various motions, clearances and forces as the seals are subjected to different aircraft maneuvers: Windmilling restart; start-ground idle; ground idle-takeoff; takeoff-burst chop, etc. Results of this model show that the seal closes appropriately and does not ram into the rotor for all of the conditions analyzed. The rig upgrade design for testing Aspirating Seals has been completed. Long lead-time items (forgings, etc.) have been ordered.

Bharat Bagepalli↗
NASA NTRS2020

AERoBOND Project Summary

Under NASA’s Convergent Aeronautics Solutions (CAS) project, the Adhesive-Free Bonding of Complex Composites (AERoBOND) project investigated off-stoichiometric epoxy polymers for fast, reliable assembly of epoxy matrix composite structures. The project goal was to demonstrate feasibility of the AERoBOND joining method by demonstrating mechanical properties greater than 80% of conventional co-cured materials while reducing structure weight by 1%. The project consisted of three convergent research areas: material and process development, systems analysis, and material and process modeling. Material and process development was the largest component of AERoBOND with approximately 6 FTE and 1WYE of support to formulate and characterize new resins, prepare carbon fiber prepregs, fabricate laminates, measure mechanical properties, analyze failure results, and select material and process improvements. The systems analysis activity estimated the potential reduction in part count and aircraft weight by comparing models of composite wing boxes with no fasteners (co-cured structure), fasteners in major joints (co-cured stringers), and fasteners in all joints. The materials and process modeling activity included a molecular model of the AERoBOND materials system to predict mechanical properties of resins with offset stoichiometry and a process model to predict the effect of resin formulation and processing conditions on the extent of mixing and degree of cure in a finished joint. As the number of airline passenger trips doubles in the next 20 years (IATA/Tourism Economics Air Passenger Forecasts, April 2019), the increased demand for new commercial aircraft is now the single greatest technical challenge to the airframe manufacturing industry. To meet efficiency requirements, new aircraft must be fabricated primarily from high performance structural composites, but manufacturing processes are inherently slow with the largest bottleneck attributed to assembly and installation of fasteners (NASA/TM–2019-220428). Manufactures of commercial transport aircraft are compelled to install more than 100,000 redundant fasteners into bonded joints to prevent failures due to unpredictable weak bonds. In structural adhesive bonds, the interface between adherend and adhesive is nearly two-dimensional making it susceptible to minute quantities of contamination, which can cause weak bonds. Currently, bond strength assessment is only possible through destructive testing (i.e., breaking the joint). For these reasons, regulatory organizations such as the Federal Aviation Administration (FAA) often require redundant load paths in secondary-bonded, primary-structures to alleviate concerns with bond performance. The AERoBOND process enables reflow of matrix resin during assembly to eliminate the material discontinuity at the interface, thereby eliminating the dependence of mechanical performance on interfacial adhesion. The AERoBOND joint is equivalent to the interlaminar region obtained during a co-cure process, so joint performance depends on the cohesive properties of the matrix resin. Conventional co-cured structures, although too costly and complex for large-scale manufacturing, are trusted by manufacturers and regulators, and are certified for flight with few or no redundant fasteners.Systems analysis performed on a composite wing model at the scale of a single-aisle commercial transport aircraft indicated that >20,000 redundant fasteners per wing could be eliminated by implementing the AERoBOND joining method. A total weight reduction of 15% was predicted in a wing box by eliminating fasteners and thinning components that must no longer support localized fastener loads and accommodate fastener dimensions. Interlaminar shear fracture toughness measured by the end-notched flexure test was greater than 1 kJ/m2 (nearly 140% of the co-cured benchmark property), which is greatly in excess of the project goals for mechanical properties. Testing was planned to measure interlaminar tensile fracture toughness as well as interlaminar tensile and shear strengths using the same AERoBOND configuration, but was delayed due to closure of LaRC facilities during the COVID-19 pandemic. The AERoBOND process model is partially validated and available for experimental use. It allows the user to input AERoBOND process parameters such as material composition, laminate configuration, and cure cycle to predict the final cure state of the AERoBOND joint. A preliminary, multi-scale material model was developed to predict AERoBOND joint mechanical properties (stiffness and strength) based on the cure state of the joint provided by the process model. The timing for transition of this technology within NASA is excellent as NASA initiates new enduring projects to address composites manufacturing rate challenges. AERoBOND technology is well suited to AAVP/AATT objectives for rapid manufacturing of a composite wing. A minimal effort (1 FTE/$15k procurement/0 WYE) is proposed in FY21 to continue a minor mechanical testing effort and maintain a SAA with ASX composites to develop commercial quality prepreg material. An RFI with the composites industry is suggested to quantify the technology gap between the current TRL and the TRL needed for transition to industry. A moderate effort [3-4 FTE/$150k/1 WYE (~$115k)] is proposed in FY22 for the “high rate composites manufacturing” project currently in planning. The partnership with ASX Composites will be expanded to produce material for sub-element/element-scale “panel-off” activities. Industry partnerships with airframe manufacturers is an expected component to explore damage tolerance and environmental stability. Further development of multi-scale modeling tools (process model, meso-scale model, and molecular model) is planned to enhance and deliver tools for rapid manufacturing infusion.

Frank Louis Palmieri↗
NASA NTRS1977

Finite element analysis of low speed viscous and inviscid aerodynamic flows

A weak interaction solution algorithm was established for aerodynamic flow about an isolated airfoil. Finite element numerical methodology was applied to solution of each of differential equations governing potential flow, and viscous and turbulent boundary layer and wake flow downstream of the sharp trailing edge. The algorithm accounts for computed viscous displacement effects on the potential flow. Closure for turbulence was accomplished using both first and second order models. The COMOC finite element fluid mechanics computer program was modified to solve the identified equation systems for two dimensional flows. A numerical program was completed to determine factors affecting solution accuracy, convergence and stability for the combined potential, boundary layer, and parabolic Navier-Stokes equation systems. Good accuracy and convergence are demonstrated. Each solution is obtained within the identical finite element framework of COMOC.

Baker, A. J.↗
NASA NTRS1980

QSO evolution and the intergalactic medium

A model is developed of a cold, photoionized intergalactic medium in which only discrete, evolving sources (quasars) produce the observed extragalactic hard X-ray background. Constraints on the model are imposed by observed optical depth limits to various redshifts, far UV and visible flux measurements, X-ray background and source counts in the region 2-6 keV, and it makes use of results of the Palomar Bright Quasar Survey and the Braccesi Quasar survey to obtain expressions for the quasar evolution function, absolute luminosity, present number density and switch-on redshift. Results of computer modeling reveal an intergalactic medium with a density with respect to closure less than 0.35, temperature of 16,000 K and fractional ion abundances of approximately 0.000001 for H I and 0.001 for He II, under the assumption of a Hubble constant of 50 km/sec per Mpc. In addition, spectral indices of -0.4 and -1.3 are obtained for the regions of the emission spectrum above and below 4 x 10 to the 18th Hz.

Sherman, R. D.↗
NASA NTRS1988

Analysis of crack propagation in roller bearings using the boundary integral equation method - A mixed-mode loading problem

Crack propagation in a rotating inner raceway of a high-speed roller bearing is analyzed using the boundary integral method. The model consists of an edge plate under plane strain condition upon which varying Hertzian stress fields are superimposed. A multidomain boundary integral equation using quadratic elements was written to determine the stress intensity factors KI and KII at the crack tip for various roller positions. The multidomain formulation allows the two faces of the crack to be modeled in two different subregions, making it possible to analyze crack closure when the roller is positioned on or close to the crack line. KI and KII stress intensity factors along any direction were computed. These calculations permit determination of crack growth direction along which the average KI times the alternating KI is maximum.

Ghosn, L. J.↗
NASA NTRS1993

Parameterization of sub-grid scale convection

The following topics are discussed: an overview of the cumulus parameterization problem; interactions between explicit and implicit processes in mesoscale models; effects of model grid size on the cumulus parameterization problem; parameterizing convective effects on momentum fields in mesoscale models; differences between slantwise and vertical cumulus parameterization; experiments with different closure hypotheses; and coupling cumulus parameterizations to boundary layer, stable cloud, and radiation schemes.

Frank, William↗
NASA NTRSDate not supplied

Coupling Online Conductivity with Offline Ion Chromatography Measurements using the Particle-into-Liquid Sampler

A particle-into-liquid sampler has been combined with a flow through conductivity cell to provide a continuous, non-destructive, online measurement in support of offline ion chromatography analysis. The conductivity measurement provides a rapid assessment of the total ion concentration augmenting the slower batch data from the offline analysis and is developed primarily to assist airborne measurements, where fast time response is essential. A model of the conductivity was developed for measured ions and excellent closure is derived for laboratory-generated aerosols. The PILS-conductivity measurement was extensively tested during the NASA Cloud, Aerosol and Monsoon Processes: Philippines Experiment (CAMP2Ex) across nineteen research flights and the conductivity data was found to augment the temporal capability of the PILS instrument for assessing ionic species, allowing sub-minute variability to be resolved. Through the sampling of a diverse range of ambient aerosol, including biomass burning, fresh and aged urban pollution, the conductivity measurement offered additional useful information to untangle the composition of complex aerosol mixtures, specifically for assessing acidic aerosol conditions.

Ewan Colin Crosbie↗
NASA NTRS2018

Preliminary Evaluation of the Simplified Higher Order Closure (SHOC) in the NASA Goddard Earth Observing System (GEOS)

Here we evaluate alternative planetary boundary layer (PBL) schemes in the NASA Goddard Earth Observing System (GEOS) model. The model currently calculates turbulent diffusivities with a local Richardson number-based scheme for stable layers, and a K-profile scheme with surface- and cloud-driven plumes based on Lock et al. (2000) for unstable layers. In this study, these are replaced with the TKE-based Simplified Higher Order Closure (SHOC). SHOC combines a 1.5-order turbulence closure with a double-Gaussian PDF serving as a subgrid condensation scheme. The new diffusivity calculations are paired with two new shallow mass flux schemes. First, the University of Washington shallow cumulus scheme of Park and Bretherton (2009), a buoyancy-sorting entraining-detraining bulk plume with overshooting entrainment and a TKE/CIN-based closure. Second, an Eddy Diffusivity Mass Flux (EDMF) scheme based on the multiple plume model of Suselj et al. (2013) with stochastic lateral entrainment. Global climatologies of low level temperature, humidity, and wind are evaluated against MERRA-2 and ERA-Interim reanalyses, and significant reductions in bias are seen with both mass flux schemes. The diurnal cycle over the continental United States is evaluated using MERRA-2 and an in situ network of surface Lidar profilers.

Arnold, Nathan↗