Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Turbulence Stabilization”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

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

At least 1,099 records · Page 61

Receptivity of flat-plate boundary layer in a non-uniform free stream (vorticity normal to the plate)

Recent progress in both the linear and nonlinear aspects of stability theory has highlighted the importance of the receptivity problem. One of the most unclear aspects of receptivity study is the receptivity of boundary-layer flow normal to vortical disturbances. Some experimental and theoretical results permit the proposition that quasi-steady outer-flow vortical disturbances may trigger by-pass transition. In present work such interaction is investigated for vorticity normal to a leading edge. The interest in these types of vortical disturbances arise from theoretical work, where it was shown that small sinusoidal variations of upstream velocity along the spanwise direction can produce significant variations in the boundary-layer profile. In the experimental part of this work, such non-uniform flow was created and the laminar-turbulent transition in this flow was investigated. The experiment was carried out in a low-turbulence direct-flow wind tunnel T-361 at the Central Aerohydrodynamic Institute (TsAGI). The non-uniform flow was produced by laminar or turbulent wakes behind a wire placed normal to the plate upstream of the leading edge. The theoretical part of the work is devoted to studying the unstable disturbance evolution in a boundary layer with strongly non-uniform velocity profiles similar to that produced by outer-flow vorticity. Specifically, the Tollmien-Schlichting wave development in the boundary layer flow with spanwise variations of velocity is investigated.

Kogan, M. N.↗

Observations of airflow around a supertall curved building and its impact on temperature and humidity in Houston's urban center

Characterization of realistically shaped skyscrapers embedded in non-uniform neighbourhoods experiencing intricate weather patterns remains inadequately investigated. Aiming to close this gap, the Center for Multiscale Applied Sensing team deployed its mobile observatory in the street canyons around the curved Wells Fargo Plaza skyscraper in downtown Houston, TX. Three deployments allowed airflow observations under different inflow wind and thermodynamic stability conditions. Doppler lidar measurements reveal that when inflow hits the curved wall of the skyscraper, perpendicular canyons experience similar vortex configurations creating two windward and two leeward circulations. Windsond measurements support that buoyancy within the deep street canyons can generate thermal updrafts as strong as 2 m s -1 which is sufficient to overturn the mechanical downwash under gentle wind conditions. Canyons experiencing venting during the daytime were observed to be more thermodynamically stable at night while thermodynamically stable canyons during the day were observed to be more thermodynamically unstable at night owing to the accumulation of heat near street-level. Fourier decomposition of the vertical velocity measurements shows that in all cases flow exhibited high Reynolds numbers and was composed of turbulent eddies of predominantly 6 and 15 min periods. Here, this observational dataset provides insights to assess wind load, pedestrian comfort, urban air mobility, and natural ventilation and may be used as a benchmark for numerical model and wind tunnel studies attempting to represent realistically complex urban conditions.

54 ENVIRONMENTAL SCIENCES↗

A Statistical Evaluation of WRF-LES Trace Gas Dispersion Using Project Prairie Grass Measurements

In recent years, new measurement systems have been deployed to monitor and quantify methane emissions from the natural gas sector. Large-eddy simulation (LES) has complemented measurement campaigns by serving as a controlled environment in which to study plume dynamics and sampling strategies. However, with few comparisons with controlled-release experiments, the accuracy of LES for modeling natural gas emissions is poorly characterized. In this paper, we evaluate LES from the Weather Research and Forecasting (WRF) Model against Project Prairie Grass campaign measurements and surface layer similarity theory. Using WRF-LES, we simulate continuous emissions from 30 near-surface trace gas sources in two stability regimes: strong convection and weak convection. We examine the impact of grid resolutions ranging from 6.25 to 52 m in the horizontal dimension on model results. We evaluate performance in a statistical framework, calculating fractional bias and conducting Welch’s t tests. WRF-LES accurately simulates observed surface concentrations at 100 m and beyond under strong convection; simulated concentrations pass t tests in this region irrespective of grid resolution. However, in weakly convective conditions with strong winds, WRF-LES substantially overpredicts concentrations—the magnitude of fractional bias often exceeds 30%, and all but one t test fails. The good performance of WRF-LES under strong convection correlates with agreement with local free convection theory and a minimal amount of parameterized turbulent kinetic energy. The poor performance under weak convection corresponds to misalignment with Monin–Obukhov similarity theory and a significant amount of parameterized turbulent kinetic energy.

