Temperature Boundary Research Articles

Numerical Study on Composite Multilayer Insulation Material for Liquid Hydrogen Storage

This study investigated the heat transfer characteristics of composite multilayer insulation (MLI) materials used for the storage and transport of liquid hydrogen at cryogenic temperatures. This research focused on analyzing the effects of thermal boundary temperature, total layer count, and vacuum level on the heat flux through the insulation material. Based on the layer-by-layer model, a heat transfer model of composite MLI was constructed. This research introduces a novel method for analyzing the heat transfer properties of composite MLI in the liquid hydrogen temperature range. Results indicate that heat flux increases with higher thermal boundary temperatures, with the MLI layers near the cold boundary playing a critical role in overall insulation performance. Additionally, numerical analysis was conducted to examine the impact of different material combinations and variations in vacuum level on heat transfer characteristics. Findings reveal that adding spray-on foam insulation reduces heat flux by 20.76% compared to using MLI alone. Furthermore, increasing the total number of MLI layers effectively mitigates heat flux increase, achieving an optimal heat flux of 0.5377 W/m2 with a total of 50 layers.

A fully coupled multiscale phase-change model at the porous interface for transpiration cooling: coupling dynamics pore-scale networks to continuum-scale free flow

Being a promising method for thermal protection system, transpiration cooling has received wide interest recently. Numerical simulation for transpiration cooling, however, has been limited due to the mis-matching between the pore-scale two-phase flow and external high-temperature aerothermodynamic environment, which induced by conventional decoupled or iterative one-way coupled simplification methods. In this work, a fully coupled continuum-scale and pore-scale model is established for transient transpiration cooling at the interface between boundary layer flow and porous medium through coupling Computational Fluid Dynamics (CFD) and Pore-Network Model (PNM), termed as the multiscale CFD-PNM coupled method. The coupled method allows to capture detailed displacement and phase change of two-phase flow at pore-scale, revealing the strong interaction of the water vapor with the external free flow within a high temperature boundary layer. After successfully validating the new coupled model by comparing with the Two-Phase Mixture Model (TPMM) solution, a number of cases mimicking the cooling at the interface of typical blunt bodies are simulated. The results show that the multiscale CFD-PNM coupled method can not only provide the thermal protection effect prediction, but also reveal many critical features that beyond the reach of continuum-scale studies. Some pore-scale phenomena that are important to the overall transpiration cooling effects are revealed, including transient phase change and composition variation of water vapor, imbibition and drainage of two-phase flow related to pore-scale capillary thresholds and applied boundary pressures, as well as the two-way mass transfer at the interface, such as the invasion of external hot air.

Numerical Investigation of Flow Inside a Channel with Elastic Vortex Generator and Elastic Wall for Heat Transfer Enhancement

In the present investigation, a detailed numerical investigation of the flow and heat transfer characteristics of a channel with an elastic fin (vortex generator) and an elastic wall has been carried out using finite element method. The Fluid-Structure Interaction (FSI) model is used to capture the interaction between the fluid and the solid structure. A sinusoidal time dependent velocity profile has been imposed at the inlet of the channel and the right half of the upper wall of the channel is heated and exposed to constant temperature boundary condition. Due to the sinusoidal velocity profile at the inlet, the elastic fin oscillates periodically and act as a vortex generator, which causes more turbulence in the flow. The obtained results showed that the Nusselt number over the heated wall is affected by the position of the flexible fin, height of flexible fin and elasticity modulus of elastic fin. Moreover, due to the elasticity of the elastic wall and sinusoidal behavior of the inlet velocity, the elastic wall oscillates periodically upward and downward. The Nusselt number values over the heated wall are increased with decrease of the elastic modulus value of the elastic wall. However, the decrease in elastic modulus value of the elastic wall contributes to an increase in the pressure drop inside the channel. It should be added that the interplay between the fluid motion and the deformable structures leads to enhanced turbulence, as the flexible fin and elastic wall introduce additional disturbances and fluctuations into the flow regime. Consequently, this heightened turbulence level has profound implications for heat transfer processes within the system.

Effect of Water-based Alumina-Copper MHD Hybrid Nanofluid on a Power-Law Form Stretching/Shrinking Sheet with Joule Heating and Slip condition: Dual Solutions Study

