210 likes | 301 Views
Modelos Hidrodinâmicos. Aula 4 Equations for 3D and 2D Hydrodynamic Models. Parameters and Boundary and Initial Conditions. Mass conservation. If P is the volumic mass , that has no Sources or Sinkes and has no diffusion because the net movement of molecules is the velocity ….
E N D
ModelosHidrodinâmicos Aula 4 Equations for 3D and 2D Hydrodynamic Models. Parameters and Boundary and Initial Conditions
Mass conservation If P isthevolumicmass, thathas no SourcesorSinkesandhas no diffusionbecausethe net movementofmoleculesisthevelocity…. Ifincompressible:
Momentum Conservation Se P for a quantidade de movimento por unidade de volume: Sourcesandsinks are Pressure forces (gravitationalis zero becausewe are interestedonlyon horizontal momentum)
ShallowWaterEquations • HydrostaticPressure (vertical acelerationnegligeable). • If “z” isthe vertical axis, pointingupwards:
Using the Leibnitz rule these equations can be integrated on vertical to obtain the equations of a 2D model.
The 2D case TheAccumulation rate = flows in – flows out
1D Case TheAccumulation rate = flows in – flows out
Momentum: 1D Case Ls Horizontal diffusion is negligible compared to vertical diffusion A Lb: wet perimeter
The 1D Spatial Grid Qi-1 zi-1 Qi zi Qi+1
Discretization A staggeredgridisconvenient. Temporal discretization can beexplicit, implicitou Crank-Nicholson
Stability • Explicit (1D): • Implicit: Incondicionally stable • Explicit 2D:
Boundary Conditions z0 Q1 z1 Q2 z2 • One can impose Free Surface levels and compute discharges or vice versa. • On sea side level is easier to know (tide) and on the land side river discharge use to be easier.
Other boundary conditions • Bathymetry! • Surface shear stress, • Diffusive fluxes, • Advective fluxes.
Initial conditions • Discharges/velocities, • Levels. • The good thing is that dissipative systems have low memory. Approximate initial conditions can be used. Usually zero velocity and horizontal free surface.
Parameters • Friction coefficient, • Diffusion coefficient. • Surface friction coefficient if flux in not known (e.g. from a meteorological model).