17 WIND ENERGY↗

Convergence Acceleration of Runge-Kutta Schemes for Solving the Navier-Stokes Equations

The convergence of a Runge-Kutta (RK) scheme with multigrid is accelerated by preconditioning with a fully implicit operator. With the extended stability of the Runge-Kutta scheme, CFL numbers as high as 1000 can be used. The implicit preconditioner addresses the stiffness in the discrete equations associated with stretched meshes. This RK/implicit scheme is used as a smoother for multigrid. Fourier analysis is applied to determine damping properties. Numerical dissipation operators based on the Roe scheme, a matrix dissipation, and the CUSP scheme are considered in evaluating the RK/implicit scheme. In addition, the effect of the number of RK stages is examined. Both the numerical and computational efficiency of the scheme with the different dissipation operators are discussed. The RK/implicit scheme is used to solve the two-dimensional (2-D) and three-dimensional (3-D) compressible, Reynolds-averaged Navier-Stokes equations. Turbulent flows over an airfoil and wing at subsonic and transonic conditions are computed. The effects of the cell aspect ratio on convergence are investigated for Reynolds numbers between 5:7 x 10(exp 6) and 100 x 10(exp 6). It is demonstrated that the implicit preconditioner can reduce the computational time of a well-tuned standard RK scheme by a factor between four and ten.

Swanson, Roy C., Jr.↗

A Flight Investigation of the Effects of Varied Lateral Damping on the Effectiveness of a Fighter Airplane as a Gun Platform

Flight tests were conducted at three conditions of damping and two atmospheric turbulence levels with fixed-reticle, gyro computing, and telescopic gunsights. Results of simulated strafing runs indicate that the gun-line dispersion could be expected to be decreased about 7 percent by increased lateral damping and to be increased about 85 percent by decreased damping and about 40 percent by rough air of the type encountered. Use of the telescopic sight indicated a 20-percent decrease in gun-line dispersion.

OPERATING PROBLEMS, PHYSIOLOGICAL↗

A Nonlocal Magneto-curvature Instability in a Differentially Rotating Disk

A global mode is shown to be unstable to nonaxisymmetric perturbations in a differentially rotating Keplerian disk containing either vertical or azimuthal magnetic fields. In an unstratified cylindrical disk model, using both global eigenvalue stability analysis and linear global initial-value simulations, it is demonstrated that this instability dominates at strong magnetic fields where local standard magnetorotational instability (MRI) becomes stable. Unlike the standard MRI mode, which is concentrated in the high flow shear region, these distinct global modes (with low azimuthal mode numbers) are extended in the global domain and are Alfvén-continuum-driven unstable modes. As its mode structure and relative dominance over MRI are inherently determined by the global spatial curvature as well as the flow shear in the presence of a magnetic field, we call it the magneto-curvature (magneto-spatial-curvature) instability. Consistent with the linear analysis, as the field strength is increased in the nonlinear simulations, a transition from MRI-driven turbulence to a state dominated by global nonaxisymmetric modes is obtained. This global instability could therefore be a source of nonlinear transport in accretion disks at a higher magnetic field than predicted by local models.

79 ASTRONOMY AND ASTROPHYSICS↗

Modeling of Wake-vortex Aircraft Encounters

There are more people passing through the world's airports today than at any other time in history. With this increase in civil transport, airports are becoming capacity limited. In order to increase capacity and thus meet the demands of the flying public, the number of runways and number of flights per runway must be increased. In response to the demand, the National Aeronautics and Space Administration (NASA), in conjunction with the Federal Aviation Administration (FAA), airport operators, and the airline industry are taking steps to increase airport capacity without jeopardizing safety. Increasing the production per runway increases the likelihood that an aircraft will encounter the trailing wake-vortex of another aircraft. The hazard of a wake-vortex encounter is that heavy load aircraft can produce high intensity wake turbulence, through the development of its wing-tip vortices. A smaller aircraft following in the wake of the heavy load aircraft will experience redistribution of its aerodynamic load. This creates a safety hazard for the smaller aircraft. Understanding this load redistribution is of great importance, particularly during landing and take-off. In this research wake-vortex effects on an encountering 10% scale model of the B737-100 aircraft are modeled using both strip theory and vortex-lattice modeling methods. The models are then compared to wind tunnel data that was taken in the 30ft x 60ft wind tunnel at NASA Langley Research Center (LaRC). Comparisons are made to determine if the models will have acceptable accuracy when parts of the geometry are removed, such as the horizontal stabilizer and the vertical tail. A sensitivity analysis was also performed to observe how accurately the models could match the experimental data if there was a 10% error in the circulation strength. It was determined that both models show accurate results when the wing, horizontal stabilizer, and vertical tail were a part of the geometry. When the horizontal stabilizer and vertical tail were removed there were difficulties modeling the sideforce coefficient and pitching moment. With the removal of only the vertical tail unacceptable errors occurred when modeling the sideforce coefficient and yawing moment. Lift could not be modeled with either the full geometry or the reduced geometry attempts.