The application of hybrid nanofluid is now being employed to augment the efficiency of heat transfer rates. A numerical study was conducted to investigate the flow characteristics of water-based-alumina copper hybrid nanofluids towards a power-law form stretching/shrinking sheet. This study also considered the influence of magnetic, Joule heating, and thermal slip parameters. This study is significant because it advances our understanding of hybrid nanofluids in the presence of magnetic fields, power-law form stretching/shrinking sheet, and heat transfer mechanisms, providing valuable insights for optimizing and innovating thermal management systems in various industrial applications such as polymers, biological fluids, and manufacturing processes like extrusion, plastic and metal forming, and coating processes. The main objective of this study is to examine the impact of specific attributes, including suction and thermal slip parameters on temperature and velocity profiles. In addition, this exploration examined the reduced skin friction and reduced heat transfer in relation to the solid volume fraction copper and magnetic effects on shrinkage sheet and thermal slip parameter on suction effect. To facilitate the conversion of a nonlinear partial differential equation into a collection of ordinary differential equations, it is necessary to incorporate suitable similarity variables into the transformation procedure. The MATLAB bvp4c solver application is utilized in the conclusion process to solve ordinary differential equations. No solution was found in the sort of when , and . As the intensity of the Eckert number increases, the temperature profile and boundary layer thickness also increase. The reduced heat transfer rate upsurged in both solutions for solid volume fraction copper for shrinking sheet, while the opposite actions can be noticed in both solutions for thermal slip parameter for suction effect. Finally, the study conducted an analysis to identify two distinct solutions for shrinking sheet and suction zone, while considering different parameter values for the copper volume fractions, magnetic and thermal slip condition effect.

The Spatiotemporal Dynamics of Temperature Variability Across Mts. Qinling: A Comparative Study from 1971 to 2022

Analyzing the spatiotemporal patterns of atmospheric temperature in sensitive areas is critically important for understanding the broader implications of global climate change, which remains a prominent topic in geosciences. It also plays a crucial role in advancing sustainable development. This study utilized daily minimum, maximum, and mean temperature data from twelve meteorological stations across the South and North Mts. Qinling (Qinling Mountains). Employing trend analysis, the Mann–Kendall mutation test, and Morlet wavelet analysis, we explored the predominant temperature trends and characteristics from 1971 to 2022. Our findings revealed consistent inter-annual warming trends in both regions, with more rapid temperature increases in the North compared to the South. Notably, significant shifts occurred in 2003 for both mean and minimum temperatures in the North, while the maximum and minimum temperature values were recorded in the 2010s and 1980s, respectively. Both regions exhibited a primary temperature fluctuation cycle of 28 years. Seasonally, the strongest warming effects appeared in spring, with the weakest in autumn, and moderate effects in winter and summer, indicating that spring contributes most significantly to regional warming. Monthly analysis showed positive temperature trends across all months, with higher rates in the North. The weakening temperature boundary effect of the Mts. Qinling suggested a weakening North–South division, particularly highlighted by the northward shift of the 1 °C isotherm curve for the coldest month, moving away from the previously observed 0 °C isotherm. This northward shift highlights the differential warming rates between the northern and southern regions. Overall, the analysis confirms a robust warming trend, with notable fluctuations in January’s temperatures since 1998, suggesting the Mts. Qinling’s emerging role as a climatic divider in the Chinese Mainland. This introduces new challenges for regional ecosystems, agricultural production, and water resource management, highlighting the pressing need to advance regional sustainable development in the face of climate change.

Multi-objective optimization of porous fins in closed cavity with heat source of sinusoidal periodic temperature boundary condition

Multi-objective optimization of porous fins in closed cavity with heat source of sinusoidal periodic temperature boundary condition

Bioinspired Hydro- and Hydrothermally Responsive Tubular Soft Actuators.

Soft actuators made of thermoresponsive polymers have great potential for intelligent robotics and biomedical devices due to their reversible deformation capability in response to temperature fluctuations. However, they are constrained by a predefined phase transition temperature, limited directional deformation, and nonbiocompatible formulations, thereby restricting their practical utility. Herein a new biomimicry approach is presented to overcome these limitations by developing hydro- and hydrothermally responsive soft actuators made of biocompatible and pliable materials i.e. cotton yarn and polyurethane. We mimic the tubular shape of elephant trunks with their unique muscle orientation by embedding a helical cotton yarn within a hydrophilic polyurethane tube, followed by targeted surface patterning. Unlike the narrow-range shape morphing across the phase transition temperature boundary of typical thermoresponsive hydrogel actuators, we harness hydrothermal stiffness variations in polyurethane to obtain consistent morphing capabilities over a much wider temperature range. The developed actuators can perform versatile activities such as linear, bending, curvilinear, and rotating movements, overcoming the unidirectional motion limitations of conventional soft actuators. The cell viability assay on the building block materials also confirms the high biocompatibility of the actuators. The reported facile fabrication strategy provides new insights for designing complex yet free-standing soft actuators from readily available supple materials.

