Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Boundary condition”

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,153 records · Page 64
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Black Holes and de Sitter Space as Time Mirrors

It is usually assumed that matter disappears together with the spacetime at the center of a Schwarzschild black hole (BH). Here, we find that if we impose a boundary condition that the field does not disappear at the BH center (that is, field flux into the singularity vanishes), the BH acts as a time mirror that totally reflects the infalling light and matter outside the BH. Namely, the reflected field propagates backward in time, passes the event horizon and moves away from the BH. In this case, a BH can be used as a time machine that allows us to send a signal into the past. We also show that de Sitter spacetime acts as a time mirror provided particles do not disappear from the spacetime at r=∞.

Svidzinsky, Anatoly (ORCID:0000000247756767)↗
DOE OSTI2020

An Assessment of the Influence of Uncertainty in Temporally Evolving Streamflow Forecasts on Riverine Inundation Modeling

Continental-scale river forecasting platforms forecast streamflow at reaches that can be used as boundary conditions to drive a local-scale flood inundation model. Uncertainty accumulated during this process stems not only from any part of the forecasting chain but can also be caused by the daily variations in weather forcing that keeps evolving as the event advances. This work aims to examine the influence of the evolving forecast streamflow on predicting the maximum inundation for extreme floods. A diagnostic case study was made on the basis of a hindcast of Hurricane Matthew striking the eastern U.S. in 2016. The U.S. National Water Model was one-way coupled to a hydrodynamic inundation model through a developed automated workflow. Although the river forcing has significantly mismatched hydrographs versus observations, the simulated peak water surface elevations and maximum extents were validated to be comparable with the observations, which indicates that the inundation model may not be sensitive to the inherited uncertainty from the weather forcing. Moreover, the uncertainty of the forecast streamflow time series caused only one order of magnitude fewer variations in inundation prediction; this dampening effect may become clearer for extreme events with large areas inundated. In addition, the forecast total volume of stream discharge appears to be an important metric for assessing the performance of river forcing for inundation mapping, as a linear correlation between the total volume and the accuracy of the predicted peak water surface elevation and maximum extent was found, with the coefficients of determination all above 0.8. Extra best-practice experience of running similar operational tasks demonstrated the tradeoff between the cost and accuracy gain.

flood, national water model, RIFT↗
DOE OSTI2020

Full Transport General Relativistic Radiation Magnetohydrodynamics for Nucleosynthesis in Collapsars

We model a compact black hole-accretion disk system in the collapsar scenario with full transport, frequency dependent, general relativistic radiation magnetohydrodynamics. We examine whether or not winds from a collapsar disk can undergo rapid neutron capture (r-process) nucleosynthesis and significantly contribute to solar r-process abundances. We find the inclusion of accurate transport has significant effects on outflows, raising the electron fraction above ${Y}_{{\rm{e}}}\sim 0.3$ and preventing third-peak r-process material from being synthesized. We analyze the time evolution of neutrino processes and electron fraction in the disk and present a simple one-dimensional model for the vertical structure that emerges. Finally, we compare our simulation to semi-analytic expectations and argue that accurate neutrino transport and realistic initial and boundary conditions are required to capture the dynamics and nucleosynthetic outcome of a collapsar.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Parametric Decay of Alfvénic Wave Packets in Nonperiodic Low-beta Plasmas

The parametric decay of finite-size Alfvén waves in nonperiodic low-beta plasmas is investigated using one-dimensional (1D) hybrid simulations. Compared with the usual small periodic system, a wave packet in a large system under the absorption boundary condition shows different decay dynamics, including reduced energy transfer, localized density cavitation, and ion heating. The resulting Alfvén wave dynamics are influenced by several factors relating to this instability, including the growth rate, central wave frequency, and unstable bandwidth. A final steady state of the wave packet may be achieved when the instability does not have enough time to develop within the residual packet, and the packet size shows well-defined scaling dependencies on the growth rate, wave amplitude, and plasma beta. Under the proper conditions, enhanced secondary decay can also be excited in the form of a narrow, amplified wave packet. These results may help to interpret laboratory and spacecraft observations of Alfvén waves, and to refine our understanding of the associated energy transport and ion heating.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