Smith, Sonya T.↗

Dynamics and Structure of Weakly-Strained Flames In Normal- and Micro-Gravity

Strained laminar flames have been systematically studied, as the understanding of their structure and dynamic behavior is of relevance to turbulent combustion. Most of these studies have been conducted in opposed-jet, stagnation-type flow configurations. Flame studies in stagnation flows also allow for the determination of fundamental flame properties such as laminar flame speeds and extinction states that can be also conveniently modeled in detail. Studies at high strain rates are important in quantifying and understanding the response of vigorously-burning flames under conditions in which the transport time scales become comparable to the chemical time scales. Studies of weakly-strained flames can be of particular interest for all stoichiometries. For example, the laminar flame speeds for any equivalence ratio, phi, can be accurately determined by using the counterflow technique only if measurements are obtained at very low strain rates. Furthermore, near-limit flames can be only stabilized by weak strain rates. Previous studies have shown that weakly-burning flames are particularly sensitive to chain mechanisms, thermal radiation, and unsteadiness. The stabilization and study of weakly-strained flames is complicated by the presence of buoyancy that can render the flames unstable to the point of extinction. Such instabilities are caused either by the induced natural convection or the fact that higher density fluid can find itself on top of a lower density fluid. Thus, the use of microgravity (mu-g) becomes essential in order to provide meaningful insight into this important combustion regime. In view of the foregoing considerations, the main objectives of the program are to: (1) Experimentally determine the laminar flame speed at near-zero strain rates; (2) Experimentally determine the extinction limits of near-limit flames; (3) Experimentally determine the response of near-limit flames to unsteadiness and heat loss; (4) Introduce Digital Particle Image Velocimetry (DPIV) in microgravity; (5) Conduct detailed numerical simulations of the experiments; and (6) Provide physical insight into the underlying physico-chemical mechanisms.

Zhang, Hai↗

A flexible gyro-fluid system of equations

Gyro-fluid equations are velocity space moments of the gyrokinetic equations. Special gyro-Landau-fluid closures have been developed that include the damping due to kinetic resonances by fitting to the collisionless local plasma response functions. This damping allows for accurate linear eigenmodes to be computed with a relatively low number of velocity space moments compared to the number of velocity quadrature points in gyrokinetic codes. However, none of the published gyro-Landau-fluid closure schemes considers the Onsager symmetries of the resulting quasi-linear fluxes as a constraint. Onsager symmetry guarantees that the matrix of diffusivities is positive definite, an important property for the numerical stability of a transport solver. A two-parameter real closure for improving the accuracy of low-resolution gyro-fluid equations, which preserves the Onsager symmetry and allows higher velocity space moments, is presented in this paper. The new linear gyro-fluid system (GFS) is used to extend the TGLF quasi-linear transport model so that it can compute the energy and momentum fluxes due to parallel magnetic fluctuations, completing the transport matrix. The GFS equations do not use a bounce average approximation. The GFS equations are fully electromagnetic with general flux surface magnetic geometry, pitch angle scattering for electron collisions, and subsonic equilibrium toroidal rotation. Using GFS eigenmodes in the quasi-linear TGLF model will be shown to yield a more accurate match to fluxes computed by CGYRO turbulence simulations. In conclusion, prospects for future applications of a quasi-linear theory to new plasma transport regimes and magnetic confinement devices in addition to tokamaks are opened by the flexibility of the GFS eigensolver.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Simulations of Boundary-Layer Transition