Anomalous heat transfer enhancement in separated flow over a zigzag-shaped dense package of inclined grooves in a channel wall at different temperature boundary conditions

Rapid development of the anomalous enhancement of separated turbulent Re = 6000 air flow and heat transfer in an in-line single-row package of 31 inclined grooves, 0.2 in dimensionless depth, in a singled-out longitudinal region of the wall of a narrow channel is studied. It is due to the interference of vortex wakes behind the grooves and the acceleration in the channel flow core with the formation of a zone of ultrahigh longitudinal velocity. The wave-shaped parameter characteristics are stabilized in the region of approximately 15th groove, whereupon the oscillation amplitudes are moderately reduced. The return flows in the grooves are enhanced with distance from the entry section, the minimum negative friction diminishing from −2 to −4. The total relative heat removal from the structured region increases atq= const by a factor of approximately 2.75 and by the factor of two atT= const with increase in the relative hydraulic losses by the factor of 1.7, as compared with the case of a plane–parallel channel. The relative heat removal from the surface bounded by the contour of the 20th inclined groove amounts to 3.7 (q= const) with increase in the hydraulic losses by the factor of 2.2. An increase in the local maximum of the longitudinal velocity up to a factor of 1.5, as compared with the mean-mass velocity, can be observable.

Method and Verification of Liquid Cooling Heat Dissipation Based on Internal Heat Source of Airborne Long-Focus Aerial Camera

The traditional passive heat dissipation method has low heat dissipation efficiency, which is not suitable for the heat dissipation of the concentrated heat source inside the long-focal aerial camera, resulting in temperature level changes and temperature gradients in the optical system near the heat source, which seriously affect the imaging performance of the aerial camera. To solve this problem, an active heat dissipation method of liquid cooling cycle is proposed in this paper. To improve the solving efficiency and ensure simulation accuracy, a dynamic boundary information transfer method based on grid area weighting is proposed. The thermal simulation results show that the liquid cooling method reduces the heat source temperature by 70.12%, and the boundary temperature transfer error is 0.015%. The accuracy of thermal simulation is verified by thermal test, and the simulation error is less than 6.44%. In addition, the performance of the optical system is further analyzed, and the results show that the MTF of the optical system is increased from 0.077 to 0.194 under the proposed active liquid cooling cycle heat dissipation method.

GTP-24-1666 - Temperature Margin Calculation During Finite Element Creep Simulations (GT2024-121776)

Abstract Creep deformation in turbomachinery applications is a highly non-linear phenomenon where +/- 15°C may halve or double the average rupture life. Even in well-controlled laboratory conditions, this stochastic nature means that identical tests may differ by a factor of two. Analysts using implicit finite element user creep subroutines may not appreciate the uncertainty of their solution, or the design changes required to account for uncertainty in either the boundary conditions or the solution itself. Designing to a maximum strain limit does not consider the sudden acceleration of tertiary creep rates. Applying a time factor may extend the simulation time by multiple factors. This paper presents a methodology where two estimations of creep damage are made in parallel to the creep simulation at actual operating temperatures. These two estimations consider the predicted change in stress relaxation, at higher or lower temperatures, to estimate the temperature change required to reach the onset of creep, at any given node in the model, and at any given time in the analysis. This margin calculation is expanded to consider the uncertainty in creep prediction. A time-based scatter factor is incorporated into the temperature margin calculation to provide a minimum temperature margin that includes model uncertainty. The analyst can also consider the uncertainty of the temperature prediction and boundary conditions to produce a robust creep prediction in a real-world simulation. The methodology is validated through the FEA of example cases and applied to a creep-limited 2nd stage turbine blade.

Evaluation of Boundary Layer Characteristics at Mount Si’e Based on UAV and Lidar Data