SAGAbg. I. A Near-unity Mass-loading Factor in Low-mass Galaxies via Their Low-redshift Evolution in Stellar Mass, Oxygen Abundance, and Star Formation Rate

Measuring the relation between star formation and galactic winds is observationally difficult. In this work we make an indirect measurement of the mass-loading factor (the ratio between the mass outflow rate and star formation rate) in low-mass galaxies using a differential approach to modeling the low-redshift evolution of the star-forming main sequence and mass–metallicity relation. We use Satellites Around Galactic Analogs (SAGA) background galaxies, i.e., spectra observed by the SAGA Survey that are not associated with the main SAGA host galaxies, to construct a sample of 11,925 spectroscopically confirmed low-mass galaxies from 0.01 ≲ z ≤ 0.21 and measure auroral line metallicities for 120 galaxies. The crux of the method is to use the lowest-redshift galaxies as the boundary condition of our model, and to infer a mass-loading factor for the sample by comparing the expected evolution of the low-redshift reference sample in stellar mass, gas-phase metallicity, and star formation rate against the observed properties of the sample at higher redshift. We infer a mass-loading factor of η m = $0.92^{+1.76}_{-0.74}$, which is in line with direct measurements of the mass-loading factor from the literature despite the drastically different sets of assumptions needed for each approach. While our estimate of the mass-loading factor is in good agreement with recent galaxy simulations that focus on resolving the dynamics of the interstellar medium, it is smaller by over an order of magnitude than the mass-loading factor produced by many contemporary cosmological simulations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2025

Transition to Petschek Reconnection in Subrelativistic Pair Plasmas: Implications for Particle Acceleration

While relativistic magnetic reconnection in pair plasmas has emerged in recent years as a candidate for the origin of radiation from extreme astrophysical environments, the corresponding subrelativistic pair-plasma regime has remained less explored, leaving open the question of how relativistic physics affects reconnection. In this paper, we investigate the differences between these regimes by contrasting two-dimensional particle-in-cell simulations of reconnection in pair plasmas with relativistic magnetization (σ ≫ 1) and subrelativistic magnetization (σ < 1). By utilizing unprecedentedly large domain sizes and outflow boundary conditions, we demonstrate that lowering the magnetization results in a change in the reconnection geometry from a plasmoid chain to a Petschek geometry, where laminar exhausts bounded by slow-mode shocks emanate from a single diffusion region. We attribute this change to the reduced plasmoid production rate in the low-σ case: When the secondary tearing rate is sufficiently low, plasmoids are too few in number to prevent the system from relaxing into a stable Petschek configuration. This geometric change also affects particle energization: We show that while high-σ plasmoid chains generate power-law energy spectra, low-σ Petschek exhausts merely heat incoming plasma and yield negligible nonthermal acceleration. These results have implications for predicting the global current sheet geometry and the resulting energy spectra in a variety of systems.

Plasma astrophysics↗
DOE OSTI2021

Outward Migration of Super-Jupiters

We report recent simulations show that giant planets of about 1 M J migrate inward at a rate that differs from the type II prediction. Here we show that at higher masses, planets migrate outward. Our result differs from previous ones because of our longer simulation times, lower viscosity, and boundary conditions that allow the disk to reach a viscous steady state. We show that, for planets on circular orbits, the transition from inward to outward migration coincides with the known transition from circular to eccentric disks that occurs for planets more massive than a few Jupiters. In an eccentric disk, the torque on the outer disk weakens due to two effects: the planet launches weaker waves, and those waves travel further before damping. As a result, the torque on the inner disk dominates, and the planet pushes itself outward. Our results suggest that the many super-Jupiters observed by direct imaging at large distances from the star may have gotten there by outward migration.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

The Neutronic Engine: A Platform for Operando Neutron Diffraction in Internal Combustion Engines