For incompressible benchmark flows, we have demonstrated the capability of the parabolized stability equations (PSE) to simulate the transition process in excellent agreement with microscopic experiments and direct Navier-Stokes simulations at modest computational cost. Encouraged by these results, we have developed the PSE methodology of three-dimensional boundary-layers in general curvilinear coordinates for the range from low to hypersonic speeds, and for both linear and nonlinear problems. For given initial and boundary conditions, the approach permits simulations from receptivity through linear and secondary instabilities into the late stages of transition where significant changes in skin friction and heat transfer coefficients occur. We have performed transition simulations for a variety of two- and three-dimensional similarity solutions and for realistic flows over swept wings at subsonic and supersonic speeds, the pressure ans suction side of turbine blades at low and medium turbulence levels, and over a blunt cone at Mach number Ma = 8. We present selected results for different transition mechanisms with emphasis on the late stage of transition and the evolution of wall-shear stress and heat transfer.

Herbert, Thorwald↗

Mid-Fidelity Computational Fluid Dynamics Analysis of the Elytron 4S UAV Concept

The Elytron 4S Unmanned Aerial Vehicle (UAV) Concept was developed to combine the advantages of fixed- and rotary-wing technology. The 4S Concept is a box-wing configuration with rotors mounted on a centrally located tiltwing. The UAV is intended to be capable of both conventional takeoff and landing (CTOL) and vertical takeoff and landing (VTOL), and is envisioned to excel in UAV performance because of the combined efficiency of fixed-wing aircraft and the hover and VTOL capabilities of regular drones or quadcopters. A mid-fidelity Unsteady Reynolds-Averaged Navier-Stokes (URANS) approach using Rotorcraft CFD (RotCFD) is performed to analyze and characterize the performance of the aircraft. The flow field is coupled with a rotor model based on blade-element momentum theory to model the 4S UAV rotors. Turbulence is modeled using a realizable k- turbulence model with special wall function. The code is used to generate aerodynamic forces and moments on the body at cruise conditions, and during VTOL. The results and their uncertainties are characterized, and an angle- of-attack and sideslip sweep are computed, both with and without rotors on. Simulations are compared with the wind tunnel tests in the 7- by 10-Foot U.S. Army Wind Tunnel at NASA Ames Research Center, performed in 2017. Results show promising comparison with experimental data, despite a late change in rotor size and rudder size of the physical model that cause the expected deviations from the simulation. A slight change in the net thrust value, when rotors are modeled, is observed because of the rotor diameter increase on the physical model. A noticeable difference in the directional stability was observed because of the increased rudder surface and added strakes. These changes were implemented to improve on the design as simulated, which is observed in the results. The simulation results paved the way to the first successful flight of the UAV Concept.

Koning, Witold J. K.↗

Design Exploration of a Transonic Cruise Slotted Airfoil

A knowledge-based aerodynamic design method, CDISC, has been leveraged for the computational design exploration of the cruise slotted airfoil as a drag-saving technology for single aisle, transonic transport aircraft. Aerodynamic predictions were generated using the NASA USM3D Reynolds-averaged Navier-Stokes flow solver, and laminar flow assessments were conducted using the NASA BLSTA3D boundary layer profile solver paired with the LASTRAC stability analysis and transition prediction software. The CDISC design method was used to parametrically vary several slotted-airfoil design variables to understand their impact on aerodynamic performance. In addition to establishing best practices for slotted airfoil design, it was observed that the historically noted skin-friction penalty for cruise slotted airfoils was attributed to the low-Reynolds-number boundary layer of the flap. At cruise, this skin-friction penalty negated any potential decrease in pressure drag enabled by the cruise slotted airfoil architecture. This observation led to the development of an Aft-Laminar Multi-element Airfoil concept that uses airfoil shaping to promote natural laminar flow on the flap. At cruise conditions, the concept is predicted to offer a 2.7% reduction in sectional drag relative to a fully turbulent supercritical airfoil. Off-design analyses showed that laminar flow could be maintained with near-cruise variations in angle of attack and Mach number, resulting in sustained aerodynamic efficiency improvements and a delay in the drag rise Mach number. A low-speed, high-lift analysis at takeoff conditions predicted a maximum lift coefficient of 2.4 could be achieved by use of a variable-camber leading edge on the main element and simple deflection of the flap. These computationally-predicted benefits have motivated future research toward an aft-laminar cruise slotted wing design as an incremental step toward the application of natural laminar flow technology on transonic commercial transports.

Brett R Hiller↗

WEST L-mode record long pulses guided by predictions using Integrated Modeling