The atmospheric boundary layer is a crucial transitional region connecting the surface with the free atmosphere, playing a bridging role in land-sea-air interactions and the interactions between different atmospheric layers. This study utilizes rotary-wing UAVs, high-resolution lidar, and WRF simulation data to analyze the vertical distribution characteristics of temperature, humidity, wind speed, and wind direction boundary layer over the Mount Si’e region in 4–6 April 2024. The results indicate that the boundary layer temperature decreases with increasing altitude, reaching up to 18°C, while humidity decreases with height, dropping to as low as 35%. Daytime wind speeds range from 4 to 8 m/s, decreasing to 2 to 4 m/s at night. The boundary layer height can reach up to 900 m during the day and drops to 100–200 m at night, showing distinct diurnal variation characteristics. UAV observations are in good agreement with lidar and WRF simulation results, highlighting the application value of UAVs in high temporal and spatial resolution boundary layer studies. The study also reveals the significant impact of complex terrain on boundary layer characteristics, providing scientific insights into the dynamic and thermal processes of the boundary layer and offering reference value for improving regional weather forecasting and numerical simulations.

Validation of an Enhanced Drinking Water Temperature Model during Distribution

Drinking water temperatures are expected to increase in the Netherlands due to climate change and the installation of district heating networks as part of the energy transition. To determine effective measures to prevent undesirable temperature increases in drinking water, a model was developed. This model describes the temperature in the drinking water distribution network as a result of the transfer of heat from the climate and above and underground heat sources through the soil. The model consists of two coupled applications. The extended soil temperature model (STM+) describes the soil temperatures using a two-dimensional finite element method that includes a drinking water pipe and two hot water pipes coupled with a micrometeorology model. The extended water temperature model (WTM+) describes the drinking water temperature as a function of the surrounding soil temperature (the boundary temperature resulting from the STM+), the thermal sphere of influence where the drinking water temperature influences the soil temperature, and the hydraulics in the drinking water network. Both models are validated with field measurements. This study describes the WTM+. Previous models did not consider the cooling effect of the drinking water on the surrounding soil, which led to an overestimation of the boundary temperature and how quickly the drinking water temperature reaches this boundary temperature. The field measurements show the improved accuracy of the WTM+ when considering one to two times the radius of the drinking water pipe as the thermal sphere of influence around the pipe.

Analysis of variables in finite element analysis for multi-pass welding in nuclear power plant components, part 1: Comprehensive study with experimental validation using scientific weld specimen

