scholarly journals Accounting for Unresolved Spatial Variability in Large Scale Models: Development and Evaluation of a Statistical Cloud Parameterization with Prognostic Higher Order Moments

2011 ◽  
Author(s):  
Robert Pincus
2010 ◽  
Vol 67 (3) ◽  
pp. 834-850 ◽  
Author(s):  
Cara-Lyn Lappen ◽  
David Randall ◽  
Takanobu Yamaguchi

Abstract In 2001, the authors presented a higher-order mass-flux model called “assumed distributions with higher-order closure” (ADHOC 1), which represents the large eddies of the planetary boundary layer (PBL) in terms of an assumed joint distribution of the vertical velocity and scalars. In a subsequent version (ADHOC 2) the authors incorporated vertical momentum fluxes and second moments involving pressure perturbations into the framework. These versions of ADHOC, as well as all other higher-order closure models, are not suitable for use in large-scale models because of the high vertical and temporal resolution that is required. This high resolution is needed mainly because higher-order closure (HOC) models must resolve discontinuities at the PBL top, which can occur anywhere on a model’s Eulerian vertical grid. This paper reports the development of ADHOC 3, in which the computational cost of the model is reduced by introducing the PBL depth as an explicit prognostic variable. ADHOC 3 uses a stretched vertical coordinate that is attached to the PBL top. The discontinuous jumps at the PBL top are “hidden” in the layer edge that represents the PBL top. This new HOC model can use much coarser vertical resolution and a longer time step and is thus suitable for use in large-scale models. To predict the PBL depth, an entrainment parameterization is needed. In the development of the model, the authors have been led to a new view of the old problem of entrainment parameterization. The relatively detailed information available in the HOC model is used to parameterize the entrainment rate. The present approach thus borrows ideas from mixed-layer modeling to create a new, more economical type of HOC model that is better suited for use as a parameterization in large-scale models.


Energies ◽  
2021 ◽  
Vol 14 (19) ◽  
pp. 6318
Author(s):  
Dan Gabriel Cacuci

This work aims at underscoring the need for the accurate quantification of the sensitivities (i.e., functional derivatives) of the results (a.k.a. “responses”) produced by large-scale computational models with respect to the models’ parameters, which are seldom known perfectly in practice. The large impact that can arise from sensitivities of order higher than first has been highlighted by the results of a third-order sensitivity and uncertainty analysis of an OECD/NEA reactor physics benchmark, which will be briefly reviewed in this work to underscore that neglecting the higher-order sensitivities causes substantial errors in predicting the expectation and variance of model responses. The importance of accurately computing the higher-order sensitivities is further highlighted in this work by presenting a text-book analytical example from the field of neutron transport, which impresses the need for the accurate quantification of higher-order response sensitivities by demonstrating that their neglect would lead to substantial errors in predicting the moments (expectation, variance, skewness, kurtosis) of the model response’s distribution in the phase space of model parameters. The incorporation of response sensitivities in methodologies for uncertainty quantification, data adjustment and predictive modeling currently available for nuclear engineering systems is also reviewed. The fundamental conclusion highlighted by this work is that confidence intervals and tolerance limits on results predicted by models that only employ first-order sensitivities are likely to provide a false sense of confidence, unless such models also demonstrate quantitatively that the second- and higher-order sensitivities provide negligibly small contributions to the respective tolerance limits and confidence intervals. The high-order response sensitivities to parameters underlying large-scale models can be computed most accurately and most efficiently by employing the high-order comprehensive adjoint sensitivity analysis methodology, which overcomes the curse of dimensionality that hampers other methods when applied to large-scale models involving many parameters.


2021 ◽  
Vol 502 (3) ◽  
pp. 3976-3992
Author(s):  
Mónica Hernández-Sánchez ◽  
Francisco-Shu Kitaura ◽  
Metin Ata ◽  
Claudio Dalla Vecchia

ABSTRACT We investigate higher order symplectic integration strategies within Bayesian cosmic density field reconstruction methods. In particular, we study the fourth-order discretization of Hamiltonian equations of motion (EoM). This is achieved by recursively applying the basic second-order leap-frog scheme (considering the single evaluation of the EoM) in a combination of even numbers of forward time integration steps with a single intermediate backward step. This largely reduces the number of evaluations and random gradient computations, as required in the usual second-order case for high-dimensional cases. We restrict this study to the lognormal-Poisson model, applied to a full volume halo catalogue in real space on a cubical mesh of 1250 h−1 Mpc side and 2563 cells. Hence, we neglect selection effects, redshift space distortions, and displacements. We note that those observational and cosmic evolution effects can be accounted for in subsequent Gibbs-sampling steps within the COSMIC BIRTH algorithm. We find that going from the usual second to fourth order in the leap-frog scheme shortens the burn-in phase by a factor of at least ∼30. This implies that 75–90 independent samples are obtained while the fastest second-order method converges. After convergence, the correlation lengths indicate an improvement factor of about 3.0 fewer gradient computations for meshes of 2563 cells. In the considered cosmological scenario, the traditional leap-frog scheme turns out to outperform higher order integration schemes only when considering lower dimensional problems, e.g. meshes with 643 cells. This gain in computational efficiency can help to go towards a full Bayesian analysis of the cosmological large-scale structure for upcoming galaxy surveys.


Author(s):  
D. Keith Walters ◽  
Greg W. Burgreen ◽  
Robert L. Hester ◽  
David S. Thompson ◽  
David M. Lavallee ◽  
...  

Computational fluid dynamics (CFD) simulations were performed for unsteady periodic breathing conditions, using large-scale models of the human lung airway. The computational domain included fully coupled representations of the orotracheal region and large conducting zone up to generation four (G4) obtained from patient-specific CT data, and the small conducting zone (to G16) obtained from a stochastically generated airway tree with statistically realistic geometrical characteristics. A reduced-order geometry was used, in which several airway branches in each generation were truncated, and only select flow paths were retained to G16. The inlet and outlet flow boundaries corresponded to the oronasal opening (superior), the inlet/outlet planes in terminal bronchioles (distal), and the unresolved airway boundaries arising from the truncation procedure (intermediate). The cyclic flow was specified according to the predicted ventilation patterns for a healthy adult male at three different activity levels, supplied by the whole-body modeling software HumMod. The CFD simulations were performed using Ansys FLUENT. The mass flow distribution at the distal boundaries was prescribed using a previously documented methodology, in which the percentage of the total flow for each boundary was first determined from a steady-state simulation with an applied flow rate equal to the average during the inhalation phase of the breathing cycle. The distal pressure boundary conditions for the steady-state simulation were set using a stochastic coupling procedure to ensure physiologically realistic flow conditions. The results show that: 1) physiologically realistic flow is obtained in the model, in terms of cyclic mass conservation and approximately uniform pressure distribution in the distal airways; 2) the predicted alveolar pressure is in good agreement with previously documented values; and 3) the use of reduced-order geometry modeling allows accurate and efficient simulation of large-scale breathing lung flow, provided care is taken to use a physiologically realistic geometry and to properly address the unsteady boundary conditions.


2017 ◽  
Vol 50 (1) ◽  
pp. 3287-3293 ◽  
Author(s):  
Erik Frisk ◽  
Mattias Krysander ◽  
Daniel Jung

Sign in / Sign up

Export Citation Format

Share Document