A new record was set on the WEST Tokamak, designed to operate long duration plasmas in a tungsten (W) environment, with an injected energy of 1.15 GJ and a plasma duration 364s. Scenario development was supported by integrated modeling using the High Fidelity Plasma Simulator (HFPS), the European IMAS-coupled version of JETTO/JINTRAC, which integrates physics-driven modules into a unified framework. In particular, a reduced model for Lower-Hybrid heating and Current-Drive (LHCD) and the quasi-linear turbulent transport model TGLF are crucial for long pulses predictions up to the Last Closed Flux Surface (LCFS). Using this workflow, a 100 s reference discharge was modeled and plasma kinetic profiles and loop voltage were quantitatively well matched. In preparation for the recent long duration experiments, non-inductive current-drive actuators (I P , n e , P LHCD ) were varied to determine the operational domain going towards fully non-inductive discharges. In particular, decreasing the plasma current is shown to ease the access to such conditions, with a careful monitoring of (n e , P LHCD ) to avoid machine limitations. In addition, post-prediction experiments conducted within the investigated parameter range validated the predicted dependencies and were shown to be in quantitative agreement. Exploratory work on the use of ECCD for MHD stability purpose is also introduced.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Performance of direct injected propane and gasoline in a high stroke-to-bore ratio SI engine: Pathways to diesel efficiency parity with ultra low soot

This work explores pathways to achieve diesel-like, high-efficiency combustion with stoichiometric 3-way catalyst compatible combustion in a single-cylinder spark ignition (SI) research engine. A unique high stroke-to-bore engine design (1.5:1) with cooled exhaust gas recirculation (EGR) and high compression ratio ( r c ) was used to improve engine efficiency by up to 30% compared with a production turbocharged gasoline direct injection spark ignition engine. Engine experiments were conducted with both 91 RON E10 gasoline and liquified petroleum gas (LPG) (i.e. autogas) and were compared to legacy gasoline data on the production engine. Geometric compression ratio ( r c ) of 13.3:1 was used for both fuels with additional experiments at 16.8:1 for LPG only. Measurements of exhaust soot particle size and number concentrations were made with both fuels. Significant reduction in soot particles across the whole particle size range were achieved with LPG due to the elimination of in-cylinder liquid films. The effects of EGR, late intake valve closing (IVC) and fuel characteristics were investigated through their effects on efficiency, combustion stability and soot production. Results of 47% gross thermal efficiency, and 45% net thermal efficiency at stoichiometric engine operation, at up to 17 bar IMEP and 2000 r/min with 16.8:1 r c were achieved with LPG. Estimated brake efficiency values were compared to a contemporary medium duty diesel engine illustrating the benefits of the chosen path for achieving diesel efficiency parity.

33 ADVANCED PROPULSION SYSTEMS↗

Stability analysis of WEST L-mode discharges with improved confinement from boron powder injection

WEST L-mode plasmas with dominant electron heating and no core torque source have observed improvements in confinement during boron (B) powder injection. These results are reminiscent of previous powder injection experiments on other devices and gaseous impurity seeding experiments on WEST. During powder injection, the stored energy increased up to 25% due to enhanced ion and electron heat and particle confinement. The improvements in confinement were not indicative of an L-H transition. To identify the dominant mechanisms and the causality chain behind these improvements in confinement, we employ interpretative modeling using METIS, predictive integrated modeling using a high-fidelity plasma simulator (HFPS), and stand-alone gyrokinetic simulations using quasi-linear gyrokinetic code. Interpretative modeling with METIS allowed for the estimation of missing data while maintaining good overall consistency with experiment. These results provided the initial conditions for fully predictive flux driven simulations using the HFPS. From these simulations, quasi-linear gyrokinetic analysis was performed at ρ = 0.5 and ρ = 0.65. Collisionality was found to be a strong candidate for the turbulence suppression mechanism at ρ = 0.5, while a combination of collisionality and the T e /T i ratio was found to be the likely mechanism at ρ = 0.65. The results additionally suggested that increased Z eff (through main ion dilution) could play a role in the improved confinement, but this could not be confirmed due to a lack of experimental measurements. The modeling framework established here can now be used to evaluate and exploit a variety of future powder injection experiments.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Ship‐Based Lidar Evaluation of Southern Ocean Low Clouds in the Storm‐Resolving General Circulation Model ICON and the ERA5 and MERRA‐2 Reanalyses