The influence of various parameters in welding simulations using finite element analysis (FEA) was investigated in this study. Welding experiments and finite element simulations were conducted to analyze residual stress and post-welding distortion. It was a specially designed V-groove specimen, filled with bead passes using Gas Tungsten Arc Welding (GTAW), that facilitated the measurement of distorted shapes and residual stress. Advanced measurement techniques, including three-dimensional laser scanning and the contour method, were used to capture detailed data on the welded specimen. Three different heat source models (body heat flux, temperature boundary, and Goldak's double-ellipsoidal model) were employed in the simulations, and their effect on residual stress distributions was analyzed. Additionally, the effects of hardening models, bead deposition methods, bead shape modeling, welding sequences, and lumping of weld passes on simulation were systematically analyzed. The findings suggest that a combined approach using different modeling techniques and careful calibration of parameters can enhance the reliability of simulations in predicting weld residual stress and distortions.

Heat transfer characteristics of mixed convection inside a differentially heated square cavity containing an oscillating porous cylinder

The current study presents mixed convective flow inside a square chamber holding a centrally placed and thermally conductive oscillating porous circular cylinder. The left boundary's temperature is kept larger than that of the right edge, and the horizontal edges are preserved at adiabatic settings. The circumferential speed of the cylinder is sinusoidal and oscillating in nature. The fluid region within the chamber is modeled employing 2D Navier-Stokes and heat energy equations. Furthermore, the fluid circulation and heat transmission within the porous cylinder are modeled using the Darcy-Brinkman-Forchheimer formulation. The leading equations are discretized utilizing the Galerkin finite element technique. The parametric study is undertaken considering three distinct diameters of the porous cylinder and three distinct oscillation frequencies. The instant Nusselt number is evaluated along the heated wall, which varies in an oscillatory pattern owing to the repeated contraction and enlargement of the thermal boundary layer. The Nusselt number is averaged over time once the value becomes statistically stationary. The study is conducted within a mixed convection region with Reynolds (Re = 100), Richardson (0.1 ≤ Ri ≤ 10), and Grashof (103 ≤ Gr ≤ 105) numbers. Upon thorough examination, it becomes clear that the system's thermal performance shows promising improvement with the largest cylinder diameter and the lowest oscillation frequency. Specifically, the average Nusselt number shows a maximum improvement of 21.50 % at the largest cylinder diameter.

Thermal optimization of PCM‐based heat sink using fins: A combination of CFD, genetic algorithms, and neural networks

AbstractSustainable development requires focusing on renewable energy, such as solar power, due to the limited availability of fossil fuels. Solar energy, available only during daylight, can be stored in thermal batteries using phase change materials (PCM). However, PCM has low thermal conductivity, slowing its melting process. Fins have been added to enhance heat diffusion, but their optimal design requires further study. This research examines an annular chamber with fixed inner and outer wall temperatures and explores fin configurations. Two series of evenly spaced fins (four and five) were analyzed, with melting times calculated using CFD simulations. Larger fins reduce melting time but limit PCM storage. Effectiveness, defined as the PCM space over melting time, was used to evaluate performance. Higher boundary temperatures decreased melting time but did not always increase effectiveness. The study also used artificial neural networks (ANN) combined with a genetic algorithm to predict optimal conditions. The configuration with five fins, a 90°C boundary temperature, a fin length of 45 mm, and a thickness of 5.4 mm achieved the highest effectiveness of 4.7, showing that smaller fins at lower temperatures are more efficient.

Kinetic Monte Carlo modelling of nano-oxide precipitation and its associated stability under neutron irradiation for the Fe-Ti-Y-O system

While developing nuclear materials, predicting their behavior under long-term irradiation regimes spanning decades poses a significant challenge. We developed a novel Kinetic Monte Carlo (KMC) model to explore the precipitation behavior of Y-Ti-O oxides along grain boundaries within nanostructured ferritic alloys (NFA). This model also assessed the response of the oxides to neutron irradiation, even up simulated radiation damage levels in the desired long dpa range for reactor components. Our simulations investigated how temperature and grain boundary sinks influenced the oxide characteristics of a 12YWT-like alloy during heat treatments at 1023, 1123, and 1223 K. The oxide characteristics observed in our simulations were in good agreement with existing literature. Furthermore, the impact of grain boundaries on precipitation was found to be minimal. The resulting oxide configurations and positions were used in subsequent simulations that exposed them to simulated neutron irradiation to a total accumulated dose of 8 dpa at three temperatures: 673, 773, and 873 K, and at dose rates of 10−3, 10−4, and 10−5 dpa/s. This demonstrated the expected inverse relationship between oxide size and dose rate. In a long-term irradiation simulation at 873 K and 10−3 dpa/s was taken out to 66 dpa and found the oxides in the vicinity of the grain boundary were more susceptible to dissolution. Additionally, we conducted irradiation simulations of a 14YWT-like alloy to reproduce findings from neutron irradiation experiments. The larger oxides in the 14YWT-like alloy did not dissolve and displayed stability similar to the experimental results.

Evidence of a Forming Nucleus in the Fourcade–Figueroa Galaxy

We analyze data from the IRAS, Wide-field Infrared Survey Explorer, and Planck satellites, revealing an unresolved dust condensation at the center of the Fourcade–Figueroa galaxy (ESO270-G017), which may correspond to a forming nucleus. We model the condensation’s continuum spectrum in the spectral range from 3 to 1300 μm using the DUSTY code. The best-fit model, based on the chi-square test, indicates that the condensation is a shell with an outer temperature of T out ≈ 12 K and an inner boundary temperature of T i ≈ 500 K. The shell’s outer radius is r o = 86.2 pc, and the inner cavity radius is r i = 0.082 pc. The condensation produces an extinction A V = 50 mag, and its luminosity is L c = 1.08 × 1034 W, which would correspond to a burst of massive star formation approximately similar to the central 5 pc of R 136 in the LMC and NGC 3603, the ionizing cluster of a giant Carina arm H ii region. The comparison with normal, luminous, and ultraluminous IR galaxies leads us to consider this obscured nucleus as the nearest and weakest object of this category.

Response of Convectively Coupled Kelvin Waves to Surface Temperature Forcing in Aquaplanet Simulations

AbstractThis study investigates changes in the propagation and maintenance of convectively coupled Kelvin waves (KWs) in response to surface warming. We use a set of three aquaplanet simulations made with the Community Atmospheric Model version 6 by varying the sea surface temperature boundary conditions to represent the current climate as well as warmer (+4 K) and cooler (−4 K) climates. Results show that KWs accelerate at the rate of about 7.1%/K and their amplitudes decrease by 4.7%/K. The dampening of KWs with warming is found to be associated with a weakening of the internal thermodynamic feedback between diabatic heating and temperature anomalies that generates KW eddy available potential energy (EAPE). The phase speed of KWs closely matches that of the second baroclinic mode KW in −4 K, while the phase speed of KWs is approximately that of the first baroclinic mode KW in +4 K. Meanwhile, the coupling between the two baroclinic modes weakens with warming. We hypothesize that in −4 K, as the first and second baroclinic modes are strongly coupled, KWs destabilize by positive EAPE generation within the second baroclinic mode and propagate more slowly, following the second baroclinic mode KW phase speed. In +4 K, as the first and second baroclinic modes decouple, KWs are damped by negative EAPE generation within the first baroclinic mode and propagate faster, following the first baroclinic mode KW phase speed.

Dynamic analysis of carbon fiber-wound reinforced composite truncated conical shell under thermal conditions: Theory and experimental validation

This work aims to investigate the dynamic behavior of a carbon fiber-wound reinforced composite (FW-RC) truncated conical shell structure under various operating conditions, including transverse excitation and winding pre-stress in thermal environment. In particular, a dynamic model for the FW-RC truncated conical shell is established, taking into account the first-order shear deformation theory and the von Karman geometric relationship. By utilizing Hamilton's principle and Galerkin method, the dynamic equations for the FW-RC truncated conical shells are derived. The natural frequencies of the system are afterward determined, and a vibration test bench is constructed specifically for the truncated conical shells. The validity of the theoretical model is verified by measuring the natural frequency of the specimen. In addition, the effects of various factors such as temperature field, pre-stress, boundary conditions, geometric parameters, and fiber winding parameters on the inherent characteristics and dynamic behavior of the system are further studied. By delving into these aspects, a comprehensive understanding of the dynamic behavior of FW-RC truncated conical shell structures can be achieved, which is of great significance for their design and application in practical engineering scenarios.

Sensitivity of horizontal resolution and land surface model in operational WRF forecast for Online Nuclear Emergency Response System (ONERS)

Accurate Meteorological forecasts are crucial for the assessment of plume dispersion and dose prediction in nuclear power plant (NPP) sites. In this work the forecast sensitivity of the Weather Research and Forecasting (WRF) model is tested by running a series of forecast simulations for horizontal resolution, and land surface models (LSM) in the context of Online Nuclear Emergency Response System (ONERS) for Indian NPP sites. 72 h forecast simulations are made for three seasons viz. summer, southeast and northeast monsoon using the Global Forecast data. Three simulation experiments, namely 2 km-NOAH, 3 km-NOAH and 3 km-NOAHMP are conducted using two different nested domain configurations (18–6–2 km and 9–3 km) and two LSM schemes (NOAH and NOAH-MP) and tested at four different NPP sites. Forecast comparison of surface winds, relative humidity, temperature, heat fluxes and planetary boundary layer heights with data from meteorological tower, radiosonde and the Modern-Era Retrospective analysis for Research and Applications 2 (MERRA-2) shows 3 km-NOAH is equally capable in predicting surface parameters as well as vertical profiles compared to 2 km-NOAH with marginal differences. 3 km-NOAHMP shows less mean bias and better correlation for boundary layer height and heat fluxes. Comparison of spatial flow-field with 5th generation European Centre for Medium-Range Weather Forecasts (ECMWF) reanalysis (ERA5) data shows synoptic scale seasonal winds, sea level pressure systems and temperature hot-spots are better captured by 3 km-NOAHMP compared to 6 km coarse domain in the 18–6–2 km configuration. The daily accumulated rainfall by all simulations is overestimated compared to ERA5 data. The predictions by 3 km-NOAHMP better agree with Integrated Multi-satellite Retrievals for Global Precipitation Measurement (GPM-IMERGE) data whereas 2 km-NOAH predicts delayed rainfall occurrence. Dispersion simulations of hypothetical plume release from a coastal NPP site with all three forecasts properly show the influence of local scale diurnal land-sea breeze and seasonal winds on the plume movement. Therefore the 9–3 km domain with NOAHMP LSM is found to be a suitable choice for operational weather forecast in ONERS for Indian NPP sites.