Neutron diffraction is a powerful tool for noninvasive and nondestructive characterization of materials and can be applied even in large devices such as internal combustion engines thanks to neutrons’ exceptional ability to penetrate many materials. While proof-of-concept experiments have shown the ability to measure spatially and temporally resolved lattice strains in a small aluminum engine on a timescale of minutes over a limited spatial region, extending this capability to timescales on the order of a crank angle degree over the full volume of the combustion chamber requires careful design and optimization of the engine structure to minimize attenuation of the incident and diffracted neutrons to maximize count rates. Here, we present the design of a “neutronic engine,” which is analogous to an optical engine in that the materials and external geometry of a typical automotive engine have been optimized to maximize access of the diagnostic while maintaining the internal combustion chamber geometry and operability of the engine. The high transparency of aluminum to neutrons makes it the ideal window material for neutron diagnostics, which allows the neutronic engine to be a truly all-metal engine with the same load and boundary condition capabilities of a modern downsized passenger car engine. The neutronic engine will enable 3D and time-resolved measurements of strain, stress, and temperature fields as well as phase transformation, texture, and microstructure throughout the metal components of the combustion chamber.

Wissink, Martin↗
DOE OSTI2021

Three-Dimensional CFD Investigation of Pre-Spark Heat Release in a Boosted SI Engine

Low-temperature heat release (LTHR) in spark-ignited internal combustion engines is a critical step toward the occurrence of auto-ignition, which can lead to an undesirable phenomenon known as engine knock. As such, correct predictions of LTHR are of utmost importance to improve the understanding of knock and enable techniques aimed at controlling it. While LTHR is typically obscured by the deflagration following the spark ignition, extremely late ignition timings can lead to LTHR occurrence prior to the spark, i.e., pre-spark heat release (PSHR). In this research, PSHR in a boosted direct-injection SI engine was numerically investigated using three-dimensional computational fluid dynamics (CFD). A hybrid approach was used, based on the G-equation model for representing the turbulent flame front and the multi-zone well-stirred reactor model for tracking the chemical reactions within the unburnt region. A recently developed best practice was also employed which keeps the well-stirred reactor model active throughout the entire simulation. This allowed for correct predictions of the previous cycle trapped residuals which have a considerable effect on the onset of LTHR. Multi-cycle simulations were conducted using Co-Optima alkylate and E30 fuels. The predicted in-cylinder pressure and heat release rate agreed well with the experimental data and served as validation for the CFD model. Following the initial validation, the dynamics of PSHR was discussed and a series of parameters of interest were assessed. First, the effect of exhaust valve temperature was investigated, qualitatively highlighting the importance of boundary conditions uncertainty. Further analyses were carried out on the effects of fuel properties, including laminar flame speed (LFS) and heat of vaporization (HOV). The results indicate that PSHR phasing is slightly advanced with lower LFS as more trapped unburnt fuel is made available for the next cycle. A similar trend was observed with lower HOV as less intense spray cooling led to higher mixture temperatures, i.e., higher mixture reactivity. Finally, a comparison of Co-Optima alkylate and E30 fuels was made using the pressure-temperature trajectory framework. It was shown that the differences between Co-Optima alkylate and E30 in PSHR tendency are correlated with both HOV and chemical effects.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

Effects of Port Angle on Scavenging of an Opposed Piston Two-Stroke Engine

Opposed-piston 2-stroke (OP-2S) engines have the potential to achieve higher thermal efficiency than a typical diesel engine. However, the uniflow scavenging process is difficult to control over a wide range of speeds and loads. Scavenging performance is highly sensitive to pressure dynamics, port timings, and port design. This study proposes an analysis of the effects of port vane angle on the scavenging performance of an opposed-piston 2-stroke engine via simulation. A CFD model of a three-cylinder opposed-piston 2-stroke was developed and validated against experimental data collected by Achates Power Inc. One of the three cylinders was then isolated in a new model and simulated using cycle-averaged and cylinder-averaged initial/ boundary conditions. This isolated cylinder model was used to efficiently sweep port angles from 12 degrees to 29 degrees at different pressure ratios. Results indicate that scavenging performance is correlated with the bulk swirl ratio generated by these port angles. Scavenging performance is also sensitive to the pressure ratio across the engine. It was concluded that, for a given pressure ratio, the smallest port angle produces the best scavenging; however, previous work indicates lower bulk swirl ratio during combustion can have an impact on emissions formation.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2022

An Investigation into the Effects of Swirl on the Performance and Emissions of an Opposed-Piston Two-Stroke Engine using Large Eddy Simulations