Global storm resolving models (GSRMs) represent the next generation of global climate models. One of them is a 5-km Icosahedral Nonhydrostatic Weather and Climate Model (ICON). Its high resolution means that parameterizations of convection and clouds, including subgrid-scale clouds, are omitted, relying on explicit simulation but necessarily utilizing microphysics and turbulence parameterizations. Standard-resolution (10–100 km) models, which use convection and cloud parameterizations, have substantial cloud biases over the Southern Ocean (SO), adversely affecting radiation and sea surface temperature. The SO is dominated by low clouds, which cannot be observed accurately from space due to overlapping clouds, attenuation, and ground clutter. We evaluated SO clouds in ICON and the ERA5 and MERRA-2 reanalyzes using approximately 2400 days of lidar observations and 2300 radiosonde profiles from 31 voyages and a Macquarie Island station during 2010–2021, compared to the model and reanalyzes using a ground-based lidar simulator. We found that ICON and the reanalyzes underestimate the total cloud fraction by about 10% and 20%, respectively. ICON and ERA5 overestimate the cloud occurrence peak at about 500 m, associated with underestimated lower tropospheric stability and overestimated lifting condensation level. The reanalyzes strongly underestimate fog and very low-level clouds, and MERRA-2 underestimates cloud occurrence at almost all heights. Outgoing shortwave radiation is overestimated in MERRA-2, implying a “too few, too bright” cloud problem. SO cloud and fog biases are a substantial issue in the analyzed model and reanalyzes and result in shortwave and longwave radiation biases.

Kuma, Peter [Stockholm Univ. (Sweden); Univ. of Ca↗

Piloted Evaluation of Modernized Limited Authority Control Laws in the NASA-Ames Vertical Motion Simulator (VMS)

The UH-60 BLACK HAWK was designed in the 1970s, when the US Army primarily operated during the day in good visual conditions. Subsequently, the introduction of night-vision goggles increased the BLACK HAWK'S mission effectiveness, but the accident rate also increased. The increased accident rate is strongly tied to increased pilot workload as a result of a degradation in visual cues. Over twenty years of research in helicopter flight control and handling qualities has shown that these degraded handling qualities can be recovered by modifying the response type of the helicopter in low speed flight. Sikorsky Aircraft Corporation initiated a project under the National Rotorcraft Technology Center (NRTC) to develop modern flight control laws while utilizing the existing partial authority Stability Augmentation System (SAS) of the BLACK HAWK. This effort resulted in a set of Modernized Control Laws (MCLAWS) that incorporate rate command and attitude command response types. Sikorsky and the US Army Aeroflightdynamics Directorate (AFDD) conducted a piloted simulation on the NASA-Ames Vertical h4otion Simulator, to assess potential handling qualities and to reduce the risk of subsequent implementation and flight test of these modern control laws on AFDD's EH-60L helicopter. The simulation showed that Attitude Command Attitude Hold control laws in pitch and roll improve handling qualities in the low speed flight regime. These improvements are consistent across a range of mission task elements and for both good and degraded visual environments. The MCLAWS perform better than the baseline UH-60A control laws in the presence of wind and turbulence. Finally, while the improved handling qualities in the pitch and roll axis allow the pilot to pay more attention to the vertical axis and hence altitude performance also improves, it is clear from pilot comments and altitude excursions that the addition of an Altitude Hold function would further reduce workload and improve overall handling qualities of the aircraft.

Sahasrabudhe, Vineet↗

Memory-based parameterization with differentiable solver: Application to Lorenz ’96

Physical parameterizations (or closures) are used as representations of unresolved subgrid processes within weather and global climate models or coarse-scale turbulent models, whose resolutions are too coarse to resolve small-scale processes. These parameterizations are typically grounded on physically based, yet empirical, representations of the underlying small-scale processes. Machine learning-based parameterizations have recently been proposed as an alternative solution and have shown great promise to reduce uncertainties associated with the parameterization of small-scale processes. Yet, those approaches still show some important mismatches that are often attributed to the stochasticity of the considered process. This stochasticity can be due to coarse temporal resolution, unresolved variables, or simply to the inherent chaotic nature of the process. To address these issues, we propose a new type of parameterization (closure), which is built using memory-based neural networks, to account for the non-instantaneous response of the closure and to enhance its stability and prediction accuracy. We apply the proposed memory-based parameterization, with differentiable solver, to the Lorenz ’96 model in the presence of a coarse temporal resolution and show its capacity to predict skillful forecasts over a long time horizon of the resolved variables compared to instantaneous parameterizations. This approach paves the way for the use of memory-based parameterizations for closure problems.

Mathematics↗