Two dynamical core formulation flaws exposed by a. baroclinic instability test case
|
|
- Hector Walters
- 5 years ago
- Views:
Transcription
1 Two dynamical core formulation flaws exposed by a baroclinic instability test case David L. Williamson, Jerry G. Olson National Center for Atmospheric Research, Boulder, CO, USA Christiane Jablonowski Department of Atmospheric, Oceanic and Space Sciences, University of Michigan, Ann Arbor, MI, USA March 6, 2008 Submitted to Mon. Wea. Rev. Corresponding author address: National Center for Atmospheric Research (NCAR), P.O. Box 3000, Boulder CO, , USA, wmson@ucar.edu The National Center for Atmospheric Research is sponsored by the National Science Foundation.
2 ABSTRACT Two flaws in the semi-lagrangian algorithm originally implemented as an optional dynamical core in the NCAR Community Atmosphere Model (CAM3.1) are exposed by a steady-state and baroclinic instability test case. Remedies are demonstrated and have been incorporated in the dynamical core. One consequence of the first flaw is an erroneous damping of the speed of a zonally uniform zonal wind undergoing advection by a zonally uniform zonal flow field. It results from projecting the transported vector wind expressed in unit vectors at the arrival point to the surface of the sphere and is eliminated by rotating the vector to be parallel to the surface. The second flaw is the formulation of an a posteriori energy fixer which, although small, systematically affects the temperature field and leads to an incorrect evolution of the growing baroclinic wave. That fixer restores the total energy each time step by changing the provisional forecasted temperature at proportionally to the magnitude of the temperature change that time step. Two other fixers are introduced that do not exhibit the flaw. One changes the provisional temperature everywhere by an additive constant, and the other changes it proportionally by a multiplicative constant. 1
3 1. Introduction Testing the dynamical core of an atmospheric General Circulation Model (GCM) in isolation from the physical parameterization package is an essential stepping stone during the model development and evaluation phase. The dynamical core tests have the potential to assess not only the diffusive characteristics and convergence of the numerical schemes but also can reveal flaws in the dynamical core formulation. Such flaws might not be obvious in weather or climate simulations if physical parameterizations or model components, like the ocean or ice, compensate or mask the inaccuracies. The term dynamical core generally refers to the numerical approximation to the resolved fluid flow component of an atmospheric model. Most commonly, the hydrostatic primitive equations, the non-hydrostatic shallow-atmosphere, and non-hydrostatic deep-atmosphere equation sets are used. The discussion here concentrates on the assessment of a hydrostatic dynamical core. In particular, we discuss two flaws in an early version of the semi-lagrangian spectral transform dynamical core that is part of the Community Atmosphere Model CAM3.1 developed at the National Center for Atmospheric Research (NCAR). The flaws have been detected with the help of two recently developed steady-state and baroclinic wave test cases (Jablonowski and Williamson 2006a,b) and have since been eliminated in the standard release of CAM3.1 (Collins et al. 2004, 2006). The first flaw resulted in a damping of the speed of a zonally uniform zonal wind undergoing advection by a zonally uniform zonal flow field. As will be illustrated in Section 3a, it was caused by the method used to obtain the vector wind at the arrival point which systematically shortened the vector every time step. The second flaw was the formulation of an a posteriori total energy fixer which, although small, systematically affected the temperature field and led to an incorrect evolution of the growing baroclinic wave. This 2
4 paper serves as a precaution for other modeling groups using similar algorithms. It is organized as follows. Section 2 briefly reviews the design and application of the steady-state and baroclinic wave test cases. Section 3 discusses the flaw in the vector wind and a posteriori energy fixer. A summary is provided in section Review of the dynamical core test cases The idealized steady-state and baroclinic instability test cases were suggested by Jablonowski (2004) and recently applied to four very different dynamical cores by Jablonowski and Williamson (2006a,b). Among them were the three dynamical cores that are part of the NCAR CAM3.1 modeling framework. These are the spectral transform Eulerian (EUL) and semi-lagrangian (SLD) models (Collins et al. 2004, 2006) as well as the Finite Volume (FV) dynamical core developed by Lin (2004). In addition, the icosahedral finite-difference model GME of the German Weather Service (DWD) was tested (Majewski et al. 2002). These hydrostatic dynamical cores represent a broad range of numerical approaches and, at very high resolutions, provide independent reference solutions. These are used to assess whether new model runs fall within the uncertainty range of the reference solutions. This is the primary measure for detecting flaws in the model formulation. In addition, the steady-state model runs are compared against the initial state which is an analytic solution to the primitive equations. The steady-state and baroclinic wave test cases have been developed for dry dynamical cores in spherical geometry with pressure-based vertical coordinates. The test cases are based on an analytically defined zonally symmetric zonal flow with two jets in the midlatitudes. In addition, a smooth surface geopotential, a constant surface pressure and a quasi-realistic meridional and 3
5 vertical temperature distribution are prescribed. The meridional wind is set to zero. All analytical expressions and further details can be found in Jablonowski and Williamson (2006a) (referred to as JW hereafter). These initial conditions are a balanced steady-state solution of the primitive equations and satisfy static, inertial and symmetric stability properties. However, they are unstable with respect to baroclinic or barotropic instability mechanisms. A baroclinic wave can be triggered if the initial conditions are overlaid with a perturbation. In particular, we select a small-amplitude but rather wide Gaussian hill perturbation that overlays the zonal wind field at all vertical levels. This triggers the evolution of a baroclinic wave over the course of ten days with explosive cyclogenesis around day Exposed flaws in the semi-lagrangian dynamical core a. Semi-Lagrangian advection of the vector wind The SLD dynamical core is initialized with the balanced initial conditions without a perturbation. This evaluates the model s ability to maintain the steady-state solution and serves as a stringent debugging tool. The model is run for 30 days at a fixed T42L26 (triangular truncation) resolution with varying time steps. This resolution corresponds to a horizontal grid spacing of approximately in the underlying Gaussian grid and utilizes 26 hybrid levels in the vertical direction. The details of the vertical level placement can be found in Jablonowski and Williamson (2006b). No explicit diffusion is used. For these model runs error norms can be directly computed since the analytical solution is known. We observed that the magnitude of the zonally uniform zonal wind in the SLD run decreased 4
6 with time whereas both the EUL and FV dynamical cores maintained the zonal wind magnitude with very minor variations (JW). Figure 1 shows the time series of the global root-mean square l 2 error norms of the zonally averaged zonal wind for the SLD T42L26 simulations. The definition of the l 2 error norm is given in JW (Eq. (15)). The figure depicts the global l 2 error norms for the three time steps t = 900 s, t = 1800 s and t = 3600 s. Both SLD versions with and without the flaw in the vector wind are shown. In the flawed simulations, the error is greater the longer the time step. This contradicts the general perception that the damping arising from semi-lagrangian interpolation is smaller with longer time steps because fewer interpolations are performed for a fixed elapsed time. The error in Fig. 1 is not caused by the damping from the interpolation. Since the flow is zonal and zonally uniform, the departure point is on the same latitude circle as the arrival point and the interpolation is between points in longitude which have the same value. The interpolation scheme should maintain a constant field in that case. The damping is a result of the formulation of the vector momentum equation. In the original CAM3.1 formulation, Williamson and Olson (1994) developed the momentum equation approximations in vector form following Ritchie (1988, 1991) and Bates et al. (1990, 1993). In general, the vector wind which is parallel to the spherical surface at the departure point is no longer parallel when transported to the arrival point. This is schematically drawn in Fig. 2 that shows the wind vector at the departure (1) and arrival points (2-4). The wind vectors at the arrival point depict the transported vector (2) that can either be projected onto the spherical surface (3) or rotated to be parallel to the surface (4). To calculate the vector components needed in a practical implementation, the original CAM3.1 formulation followed Bates et al. (1990) in which the equations for the individual vector components are obtained by relating the unit vectors at the departure points to those at the arrival points and ignoring the vertical components (Eq. (10) of Williamson and Olson 5
7 (1994)). This essentially projects the transported vector onto the spherical surface and thereby shortens it. Thus the zonal wind decreases with time in this test case. With the longer time step, the wind at the arrival point is less parallel to the surface and thus shortened more by the projection to the surface. In contrast, Ritchie (1988, 1991) rotated the vector to be parallel to the spherical surface, maintaining its length. Solutions with this rotated approach are also shown in Fig. 1 for the same three time steps. All errors from the rotated vector cases are clearly reduced and no longer exhibit the dependence on the time step length. They graphically overlay each other and now match those from the Eulerian dynamical core seen in Fig. 4a of JW. This revised formulation was implemented into CAM3.1 (Eq. (3.298), Collins et al. (2004)) and used in the experiments reported in JW. The semi-lagrangian curves presented here with the revised formulation are not the same as presented in Fig. 4b of JW. Those errors were calculated with the model configured as it would be for an actual GCM climate simulation or weather forecast. Here we configured the SLD model to isolate the vector wind problem. In particular here we did not revert to a geodesic trajectory calculation near the poles (Williamson and Rasch 1989) since the flow is zonal. We also did not include a decentering mechanism to minimize orographic resonance (Collins et al. 2004). Both of those approximations contribute to the l 2 error growth seen in Fig. 4a of JW, with the decentering mechanism being the dominant source. b. Total energy fixer The second test imposes a well resolved Gaussian-hill perturbation on the balanced zonal flow. This perturbation grows with time and evolves into a baroclinic wave. We run the test for 30 days with both the spectral transform Eulerian and semi-lagrangian dynamical cores at T170 6
8 triangular truncation and 26 hybrid levels. This resolution corresponds to a horizontal grid spacing of approximately in the underlying Gaussian grid. The time steps are t = 300 s for EUL and t = 900 s for SLD. No explicit diffusion is used in the SLD simulations whereas in EUL second- and fourth-order diffusion operators are applied to the prognostic variables (JW). The SLD simulations include three variants of an a posteriori total energy fixer (FIXER 1-3) plus a simulation without the energy fix. The details of the three energy fixers are explained below. In short, they artificially restore the global conservation of total energy in non-energy-conserving simulations. This is done by adjusting the temperature field after each time step. The different adjustment techniques are labeled as FIXER 1-3. Figure 3 shows the root-mean square l 2 difference of surface pressure between the four semi-lagrangian integrations and the Eulerian run. The definition of the l 2 difference can be found in JW (Eq. (16)). The label FIXER 1 (short dashed line) denotes the original implementation. Note that FIXER 1 lies outside the stippled region which, at this resolution, indicates a flaw in the SLD model simulation. As argued in JW the stippled region defines the uncertainty in the baroclinic wave reference solutions obtained from a set of four very different dynamical cores applied at high horizontal resolution. One of these reference solutions is that from a corrected SLD (FIXER 2) model run. However, differences between the reference solutions with different dynamical cores are introduced by different truncation errors. Although initially these errors might be very small at high resolutions, they do grow with time because the basic state is baroclinically unstable. This error growth limits how well the reference solutions are actually known. The stippled region in Fig. 3 is bounded by the envelope of differences between all combinations of the high resolution solutions (JW). If the differences of one model against one of those reference solutions fall within the stippled region then the solution is captured to within the uncertainty. Differences outside the stippled region indicate that the solution is not as accurate 7
9 as any of those reference solutions. This analysis technique is also explained in more detail in JW. Clearly, the original semi-lagrangian implementation with FIXER 1 does not produce an accurate evolution of the baroclinic wave. This can further be seen in Fig. 4 which plots the surface pressure at day 10 from three T170L26 SLD solutions and the T170L26 EUL run. Even the resolved scales, well away from the truncation limit, are incorrect in the semi-lagrangian FIXER 1 solution. Plots for T340L26 SLD with FIXER 1 (not shown) look very similar to those of T170L26 SLD with FIXER 1 indicating that the flawed SLD simulation is converging to a different solution in comparison to models like EUL, FV, GME and SLD (FIXER 2) at equivalent resolutions. The latter are shown at day 10 in Jablonowski and Williamson (2006b). After much searching we isolated the source of the error in the semi-lagrangian integration to an a posteriori total energy fixer that was included in that dynamical core. For coupled atmosphereocean climate change applications, the atmospheric component of the climate model must conserve energy to tenths of W m 2 or better (Boville 2000). A greater imbalance could produce a drift of the deep ocean in a coupled system which is large enough to imply a non-equilibrium solution. Since the models have finite resolution, a dissipation mechanism must be included to control the build-up of kinetic energy at the smallest resolved scales. Some models, such as the EUL dynamical core, include horizontal diffusion to control the energy at the smallest scales. EUL also includes a frictional heating term corresponding to the explicit momentum diffusion to make the damping of the momentum conservative (Collins et al. 2004). This frictional heating process appears to be a reasonable approach that provides for a greater degree of global energy conservation. However, it is somewhat arbitrary and does not capture the true energetics of the system. In addition, other terms like the diffusion on temperature still lead to a lack of conservation in the tenths of W m 2 range. Other models, such as the FV dynamical core, control the smallest scales 8
10 via monotonicity constraints in the numerical scheme and also utilize a 2 divergence damping mechanism. In semi-lagrangian schemes the interpolants control the energy at the smallest scales to some extent, but explicit diffusion is often applied as well. In our experience, the energy loss associated with these damping mechanisms described above in all three CAM3.1 dynamical cores is around 2 W m 2 at resolutions normally applied to climate simulations (1 to 2 range). The energy loss due to the inherent damping by the numerical methods is not explicitly known. Therefore, an a posteriori energy fixer is generally applied at every time step in all CAM3.1 dynamical cores. However, the fixer was not applied to the EUL reference solution used to calculate the differences in Fig. 3. In the continuous primitive equations total energy is conserved if the rate of change ηs ( v 2 ) p ( t A η top 2 + c p s pt dη da = Φ s η A t Φ p ) top top da (1) t is obtained (Laprise and Girard 1990). This equation is valid in the absence of diabatic and frictional effects for hybrid vertical coordinate systems as used here. Φ s, p s and Φ top, p top are the geopotential and pressure at the surface and the model top, c p is the specific heat of dry air at constant pressure and v = (u, v) T stands for the horizontal velocity vector with the zonal and meridional wind components u and v. Furthermore, T symbolizes the temperature, p is the pressure, and t denotes the time. The integrals span the 3D and 2D domains where A symbolizes the horizontal area of the sphere and η stands for the hybrid vertical coordinate. The vertical integral is bounded by the value η s at the surface and η top at the model top. Here, η s is identical to unity and η top is set p top /p 0 with p 0 = 1000 hpa. Note that η top is non-zero for constant p top > 0 hpa which is the case here. A constant pressure at the model top ensures the global conservation of 9
11 total energy in the continuous equations and simplifies the 2D integral. Eq. (1) then becomes { t A [ ( ) ] } 1 ηs ( v 2 ) p g η top 2 + c pt η dη + Φ s p s da = 0. (2) Here we divided Eq. (1) by the gravitational constant g to recover energy units (Kasahara 1974). This expression is equivalent to t E = 0 where E symbolizes the global integral of the total energy as shown by the term in the curly bracket in Eq. (2). In the semi-discrete system with p/ η p/ η and dη η, the domain-integrated total energy E is given by E = A 1 g [( K k=1 ( vk c pt k ) p k ) + Φ s p s ]da. (3) The summation index k indicates the vertical index of a full model level with the maximum level number K near the surface. The pressure difference p k is defined as p k = p k+1/2 p k 1/2 = p 0 A k + p s B k (4) with A k = A k+1/2 A k 1/2 and B k = B k+1/2 B k 1/2. The discrete positions of the hybrid coefficients A k+1/2 and B k+1/2 at the model interface levels are listed in Jablonowski and Williamson (2006b). η k is given by η k = η k+1/2 η k 1/2 = A k + B k. ( ) Let ˆT +, ˆv +, ˆp + s denote the temperature, horizontal wind vector and surface pressure at the end of a time step and (T, v, p s ) denote the values at the beginning of the time step. In energyconserving model formulations the residual RES = Ê+ E (5) 10
12 would be zero. Ê + and E are given by Ê + = E = A A 1 g 1 g { [ K k=1 { [ K k=1 ( (ˆv + k )2 2 + c p ˆT + k ( (v k )2 + c p T k 2 ) ) (p 0 A k + ˆp + s B k ) (p 0 A k + p s B k ) ] ] + Φ sˆp + s + Φ s p s } } da, (6) da. (7) Note that diabatic energy sources and sinks need to be included in Eq. (7) if simulations with physical parameterizations are considered (Collins et al. 2004). This is not the case here. In general, RES is not zero due to implicit and explicit diffusion processes in the dynamical ( ) cores. Therefore, modifications can be made to the provisional forecast values ˆT +, ˆv +, ˆp + s. This adjustment yields updated values (T +, v +, p + s ) which, if substituted for the provisional values in Eq. (5), yield a zero residual. This is the underlying concept of the energy fixer. The form of the energy fixer first used with the semi-lagrangian model (FIXER 1) only modified the temperature. The modification was proportional to the magnitude of the local change in T at that time step and given by (FIXER 1) T + (λ, ϕ, η) = ˆT + (λ, ϕ, η) + β 1 ˆT + (λ, ϕ, η) T (λ, ϕ, η) (8) where λ and ϕ denote the longitudinal and latitudinal position in the spherical grid. This adjustment follows the spirit of the water vapor fixer developed by Rasch and Williamson (1991) and Williamson and Rasch (1994) for predecessors of the CAM3 model generation. The constant β 1 is determined by replacing ˆT + with T + in Eq. (6) and setting RES = 0 in Eq. (5). Eq. (8) is then substituted for T + in Eq. (6), Eqs. (6) and (7) are plugged into Eq. (5) which is solved for β 1. After determining that FIXER 1 of Eq. (8) was the source of the error in the semi-lagrangian integrations two other fixers were tried, both of which also only modify the temperature field. 11
13 FIXER 2 is given by T + (λ, ϕ, η) = ˆT + (λ, ϕ, η) + β 2 (9) and FIXER 3 is formulated as T + (λ, ϕ, η) = (1 + β 3 ) ˆT + (λ, ϕ, η). (10) The future wind and surface pressure fields are set to v + = ˆv + and p + s = M ˆp + s where M symbolizes a mass fixer for non-mass-conserving model implementations. The details of the mass fixer for EUL and SLD are provided in Collins et al. (2004). They are unimportant for the discussion here. In short, the mass fixer ensures that the global integral of the surface pressure and thereby the total mass is conserved. Note that M is not needed for the mass-conserving model FV. The energy fixer FIXER 2 changes the provisional temperature by a constant, whereas FIXER 3 changes it proportionally. Both constants β 2 and β 3 are determined as described above for FIXER 1. This yields the equation for β 2 ge { [ K ( (ˆv + ) k )2 + A k=1 + c 2 p ˆT k (p0 A k + M ˆp + s B k ) ] } + Φ s M ˆp + s da β 2 = K A k=1 c. (11) p(p 0 A k + M ˆp + s B k ) da FIXER 2 is adopted in the released version of CAM3.1 and used operationally for climate applications. Figure 3 confirms that the SLD FIXER 2 and FIXER 3 simulations correct the problem with FIXER 1. The figure shows the difference curves for the T170L26 Eulerian simulation versus the T170L26 semi-lagrangian runs with FIXER 2, FIXER 3, and with no fixer. These three curves all lie on top of each other and fall with the stippled region of the plot. Thus all three SLD versions calculate the growth of the imposed perturbation to within the uncertainty of the reference solutions. The bottom two panels of Figure 4 depict the surface pressure at day 10 from the solution with no fixer and with FIXER 2. The two are visually indistinguishable from each other and 12
14 both are very similar to the T170L26 Eulerian reference solution. The solution with FIXER 3 (not shown) is also visually indistinguishable from the two other semi-lagrangian solutions. In fact, the l 2 difference between the solutions with FIXER 2 and FIXER 3 is hpa at day 10, which is well below the differences of any of the three SLD simulations with the Eulerian solution. 4. Conclusions The baroclinic instability test case developed by Jablonowski (2004) exposed two flaws in the semi-lagrangian algorithms originally implemented as an optional dynamical core in NCAR s CAM3.1 modeling framework. The first flaw resulted in a damping of the speed of a zonally uniform zonal wind undergoing advection by a zonally uniform zonal flow field. It was caused by projecting the transported wind expressed in unit vectors at the arrival point to the surface of the sphere following Bates et al. (1990, 1993), thus systematically shortening the vector every time step. It is eliminated by rotating the vector to be parallel to the surface following Ritchie (1988, 1991), maintaining its length. The second flaw was the formulation of an a posteriori energy fixer which, although small, systematically affected the temperature field and lead to an incorrect evolution of the baroclinicly growing wave. That fixer changed the provisional temperature forecasted at the end of the time step proportionally to the magnitude of the temperature change that time step. Two other fixers were introduced that do not exhibit the flaw. The dynamical core without an energy fixer also does not exhibit the flaw. The two fixers change the provisional temperature by a constant, or proportionally to the temperature. The former has been adopted by CAM3.1 for all dynamical cores (Collins et al. 2004) and was used in the experiments reported in Jablonowski and Williamson (2006a). 13
15 Acknowledgments. This work was partially supported by the Office of Science (BER), U.S. Department of Energy, Cooperative Agreement No. DE-FC02-97ER
16 REFERENCES Bates, J. R., S. Moorthi, and R. W. Higgins, 1993: A global multilevel atmospheric model using a vector semi-lagrangian finite difference scheme. Part I: Adiabatic formulation. Mon. Wea. Rev., 121, Bates, J. R., F. H. M. Semazzi, R. W. Higgins, and S. R. M. Barros, 1990: Integration of the shallow water equations on the sphere using a vector semi-lagrangian scheme with a multigrid solver. Mon. Wea. Rev., 118, Boville, B. A., 2000: Toward a complete model of the climate system. Numerical Modeling of the Global Atmosphere in the Climate System, P. Mote and A. O Neill, Eds., Kluwer Academic Publishers, Collins, W. D., et al., 2004: Description of the NCAR Community Atmosphere Model (CAM3.0). NCAR Tech. Note NCAR/TN-464+STR, National Center for Atmospheric Research, Boulder, Colorado. 226 pp. Collins, W. D., et al., 2006: The formulation and atmospheric simulation of the Community Atmosphere Model: CAM3. J. Climate, 19, Jablonowski, C., 2004: Adaptive grids in weather and climate modeling. Ph.D. thesis, University of Michigan, Ann Arbor, MI, Department of Atmospheric, Oceanic and Space Sciences, 292 pp. Jablonowski, C. and D. L. Williamson, 2006a: A baroclinic instabilitiy test case for atmospheric model dynamical cores. Quart. J. Roy. Meteor. Soc., 132 (621C),
17 Jablonowski, C. and D. L. Williamson, 2006b: A baroclinic wave test case for dynamical cores of General Circulation Models: Model intercomparisons. NCAR Tech. Note NCAR/TN-469+STR, National Center for Atmospheric Research, Boulder, Colorado. 89 pp. Kasahara, A., 1974: Various vertical coordinate systems used for numerical weather prediction. Mon. Wea. Rev., 102, Laprise, R. and C. Girard, 1990: A spectral General Circulation Model using the piecewiseconstant finite-element representation on a hybrid vertical coordinate system. J. Climate, 3, Lin, S.-J., 2004: A vertically Lagrangian finite-volume dynamical core for global models. Mon. Wea. Rev., 132, Majewski, D., et al., 2002: The operational global icosahedral-hexagonal gridpoint model GME: Description and high-resolution tests. Mon. Wea. Rev., 130, Rasch, P. J. and D. L. Williamson, 1991: The sensitivity of a general circulation model climate to the moisture transport formulation. J. Geophy. Res., 96, Ritchie, H., 1988: Application of the semi-lagrangian method to a spectral model of the shallow water equations. Mon. Wea. Rev., 116, Ritchie, H., 1991: Application of the semi-lagrangian method to a multilevel spectral primitiveequations model. Quart. J. Roy. Met. Soc., 117, Williamson, D. L. and J. G. Olson, 1994: Climate simulations with a semi-lagrangian version of the NCAR Community Climate Model. Mon. Wea. Rev., 122,
18 Williamson, D. L. and P. J. Rasch, 1989: Two-dimensional semi-lagrangian transport with shapepreserving interpolation. Mon. Wea. Rev., 117, Williamson, D. L. and P. J. Rasch, 1994: Water vapor transport in the NCAR CCM2. Tellus A, 46,
19 List of Figures Fig. 1 l 2 norm ū ū t=0 in m s 1 for T42L26 semi-lagrangian integrations with time steps t = 900 s, t = 1800 s and t = 3600 s. Top three lines are from the model with the flaw in the vector transport formulation, bottom three (which have very similar values and thus appear to be a single line) from the model with the modified (rotated) vector wind. Fig. 2 Sketch of the vector winds. (1) Vector at the departure point, (2) transported vector at the arrival point, (3) vector projected onto the sphere (shortened), (4) rotated vector parallel to the sphere (length maintained). Fig. 3 l 2 norm of the surface pressure differences (in hpa) between various versions of the semi-lagrangian model and the Eulerian spectral, all at T170L26 resolution. All curves except FIXER 1 lie on top of each other. Fig. 4 Surface pressure at day 10 from various versions of the semi-lagrangian model and the Eulerian spectral, all at T170L26 resolution. Contour interval is 7.5 hpa with the 980 hpa and 1010 hpa contours thicker. 18
20 t = 3600 s t = 1800 s t = 900 s runs with rotated vector all three time steps l 2 (ū -ū 0 ) Days FIG. 1. l 2 norm ū ū t=0 in m s 1 for T42L26 semi-lagrangian integrations with time steps t = 900 s, t = 1800 s and t = 3600 s. Top three lines are from the model with the flaw in the vector transport formulation, bottom three (which have very similar values and thus appear to be a single line) from the model with the modified (rotated) vector wind. 19
21 Departure point Arrival point FIG. 2. Sketch of the vector winds. (1) Vector at the departure point, (2) transported vector at the arrival point, (3) vector projected onto the sphere (shortened), (4) rotated vector parallel to the sphere (length maintained). 20
22 FIG. 3. l 2 norm of the surface pressure differences (in hpa) between various versions of the semi- Lagrangian model and the Eulerian spectral, all at T170L26 resolution. All curves except FIXER 1 lie on top of each other. 21
23 FIG. 4. Surface pressure at day 10 from various versions of the semi-lagrangian model and the Eulerian spectral, all at T170L26 resolution. Contour interval is 7.5 hpa with the 980 hpa and hpa contours thicker.
A Test Suite for GCMs: An Intercomparison of 11 Dynamical Cores
A Test Suite for GCMs: An Intercomparison of 11 Dynamical Cores Christiane Jablonowski 1, Peter Lauritzen 2, Mark 3 Taylor, Ram Nair 2 1 University of Michigan, 2 National Center for Atmospheric Research,
More informationA mass-conservative version of the semi- Lagrangian semi-implicit HIRLAM using Lagrangian vertical coordinates
A mass-conservative version of the semi- Lagrangian semi-implicit HIRLAM using Lagrangian vertical coordinates Peter Hjort Lauritzen Atmospheric Modeling & Predictability Section National Center for Atmospheric
More informationIFS DOCUMENTATION CY28r1
IFS Documentation Cycle CY28r1 IFS DOCUMENTATION CY28r1 III: DYNAMICS AND NUMERICAL PROCEDURES (Operational implementation 9 March 2004) Table of contents Chapter 1. Introduction Chapter 2. Basic equations
More informationA Spectral Element Eulerian-Lagrangian Atmospheric Model (SEELAM)
A Spectral Element Eulerian-Lagrangian Atmospheric Model SEELAM Francis X. Giraldo Naval Research Laboratory Monterey, CA 93943 E-mail: giraldo@nrlmry.navy.mil A new dynamical core for numerical weather
More informationAn explicit and conservative remapping strategy for semi-lagrangian advection
An explicit and conservative remapping strategy for semi-lagrangian advection Sebastian Reich Universität Potsdam, Potsdam, Germany January 17, 2007 Abstract A conservative semi-lagrangian advection scheme
More informationIFS Documentation Cycle CY25r1 IFS DOCUMENTATION PROCEDURES (CY25R1) (Operational implementation 9 April 2002) Edited by Peter W.
IFS Documentation Cycle CY25r IFS DOCUMENTATION Part III: DYNAMICS AND NUMERICAL PROCEDURES (CY25R) (Operational implementation 9 April 2002) Edited by Peter W. White (Text written and updated by members
More informationAdvanced Numerical Methods for Numerical Weather Prediction
Advanced Numerical Methods for Numerical Weather Prediction Francis X. Giraldo Naval Research Laboratory Monterey, CA 93943-5502 phone: (831) 656-4882 fax: (831) 656-4769 e-mail: giraldo@nrlmry.navy.mil
More informationA test case for the inviscid shallow-water equations on the sphere
Quarterly Journal of the Royal Meteorological Society Q. J. R. Meteorol. Soc. 142: 488 495, January 216 A DOI:1.12/qj.2667 A test case for the inviscid shallow-water equations on the sphere R. K. Scott,
More informationBlock-structured Adaptive Meshes and Reduced Grids for Atmospheric General Circulation Models
Block-structured Adaptive Meshes and Reduced Grids for Atmospheric General Circulation Models B Y CHRISTIANE JABLONOWSKI 1, ROBERT C. OEHMKE 2, QUENTIN F. STOUT 1 1 University of Michigan, 2455 Hayward
More informationATM 298, Spring 2013 Lecture 4 Numerical Methods: Horizontal DiscreDzaDons April 10, Paul A. Ullrich (HH 251)
ATM 298, Spring 2013 Lecture 4 Numerical Methods: Horizontal DiscreDzaDons April 10, 2013 Paul A. Ullrich (HH 251) paullrich@ucdavis.edu Outline 1. Introduction / Motivation 2. Finite Difference Methods
More informationAdvanced Numerical Methods for Numerical Weather Prediction
Advanced Numerical Methods for Numerical Weather Prediction Francis X. Giraldo Naval Research Laboratory Monterey, CA 93943-5502 phone: (831) 656-4882 fax: (831) 656-4769 e-mail: giraldo@nrlmry.navy.mil
More informationPART III: DYNAMICS AND NUMERICAL PROCEDURES
IFS DOCUMENTATION Cy40r1 Operational implementation 22 November 2013 PART III: DYNAMICS AND NUMERICAL PROCEDURES Table of contents Chapter 1 Chapter 2 Chapter 3 References Introduction Basic equations
More informationNumerically Converged Solutions of the Global Primitive Equations for Testing the Dynamical Core of Atmospheric GCMs
NOVEMBER 4 POLVANI ET AL. 539 Numerically Converged Solutions of the Global Primitive Equations for Testing the Dynamical Core of Atmospheric GCMs L. M. POLVANI Department of Applied Physics and Applied
More informationAn explicit and conservative remapping strategy for semi-lagrangian advection
ATMOSPHERIC SCIENCE LETTERS Atmos. Sci. Let. 8: 58 63 (2007) Published online 22 May 2007 in Wiley InterScience (www.interscience.wiley.com).151 An explicit and conservative remapping strategy for semi-lagrangian
More informationFluent User Services Center
Solver Settings 5-1 Using the Solver Setting Solver Parameters Convergence Definition Monitoring Stability Accelerating Convergence Accuracy Grid Independence Adaption Appendix: Background Finite Volume
More informationGlobal Numerical Weather Predictions and the semi-lagrangian semi-implicit dynamical 1 / 27. ECMWF forecast model
Global Numerical Weather Predictions and the semi-lagrangian semi-implicit dynamical core of the ECMWF forecast model Michail Diamantakis ECMWF Symposium on Modelling and Computations in Atmospheric Sciences
More informationAdvanced Numerical Methods for NWP Models
Advanced Numerical Methods for NWP Models Francis X. Giraldo Naval Postgraduate School Monterey, CA 93943-5502 phone: (831) 656-2293 fax: (831) 656-2355 e-mail: fxgirald@nps.edu Award #: N0001406WX20257
More informationGlobal non-hydrostatic modelling using Voronoi meshes: The MPAS model
Global non-hydrostatic modelling using Voronoi meshes: The MPAS model William C. Skamarock, Joseph B. Klemp, Michael Duda, Laura Fowler and Sang-Hun Park National Center for Atmospheric Research Boulder,
More informationExtending scalability of the community atmosphere model
Journal of Physics: Conference Series Extending scalability of the community atmosphere model To cite this article: A Mirin and P Worley 2007 J. Phys.: Conf. Ser. 78 012082 Recent citations - Evaluation
More informationMulti-Domain Pattern. I. Problem. II. Driving Forces. III. Solution
Multi-Domain Pattern I. Problem The problem represents computations characterized by an underlying system of mathematical equations, often simulating behaviors of physical objects through discrete time
More informationNumerical Methods for (Time-Dependent) HJ PDEs
Numerical Methods for (Time-Dependent) HJ PDEs Ian Mitchell Department of Computer Science The University of British Columbia research supported by National Science and Engineering Research Council of
More informationTable of contents for: Waves and Mean Flows by Oliver Bühler Cambridge University Press 2009 Monographs on Mechanics. Contents.
Table of contents for: Waves and Mean Flows by Oliver Bühler Cambridge University Press 2009 Monographs on Mechanics. Preface page 2 Part I Fluid Dynamics and Waves 7 1 Elements of fluid dynamics 9 1.1
More informationTrajectory Calculations for Spherical Geodesic Grids in Cartesian Space
See discussions, stats, author profiles for this publication at: http://www.researchgate.net/publication/4961110 Trajectory Calculations for Spherical Geodesic Grids in Cartesian Space ARTICLE in MONTHLY
More informationPART III: DYNAMICS AND NUMERICAL PROCEDURES
IFS DOCUMENTATION Cy43r1 Operational implementation 22 Nov 2016 PART III: DYNAMICS AND NUMERICAL PROCEDURES Table of contents Chapter 1 Chapter 2 Chapter 3 References Introduction Basic equations and discretisation
More informationLong time integrations of a convective PDE on the sphere by RBF collocation
Long time integrations of a convective PDE on the sphere by RBF collocation Bengt Fornberg and Natasha Flyer University of Colorado NCAR Department of Applied Mathematics Institute for Mathematics Applied
More informationSemi-Lagrangian Advection. The Basic Idea. x 1.
3.2.6. Semi-Lagrangian Advection We have studied the Eulerian leapfrog scheme and found it to be conditionally stable. The criterion for stability was the CFL condition µ c t x 1. For high spatial resolution
More informationContinued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data
Continued Investigation of Small-Scale Air-Sea Coupled Dynamics Using CBLAST Data Dick K.P. Yue Center for Ocean Engineering Department of Mechanical Engineering Massachusetts Institute of Technology Cambridge,
More informationM.Sc. in Computational Science. Fundamentals of Atmospheric Modelling
M.Sc. in Computational Science Fundamentals of Atmospheric Modelling Peter Lynch, Met Éireann Mathematical Computation Laboratory (Opp. Room 30) Dept. of Maths. Physics, UCD, Belfield. January April, 2004.
More informationA Multiscale Nested Modeling Framework to Simulate the Interaction of Surface Gravity Waves with Nonlinear Internal Gravity Waves
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. A Multiscale Nested Modeling Framework to Simulate the Interaction of Surface Gravity Waves with Nonlinear Internal Gravity
More informationChapter 6. Semi-Lagrangian Methods
Chapter 6. Semi-Lagrangian Methods References: Durran Chapter 6. Review article by Staniford and Cote (1991) MWR, 119, 2206-2223. 6.1. Introduction Semi-Lagrangian (S-L for short) methods, also called
More informationThe semi-lagrangian technique: current status and future 4 September developments / 42
The semi-lagrangian technique: current status and future developments Michail Diamantakis ECMWF 4 September 2013 (thanks to: Adrian Simmons, Niels Borman and Richard Forbes) The semi-lagrangian technique:
More informationUniversity of Reading. Discretization On Non-uniform Meshes: Tests solving shallow-water equations
University of Reading School of Mathematics, Meteorology and Physics Discretization On Non-uniform Meshes: Tests solving shallow-water equations By FAWZI B ALBUSAIDI August 2008 This dissertation is submitted
More informationThe Evolution of Dynamical Cores for Global Atmospheric Models
Journal of the Meteorological Society of Japan, Vol. 85B, pp. 241--269, 2007 241 The Evolution of Dynamical Cores for Global Atmospheric Models David L. WILLIAMSON National Center for Atmospheric Research*,
More informationIntroduction to Computational Fluid Dynamics Mech 122 D. Fabris, K. Lynch, D. Rich
Introduction to Computational Fluid Dynamics Mech 122 D. Fabris, K. Lynch, D. Rich 1 Computational Fluid dynamics Computational fluid dynamics (CFD) is the analysis of systems involving fluid flow, heat
More informationIntroduction to C omputational F luid Dynamics. D. Murrin
Introduction to C omputational F luid Dynamics D. Murrin Computational fluid dynamics (CFD) is the science of predicting fluid flow, heat transfer, mass transfer, chemical reactions, and related phenomena
More informationLagrangian methods and Smoothed Particle Hydrodynamics (SPH) Computation in Astrophysics Seminar (Spring 2006) L. J. Dursi
Lagrangian methods and Smoothed Particle Hydrodynamics (SPH) Eulerian Grid Methods The methods covered so far in this course use an Eulerian grid: Prescribed coordinates In `lab frame' Fluid elements flow
More informationCFD MODELING FOR PNEUMATIC CONVEYING
CFD MODELING FOR PNEUMATIC CONVEYING Arvind Kumar 1, D.R. Kaushal 2, Navneet Kumar 3 1 Associate Professor YMCAUST, Faridabad 2 Associate Professor, IIT, Delhi 3 Research Scholar IIT, Delhi e-mail: arvindeem@yahoo.co.in
More informationMODELED MESOSCALE GRAVITY WAVES: CONTINUOUS SPECTRUM AND ENERGY CASCADE Chungu Lu 1, Steven Koch, Fuqing Zhang, and Ning Wang
P6.5 MODELED MESOSCALE GRAVITY WAVES: CONTINUOUS SPECTRUM AND ENERGY CASCADE Chungu Lu 1, Steven Koch, Fuqing Zhang, and Ning Wang NOAA Research-Forecast Systems Laboratory, Boulder, Colorado [In collaboration
More informationAccuracy analysis of mimetic finite volume operators on geodesic grids and a consistent alternative
Accuracy analysis of mimetic finite volume operators on geodesic grids and a consistent alternative Pedro S. Peixoto College of Engineering, Mathematics and Physical Sciences - University of Exeter Instituto
More informationInterpolation of gridded data using Geostrophic Balance. Lennon O Naraigh *
Interpolation of gridded data using Geostrophic Balance Lennon O Naraigh * Date: 6 th September 2004 Abstract: In this report, a method of interpolating data from the grid of a global climate model (GCM)
More informationRecent research for dynamical cores of nonhydrostatic, deep-atmosphere, unified models
Recent research for dynamical cores of nonhydrostatic, deep-atmosphere, unified models Andrew Staniforth and Nigel Wood Met Office, FitzRoy Road, Exeter EX1 3PB, United Kingdom Andrew.Staniforth@metoffice.gov.uk,
More informationSurface effects in QG dynamics
Surface effects in QG dynamics Ben Harvey b.j.harvey@reading.ac.uk Supervisor: Maarten Ambaum 2/11/2010 Introduction Surface Outline QG dynamics Surface QG (temperature) 2-d Euler (vorticity) 2/11/2010
More informationThe CSLaM transport scheme & new numerical mixing diagnostics
The CSLaM transport scheme & new numerical mixing diagnostics CESM Ocean Model Working Group Meeting 14 15 December, 2011 Peter Hjort Lauritzen Outline: Two parts Briefly introduce the CSLaM scheme (Conservative
More informationA Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere
A Semi-Lagrangian Discontinuous Galerkin (SLDG) Conservative Transport Scheme on the Cubed-Sphere Ram Nair Computational and Information Systems Laboratory (CISL) National Center for Atmospheric Research
More informationCS205b/CME306. Lecture 9
CS205b/CME306 Lecture 9 1 Convection Supplementary Reading: Osher and Fedkiw, Sections 3.3 and 3.5; Leveque, Sections 6.7, 8.3, 10.2, 10.4. For a reference on Newton polynomial interpolation via divided
More informationGlobal Stokes Drift and Climate Wave Modeling
Global Stokes Drift and Climate Wave Modeling Adrean Webb University of Colorado, Boulder Department of Applied Mathematics February 20, 2012 In Collaboration with: Research funded by: Baylor Fox-Kemper,
More informationA 3-D Finite-Volume Nonhydrostatic Icosahedral Model (NIM) Jin Lee
A 3-D Finite-Volume Nonhydrostatic Icosahedral Model (NIM) Jin Lee Earth System Research Laboratory(ESRL) Director Dr. A.E. (Sandy) MacDonald GFDLNSSLARLAOMLGLERLPMEL Aeronomy Lab. Climate Diagnostic center
More informationAdaptive Grids for Weather and Climate Models
Adaptive Grids for Weather and Climate Models C. Jablonowski National Center for Atmospheric Research, 185 Table Mesa Drive, Boulder, CO 835, USA cjablono@ucar.edu M. Herzog 1, J. E. Penner 1, R. C. Oehmke
More informationA finite-volume integration method for computing pressure gradient force in general vertical coordinates
Q. J. R. Meteorol. Soc. (1997), 123, pp. 1749-1 762 A finite-volume integration method for computing pressure gradient force in general vertical coordinates By SHIAN-JIANN LIN* University of Maryland,
More informationRegional Cooperation for Limited Area Modeling in Central Europe. Dynamics in LACE. Petra Smolíková thanks to many colleagues
Dynamics in LACE Petra Smolíková thanks to many colleagues Outline 1 1. VFE NH Design of vertical finite elements scheme for NH version of the model 2. SL scheme Application of ENO technique to SL interpolations
More informationA Polygon-based Line Integral Method for Calculating Vorticity, Divergence, and Deformation from Non-uniform Observations
2.7 olygon-based Line Integral Method for Calculating Vorticity, Divergence, and Deformation from Non-uniform Observations Charles N. Helms and Robert E. Hart Florida State University, Tallahassee, Florida
More informationMulti-Mesh CFD. Chris Roy Chip Jackson (1 st year PhD student) Aerospace and Ocean Engineering Department Virginia Tech
Multi-Mesh CFD Chris Roy Chip Jackson (1 st year PhD student) Aerospace and Ocean Engineering Department Virginia Tech cjroy@vt.edu May 21, 2014 CCAS Program Review, Columbus, OH 1 Motivation Automated
More informationAbout the SPEEDY model (from Miyoshi PhD Thesis):
SPEEDY EXPERIMENTS. About the SPEEDY model (from Miyoshi PhD Thesis): The SPEEDY model (Molteni 2003) is a recently developed atmospheric general circulation model (AGCM) with a spectral primitive-equation
More informationIMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL RESERVOIR SIMULATIONS USING THE CIP SCHEME WITH THIRD-ORDER ACCURACY
PROCEEDINGS, Thirty-Seventh Workshop on Geothermal Reservoir Engineering Stanford University, Stanford, California, January 30 - February 1, 2012 SGP-TR-194 IMPROVING THE NUMERICAL ACCURACY OF HYDROTHERMAL
More informationImpact of Consistent Semi-Lagrangian Trajectory Calculations on Numerical Weather Prediction Performance
OCTOBER 2017 H U S A I N A N D G I R A R D 4127 Impact of Consistent Semi-Lagrangian Trajectory Calculations on Numerical Weather Prediction Performance SYED ZAHID HUSAIN AND CLAUDE GIRARD Atmospheric
More informationIntroduction to ANSYS CFX
Workshop 03 Fluid flow around the NACA0012 Airfoil 16.0 Release Introduction to ANSYS CFX 2015 ANSYS, Inc. March 13, 2015 1 Release 16.0 Workshop Description: The flow simulated is an external aerodynamics
More informationThe Global Nonhydrostatic Model NICAM: numerics and performance of its dynamical core
The Global Nonhydrostatic Model NICAM: numerics and performance of its dynamical core H. Tomita and M. Satoh Frontier Research Center for Global Change 3173-25 Showamachi, Kanazawa-ku, Yokohama-city, Kanagawa,236-0001,
More informationSensitivity of resolved and parameterized surface drag to changes in resolution and parameterization
Sensitivity of resolved and parameterized surface drag to changes in resolution and parameterization Annelize van Niekerk Ted Shepherd With thanks to: Simon Vosper, Stuart Webster, Andy Elvidge, Irina
More informationCoastal Ocean Modeling & Dynamics
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Coastal Ocean Modeling & Dynamics Roger M. Samelson College of Oceanic and Atmospheric Sciences Oregon State University
More informationBaylor Fox-Kemper (Brown Geo.) CARTHE Spring 2014 All-Hands Meeting
Frontogenesis in the Presence of Stokes Forces Baylor Fox-Kemper (Brown Geo.) with Jim McWilliams (UCLA), Nobu Suzuki (Brown), and Qing Li (Brown) Expanding on past work with: Peter Hamlington (CU-Boulder),
More informationCGT 581 G Fluids. Overview. Some terms. Some terms
CGT 581 G Fluids Bedřich Beneš, Ph.D. Purdue University Department of Computer Graphics Technology Overview Some terms Incompressible Navier-Stokes Boundary conditions Lagrange vs. Euler Eulerian approaches
More informationProgress on Advanced Dynamical Cores for the Community Atmosphere Model. June 2010
Progress on Advanced Dynamical Cores for the Community Atmosphere Model June 2010 Art Mirin, P. O. Box 808, Livermore, CA 94551 This work performed under the auspices of the U.S. Department of Energy by
More informationNew formulations of the semi-lagrangian method for Vlasov-type equations
New formulations of the semi-lagrangian method for Vlasov-type equations Eric Sonnendrücker IRMA Université Louis Pasteur, Strasbourg projet CALVI INRIA Nancy Grand Est 17 September 2008 In collaboration
More informationSystem Design for Visualizing Scientific Computations
25 Chapter 2 System Design for Visualizing Scientific Computations In Section 1.1 we defined five broad goals for scientific visualization. Specifically, we seek visualization techniques that 1. Can be
More information1.2 Numerical Solutions of Flow Problems
1.2 Numerical Solutions of Flow Problems DIFFERENTIAL EQUATIONS OF MOTION FOR A SIMPLIFIED FLOW PROBLEM Continuity equation for incompressible flow: 0 Momentum (Navier-Stokes) equations for a Newtonian
More informationThe Ventilated Ocean: What controls ocean stratification and overturning circulation and Lagrangian ocean modeling. Alexey Fedorov
The Ventilated Ocean: What controls ocean stratification and overturning circulation and Lagrangian ocean modeling Alexey Fedorov Many thanks to Patrick Haertel Yale University July 2011 1 What controls
More informationHybrid Scheme for Three-Dimensional Advective Transport
JUNE 1999 ZUBOV ET AL. 1335 Hybrid Scheme for Three-Dimensional Advective Transport VLADIMIR A. ZUBOV Main Geophysical Observatory, St. Petersburg, Russia EUGENE V. ROZANOV AND MICHAEL E. SCHLESINGER Climate
More informationQuantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Quantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement Lian Shen Department of Mechanical Engineering
More informationExample 24 Spring-back
Example 24 Spring-back Summary The spring-back simulation of sheet metal bent into a hat-shape is studied. The problem is one of the famous tests from the Numisheet 93. As spring-back is generally a quasi-static
More informationQuantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Quantifying the Dynamic Ocean Surface Using Underwater Radiometric Measurement Lian Shen Department of Mechanical Engineering
More informationA spectral element shallow water model on spherical geodesic grids
INTERNATIONAL JOURNAL FOR NUMERICAL METHODS IN FLUIDS Int. J. Numer. Meth. Fluids 2001; 35: 869 901 A spectral element shallow water model on spherical geodesic grids Francis X. Giraldo*,1 Naal Research
More informationThe diabatic contour advective. semi-lagrangian algorithms for the spherical. shallow water equations
The diabatic contour advective semi-lagrangian algorithms for the spherical shallow water equations Ali R. Mohebalhojeh Institute of Geophysics, University of Tehran, Tehran, Iran David G. Dritschel School
More informationAn inherently mass-conservative semi-implicit semi-lagrangian model. Peter Hjort Lauritzen
An inherently mass-conservative semi-implicit semi-lagrangian model by Peter Hjort Lauritzen ph.d. thesis Department of Geophysics University of Copenhagen Denmark September, 2005 Abstract A locally mass-conservative
More informationDependence of aqua-planet simulations on time step
Q. J. R. Meteorol. Soc. (2003), 129, pp. 2049 2064 doi: 10.1256/qj.02.62 Dependence of aqua-planet simulations on time step By DAVID L. WILLIAMSON and JERRY G. OLSON National Center for Atmospheric Research,
More informationApplication of Finite Volume Method for Structural Analysis
Application of Finite Volume Method for Structural Analysis Saeed-Reza Sabbagh-Yazdi and Milad Bayatlou Associate Professor, Civil Engineering Department of KNToosi University of Technology, PostGraduate
More informationUsing a Single Rotating Reference Frame
Tutorial 9. Using a Single Rotating Reference Frame Introduction This tutorial considers the flow within a 2D, axisymmetric, co-rotating disk cavity system. Understanding the behavior of such flows is
More informationDriven Cavity Example
BMAppendixI.qxd 11/14/12 6:55 PM Page I-1 I CFD Driven Cavity Example I.1 Problem One of the classic benchmarks in CFD is the driven cavity problem. Consider steady, incompressible, viscous flow in a square
More informationVI Workshop Brasileiro de Micrometeorologia
Validation of a statistic algorithm applied to LES model Eduardo Bárbaro, Amauri Oliveira, Jacyra Soares November 2009 Index Objective 1 Objective 2 3 Vertical Profiles Flow properties 4 Objective 1 The
More informationAn explicit flux-form semi-lagrangian shallow-water model on the sphere
Q. J. R. Meteorol. Soc. (1997), 123, pp. 2477-2498 An explicit flux-form semi-lagrangian shallow-water model on the sphere By SHIAN-JIA" LI"* and RICHARD B. ROOD2 'Goddard Space Flight Center and University
More information21. Efficient and fast numerical methods to compute fluid flows in the geophysical β plane
12th International Conference on Domain Decomposition Methods Editors: Tony Chan, Takashi Kako, Hideo Kawarada, Olivier Pironneau, c 2001 DDM.org 21. Efficient and fast numerical methods to compute fluid
More informationCalculate a solution using the pressure-based coupled solver.
Tutorial 19. Modeling Cavitation Introduction This tutorial examines the pressure-driven cavitating flow of water through a sharpedged orifice. This is a typical configuration in fuel injectors, and brings
More informationLab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders
Lab 9: FLUENT: Transient Natural Convection Between Concentric Cylinders Objective: The objective of this laboratory is to introduce how to use FLUENT to solve both transient and natural convection problems.
More informationA Direct Simulation-Based Study of Radiance in a Dynamic Ocean
A Direct Simulation-Based Study of Radiance in a Dynamic Ocean Lian Shen Department of Civil Engineering Johns Hopkins University Baltimore, MD 21218 phone: (410) 516-5033 fax: (410) 516-7473 email: LianShen@jhu.edu
More informationQUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD
QUASI-3D SOLVER OF MEANDERING RIVER FLOWS BY CIP-SOROBAN SCHEME IN CYLINDRICAL COORDINATES WITH SUPPORT OF BOUNDARY FITTED COORDINATE METHOD Keisuke Yoshida, Tadaharu Ishikawa Dr. Eng., Tokyo Institute
More informationICON Performance Benchmark and Profiling. March 2012
ICON Performance Benchmark and Profiling March 2012 Note The following research was performed under the HPC Advisory Council activities Participating vendors: Intel, Dell, Mellanox Compute resource - HPC
More informationFlow Structures of Jupiter s Great Red Spot Extracted by Using Optical Flow Method
Flow Structures of Jupiter s Great Red Spot Extracted by Using Optical Flow Method Tianshu Liu & Bo Wang Department of Mechanical & Aeronautical Engineering Western Michigan University, Kalamazoo, MI 49008
More informationAn added mass partitioned algorithm for rigid bodies and incompressible flows
An added mass partitioned algorithm for rigid bodies and incompressible flows Jeff Banks Rensselaer Polytechnic Institute Overset Grid Symposium Mukilteo, WA October 19, 216 Collaborators Bill Henshaw,
More informationPartial Differential Equations
Simulation in Computer Graphics Partial Differential Equations Matthias Teschner Computer Science Department University of Freiburg Motivation various dynamic effects and physical processes are described
More informationExtension of NHWAVE to Couple LAMMPS for Modeling Wave Interactions with Arctic Ice Floes
DISTRIBUTION STATEMENT A. Approved for public release; distribution is unlimited. Extension of NHWAVE to Couple LAMMPS for Modeling Wave Interactions with Arctic Ice Floes Fengyan Shi and James T. Kirby
More informationcuibm A GPU Accelerated Immersed Boundary Method
cuibm A GPU Accelerated Immersed Boundary Method S. K. Layton, A. Krishnan and L. A. Barba Corresponding author: labarba@bu.edu Department of Mechanical Engineering, Boston University, Boston, MA, 225,
More informationStudies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization
Studies of the Continuous and Discrete Adjoint Approaches to Viscous Automatic Aerodynamic Shape Optimization Siva Nadarajah Antony Jameson Stanford University 15th AIAA Computational Fluid Dynamics Conference
More informationCompressible Flow in a Nozzle
SPC 407 Supersonic & Hypersonic Fluid Dynamics Ansys Fluent Tutorial 1 Compressible Flow in a Nozzle Ahmed M Nagib Elmekawy, PhD, P.E. Problem Specification Consider air flowing at high-speed through a
More informationSupplemental Material Deep Fluids: A Generative Network for Parameterized Fluid Simulations
Supplemental Material Deep Fluids: A Generative Network for Parameterized Fluid Simulations 1. Extended Results 1.1. 2-D Smoke Plume Additional results for the 2-D smoke plume example are shown in Figures
More informationRK-Time Integration and High-Order Spatial Discretization a New Dynamical Core for the LMK
RK-Time Integration and High-Order Spatial Discretization a New Dynamical Core for the LMK Jochen Förstner Deutscher Wetterdienst, Kaiserleistraße 35, 63067 Offenbach am Main jochen.foerstner@dwd.de COSMO
More informationUsing the Eulerian Multiphase Model for Granular Flow
Tutorial 21. Using the Eulerian Multiphase Model for Granular Flow Introduction Mixing tanks are used to maintain solid particles or droplets of heavy fluids in suspension. Mixing may be required to enhance
More informationInherently Conservative Nonpolynomial-Based Remapping Schemes: Application to Semi-Lagrangian Transport
5044 M O N T H L Y W E A T H E R R E V I E W VOLUME 136 Inherently Conservative Nonpolynomial-Based Remapping Schemes: Application to Semi-Lagrangian Transport MATTHEW R. NORMAN Department of Marine, Earth,
More informationOptimization of a two-link Robotic Manipulator
Optimization of a two-link Robotic Manipulator Zachary Renwick, Yalım Yıldırım April 22, 2016 Abstract Although robots are used in many processes in research and industry, they are generally not customized
More informationEssay 1: Dimensional Analysis of Models and Data Sets: Similarity Solutions
Table of Contents Essay 1: Dimensional Analysis of Models and Data Sets: Similarity Solutions and Scaling Analysis 1 About dimensional analysis 4 1.1 Thegoalandtheplan... 4 1.2 Aboutthisessay... 5 2 Models
More informationJ. Vira, M. Sofiev SILAM winter school, February 2013, FMI
Numerical aspects of the advection-diffusion equation J. Vira, M. Sofiev SILAM winter school, February 2013, FMI Outline Intro Some common requirements for numerical transport schemes Lagrangian approach
More informationTutorial 1. Introduction to Using FLUENT: Fluid Flow and Heat Transfer in a Mixing Elbow
Tutorial 1. Introduction to Using FLUENT: Fluid Flow and Heat Transfer in a Mixing Elbow Introduction This tutorial illustrates the setup and solution of the two-dimensional turbulent fluid flow and heat
More information