Opposed-piston two-stroke (OP-2S) engines have the potential to achieve higher thermal efficiency than a conventional four-stroke diesel engine. However, the uniflow scavenging process is difficult to control over a wider range of speed and loads due to its sensitivity to pressure dynamics, port timings, and port design. Specifically, the angle of the intake ports can be used to generate swirl which has implications for open and closed cycle effects. This study proposes an analysis of the effects of port angle on the in-cylinder flow distribution and combustion performance of an OP-2S using computational fluid dynamics engine. Large Eddy Simulation (LES) was used to model turbulence given its ability to predict in-cylinder mixing and cyclic variability. A three-cylinder model was validated to experimental data collected by Achates Power and the grid was verified using an LES quality approach from the literature. The model was used to simulate port angles from 12 to 29 degrees at constant pressure and temperature boundary conditions. Results indicated that the higher bulk swirl ratio generated by larger port angles tends to trap more internal residuals. This effect on the scavenging performance, combined with the larger trapped swirl ratio, also has a significant impact on the combustion performance in a two-stroke engine. It was concluded that there exists a tradeoff of efficiency and emissions that must be considered when increasing the port angle in a uniflow two-stroke engine. 1. Introduction Two-stroke engines were developed in the 1800s and are commonly used in practice today for lightweight power applications such as motorcycles or handheld power tools, as well as large-bore low speed engines for stationary power and marine applications. Two-stroke engines provide the distinct advantage of higher torque/power output but present challenges due to its coupled intake and exhaust process. Four-stroke engines have the benefit of using two full piston strokes to induct intake charge and to expel exhaust gases. Controlling the trapped charge composition in a twostroke engine via a scavenging process is crucial for reliable and efficient engine operation.

33 ADVANCED PROPULSION SYSTEMS↗
DOE OSTI2023

Data-driven, structure-preserving approximations to entropy-based moment closures for kinetic equations

In this study, we present a data-driven approach for approximating entropy-based closures of moment systems from kinetic equations. The proposed closure learns the entropy function by fitting the map between the moments and the entropy of the moment system, and thus does not depend on the spacetime discretization of the moment system or specific problem configurations such as initial and boundary conditions. With convex and C 2 approximations, this data-driven closure inherits several structural properties from entropy-based closures, such as entropy dissipation, hyperbolicity, and H-Theorem. We construct convex approximations to the Maxwell–Boltzmann entropy using convex splines and neural networks, test them on the plane source benchmark problem for linear transport in slab geometry, and compare the results to the standard, entropy-based systems which solve a convex optimization problem to find the closure. Numerical results indicate that these data-driven closures provide accurate solutions in much less computation time than that required by the optimization routine.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Development of a hardware-in-the-loop testbed for laboratory performance verification of flexible building equipment in typical commercial buildings

The goals of reducing energy costs, shifting electricity peaks, increasing the use of renewable energy, and enhancing the stability of the electric grid can be met in part by fully exploiting the energy flexibility potential of buildings and building equipment. The development of strategies that exploit these flexibilities could be facilitated by publicly available high-resolution datasets illustrating how control of HVAC systems in commercial buildings can be used in different climate zones to shape the energy use profile of a building for grid needs. This article presents the development and integration of a Hardware-In-the-Loop Flexible load Testbed (HILFT) that integrates physical HVAC systems with a simulated building model and simulated occupants with the goal of generating datasets to verify load flexibility of typical commercial buildings. Compared to simulation-only experiments, the hardware-in-the-loop approach captures the dynamics of the physical systems while also allowing efficient testing of various boundary conditions. The HILFT integration in this article is achieved through the co-simulation among various software environments including LabVIEW, MATLAB, and EnergyPlus. Although theoretically viable, such integration has encountered many real-world challenges, such as: 1) how to design the overall data infrastructure to ensure effective, robust, and efficient integration; 2) how to avoid closed-loop hunting between simulated and emulated variables; 3) how to quantify system response times and minimize system delays; and 4) how to assess the overall integration quality. Lessons-learned using the examples of an AHU-VAV system, an air-source heat pump system, and a water-source heat pump system are presented.

Chen, zhelun↗
DOE OSTI2025

Infinite temperature at zero energy

We construct a family of static, geometrically local Hamiltonians that inherit eigenstate properties of periodically-driven (Floquet) systems. Our construction is a variation of the Feynman-Kitaev clock -- a well-known mapping between quantum circuits and local Hamiltonians -- where the clock register is given periodic boundary conditions. Assuming the eigenstate thermalization hypothesis (ETH) holds for the input circuit, our construction yields Hamiltonians whose eigenstates have properties characteristic of infinite temperature, like volume-law entanglement entropy, across the whole spectrum -- including the ground state. We then construct a family of exactly solvable Floquet quantum circuits whose eigenstates are shown to obey the ETH at infinite temperature. Combining the two constructions yields a new family of local Hamiltonians with provably volume-law-entangled ground states, and the first such construction where the volume law holds for all contiguous subsystems.

FOS: Physical sciences↗
DOE OSTI2026

Quantum fate of the Choptuik naked singularity

Classical critical collapse provides a dynamical route from smooth initial data to a naked singularity, representing a sharper violation of predictability than ordinary black hole singularities. We argue that this distinction is erased by quantum backreaction. Building on the semiclassical interior analysis, where quantum self-energy of the collapsing matter generates a universal growing mode and a finite mass gap, we study the exterior naked singularity region that determines global visibility in the Einstein-scalar system. We analyze controlled exterior models in both $2+1$ and $3+1$ dimensions. In the former, smooth matching and physical boundary conditions analytically select a vacuum polarization state, whose backreaction cloaks the classically naked region by a quantum trapped branch. In the latter, numerical horizon tracing shows that near a quantum-shifted threshold the exterior develops finite-mass marginally trapped surfaces rather than a zero-mass naked endpoint. These results suggest a global quantum picture in which the Choptuik naked singularity shares the fate of an ordinary black hole singularity: quantum effects push the putative Cauchy horizon behind a quantum-generated horizon, thereby reducing the loss of predictability to the standard black hole evaporation problem.

FOS: Physical sciences↗
DOE OSTI2024

Aerosol-induced closure of marine cloud cells: enhanced effects in the presence of precipitation

Abstract. The Weather Research Forecasting (WRF) version 4.3 model is configured within a Lagrangian framework to quantify the impact of aerosols on evolving cloud fields. Kilometer-scale simulations utilizing meteorological boundary conditions are based on 10 case study days offering diverse meteorology during the Aerosol and Cloud Experiments in the Eastern North Atlantic (ACE-ENA). Measurements from aircraft, the ground-based Atmosphere Radiation Measurement (ARM) site at Graciosa Island in the Azores, and A-Train and geostationary satellites are utilized for validation, demonstrating good agreement with the WRF-simulated cloud and aerosol properties. Higher aerosol concentration leads to suppressed drizzle and increased cloud water content in all case study days. These changes lead to larger radiative cooling rates at cloud top, enhanced vertical velocity variance, and increased vertical and horizontal wind speed near the base of the lower-tropospheric inversion. As a result, marine cloud cell area expands, narrowing the gap between shallow clouds and increasing cloud optical thickness, liquid water content, and the top-of-atmosphere outgoing shortwave flux. While similar aerosol effects are observed in lightly to non-raining clouds, they tend to be smaller by comparison. These simulations show a relationship between cloud cell area expansion and the radiative adjustments caused by liquid water path and cloud fraction changes. The adjustments are positive and scale as 74 % and 51 %, respectively, relative to the Twomey effect. While higher-resolution large-eddy simulations may provide improved representation of cloud-top mixing processes, these results emphasize the importance of addressing mesoscale cloud-state transitions in the quantification of aerosol radiative forcing that cannot be attained from traditional climate models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

SICOPOLIS-AD v1: an open-source adjoint modeling framework for ice sheet simulation enabled by the algorithmic differentiation tool OpenAD

We present a new capability of the ice sheet model SICOPOLIS that enables flexible adjoint code generation via source transformation using the open-source algorithmic differentiation (AD) tool OpenAD.The adjoint code enables efficient calculation of the sensitivities of a scalar-valued objective function or quantity of interest (QoI) to a range of important, often spatially varying and uncertain model input variables, including initial and boundary conditions, as well as model parameters.Compared to earlier work on the adjoint code generation of SICOPOLIS, our work makes several important advances:(i) it is embedded within the up-to-date trunk of the SICOPOLIS repository – accounting for 1.5 decades of code development and improvements – and is readily available to the wider community;(ii) the AD tool used, OpenAD, is an open-source tool;(iii) the adjoint code developed is applicable to both Greenland and Antarctica, including grounded ice as well as floating ice shelves, with an extended choice of thermodynamical representations.A number of code refactorization steps were required. They are discussed in detail in an Appendix as they hold lessons for the application of AD to legacy codes at large.As an example application, we examine the sensitivity of the total Antarctic Ice Sheet volume to changes in initial ice thickness, austral summer precipitation, and basal and surface temperatures across the ice sheet.Simulations of Antarctica with floating ice shelves show that over 100 years of simulation the sensitivity of total ice sheet volume to the initial ice thickness and precipitation is almost uniformly positive, while the sensitivities to surface and basal temperature are almost uniformly negative. Sensitivity to austral summer precipitation is largest on floating ice shelves from Queen Maud to Queen Mary Land. The largest sensitivity to initial ice thickness is at outlet glaciers around Antarctica. Comparison between total ice sheet volume sensitivities to surface and basal temperature shows that surface temperature sensitivities are higher broadly across the floating ice shelves, while basal temperature sensitivities are highest at the grounding lines of floating ice shelves and outlet glaciers. A uniformly perturbed region of East Antarctica reveals that, among the four control variables tested here, total ice sheet volume is the most sensitive to variations in austral summer precipitation as formulated in SICOPOLIS.Comparison between adjoint- and finite-difference-derived sensitivities shows good agreement, lending confidence that the AD tool is producing correct adjoint code.The new modeling infrastructure is freely available at http://www.sicopolis.net (last access: 2 April 2020) under the development trunk.

58 GEOSCIENCES↗
DOE OSTI2021

Effects of spatial resolution on WRF v3.8.1 simulated meteorology over the central Himalaya

The sensitive ecosystem of the central Himalayan (CH) region, which is experiencing enhanced stress from anthropogenic forcing, requires adequate atmospheric observations and an improved representation of the Himalaya in the models. However, the accuracy of atmospheric models remains limited in this region due to highly complex mountainous topography. This article delineates the effects of spatial resolution on the modeled meteorology and dynamics over the CH by utilizing the Weather Research and Forecasting (WRF) model extensively evaluated against the Ganges Valley Aerosol Experiment (GVAX) observations during the summer monsoon. The WRF simulation is performed over a domain (d01) encompassing northern India at 15 km x 15 km resolution and two nests (d02 at 5 km x 5 km and d03 at 1 km x 1 km) centered over the CH, with boundary conditions from the respective parent domains. WRF simulations reveal higher variability in meteorology, e.g., relative humidity (RH = 70.3%-96.1%) and wind speed (WS = 1.1-4.2m s(-1)), compared to the ERA-Interim reanalysis (RH D 80.0%-85.0%, WS D 1.2-2.3m s(-1)) over northern India owing to the higher resolution. WRF-simulated temporal evolution of meteorological variables is found to agree with balloon-borne measurements, with stronger correlations aloft (r = 0.44-0.92) than those in the lower troposphere (r = 0.18-0.48). The model overestimates temperature (warm bias by 2.8 degrees C) and underestimates RH (dry bias by 6.4 %) at the surface in d01. Model results show a significant improvement in d03 (P = 827.6 hPa, T = 19.8 degrees C, RH = 92.3 %), closer to the GVAX observations (P = 801.4 hPa, T = 19.5 degrees C, RH = 94.7 %). Interpolating the output from the coarser domains (d01, d02) to the altitude of the station reduces the biases in pressure and temperature; however, it suppresses the diurnal variations, highlighting the importance of well-resolved terrain (d03). Temporal variations in near-surface P, T, and RH are also reproduced by WRF in d03 to an extent (r>0:5). A sensitivity simulation incorporating the feedback from the nested domain demonstrates the improvement in simulated P, T, and RH over the CH. Our study shows that the WRF model setup at finer spatial resolution can significantly reduce the biases in simulated meteorology, and such an improved representation of the CH can be adopted through domain feedback into regional-scale simulations. Interestingly, WRF simulates a dominant easterly wind component at 1 km x 1 km resolution (d03), which is missing in the coarse simulations; however, the frequency of southeasterlies remains underestimated. The model simulation implementing a highresolution (3 s) topography input (SRTM) improved the prediction of wind directions; nevertheless, further improvements are required to better reproduce the observed local-scale dynamics over the CH.

Singh, Jaydeep↗